commande AIC() et ddl - modele glmer

[Raphaëlle Métras]

Bonjour,

Il y aurait-il un moyen d'extraire de facon systematique l'AIC dun modele glmer : La commande AIC() ne semble pas marcher.

Aussi comment peut-on faire pour extraire les degres de liberte du modele glmer pour ensuite calculer le parametre de dispersion (deviance du modele / degre de liberte). Je ne trouve pas cette valeur dans le summary () - en revanche la deviance est affichee.

Merci !

Raphaelle

[Aurélien Madouasse]

Bonjour,

Concernant l'AIC, il est possible de l'obtenir avec:

Code : Tout sélectionner

fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy)

summary(fm1)@AICtab[1]

En ce qui concerne les degrés de liberté, je crois qu'ils ne sont pas calculés pour les raisons indiquées ici:

https://stat.ethz.ch/pipermail/r-help/2006-May/094765.html


Aurélien

[Raphaëlle Métras]

Merci Aurelien pour la commande et le lien.
Bonne journee !