Keyword: コルモゴロフ・スミルノフ, Kolmogorov-Smirnov Test, 検定, ２標本
概要
本サンプルは２標本コルモゴロフ・スミルノフ検定（the two sample Kolmogorov-Smirnov Test) を行うサンプルプログラムです。 本サンプルは以下に示される独立した２つの標本を分析対象とし、検定統計量DとZ統計量、裾確率を算出します。
※本サンプルはnAG Toolbox for MATLAB®が提供する関数 nag_2_sample_ks_test() のExampleコードです。実行にはMATLAB®本体（他社製品）とnAG Toolbox for MATLAB®が必要です。
本サンプル及び関数の詳細情報は g08cd のマニュアルページをご参照ください。
入力データ
x = [1.5902480498365;
0.4514344715575777;
0.7425605404715657;
0.4500701410943835;
1.757489613962739;
0.09494659709066995;
0.3611426417219349;
0.8655310991929152;
0.07874046899432821;
1.15034115657093;
1.909225474130498;
1.248205642566939;
1.278136611206321;
1.376242882235565;
0.4024793749242044;
1.762011956845783;
0.6952326353936773;
0.9876414433367825;
0.6059391432983038;
0.2078601007360235;
1.564305035404043;
1.746505783532216;
0.1778870249628505;
0.1076281399262869;
0.07704165044712488;
0.04567587466526592;
0.8039773560869385;
1.241984667916896;
1.797786082814494;
0.852819977142019;
1.747622472821148;
1.912044530456294;
0.4706479344188908;
1.243486362938038;
1.279342351641843;
0.8903608734478126;
1.727202408903865;
0.2581867825828414;
0.9415224497978657;
1.74941093213823;
1.20905525205517;
0.9257769700555145;
0.9774905558012059;
0.8557640721625349;
0.2648216307181972;
0.2305151298211215;
1.214582160870085;
0.0008087845523660937;
0.04826634252569287;
1.422324123840062;
1.64461234171838;
1.242902747480194;
1.859348219383043;
0.8759548065432361;
0.9770887783336585;
1.077045946574884;
0.07040286404180876;
0.725208635893021;
1.129721212172624;
0.7794878837807858;
1.236501744194477;
0.8327862848080944;
1.932002357334167;
0.01096554652191538;
1.430490703604279;
0.9004992406397488;
0.8208346310843287;
0.7213215653615261;
1.032756476884525;
0.3378646599255028;
1.975473898973496;
1.307202749307935;
1.022550018984354;
0.5703968003457041;
1.174687331604432;
1.691927534708742;
0.09038851332572448;
0.08328616335992274;
0.3974243748119088;
0.4000804015026145;
1.734779923254609;
0.072131066247919;
1.330366567063721;
0.08900990107193633;
0.6128580701836494;
1.25978558572682;
0.3681357787827774;
1.766612133578424;
0.5447167419613813;
0.1468692016200379;
0.08565369904623682;
1.73695078139111;
1.502452715832918;
1.023663015528883;
1.548981938796881;
0.06591087486198742;
1.089678519368579;
0.7863284085754246;
0.4083539838920379;
0.1336746295053808];
y = [0.3959382482212902;
2.189059814171972;
0.5309091506354565;
1.398213120735397;
1.226993173102978;
2.092465206711672;
1.376152675038049;
0.2635476162603962;
0.6826139872981392;
1.725437958005721;
1.218220941115947;
0.2907559334793868;
0.2888496459940526;
2.237322438311175;
1.335206853771786;
1.482575912633977;
1.511554656117706;
1.49444503365132;
1.647783675961666;
1.093415807808477;
1.929516829926758;
1.771633171736251;
1.362424985201225;
0.7600483503138998;
1.474848032784191;
1.579250523709051;
1.100106550316349;
0.9988674396189126;
1.974550189113563;
0.6728895918936394;
0.6568563195351251;
1.038967939288667;
1.420212415024878;
2.156198816080763;
1.879016986416847;
1.409577351753394;
1.953857139342458;
0.3220659541919986;
2.244466171121614;
1.738164359951565;
0.4534069003591212;
0.8715218615975232;
1.596006863571959;
1.501743503682336;
2.099525595470502;
0.8538189862774034;
0.9554670591088673;
1.668907067970278;
1.517236797575171;
1.478170841410264];
ntype = int32(1);
[d, z, p, sx, sy, ifail] = g08cd(x, y, ntype)
- x には一つめの標本の観測値を与えています。
- y には二つめの標本の観測値を与えています。
- ntype には計算する検定統計量の種類を指定しています。"1" は２つの標本累積分布関数間の最大絶定偏差を計算することを意味します。
- 最後に本関数を呼び出す構文を指定しています。
出力結果
d =
0.3600
z =
0.0624
p =
2.8438e-04
sx =
array elided
sy =
array elided
ifail =
0
- d は検定統計量Dを示しています。
- z はZ統計量を示しています。
- p は裾確率を示しています。
- sx は一つ目の標本の観測値を示していますが、ここでは出力が省略されています。
- sy は二つ目の標本の観測値を示していますが、ここでは出力が省略されています。
- ifail は関数がエラーを検知しなければ"0"が出力されます。