Reprenons ici le code pour récupérer la base de la taille (et du poids) d’un ensemble de 200 personnes.
Davis=read.table( "http://socserv.socsci.mcmaster.ca/jfox/Books/Applied-Regression-2E/datasets/Davis.txt") Davis[12,c(2,3)]=Davis[12,c(3,2)] X=Davis$height n=length(X)
La première chose que l’on avait vu la semaine passée est le PP-plot,
On va ici ‘tester’ une loi normale (de même moyenne et de même variance que notre échantillon… faisons simple)
PP=data.frame(empirique=(1:n)/n, theorique=pnorm(sort(X),mean(X),sd(X)))