这工作正常:
fit.mc1 <-MCMCglmm(bull~1,random=~school,data=dt1,family="categorical",
prior=list(R=list(V=1, fix=1), G=list(G1=list(V=1, nu=0))), slice=T)
这样做也是如此:
fit.glmer <- glmer(bull~(1|school),data=dt1,family=binomial)
但是现在我正在尝试使用该软件包glmmadmb
,但这不起作用:
fit.mc12 <- glmmadmb(bull~1+(1|school), data=dt1, family="binomial",
mcmc=TRUE, mcmc.opts=mcmcControl(mcmc=50000))
它会产生错误:
Error in glmmadmb(bull~ 1 + (1 | school), data = dt1, family = "binomial", :
The function maximizer failed (couldn't find STD file)
In addition: Warning message:
running command '<snip>\cmd.exe <snip>\glmmadmb.exe" -maxfn 500 -maxph 5
-noinit -shess -mcmc 5000 -mcsave 5 -mcmult 1' had status 1