Создание функции прогнозирования в регрессии Пуассона

Question

У меня есть набор данных, в котором мне нужно выполнить регрессионный анализ Пуассона, как количество посещений врача в течение двух недель изменяется в зависимости от возрастной группы (т.е. < 30 , между 30 и 50, и> 50), пола и болезни. Проведение секса и количество заболеваний постоянны по их средним значениям.

Вот пример моих данных:

visits gender age illness 
    1 female 19  1 
    1 female 19  1 
    1 male 19  3 
    1 male 19  1 
    1 male 19  2 
    1 female 19  5 
    1 female 19  4 
    1 female 19  3 
    1 female 19  2 
    1 male 19  1 

Однако, я не знаю, как идти об этом, я не знаю, как правильно ввести эти группы. Поскольку мне нужно узнать предсказанные темпы посещений врача в течение двух недель для разных возрастных групп.

Я знаю, как ввести исходное уравнение: GLM (посещения ~ возраст + пол + болезни, данные = DoctorVisits, семья = пуассоновы)

Но я не знаю, как бы я идти о создании функции прогнозирования ,

Answer 1

Допустим, вы хотите, чтобы предсказать, мужчина 21 г/о и болезни = 3

predict(your_glm, 
newdata = data.frame(gender = "male", age = 21, illness = 3), 
type = "response") 

Здесь вы в основном создаете фрейм данных с obseravations вы хотите прогнозы внутри функции. Если у вас есть несколько наблюдений, для которых требуется предсказание, тогда может быть более целесообразным сначала создать кадр данных, а затем подать его для прогнозирования функции, просто поменяйте "newdata=" на "data=".

type = "response" даст вам прогноз в том же формате, что и в glm, иначе это будет log ods.