J'utilise un modèle de régression logistique sous la forme:
lmer(response~1+(1|site), family=binomial, REML = FALSE)
Normalement, je calculerais l'ICC à partir des variances d'interception et résiduelles, mais le résumé du modèle n'inclut pas la variance résiduelle. Comment puis-je calculer cela?
intercept_variance / (intercept_variance + pi^2/3)
- alors ne quadrillez pas la variance.Réponses:
Vous pouvez utiliser
icc()
la fonction-du paquet sjstats .Dans le fichier d'aide,
?sjstats::icc
vous trouverez une référence à la formule pour les modèles mixtes avec réponse binaire:La déviance résiduelle dans la régression logistique est fixée à
(pi ^ 2) / 3
.la source