672
Le logiciel R
13.2- note ik = μ + α i + ik o` u les ik sont des variables al´ eatoires i.i.d. suivant une
loi N(0, σ
2 ).
13.3- modele <- aov(note∼sonore)
summary(modele)
13.4- pairwise.t.test(note,sonore,p.adjust="bonf")
TukeyHSD(modele)
• ´
Etude sur l’intima-m´ edia
13.1- url <- "http://biostatisticien.eu/springeR/Intima_Media.txt"
intima.media <- read.table(url,header=TRUE,dec=",",sep=" ")
attach(intima.media)
13.2- plot(mesure∼factor(alcool))
13.3- modele.intima <- aov(mesure∼factor(alcool))
summary(modele.intima)
Il existe une diff´ erence de mesure de l’´ epaisseur de l’intima-m´ edia suivant
la consommation d’alcool (valeur-p ≤ 0.05).
13.4- plot(modele.intima)
# Homosc´ edasticit´ e :
require("car")
leveneTest(mesure,factor(alcool))
• ´
Etude chez les sportifs
13.1- T´ el´ echargez le fichier http://biostatisticien.eu/springeR/sportif.
RData dans votre r´ epertoire courant.
load("sportif.RData")
attach(sportif)
Facteur temps de pratique sportive (7 modalit´ es). score ik = μ + α i + ik o` u
les ik sont des variables al´ eatoires i.i.d. suivant une loi N(0, σ
2 ).
13.2- modele.sport <- aov(score∼temps)
summary(modele.sport)
13.3- cmat <- rbind(" : peu sportif versus moyennement sportif"=
c(1,1,-1,-1,0,0,0),
" : tr` es sportif versus peu sportif"=c(3,3,0,0,-2,-2,-2))
13.4- require("gmodels")
fit.contrast(modele.sport,temps,cmat)