{"word_count": 626, "character_count": 4815}
Toujours pour répondre à une question posée par mail ("ça quoi ça sert le test lillie sous R ? "), un court billet. Avant de parler de Lillie1, il faut revenir sur Kolmogorov-Smirnov. L'idée est de testercontreJe prends une forme paramétrique pour la fonction que l'on cherche à tester car c'est ce qui est souvent fait en pratique. Et surtout, ça permettra de mieux comprendre (je pense) l'idée du test lilllie. L'idée du test repose sur le théorème de Glivenko-Cantelli. En fait, Andreï Kolmogorov et Vladimir Smirnov ont montré que l'on pouvait encadrer la vitesse de convergence,ce qui se montre en notant que  converge en loi (si l'hypothèse H0 était la bonne) vers un pont brownien changé de temps par la fonction quantile .La loi limite est appelée loi de Kolmogorv-Smirnov, et ne dépend pas de la loi sous-jacente. Par contre elle n'est qu'asymptotique. Pour récupérer la loi à distance finie, regardons les simulations suivantes. On notera en particulier que la loi de la différence ne dépend pas de la vraie loi des observations (même si c'est cette loi que l'on teste, comme toujours lorsque l'on cherche la loi sous l'hypothèse H0)&#62; n=20&#62; ns=100000&#62; SU=SN=SE=rep(NA,ns)&#62; for(i in 1:ns){+ X=runif(n)+ SU[i]=max(c(abs(rank(X)/n-X),abs((rank(X)-1)/n-X)))+ X=rnorm(n,0,1)+ SN[i]=max(c(abs(rank(X)/n-pnorm(X,0,1)),abs((rank(X)-1)/n-pnorm(X,0,1))))+ X=rexp(n,1)+ SE[i]=max(c(abs(rank(X)/n-pexp(X,1)),abs((rank(X)-1)/n-pexp(X,1))))+ }&#62; plot(density(SN),col="red",lwd=2)&#62; lines(density(SU),col="blue",lwd=2)&#62; lines(density(SE),col="green",lwd=2)avec la densité de la loi du maximum entre la fonction de répartition théorique et la fonction de répartition empirique enrouge pour l'échantillon Gaussien, enbleu pour l'échantillon uniforme et envert pour l'échantillon exponentiel. Peu importe la loi sous-jacente, la statistique de test est la même.On peut en particulier récupérer les quantiles de cette lois pour n=20 observations (la loi du test dépend simplement du nombre d'observations).&#62; quantile(SE,.95)95%0.2937023&#62; quantile(SE,.90)90%0.264661On peut d'ailleurs confirmer ces valeurs en regardant dans les tables (ici ou là), par exemple 0.294 à 95%,Mais ce test ne sert que pour juger de l'adéquation à une loi précise (que nous avions noté  ici). En pratique, on ne connaît pas cette loi, mais on peut l'estimer, en particulier pour les lois paramétriques. On peut alors chercher à tester autre chose, commecontrequi s'apparente plus à un test d'ajustement à une famille de lois. Ce test étant différent du précédant, il n'y a pas de raison d'avoir la même loi pour la statistique de test (qui néanmoins, à  peut rester celle que nous avions considérée, à savoir la norme infinie entre les fonctions de répartition). Graphiquement par exemple, on cherchait aurparavant à tester l'ajustement de la loirouge, et maintenant on va tester l'ajustement de la loibleue,qui va s'adapter à l'échantillon, contrairement au test initial,Le code ici, dans le cas Gaussien par exemple (la loi sous-jacente n'est plus neutre dans ce cas) donne&#62; n=20&#62; ns=100000&#62; SN=SNML=rep(NA,ns)&#62; for(i in 1:ns){+ X=rnorm(n,0,1)+ SN[i]  =max(c(abs(rank(X)/n-pnorm(X,0,1)),+               abs((rank(X)-1)/n-pnorm(X,0,1))))+ SNML[i]=max(c(abs(rank(X)/n-pnorm(X,mean(X),sd(X))),+               abs((rank(X)-1)/n-pnorm(X,mean(X),sd(X)))))+ }&#62; plot(density(SN),col="blue",lty=2)&#62; lines(density(SNML),col="red",lwd=2)Les densités sont sensiblement différentes,et les quantiles aussi,&#62; quantile(SN,.95)&#38;nbsp;     95%0.2943206&#62; quantile(SNML,.95)95%0.1923295Sous R, deux fonctions distinctes sont à utiliser. Par exemple l'ajustement à une loi  donne&#62;  X=rnorm(20)&#62; mean(X)[1] -0.553641&#62; sd(X)[1] 0.8293167&#62;  ks.test(X,"pnorm", mean=0,sd=1)One-sample Kolmogorov-Smirnov testdata:  XD = 0.3096, p-value = 0.03362alternative hypothesis: two-sided(on rejette ici l'ajustement à une loi normale centrée réduite) et l'ajustement à une loi  donne&#62; lillie.test(X)Lilliefors (Kolmogorov-Smirnov) normality testdata:  XD = 0.1519, p-value = 0.2629avec une p-value de l'ordre de 26%. Notons que cette p-value est interprétable, contrairement au test suivant, qui n'est pas valide,&#62;  ks.test(X,"pnorm", mean=mean(X),sd=sd(X))One-sample Kolmogorov-Smirnov testdata:  XD = 0.1519, p-value = 0.69alternative hypothesis: two-sidedEncore une fois, la statistique de test est la même, sauf que la loi de la statistique de test sous H0 n'est pas juste dans le dernier cas. En particulier, n sera beaucoup plus exigent dans le test car on a choisi la meilleure loi normale possible (à l'aide d'un critère de maximisation de vraisemblance). Sinon je renvoie ici pour un document sur le même sujet.1 Lillie pour Hubert Lilliefors qui, le premier, a proposé ce test