我是glmer.nb和一般多级建模的新手。任何人都可以指出我在这个错误消息的正确方向？据我所知，glmer.nb应该使用MASS包从g0估算theta，然后将其用于glmer.nb模型。任何见解都会非常感激。我已经尝试更新我的所有包裹。
x1<-rbinom(750, prob=.2, size=1)
x2<-rbinom(750, prob=.15, size=1)
x3<-rbinom(750, prob=.10, size=1)
id<-seq(1,750,1)
x<-c("x1","x2","x3")
data<-data.frame(id,x1,x2,x3)
data.long <- reshape(data,idvar="id", varying=list(x),
v.names=c("x"), timevar="year",
times = c(0,1,2), direction="long")
mod <- glmer.nb(x ~ (1|id),data = data.long, na.action=na.omit,
interval = log(th) + c(-3, 3), verbose = FALSE)
Error in optimize(function(t) { : object 'th' not found