8 votos

¿Cuáles son los valores ajustados en un modelo de efectos aleatorios?

Considere la posibilidad de un lineal intercepto aleatorio modelo:

\begin{align} y_{ij} &= A_{i} + \varepsilon_{ij} \\ A_{i} &\sim N(0,\tau^2) \\ \varepsilon_{ij} &\sim N(0,\sigma^2) \end{align}

donde, $A_i$ $\varepsilon_{ij}$ son iid e independientes el uno del otro. Yo encajaría esta en R con algo como lmer(y ~ (1|id)), donde id es el índice del grupo ($i$ en la frase anterior). ¿Cuáles son los valores ajustados puedo obtener llamando fitted() sobre el objeto devuelto por lmer?

He oído que se dijo (tal vez el orador estaba hablando vagamente) que haciendo un efecto aleatorio en lugar de fijo no afecta a la estimación de punto, sólo su varianza. Bueno, si he usado un fijo de intercepción de arriba, yo esperaría OLS para dar el grupo de medios como el conjunto de valores, que no es lo que me sale cuando ejecuto fitted sobre el modelo de efectos aleatorios. También el comentario que he citado, no parece muy justo, porque en el caso de efectos fijos que sería adecuado un plazo para cada grupo, mientras que en el de efectos aleatorios caso de que usted sólo tiene 1 df, a la derecha, la varianza del efecto de grupo? Es ese derecho? ¿El efecto aleatorio de probabilidad de la ecuación de integrar a cabo la intersección de los valores de tan sólo el efecto de grupo de la varianza debe ser estimado? Si mi interpretación es correcta, y la intersección de los términos nunca se estima que, ¿cómo interpretar los valores ajustados en un modelo de efectos aleatorios?

6voto

Sean Hanley Puntos 2428

Su modelo se ajusta a una media de población de intercepción. También se ajusta la varianza de la distribución de la población. A partir de los datos y estos dos, al azar de las intercepciones se predijo para cada estudio de la unidad. (Véase mi respuesta a esta pregunta: ¿por Qué los valores estimados a partir de un Best Linear Unbiased Predictor (BLUP) difieren de un Best Linear Unbiased Estimator (AZUL)?) En su (relativamente fácil) modelo, los valores ajustados son la suma de la estimación de efectos fijos coeficiente y el de la predicción aleatoria de interceptar. Dicho de otra manera, los valores son los valores de predicción ($\hat y_i$) para las unidades de estudio basada en el modelo, conociendo tanto la estimación de los fijos y los de predicción de efectos aleatorios.

Aquí está una demostración rápida utilizando el código de mi vinculado respuesta:

cbind(fitted(re.mod3), 
      rep(coef(summary(re.mod3))[1]+ranef(re.mod3)[[1]][[1]], each=5))
#        [,1]     [,2]
# 1  13.19965 13.19965
# 2  13.19965 13.19965
# 3  13.19965 13.19965
# 4  13.19965 13.19965
# 5  13.19965 13.19965
# 6  16.31164 16.31164
# 7  16.31164 16.31164
# 8  16.31164 16.31164
# 9  16.31164 16.31164
# 10 16.31164 16.31164
# 11 17.47962 17.47962
# 12 17.47962 17.47962
# 13 17.47962 17.47962
# 14 17.47962 17.47962
# 15 17.47962 17.47962
# 16 15.49802 15.49802
# 17 15.49802 15.49802
# 18 15.49802 15.49802
# 19 15.49802 15.49802
# 20 15.49802 15.49802
# 21 13.82224 13.82224
# 22 13.82224 13.82224
# 23 13.82224 13.82224
# 24 13.82224 13.82224
# 25 13.82224 13.82224

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X