[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