錯誤glmmadmb（.....）函數最大化失敗（找不到STD文件）

Question

這工作得很好：

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 

Answer 1

嗯。任何可重複的例子......？

下面的簡單模擬的情況下，似乎工作（雖然glmmADMB與mcmc比MCMCglmm慢多 - 它實際上並沒有完呢對我來說，雖然它似乎沒有抱怨車輪繼續滾滾向前）。

對於這種簡單的情況下我懷疑glmmADMB由MCMCglmm爲主，雖然它可能是，如果你正在處理的反貝葉斯裁判有用...

nschool <- 20 
nrep <- 20 
dt1 <- expand.grid(school=LETTERS[1:nschool],rep=seq(nrep)) 
set.seed(101) 
u.school <- rnorm(nrep) 
dt1$eta <- u.school[dt1$school] 
dt1$bull <- rbinom(nrow(dt1),size=1,prob=plogis(dt1$eta)) 

library(MCMCglmm) 
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=TRUE) 

library(lme4) 
fit.glmer <- glmer(bull~(1|school),data=dt1,family=binomial) 

library(glmmADMB) 
fit.mc12 <- glmmadmb(bull~1+(1|school), data=dt1, family="binomial", 
        mcmc=TRUE, mcmc.opts=mcmcControl(mcmc=50000))