零膨胀GAM置信带

Question

我正在运行一个ziplss GAM模型，并希望提取（针对95％的置信区间）响应等级的预测（即，组合计数和概率）。

我可以很容易地使用se来获得se.fit = TRUE的值来进行链接扩展，也可以使用type = "response"来获得响应预测，但是我不确定如何做到这两个。 ..任何输入将不胜感激！ 这是一个使用the数据的ziplss模型。

library(mgcv)   
library(glmmTMB)
data(Salamanders)               
zipgam <- gam(list(count ~ spp * mined + s(site, bs = "re"),
                ~ spp * mined),
           data = Salamanders, family = ziplss, method = "REML")

# response preds
preds.response <- data.frame(Predict = predict(zipgam, type = "response"))
## preds and SEs of link-scale predictions...
preds.link.se <- predict(zipgam, type = "link", se.fit = TRUE)