[R] AD-MB problem in package:glmmADMB

coldeyes.Rhelp coldeyes.rhelp at gmail.com
Sat Mar 29 15:25:18 CET 2008


Hi here,
i recently notice the software built by otter company  know as AD-MB 
tools . i explore a little bit this tools under  the R Package for mixed 
logistics model. it is a very interesting tool. but i get several questions.

first question: i tried the R code below, and have the error , i can run 
the code for switching the setting to impSamp=0 , it means i can only 
run the code under Laplace approximation but importance sampling.  does  
anyone  have same problem here ?

second question: although i used save.dir=NULL option , but i can not 
see the folder which saves the intermedia output, i didn't find the 
nbmm.std file which save the BLUP. the strangest thing is i can see a 
temporary folder name as _glmm_ADMB_temp_dir_ when programing is running 
, but this fold mysteriously disappeared after the programing done . 
searched the whole computer desn't find any clue.

I posted my output below.

I really appreciate if you can help me with the questions.

source("http://www.mi.uib.no/~skaug/cash/lesaffre_dat.s 
<http://www.mi.uib.no/%7Eskaug/cash/lesaffre_dat.s>")
fit_glmmADMB_ip <- glmm.admb(y~treat,random=~1,group="subject",impSamp 
=100,data=lesaffre,family="binomial",link="logit",save.dir= NULL)

Initial statistics: 2 variables; iteration 0; function evaluation 0
Function value 1.5926934e+003; maximum gradient component mag 1.2500e+001
Var Value Gradient |Var Value Gradient |Var Value Gradient
1 0.0000 1.2500e+001 | 2 0.0000 6.2441e-001 |

- final statistics:
2 variables; iteration 9; function evaluation 13
Function value 1.2593e+003; maximum gradient component mag -5.7984e-005
Exit code = 1; converg criter 1.0000e-004
Var Value Gradient |Var Value Gradient |Var Value Gradient
1-56.9607 -1.2327e-005 | 2 -3.7159 -5.7984e-005 |
inner maxg = 9.681e-005
max separable g 9.681e-005
Newton raphson 1
f = 971.7651636 max g = 9.245773036e-010
max separable g 9.245775812e-010
Newton raphson 2
f = 971.7651636 max g = 2.220446049e-015
971.7651636
Error in glmm.admb(y ~ treat, random = ~1, group = "subject", impSamp = 
100, :
The function maximizer failed
In addition: Warning message:
In shell(paste(.path.package("glmmADMB"), "/admb/", file_name, ".exe", :
'cmd' execution failed with error code 1



More information about the R-help mailing list