[R] AD-model builder under R for mix logistics model, a problem.
coldeyes.Rhelp
coldeyes.rhelp at gmail.com
Sat Mar 29 13:31:52 CET 2008
the code is list below , i can run this code for laplace approximation
but not importance sampling. the code work properly if use impSamp=0 ,
but error for any other number >0 .
could you guys help me out.
thanks so much !
source("http://www.mi.uib.no/~skaug/cash/lesaffre_dat.s")
library(glmmADMB)
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