web-dev-qa-db-fra.com

Test de Kolmogorov-Smirnov en R

J'ai essayé d'utiliser le test de Kolmogorov-Smirnov pour tester la normalité d'un échantillon. Voici un petit exemple simple de ce que je fais:

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

Voici le résultat que R me donne:

        One-sample Kolmogorov-Smirnov test

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

La valeur de p est très faible alors que le test devrait accepter l'hypothèse nulle.

Je ne comprends pas pourquoi cela ne fonctionne pas.

17
gagaouthu

Comme indiqué dans le ks.testaide , vous devez donner au ks.test fonction les arguments de pnorm. Si vous ne précisez pas la variation moyenne et standard, le test est effectué sur une distribution gaussienne standard.

Ici, vous devez écrire:

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