Regresión de Poisson para datos de recuento - predicciones

Question

Probablemente se trate de un error elemental en mi comprensión o en mi implementación de R: Estoy tratando de utilizar un modelo de Poisson para hacer algunas predicciones. Los datos originales son datos de recuento discretos. Esperaría que las predicciones también fueran resultados discretos (por ejemplo, 39, 40, 41). En cambio, las predicciones incluyen decimales (41,2), lo que parece impar para predicciones de recuento / la distribución de Poisson. ¿Qué estoy haciendo mal?

Por ejemplo:

warpbreaks
breaksmodel<-glm(breaks~wool*tension, warpbreaks, family=poisson)
predict(breaksmodel,warpbreaks,type="response")

Answer 1

Las predicciones de un modelo de Poisson son el parámetro medio de una distribución de Poisson, que no está obligada a ser un número entero. (Piense en el viejo enigma de "la familia media tiene 2,4 hijos"). Resultados simulados de su modelo utilizando esas medias predichas serán, por supuesto, números enteros: para la combinación de predictores que da una media predicha de 41,2, podría obtener 45, 39, 42 o 40 como resultados simulados.