Tengo un recuento de la enfermedad ausencias antes y después de un accidente, y quiero averiguar si un accidente aumenta la enfermedad las ausencias de manera diferente en los distintos grupos.
Estoy tratando de formular un modelo de Poisson para esto, pero no estoy seguro de si lo estoy haciendo correctamente, o si debo de estar haciendo algo completamente diferente.
Algunos de mis temas tienen un accidente (una vez), y que he dividido el de datos en dos filas para esas personas, antes y después de la accidente. Algunos nunca se enfrentan a un accidente, y para que contar siempre como sano, y por lo tanto sólo tienen una fila de los datos. Yo tengo una variable ("estado") que indica si la fila preocupaciones tiempo antes o después del accidente.
El modelo que yo he llegado con:
fit <- glm(count~state*group+age, family="poisson", data=d)
Es este un enfoque correcto?
Para llevar esto más lejos, también me gustaría tener en cuenta años de los sujetos en el estudio antes y después de accidente.
Sería añadiendo +offset(log(person_years))
a las variables dependientes lograr esto?