CommencerCommencer gratuitement

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

Afficher le cours

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)
Modifier et exécuter le code