Encuentre las estimaciones de los parámetros para cada término aleatorio en un binomial GLMM (lme4).
¿Alguien sabe cómo extraer las estimaciones de parámetros de un término aleatorio cuando se utiliza la sintaxis (1 | …) en un modelo glmer (incluyendo se, t ratio y p-value)? Solo puedo acceder a la varianza promedio y la desviación estándar con la función summary.
Algunos antecedentes: Utilicé términos aleatorios de cohorte y período (ambos factorizados), donde período = cada año de encuesta y cohorte = 8 cohortes de nacimiento. Mi modelo vacío se ve así:
glmer(pid ~ edad + edad2 + (1 | cohorte) + (1| período)
davy.ai
Sí, es posible extraer las estimaciones de parámetros de los términos aleatorios cuando se utiliza la sintaxis (1 | …) en un modelo glmer. Puede utilizar la función
ranef()
para extraer los efectos aleatorios y sus respectivos errores estándar, t-ratios y p-values. Por ejemplo, utilizando su modelo, puede extraer los efectos aleatorios para los términos de cohorte y periodo de la siguiente manera:Esto le dará un data frame con los efectos aleatorios estimados para cada nivel de los factores de cohorte y periodo, junto con sus errores estándar, t-ratios y p-values. Tenga en cuenta que estas estimaciones representan las desviaciones de los efectos aleatorios desde la intercepción general, que está incluida en los efectos fijos del modelo.
También puede utilizar la función
ranef()
para graficar los efectos aleatorios y sus intervalos asociados. Por ejemplo:Esto le dará un gráfico de los efectos aleatorios estimados para cada nivel de los factores de cohorte y periodo, junto con sus intervalos de confianza del 95%.
En general, la función
ranef()
es una herramienta útil para explorar los efectos aleatorios en un modelo glmer, y puede proporcionar información valiosa sobre la variabilidad de los datos más allá de los efectos fijos del modelo.