[問題] MATLAB做JB test 對小樣本與大樣本的結果比較已回收

看板MATLAB作者 (無聊人生)時間13年前 (2010/09/18 01:45), 編輯推噓1(101)
留言2則, 2人參與, 最新討論串1/1
如題,想做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
你觀察一下jb是不是一個iby1的矩陣
09/18 10:07, 1F

09/18 17:48, , 2F
應該不是,就我的寫法是每次迴圈都會更新jb,再檢定一次
09/18 17:48, 2F
文章代碼(AID): #1CawbFoW (MATLAB)