Chapitre 12. R´ egression lin´ eaire simple et multiple
515
> modele6 <- lm(BWT∼AGE+SMOKE+AGE:SMOKE)
Astuce
Pour introduire un terme d’interaction entre deux variables X 1 et X 2
dans un mod` ele lin´ eaire, il suffit de taper X1 :X2. Notez que l’´ ecriture X1 * X2
dans une formule correspond ` a l’´ ecriture X1+X2+X1 :X2. Ainsi, l’appel pr´ ec´ edent aurait pu ˆ etre remplac´ e par :
> modele6 <- lm(BWT∼AGE * SMOKE)
La diff´ erence entre ces deux mod` eles peut ˆ etre visualis´ ee sur les deux graphiques qui suivent. Le premier pr´ esente deux droites parall` eles dont le coefficient directeur repr´ esente l’effet de l’AGE sur BWT qui est donc le mˆ eme chez les
fumeuses et les non fumeuses.
> co <- coef(modele5); a0 <- co[1] ; a1 <- co[1]+co[3]
> b <- co[2]
> fSMOKE <- as.factor(SMOKE)
> plot(BWT∼AGE, xlab="AGE en années", ylab=
+ "Poids de naissance (g)",main=expression(BWT∼"="∼
+ beta[0]+beta[1] * AGE+beta[2] * SMOKE+epsilon),
+ col = c(’blue’,’red’)[fSMOKE], pch=c(1,18)[fSMOKE])
> abline(a=a0,b,col="blue") ; abline(a=a1,b,col="red")