Statistique théorique et appliquée - Tome 2
Utilisation de R pour l'exemple 9.3.7
par Emmanuel Nowak
Il s'agit ici de comparer les hauteurs moyennes des arbres de trois types de hêtraies à l'aide du test non paramétrique de Kruskal et Wallis, ceci à titre d'exemple puisque ce problème a déjà été résolu à l'aide d'une analyse de la variance à un facteur dans les exemples 9.2.1 à 9.3.2. Il faut importer les données relatives à l'exemple 2.3.1, fournies dans le fichier 's2e02031.txt', et diviser les hauteurs par 10 afin d'obtenir des valeurs en mètres :
| s2e02031 <- read.table("C:/Dagnelie/st2donn/txt.2/s2e02031.txt",sep="\t",header=T) s2e02031$Haut <- s2e02031$Haut/10 s2e02031[1:5,] |
Types Haut 1 1 23.4 2 1 24.4 3 1 24.6 4 1 24.9 5 1 25.0 |
Le test de Kruskal-Wallis, comme l'anayse de la variance, met en évidence des différences hautement significatives :
| attach(s2e02031) kruskal.test(Haut~Types) |
Kruskal-Wallis rank sum test data: Haut by Types Kruskal-Wallis chi-squared = 9.3337, df = 2, p-value = 0.009402 |
Avant de traiter d'autres exemples, ne pas oublier de détacher le fichier de travail :
| detach(s2e02031) |
| Haut de la page |
| Autres illustrations avec R |
Dernière mise à jour : septembre 2006