Articulo de referencia

Modelo lineal generalizado jerárquico

En estadística , los modelos lineales generalizados jerárquicos extienden los modelos lineales generalizados al relajar el supuesto de que los componentes de error son independi...

En estadística , los modelos lineales generalizados jerárquicos extienden los modelos lineales generalizados al relajar el supuesto de que los componentes de error son independientes . [ 1 ] Esto permite construir modelos en situaciones donde se necesita más de un término de error y también permite dependencias entre términos de error. [ 2 ] Los componentes de error pueden estar correlacionados y no necesariamente seguir una distribución normal . Cuando hay diferentes clústeres, es decir, grupos de observaciones, las observaciones en el mismo clúster están correlacionadas. De hecho, están correlacionadas positivamente porque las observaciones en el mismo clúster comparten algunas características comunes. En esta situación, usar modelos lineales generalizados e ignorar las correlaciones puede causar problemas. [ 3 ]

Descripción general y modelo

Modelo

En un modelo jerárquico, las observaciones se agrupan en clústeres, y la distribución de una observación está determinada no solo por la estructura común entre todos los clústeres, sino también por la estructura específica del clúster al que pertenece esta observación. Por lo tanto, se introduce en el modelo un componente de efecto aleatorio, diferente para cada clúster.y{\displaystyle y}ser la respuesta,{\displaystyle u}sea ​​el efecto aleatorio,gramo{\displaystyle g}ser la función de enlace,η=incógnitaβ{\displaystyle \eta =X\beta }, yv=v(){\displaystyle v=v(u)}es alguna función estrictamente monótona de{\displaystyle u}. En un modelo lineal generalizado jerárquico, la suposición sobrey|{\displaystyle y|u}y{\displaystyle u}deben hacerse: [ 2 ]y F(θ,ϕ){\displaystyle y\mid u\sim \ f(\theta,\,\phi)}y F(α).{\displaystyle u\sim \ f_{u}(\alpha).}

El predictor lineal tiene la siguiente forma:

gramo(mi(y))=gramo(μ)=η=incógnitaβ+v{\displaystyle g(E(y))=g(\mu )=\eta =X\beta +v\,}

dóndegramo{\displaystyle g}es la función de enlace,μ=mi(y){\displaystyle \mu =E(y)}, η=incógnitaβ+v{\displaystyle \eta =X\beta +v}, yv=v(){\displaystyle v=v(u)}es una función monótona de{\displaystyle u}En este modelo lineal generalizado jerárquico, el efecto fijo se describe medianteβ{\displaystyle \beta }, que es el mismo para todas las observaciones. El componente aleatorio{\displaystyle u}no se observa y varía aleatoriamente entre los grupos. Por lo tanto,v{\displaystyle v}toma el mismo valor para las observaciones en el mismo clúster y valores diferentes para las observaciones en clústeres diferentes. [ 3 ]

Identificabilidad

La identificabilidad es un concepto en estadística . Para realizar la inferencia de parámetros, es necesario asegurarse de que se cumpla la propiedad de identificabilidad. [ 4 ] En el modelo mencionado anteriormente, la ubicación de v no es identificable, ya que

incógnitaβ+v=(incógnitaβ+a)+(va){\displaystyle X\beta +v=(X\beta +a)+(va)\,}

para constantea{\displaystyle a}. [ 2 ] Para que el modelo sea identificable, necesitamos imponer restricciones a los parámetros. La restricción se suele imponer a los efectos aleatorios, comomi(v)=0{\displaystyle E(v)=0}. [ 2 ]

Al asumir diferentes distribuciones dey{\displaystyle y\mid u} y{\displaystyle u}y utilizando diferentes funciones degramo{\displaystyle g}y 'v{\displaystyle v}, podremos obtener diferentes modelos. Además, el modelo lineal generalizado mixto (GLMM) es un caso especial del modelo lineal generalizado jerárquico. En los modelos lineales generalizados jerárquicos, las distribuciones del efecto aleatorio{\displaystyle u}no necesariamente siguen una distribución normal . Si la distribución de{\displaystyle u}es normal y la función de enlace dev{\displaystyle v}Si la función identidad es , entonces el modelo lineal generalizado jerárquico es lo mismo que GLMM. [ 2 ]

Distribuciones dey{\displaystyle y\mid u}y{\displaystyle u}También se puede elegir que sean conjugados, ya que se mantienen buenas propiedades y es más fácil para el cálculo y la interpretación. [ 2 ] Por ejemplo, si la distribución dey{\displaystyle y\mid u}es Poisson con cierta media, la distribución de{\displaystyle u}es Gamma y se utiliza el enlace logarítmico canónico, entonces llamamos al modelo modelos lineales generalizados jerárquicos conjugados de Poisson. Siy{\displaystyle y\mid u}sigue una distribución binomial con cierta media,{\displaystyle u}Si tiene la distribución beta conjugada y se utiliza la función de enlace logit canónica, entonces llamamos al modelo modelo beta conjugado. Además, el modelo lineal mixto es el modelo lineal generalizado jerárquico conjugado normal. [ 2 ]

Un resumen de los modelos comúnmente utilizados son: [ 5 ]

Ajuste de los modelos lineales generalizados jerárquicos

Los modelos lineales generalizados jerárquicos se utilizan cuando las observaciones provienen de diferentes grupos. Existen dos tipos de estimadores: estimadores de efectos fijos y estimadores de efectos aleatorios, que corresponden a parámetros en  :η=incógnitaβ{\displaystyle \eta =\mathbf {x} {\boldsymbol {\beta }}} y env(){\displaystyle \mathbf {v(u)} } , respectivamente. Hay diferentes maneras de obtener estimaciones de parámetros para un modelo lineal generalizado jerárquico. Si solo interesan los estimadores de efectos fijos, se puede utilizar el modelo promediado en la población. Si la inferencia se centra en los individuos, habrá que predecir los efectos aleatorios. [ 3 ] Hay diferentes técnicas para ajustar un modelo lineal generalizado jerárquico.

Ejemplos y aplicaciones

Los modelos lineales generalizados jerárquicos se han utilizado para resolver diferentes problemas de la vida real.

Ingeniería

Por ejemplo, este método se utilizó para analizar la fabricación de semiconductores, ya que los procesos interrelacionados forman una jerarquía compleja. [ 6 ] La fabricación de semiconductores es un proceso complejo que requiere diferentes procesos interrelacionados. [ 7 ] El modelo lineal generalizado jerárquico, que requiere datos agrupados, es capaz de abordar procesos complejos. Los ingenieros pueden utilizar este modelo para descubrir y analizar subprocesos importantes y, al mismo tiempo, evaluar la influencia de estos subprocesos en el rendimiento final. [ 6 ]

Negocio

Los problemas de investigación de mercado también pueden analizarse utilizando modelos lineales generalizados jerárquicos. Los investigadores aplicaron el modelo a consumidores dentro de países para resolver problemas en la estructura de datos anidados en la investigación de marketing internacional. [ 8 ]

Referencias

  1. Modelos lineales generalizados . Chapman and Hall/CRC. 1989. ISBN 0-412-31760-5.
  2. 1 2 3 4 5 6 7 Y. Lee; JA Nelder (1996). "Modelos lineales generalizados jerárquicos". Journal of the Royal Statistical Society, Serie B. 58 ( 4): 619– 678. doi : 10.1111/j.2517-6161.1996.tb02105.x . JSTOR 2346105 . 
  3. 1 2 3 Agresti, Alan (2002). Análisis de datos categóricos . Hoboken, Nueva Jersey: John Wiley & Sons, Inc. ISBN 0-471-36093-7.
  4. Allman, Elizabeth S. ; Matias, Catherine; Rhodes, John A. (2009). "Identificabilidad de parámetros en modelos de estructura latente con muchas variables observadas". The Annals of Statistics . 37, No. 6A (6A): 3099– 3132. arXiv : 0809.5032 . Bibcode : 2008arXiv0809.5032A . doi : 10.1214/09-AOS689 . S2CID 16738108 . 
  5. Lars Rönnegård; Xia Shen; Moudud Alam (dic. 2010). "hglm: Un paquete para ajustar modelos lineales generalizados jerárquicos". The R Journal . 2/2 .
  6. 1 2 Naveen Kumar; Christina Mastrangelo; Doug Montgomery (2011). "Modelado jerárquico mediante modelos lineales generalizados". Quality and Reliability Engineering International .
  7. Chung Kwan Shin; Sang Chan Park (2000). "Un enfoque de aprendizaje automático para la gestión del rendimiento en la fabricación de semiconductores". International Journal of Production Research . 38 (17): 4261– 4271. doi : 10.1080/00207540050205073 . S2CID 111295634 . 
  8. Burcu Tasoluk; Cornelia Dröge; Roger J. Calantone (2011). "Interpretación de interrelaciones a través de múltiples niveles en modelos HGLM: una aplicación en la investigación de marketing internacional" . International Marketing Review . 28 (1): 34– 56. doi : 10.1108/02651331111107099 .