[R-sig-ME] glmer and start values help

Adam D. I. Kramer adik-rhelp at ilovebacon.org
Thu Sep 8 19:36:24 CEST 2011


Hi,

I'm trying to muck about a bit with the start values for a glmer model.
Here's what I've set:

> ST
[[1]]
             (Intercept)
(Intercept)      -0.779

> fef
[1] 0.04519647 0.01011628
> has_pos <- glmer(has_fpos_d1 ~ has_upos_d1 + (1 | userid), data=x,
+ family='binomial', verbose=TRUE, control=list(msVerbose=TRUE,
+ maxIter=500, maxFN=5000), start=list(ST=ST, fixef=fef))
   0:     10415953.: 0.608675 0.0451965 0.0101163
   1:     10305257.:  0.00000 -0.139832 -0.361935

...as you can see, the two fixed effects startvalues are as I said they
should be at step 0, but the random effect start value is not. (Or at least,
it seems really weird that the fixef params would be spot-on the same if the
report from step 0 is actually the report of the values after the
iteration.)

The help page is pretty clear on the topic:

   If the list contains components named 'fixef' and/or 'ST', these are used
   as the starting values for those slots.

Any thoughts?

(has_fpos_d1 and has_upos_d1 are both logical.)

R is [1] "R version 2.13.1 (2011-07-08)"
lme4 is     note = {R package version 0.999375-41}

--Adam




More information about the R-sig-mixed-models mailing list