1 votos

Regresión de Poisson con desplazamiento frente a regresión de Poisson con peso

Me gustaría conocer la expresión lineal del peso y el desplazamiento en términos de regresión de poisson en glm.

por ejemplo, para la compensación

   glm( y ~ x + offset(of), data, family=poisson(link="log"))

el modelo glm anterior puede expresarse como

  log(y) = constant + beta*x + log(of).

Mi pregunta es cuál es la expresión equivalente para el siguiente modelo glm:

glm(y ~ x, weight=weight, data, family=poisson(link="log")).

Según tengo entendido, el desplazamiento y el peso son diferentes y afectan a las estimaciones de forma distinta.

4voto

Gordon Smyth Puntos 170

Creo que su pregunta puede estar motivada por el hecho de que los glms log-lineales de Poisson pueden utilizarse para modelar tanto recuentos como frecuencias. Si count es un vector de recuentos y of es un vector de desplazamientos positivos, entonces los dos ajustes siguientes son equivalentes:

fit.count <- glm(count ~ x + offset(of), family=poisson)

y

frequency <- counts / of
fit.freq <- glm(frequency ~ x, family=poisson, weights=of)

Normalmente aquí of es el tiempo total de exposición. En la primera versión, los recuentos se modelan de la forma habitual de glm, y el hecho de que los recuentos esperados sean proporcionales al tiempo de exposición se codifica en el desplazamiento.

En la segunda versión, calculamos la frecuencia de recuentos por unidad de tiempo de exposición. Esto elimina explícitamente la dependencia de los recuentos del tiempo de exposición, pero cambia la varianza. El argumento de las ponderaciones debe entonces reflejar el hecho de que las frecuencias basadas en los tiempos de exposición más largos son más precisas.

Aunque las dos variables de respuesta están en escalas diferentes, los dos ajustes glm son completamente equivalentes en todos los aspectos importantes. Darán las mismas estimaciones de coeficientes, las mismas desviaciones y los mismos valores p.

La equivalencia de los dos modelos ha sido señalada anteriormente por Alan Chalk en este foro: ¿Pueden las ponderaciones y la compensación conducir a resultados similares en la regresión de Poisson?

Este mismo truco se puede utilizar para cualquier glm con un offset y un log-link. También puede utilizarse para glms binomiales en los que la variable de respuesta puede ser el recuento binomial o la proporción de éxitos.

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