Recaler les pentes
Le dernier graphique a montré que l'évolution du taux de criminalité variait selon le comté. Cela vous indique que vous devriez inclure Year à la fois comme effet fixe et comme effet aléatoire dans votre modèle. Inclure Year de cette façon permet d'estimer une pente globale pour l'ensemble des comtés ainsi qu'une pente propre à chaque comté. La pente de l'effet fixe estime le changement des crimes majeurs pour l'ensemble des comtés du Maryland. La pente de l'effet aléatoire modélise le fait que les comtés ont des changements différents en matière de criminalité.
Cependant, l'ajustement de ce modèle produit un message d'avertissement ! Pour y remédier, faites commencer Year non pas à 2006, mais à 0. Nous vous fournissons cette nouvelle variable, Year2 (par exemple, 2006 dans Year correspond à 0 dans Year2). Parfois, lors de l'ajustement d'une régression, vous devez mettre à l'échelle ou centrer l'ordonnée à l'origine pour qu'elle commence à 0. Cela améliore la stabilité numérique du modèle.
Cette activité fait partie du cours
Modèles hiérarchiques et à effets mixtes en R
Instructions de l’exercice
- Créez un
lmer()pour prédireCrimeavecYearà la fois comme pente d'effet fixe et pente d'effet aléatoire, etCountycomme intercept d'effet aléatoire. - Créez un deuxième
lmer()pour prédireCrimeavecYear2à la fois comme pente d'effet fixe et pente d'effet aléatoire, etCountycomme intercept d'effet aléatoire.
Exercice interactif pratique
Essayez cet exercice en complétant ce code d’exemple.
# Fit the model with Year as both a fixed and random-effect
lmer(___ ~ Year + (1 + Year | ___) , data = md_crime)
# Fit the model with Year2 rather than Year
lmer(___ ~ Year2 + (1 + Year2 | ___) , data = md_crime)