我正在尝试运行一个混合效果模型，其中包括三个具有交互作用的固定效果以及随机截距和斜率 . 我想在glmmadmb中指定的模型是：
> fit_zipoiss_ambig<-glmmadmb(AmbigCount~Posn.c*mood.c*Valence.c + offset(InputAmbig) + (1+Valence.c|mood.c/Chain), data = Data, zeroInflation = TRUE, family="poisson")
首先我收到此错误消息：
Error in Droplevels(eval(parse(text = x), data)) :
all grouping variables in random effects must be factors
所以我用（例如） fPosn.c=as.factor(Data$Posn.c) 将我的所有预测变换为因子 . 然后我跑了这个模型：
> fit_zipoiss_ambig<-glmmadmb(AmbigCount~fPosn.c*fmood.c*fValence.c + offset(InputAmbig) + (1+fValence.c|fmood.c/Chain), data = Data, zeroInflation = TRUE, family="poisson")
然后我收到了这个错误：
Error in glmmadmb(AmbigCount ~ fPosn.c * fmood.c * fValence.c + offset(InputAmbig) + :
The function maximizer failed (couldn't find STD file) Troubleshooting steps include (1) run with 'save.dir' set and inspect output files; (2) change run parameters: see '?admbControl'
In addition: Warning message:
running command 'C:\Windows\system32\cmd.exe /c "C:/Program Files/R/R-3.2.2/library/glmmADMB/bin/windows64/glmmadmb.exe" -maxfn 500 -maxph 5 -noinit -shess' had status 1
我尝试按照故障排除建议进行操作，以便在模型结束时包含 , admb.opts=admbControl(shess=FALSE,noinit=FALSE)) . 现在我收到此错误：
Error in glmmadmb(AmbigCount ~ fPosn.c * fmood.c * fValence.c + offset(InputAmbig) + :
rank of X = 106 < ncol(X) = 107
我不知道这个错误意味着什么 . 我希望有人可以帮我弄清楚如何在glmmadmb中指定我的模型或者失败，其他一些包将允许我测试泊松或负二项分布 .
1 回答
没有能够自己运行它，对我跳出的是：
就您的第一条错误消息而言，它表示嵌套随机效果公式中的变量需要是因子 .
然后，在您的代码中：fPosn.c = as.factor（Data $ Posn.c）您没有在数据框中创建“fPosn.c” . 为此，您需要运行：Data $ fPosn.c = as.factor（Data $ Posn.c）