我使用glmmTMB
运行了多级模型:
ml.zip <- glmmTMB(y ~ x1*x2*x3 + offset(log(exposure)) + (1 | id),
zi = ~ 1,
weights = wts,
data = df,
family = poisson)
现在,我想构建一些自举置信区间。使用lme4
开发模型时,可以使用confint
函数实现:
confint(ml.zip, method = "boot")
但是,使用glmmTMB
构建模型时,同一函数会抛出错误:
> confint(ml.zip1, method = "boot")
Error in match.arg(method) :
'arg' should be one of “wald”, “Wald”, “profile”, “uniroot”
有没有办法使用glmmTMB
构建自举置信区间?