我在matlab中进行了GOF测试,但两个样本Kolmogorov Smirnov测试给出了相反的结果。 在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偏差水平)。
答案 0 :(得分:0)
应使用一个样本Kolmogorov Smirnov检验代替两个样本KS检验,以获得理论分布的GOF。