Je poursuis ma découverte des méthodes loglinéaires. J'aimerais à présent évaluer la significativité des effets à l'aide de tests de Wald. Comment faire quand les effets sont considérés comme redondants et n'apparaissent pas dans summary.glm(). Par exemple, soit la table
Code : Tout sélectionner
MI MM PPD Freq
1 a w 0 81
2 b w 0 159
3 a x 0 251
4 b x 0 283
5 a y 0 8
6 b y 0 20
7 a z 0 12
8 b z 0 20
9 a w 1 13
10 b w 1 74
11 a x 1 24
12 b x 1 77
13 a y 1 4
14 b y 1 8
15 a z 1 2
16 b z 1 9
Je sélectionne le modèle
Code : Tout sélectionner
Call:
glm(formula = Freq ~ MM + MI + PPD + MM:MI + MM:PPD + MI:PPD,
family = poisson, data = data.frame(tmp))
Deviance Residuals:
1 2 3 4 5 6 7
0.090599 -0.064338 0.070038 -0.065776 -0.637618 0.449021 0.006014
8 9 10 11 12 13 14
-0.004654 -0.220857 0.094817 -0.222774 0.126876 1.196424 -0.637618
15 16
-0.014673 0.006946
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 4.38437 0.10714 40.922 < 2e-16 ***
MMx 1.13666 0.12174 9.337 < 2e-16 ***
MMy -2.08765 0.31298 -6.670 2.55e-11 ***
MMz -1.90120 0.29162 -6.519 7.06e-11 ***
MIb 0.68964 0.12842 5.370 7.86e-08 ***
PPD1 -1.75878 0.19733 -8.913 < 2e-16 ***
MMx:MIb -0.56131 0.14856 -3.778 0.000158 ***
MMy:MIb -0.09273 0.37231 -0.249 0.803308
MMz:MIb -0.17603 0.35318 -0.498 0.618188
MMx:PPD1 -0.53906 0.16876 -3.194 0.001402 **
MMy:PPD1 0.18407 0.37329 0.493 0.621933
MMz:PPD1 -0.02088 0.37722 -0.055 0.955847
MIb:PPD1 0.97780 0.18646 5.244 1.57e-07 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for poisson family taken to be 1)
Null deviance: 1502.8536 on 15 degrees of freedom
Residual deviance: 2.5916 on 3 degrees of freedom
AIC: 110.18
Number of Fisher Scoring iterations: 4
Je peux facilement déduire l'effet et la probabilité de PPD0 à partir de ceux de PPD1, puisque la variable est dichotomique. Mais comment faire concrètement quand la variable n'est pas dichotomique. Par exemple dans le cas de MMw ou de MMw:MIb? J'ai fouillé dans les objets R glm et summary.glm(), mais les effets redondants sont éliminés d'office. Comment les restituer?
Merci pour toute aide.