predict.glmmTMB() 具有随机效应和 beta 误差族

时间:2021-04-14 04:26:07

标签: r predict random-effects glmmtmb

我在使用 predict.glmmTMB() 时遇到了一些麻烦,当我从我的模型中去掉随机效应时,我可以让它正常工作,但是一旦我有了它们,它就拒绝工作.

模型

MOD10_temp<-glmmTMB (Adj_Prop ~  AvgGorg + TOD
                + (1|Site_ID),
                family=beta_family(link = "logit"),
                data = temp)

然后我设置了预测组件:

new_data_temp = data.frame(AvgGorg = rep(seq(0,50,1), 2), 
                           TOD = c(rep(c('midday'),51), rep(c('Crepuscular'), 51)))

# then the prediction
predict(MOD10_temp, type = 'response', newdata = new_data_temp)

但是我总是收到这个错误:

Error in eval(predvars, data, env) : object 'Site_ID' not found

所以在这个模型中,Site_ID 是随机效应,但我不明白为什么 predict() 函数不能用具有随机效应的模型进行预测?使用 glmer() 进行预测时我没有这个问题,但在这种情况下我需要一个 beta 分布,因此 glmmTMB

非常感谢任何想法/帮助!

0 个答案:

没有答案
相关问题