2013-03-21 41 views
0

我做了一個在matlab中的GOF測試,但兩個樣本Kolmogorov Smirnov測試給了我相反的結果。 在MATLAB我用這個代碼2樣本kolmogorov-smirnov gof matlab的不同結果

file=importdata('all_subj_1file.xls','\t'); 
[ill1]=xlsread('all_subj_1file.xls','ill1','B2:E1608'); 
ill1_frac=ill1(:,2); 
[a,b]=gamfit(ill1_frac); 
y=gamrnd(a(1),a(2),1607,1); %same length gamma 
[h, p, ks2stat]=kstest2(ill1_frac, y) 

但不時我得到不同的pks2stat太(可能是因爲我產生新的gamrnd我每次運行的代碼)。
問題是如何選擇Kolmogorov-Smirnov的結果?我應該多次運行代碼並選擇一個我最喜歡的代碼?但它們相反(如Ho或H1的偏差水平)。

回答

0

爲了獲得理論分佈的GOF,應該使用一個樣本Kolmogorov Smirnov檢驗來代替兩個樣本KS檢驗。