Comprender la varianza de los efectos aleatorios en los modelos lmer ()


16

Tengo problemas para entender el resultado de mi lmer()modelo. Es un modelo simple de una variable de resultado (Apoyo) con interceptaciones estatales / efectos aleatorios estatales variables:

mlm1 <- lmer(Support ~ (1 | State))

Los resultados de summary(mlm1)son:

Linear mixed model fit by REML 
Formula: Support ~ (1 | State) 
   AIC   BIC logLik deviance REMLdev
 12088 12107  -6041    12076   12082
Random effects:
 Groups   Name        Variance  Std.Dev.
 State    (Intercept) 0.0063695 0.079809
 Residual             1.1114756 1.054265
Number of obs: 4097, groups: State, 48

Fixed effects:
            Estimate Std. Error t value
(Intercept)  0.13218    0.02159   6.123

Supongo que la variación de las intercepciones de estado variable / efectos aleatorios es 0.0063695. Pero cuando extraigo el vector de estos efectos aleatorios de estado y calculo la varianza

var(ranef(mlm1)$State)

El resultado es: 0.001800869considerablemente menor que la varianza reportada por summary().

Por lo que yo entiendo, el modelo que he especificado se puede escribir:

yi=α0+αs+ϵi, for i={1,2,...,4097}

αsN(0,σα2), for s={1,2,...,48}

αsσα2lmer()


lmer()σα2α^syiyis

Aquí hay una pregunta muy similar, con una respuesta algo diferente
Arne Jonas Warnke

Respuestas:


11

Este es un clásico anova de una vía. Una respuesta muy breve a su pregunta es que el componente de varianza se compone de dos términos.

σ^α2=E[148s=148αs2]=148s=148α^s2+148s=148var(α^s)

Entonces, el término que calculó es el primer término en la rhs (ya que los efectos aleatorios tienen media cero). El segundo término depende de si se utiliza REML de ML y la suma de los errores estándar al cuadrado de sus efectos aleatorios.


2
¡Ok lo tengo! Entonces, la suma de los SE al cuadrado de los RE - 1/48 * sum((se.ranef(mlm1)$State)^2)- es 0.004557198. La varianza de las estimaciones puntuales de las ER (obtenidas, como anteriormente, usando var(ranef(mlm1)$State)) es 0.001800869. La suma es 0.006358067, que es la varianza informada usando summary()en el lmer()modelo anterior, de 4 o 5 dígitos como mínimo. Muchas gracias @probability
nomad545

2
Para aquellos que buscan esta respuesta y el comentario de ayuda, tenga en cuenta que nomad545 también ha utilizado el armpaquete R para la se.ranef()función.
ndoogan

1
@probabilityislogic: ¿Puede proporcionar más detalles sobre cómo se calculó esa ecuación? Específicamente, ¿cómo se logró la segunda igualdad? Además, ¿no debería haber un sombrero en el alfa después de la primera igualdad?
user1357015

1
Ynorteormetrounl(1norteα0 0,Σ) dónde Σ=yonorteσmi2+σα2ZZTes la varianza "incondicional" de Y. Si haces esto (además de usar algunas manipulaciones) obtienes la igualdad anterior. La segunda igualdad sigue porquemi(αs)=0 0 (bajo el modelo) significado vunr(αs)=mi(αs2)
probabilityislogic
Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.