668
Le logiciel R
Le R
2 a (l´ eg` erement) augment´ e.
B- ´
Etude sur la r´ egression lin´ eaire multiple
• ´
Etude sur l’intima-m´ edia
On constate (tabac[is.na(paqan)] et paqan[tabac==0]) que la variable
paqan (donnant le nombre de paquets de cigarettes par ann´ ee) pr´ esente des
valeurs manquantes NA ` a chaque fois que la variable tabac prend la valeur 0 (qui
signifie ne fume pas). Il est donc naturel de remplacer les donn´ ees manquantes
de la variable paquan par des 0 :
intima.media$paqan[is.na(paqan)] <- 0
Cr´ eation de la variable IMC :
IMC <- intima.media$poids/((intima.media$taille/100)^2)
attach(intima.media)
mon.data <- data.frame(AGE,SPORT,alcool,paqan,IMC,mesure)
12.1- pairs(mon.data)
Certains nuages de points ont une allure lin´ eaire, pouvant entraˆ ıner de la
colin´ earit´ e.
12.2- modele.age <- lm(mesure∼AGE,data=mon.data)
summary(modele.age)
modele.sport <- lm(mesure∼SPORT,data=mon.data)
summary(modele.sport) ### valeur-p > 0.25
modele.alcool <- lm(mesure∼factor(alcool),data=mon.data)
summary(modele.alcool)
modele.paqan <- lm(mesure∼paqan,data=mon.data)
summary(modele.paqan)
modele.IMC <- lm(mesure∼IMC,data=mon.data)
summary(modele.IMC)
12.3- modele.age.inter <- lm(mesure∼AGE * paqan)
# Interaction non-significative au seuil de 0.25:
summary(modele.age.inter)
modele.alcool.sans.inter <- lm(mesure∼factor(alcool)+paqan,
data=mon.data)
modele.alcool.inter <- lm(mesure∼factor(alcool) * paqan,
data=mon.data)
summary(modele.alcool.inter)
anova(modele.alcool.sans.inter,modele.alcool.inter)