[R-sig-ME] Fwd: glmmadmb help
Ben Bolker
bbolker at gmail.com
Mon Sep 30 01:58:37 CEST 2013
[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
> /
>
>
More information about the R-sig-mixed-models
mailing list