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