r - R中的零膨胀泊松GLM

Question

我的数据框如下：

> 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 的代码？

谢谢你。

标签： rggplot2glm