[cc'ing to r-sig-mixed-models] I had some trouble reconstructing your results, as I had difficulty reading the file, but I did eventually succeed. I tried it with both lme4::glmer and glmmADMB::glmmadmb. glmer finished after 100 seconds (note that you also have values of "" for GENDER_; unless you know you have claims for which an indeterminate GENDER_ is acceptable, you should double-check your data cleaning process!) glmmadmb ran for 94 minutes and then quit. It *is* quite possible to do such large problems with AD Model Builder (the foundation of glmmADMB), but at present glmmadmb is built for flexibility rather than performance -- there are some tricks for speeding up random effects that haven't been applied. See e.g. the 'wildflowers' or 'owls' projects at https://groups.nceas.ucsb.edu/non-linear-modeling/projects . ===================== library(data.table) y <- as.data.frame(fread("Article.txt")) library(lme4) t1 <- system.time(tt <- glmer(Claims~GENDER_+(1|T), data=y, family="poisson", verbose=1)) ## t1: 100 seconds elapsed library(glmmADMB) y$T <- factor(y$T) t2 <- system.time(tt2 <- glmmadmb(Claims~GENDER_+(1|T), data=y, family="poisson", extra.args="-ndi 500000"))
On 13-09-29 11:50 AM, Younes Mouatassim wrote:
Hi,
I am trying to fit random effects to the number of accident using
glmmadmb. However i kieep getting two error messages:
1) tt = glmmadmb(Claims~GENDER_+MARITAL_STATUS,random=~T,data=x,family="poisson") Error in `colnames<-`(`*tmp*`, value = character(0)) : attempt to set 'colnames' on an object with less than two dimensions
2)
tt =
glmmadmb(Claims~GENDER_+MARITAL_STATUS+(1|T),data=x,zeroInflation=TRUE,family="nbinom") Error in glmmadmb(Claims ~ GENDER_ + MARITAL_STATUS + (1 | T), data = x, : The function maximizer failed (couldn't find STD file) Troubleshooting steps include (1) run with 'save.dir' set and inspect output files; (2) change run parameters: see '?admbControl' In addition: Warning message: running command 'C:\WINDOWS\system32\cmd.exe /c "C:/Program Files/R/R-3.0.1/library/glmmADMB/bin/windows32/glmmadmb.exe" -maxfn 500 -maxph 5 -noinit -shess' had status 1 I attach data.
If any one understand what s gone wrong, please give a hand. Youn?s
Cordialement *Younes MOUATASSIM* *Actuariat pricing* *Zurich Insurance * /T?l: 0//614 22 93 30 / -- Bonjour, Cordialement *Younes MOUATASSIM* *Actuariat pricing* *Zurich Insurance * /T?l: 0//614 22 93 30 /