Завершенные нулями доверительные интервалы GAM - PullRequest
Новая
       31

Завершенные нулями доверительные интервалы GAM

0 голосов
/ 28 октября 2018

Я использую модель ziplss GAM и хотел бы извлечь прогнозы (вместе с 95% доверительными интервалами) для шкалы ответов (т. Е. Комбинированные подсчеты и вероятности).

Я могу легко получить значения se, используя se.fit = TRUE для масштаба ссылки, и я могу получить предсказания ответа, используя type = "response", но я не уверен, как сделать оба ... любойвход будет оценен!Вот модель ziplss, использующая данные Salamanders. 

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)
...