4 votos

Regresión multinivel: pregunta sobre notación.

Tengo algunas dificultades en la comprensión de la notación de los modelos de regresión multinivel. Consideremos, por ejemplo, una variable intercepto y la variación de la pendiente del modelo con un solo nivel-I predictor. Tenemos:

Modelo de nivel I: $y_{ij}=\beta _{0j}+\beta _{1j}x_{ij}+\epsilon _{ij}$
Modelo de nivel II:
$\beta _{0j}=\gamma _{00}+u_{0j}$
$\beta _{1j}=\gamma _{10}+u_{1j}$
Modelo combinado: $y_{ij}=\gamma _{00}+\gamma _{10}x_{ij}+u_{0j}+u_{1j}x_{ij}+\epsilon _{ij}$

donde:
$\gamma _{00},\gamma _{10}$ son fijos

$u_{0j}\sim N(0,\sigma _{\beta_{0}}^{2})$
$u_{1j}\sim N(0,\sigma _{\beta_{1}}^{2})$
$\epsilon _{ij} \sim N(0,\sigma _{\epsilon}^{2})$

$\epsilon _{ij} \bot u_{0j}$
$\epsilon _{ij} \bot u_{1j}$

$\beta _{0j}\sim N(\gamma _{00},\sigma _{\beta_{0}}^{2})$
$\beta _{1j}\sim N(\gamma _{10},\sigma _{\beta_{1}}^{2})$
$Cov(\beta _{0j},\beta _{1j})= \rho \sigma _{\beta_{0}}\sigma _{\beta_{1}}$

$Var(y_{ij} |x_{ij})=\sigma _{\beta_{0}}^{2}+x_{ij}^{2}\sigma _{\beta_{1}}^{2}+2\rho \sigma _{\beta_{0}}\sigma _{\beta_{1}}x_{ij}+\sigma _{\epsilon}^{2}$
$Var(y_{ij} |x_{ij})=\sigma _{y_{i}}^{2}$

Ahora mi pregunta es la siguiente: ¿cuál de estas notaciones son correctas y por qué?
1) $y_{ij}\sim N(\beta _{0j}+\beta _{1j}x_{ij}, \, \sigma _{\epsilon}^{2})$
2) $y_{ij}\sim N(\beta _{0j}+\beta _{1j}x_{ij}, \, \sigma _{y_{i}}^{2})$
3) $y_{ij} | \beta _{0j},\beta _{1j}\sim N(\beta _{0j}+\beta _{1j}x_{ij}, \, \sigma _{\epsilon}^{2})$
4) $y_{ij}\sim N(\gamma _{00}+\gamma _{10}x_{ij},\,\sigma _{\epsilon}^{2})$
5) $y_{ij}\sim N(\gamma _{00}+\gamma _{10}x_{ij},\,\sigma _{y_{i}}^{2})$
El hecho es que yo no lo hago ahora, si en el nivel individual se debería considerar la posibilidad de $\beta _{0j}$ $\beta _{1j}$ fijo o no. Así que no sé si debo añadir que el plazo de la varianza entre la varianza o no.
Gracias.

2voto

JavaAndCSharp Puntos 1227

Esto es lo que haría yo para denotar este modelo jerárquico,

$$y_{ij} \,\vert\, \alpha_{j},\beta_{j} \sim \textsf{N}\left( \alpha_j+\beta_j x_{ij},\sigma^2 \right)$$ $$\alpha_j\sim \textsf{N}\left( \gamma,\tau^2 \right)$$ $$\beta_j \sim\textsf{N}\left( \kappa,\phi^2 \right)$$ donde $\gamma,\tau,\kappa,\phi $ se dan constantes. He cambiado la notación ya que creo que esto se ve más limpio, sin todos los subíndices y las libera de cualquier idea de la interdependencia. Si hay covarianza, entonces esto podría funcionar mejor, $$y_{ij}\,\vert\,\vec{\beta}\sim\textsf{N}\left( \vec{\beta}\mathbf{X},\sigma^2 \right)$$ $$\vec{\beta}\sim\textsf{MVN}\left( \vec{\gamma},\Sigma \right)$$

donde $\vec{\beta},\vec{\gamma}$ da $2\times1$ vectores, y $\Sigma$ es $2\times2$ matriz.

A menos que haya ajustado el modelo, no puede sostener el nivel superior fija, ya que depende de los niveles inferiores.

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