我的数据框如下:
> 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的代码吗?
谢谢。