R中的零膨胀泊松GLM

时间:2020-10-12 12:26:25

标签: r ggplot2 glm

我的数据框如下:

> head(dataset_presence)
        Date  Int_Bk    ui_avg     ui_sd dcalls sep2 presence  ui_lag02
5  2016-01-12 -9044.3  1.136538  4.502629     41    2        1        NA
6  2016-01-13 -8761.7 -1.841046 12.318525     17    0        1        NA
7  2016-01-14 -8930.4  3.887893 12.243056      0    1        1  1.136538
8  2016-01-15 -9267.8 10.980019  7.707563      0    0        0 -1.841046
9  2016-01-16 -9092.2 19.544113 15.408313      0    0        0  3.887893
10 2016-01-17 -9146.0  1.304933 13.310565      0    0        0 10.980019

我想对dcalls〜ui_avg运行零膨胀的GLM,然后使用ggplot显示响应曲线。 这是我用于非零膨胀GLM的代码(可以正常工作):

p7 <- ggplot(dataset_presence, aes(ui_avg, dcalls)) +
    geom_point() +
    geom_smooth(method="glm", method.args = list(family = poisson(link = "log"))) 

有人可以帮我提供用于零膨胀GLM的代码吗?

谢谢。

0 个答案:

没有答案