0

我在使用 时遇到了一些麻烦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.

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

4

0 回答 0