Dự đoán poM GLM có bù

Tôi biết đây có lẽ là một câu hỏi cơ bản ... Nhưng dường như tôi không tìm thấy câu trả lời.

Tôi đang lắp GLM với một gia đình Poisson, và sau đó cố gắng xem xét các dự đoán, tuy nhiên phần bù dường như được xem xét:

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")

Tôi nhận được các trường hợp không giá ...

Tôi cũng đã thử

model_glm=glm(cases~rhs(data$year,2003)+lhs(data$year,2003)+
offset(log(population)), data=data, subset=28:36, family=poisson())

với kết quả tương tự. Tuy nhiên, khi tôi dự đoán từ GAM, sử dụng mgcv, các dự đoán sẽ cân nhắc bù (tôi nhận được tỷ lệ).

Tôi đang thiếu một cái gì đó?

Đó là chính xác bạn có được trường hợp thay vì tỷ lệ vì bạn đang dự đoán trường hợp. Nếu bạn muốn đạt được tỷ lệ, bạn nên sử dụng phương pháp dự đoán trên một tập dữ liệu mới có tất cả các cột bằng dữ liệu nhưng cột dân số bằng 1, do đó, có log (populaton) = 0. Trong trường hợp này, bạn sẽ nhận được số lượng các trường hợp của một đơn vị dân số, tức là tỷ lệ.