Statistique théorique et appliquée - Tome 2
Utilisation de R pour l'exemple 8.2.1
par Emmanuel Nowak
Il s'agit d'estimer la longueur moyenne du corps pour la population dont proviennent les mesures. Il faut tout d'abord importer les données fournies dans le fichier 's2e07021.txt', les longueurs des corps ayant été l'objet de l'exemple 7.2.1 :
| s2e07021 <- read.table("C:/Dagnelie/st2donn/txt.2/s2e07021.txt",sep="\t",header=T) summary(s2e07021) |
Long Min. :145.0 1st Qu.:150.3 Median :157.5 Mean :157.4 3rd Qu.:162.5 Max. :176.0 |
On peut maintenant estimer la longueur moyenne, ponctuellement comme cela figure déjà ci-dessus avec une décimale :
| attach(s2e07021) mean(Long) |
[1] 157.4091 |
et par intervalle de confiance à 95%, en utilisant une loi de Student :
| t.test(Long)$conf.int |
[1] 153.6403 161.1779 attr(,"conf.level") [1] 0.95 |
On peut également tester la normalité de la longueur du corps des vaches laitières dans la population à l'aide d'un test de Shapiro-Wilk (présenté au paragraphe 3.3.3), la validité des limites de confiance de la moyenne y étant conditionnée, d'autant moins cependant que la taille de l'échantillon est grande :
| shapiro.test(Long) |
Shapiro-Wilk normality test data: s2e07021$Long W = 0.9593, p-value = 0.4748 |
Avant de traiter d'autres exemples, ne pas oublier de détacher le fichier de travail :
| detach(s2e07021) |
| Haut de la page |
| Autres illustrations avec R |
Dernière mise à jour : juillet 2006