[問題] MATLAB做JB test 對小樣本與大樣本的結果比較已回收
如題,想做JBtest(一種常態性檢定)
"當樣本數不足時(小樣本),會有過度拒絕的問題(樣本其實是常態分配)"
的結果驗証。 編碼如下:
clc
clear all
h=0;
i=10; %i為樣本數
for k=1:1000
A=normrnd(0.0513,1.1030,i,1);
sA=skewness(sampleA)
kA=kurtosis(sampleA)
jb=((sA.^2)./6+((kA-3).^2)./24).*i %jb統計量公式
if(jb<=0.0506356||jb>=7.37776) %不確定是單尾還雙尾
h=h+1;
end
end
結果做出來小樣本根本沒有這問題,想請問是那邊寫錯。
(i不管等於多少都沒啥差異出現)
另一個問題,如果是用內建函數寫
h=jbtest(A) 代替手寫的jb統計量(程式第九行),
效果也是一樣吧? (當然,我知道if 那邊也要改)
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 61.217.152.43
推
09/18 10:07, , 1F
09/18 10:07, 1F
→
09/18 17:48, , 2F
09/18 17:48, 2F