我在 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 偏差水平)。