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?
Odpowiedzi:
Prawidłowe jest otrzymywanie spraw zamiast stawek, ponieważ przewidujesz przypadki. Jeśli chcesz uzyskać stawki, powinieneś użyć metody przewidywania w nowym zestawie danych, w którym wszystkie kolumny są równe danym, ale kolumna populacji identycznie równa 1, więc log (populaton) = 0. W takim przypadku otrzymasz liczbę przypadków jednej jednostki populacji, tj. Wskaźnik.
źródło