私は MATLAB で GOF テストを作成しましたが、2 つのサンプル コルモゴロフ スミルノフ テストでは反対の結果が得られました。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 を生成するためです)。
問題は、コルモゴロフ-スミルノフの結果をどのように選択するかです。コードを何度も実行して、最も気に入ったものを選択する必要がありますか? しかし、それらは反対です (Ho または H1 レベルの偏差のように)。