web-dev-qa-db-ja.com

Rのコルモゴロフ-スミルノフ検定

Kolmogorov-Smirnovテストを使用して、サンプルの正常性をテストしようとしました。これは私がすることの小さな簡単な例です:

x <- rnorm(1e5, 1, 2)
ks.test(x, "pnorm")

Rから得られる結果は次のとおりです。

        One-sample Kolmogorov-Smirnov test

data:  x
D = 0.3427, p-value < 2.2e-16
alternative hypothesis: two-sided

P値は非常に低いのに対し、検定は帰無仮説を受け入れる必要があります。

なぜ機能しないのか分かりません。

17
gagaouthu

ks.testhelp で指摘されているように、ks.test関数にpnormの引数を与える必要があります。平均と標準のばらつきを正確に把握していない場合、テストは標準のガウス分布で行われます。

ここに書くべきです:

ks.test(x, "pnorm", 1, 2) #or ks.test(x, "pnorm", mean=1, sd=2) 
22
Pop