Je compare les modèle suivants :
mm1ws<- lmer(maisparceWS[,i+14]~ prec + prec/rep + (1|prec/rep:bloc) + (1|prec/code) + (1 | lignenord),maisparceWS)
mm2ws<- lmer(maisparceWS[,i+14]~ rep + (1|rep:bloc) + (1|code) + (1 | lignenord),maisparceWS)
mm3ws<- lmer(maisparceWS[,i+14]~ prec + prec/rep + prec/rep:bloc + (1|prec/code) + (1 | lignenord),maisparceWS)
mm4ws<- lmer(maisparceWS[,i+14]~ prec + (1|prec/rep:bloc) + (1|prec/code) + (1 | lignenord),maisparceWS)
J'obtiens des résultats différent des critères AIC et BIC selon que j'utilise la commande :
Code : Tout sélectionner
anova(mm1ws,mm2ws,mm3ws,mm4ws)Df AIC BIC logLik Chisq Chi Df Pr(>Chisq)
mm2ws 6 174.71 198.66 -81.353
mm4ws 8 152.17 184.10 -68.085 26.536 2 1.729e-06 ***
mm1ws 10 141.05 180.96 -60.523 15.125 2 0.0005197 ***
mm3ws 44 139.13 314.75 -25.564 69.918 34 0.0002776 ***
ou la commande
Code : Tout sélectionner
summary(mm1ws)@AICtab> summary(mm3ws)@AICtab
AIC BIC logLik deviance REMLdev
263.5587 439.1831 -87.77935 51.12799 175.5587
> summary(mm1ws)@AICtab
AIC BIC logLik deviance REMLdev
156.9296 196.8443 -68.4648 121.0458 136.9296
> summary(mm2ws)@AICtab
AIC BIC logLik deviance REMLdev
184.8749 208.8237 -86.43746 162.7061 172.8749
> summary(mm4ws)@AICtab
AIC BIC logLik deviance REMLdev
159.7908 191.7225 -71.8954 136.1704 143.7908
Selon la méthode je ne sélectionne pas du tout le même modèle (sur le critère AIC) :
dans le premier cas c'est le modèle mm3ws
dans le second c'est le modèle mm1ws
Quelle est la différence entre les deux méthodes ?
Quelle méthode dois-je utilisé pour la sélection de modèle?
Merci d'avance
Jérôme