Reescalar pendientes
El último gráfico mostró que los cambios en la tasa de delincuencia variaban según el condado. Esto te indica que deberías incluir Year como efecto fijo y como efecto aleatorio en tu modelo. Incluir Year de esta forma estimará una pendiente global para todos los condados y, además, una pendiente para cada condado. La pendiente del efecto fijo estima el cambio en delitos graves en el conjunto de condados de Maryland. La pendiente del efecto aleatorio modela que los condados tienen cambios distintos en la delincuencia.
¡Pero ajustar este modelo produce un mensaje de advertencia! Para solucionarlo, cambia Year para que empiece en 0 en lugar de 2006. Te proporcionamos esta nueva variable, Year2 (por ejemplo, 2006 en Year es 0 en Year2). A veces, al ajustar una regresión, necesitas escalar o centrar el intercepto para que empiece en 0. Esto mejora la estabilidad numérica del modelo.
Este ejercicio forma parte del curso
Modelos jerárquicos y de efectos mixtos en R
Instrucciones del ejercicio
- Construye un
lmer()para predecirCrimeconYearcomo pendiente tanto de efecto fijo como de efecto aleatorio yCountycomo intercepto de efecto aleatorio. - Construye un segundo
lmer()para predecirCrimeconYear2como pendiente tanto de efecto fijo como de efecto aleatorio yCountycomo intercepto de efecto aleatorio.
Ejercicio interactivo práctico
Prueba este ejercicio y completa el código de muestra.
# 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)