Wiem, że to prawdopodobnie podstawowe pytanie ... Ale nie wydaje mi się, aby znaleźć odpowiedź.
Dopasowuję GLM do rodziny Poisson, a następnie próbowałem przyjrzeć się prognozom, jednak wydaje się, że uwzględniono przesunięcie:
model_glm=glm(cases~rhs(data$year,2003)+lhs(data$year,2003),
offset=(log(population)), data=data, subset=28:36, family=poisson())
predict (model_glm, type="response")
Dostaję skrzynki, a nie stawki ...
Próbowałem też
model_glm=glm(cases~rhs(data$year,2003)+lhs(data$year,2003)+
offset(log(population)), data=data, subset=28:36, family=poisson())
z tymi samymi wynikami. Jednak gdy przewiduję z GAM, używając mgcv, przewidywania uwzględniają przesunięcie (otrzymuję stawki).
Czegoś mi brakuje?