Pente à effet aléatoire non corrélée
Dans l’exercice précédent, vous avez utilisé le paramétrage par défaut de lme4 et supposé que les pentes et les intercepts au sein de chaque groupe étaient corrélés pour les estimations d’effets aléatoires. Cependant, cette hypothèse n’est pas toujours valide, et il peut être utile de simplifier le modèle si vous rencontrez des difficultés de convergence numérique.
Construire un modèle avec des effets aléatoires non corrélés est une façon potentielle de simplifier le modèle. Par ailleurs, les modèles lmer() peuvent être difficiles à ajuster ; examiner les sorties du modèle est donc une étape utile pour diagnostiquer votre modèle. Autre possibilité : votre expertise métier peut vous amener à supposer que les effets aléatoires ne sont pas corrélés.
Pour ajuster un modèle avec une pente d’effet aléatoire non corrélée, utilisez || plutôt que | dans la syntaxe de lmer().
Le deuxième modèle que vous avez construit dans l’exercice précédent, model_b, a été chargé pour vous. Comparez les sorties de model_c avec les anciennes sorties de model_b.
La vignette lme4 lmer comporte une section sur les effets aléatoires non corrélés.
Cet exercice fait partie du cours
Modèles hiérarchiques et à effets mixtes en R
Exercice interactif pratique
Essayez cet exercice en complétant cet exemple de code.
# Include AverageAgeofMother as fixed-effect and LogTotalPop and State as uncorrelated random-effects
model_c <- lmer(BirthRate ~ ___,
county_births_data)
# Compare outputs of both models
summary(model_b)
summary(model_c)