【发布时间】:2017-02-21 09:20:27
【问题描述】:
我正在尝试运行一个混合效应模型,其中包括三个具有交互作用的固定效应以及一个随机截距和斜率。我试图在 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 中指定我的模型或失败,其他一些包将允许我测试泊松或负二项分布。
【问题讨论】:
标签: r mixed-models poisson