[R-sig-ME] Log-likelihood and mixed models in glmer

Ben Bolker bbolker at gmail.com
Mon Jan 16 18:09:53 CET 2012


Andrew Miles <rstuff.miles at ...> writes:

>  Can someone point me to a reference that will explain why, when
> using mixed models (glmer and lmer) adding explanatory variables
> decreases the log likelihood?  This makes no sense to me, as adding
> explanatory power should make the model fit the data worse.  I've
                        ^^^^^
                        never?

> attached the data I am using, which contains no missing values, and
> here are the models I am running, and the results:

  The attachment didn't make it through to the mailing list.
Could you post it somewhere (or send it to me)?

> #note, models do not fully converge, but examination of estimates
>   using verbose=T suggests they are resonable

    The fact that they didn't converge (combined with your observation)
seems like a giant warning message to me ...  hard to say more without
seeing the data (see above)

> mod.null = glmer(res.lifesat.last5 ~ 1 + (1|hhidpn) + 
> (1|hhid), data=data.nomiss, family=binomial(link="probit"))
> mod1 = glmer(res.lifesat.last5 ~ networth2.gmc + (1|hhidpn) + 
>  (1|hhid), data=data.nomiss, family=binomial(link="probit"))
> mod2 = glmer(res.lifesat.last5 ~ networth2.gmc + married + (1|hhidpn) +
>  (1|hhid), data=data.nomiss, family=binomial(link="probit"))
> mod3 = glmer(res.lifesat.last5 ~ networth2.gmc + married + depscore + 
> selfhealth + (1|hhidpn) +
> (1|hhid), data=data.nomiss, family=binomial(link="probit"))

  I have to add some text so the Gmane portal will be happy,
so let me just add that

  mod1 <- update(mod.null, . ~ . + networth2.gmc)
  mod2 <- update(mod1, . ~ . + married)
  mod3 <- update(mod2, . ~ . + depscore + selfhealth)

  would be a little bit clearer.

  Have you tried centering any continuous predictors?
 
> #note that mod2 and mod3 have lower log-likelihoods than mod1, 
>  and mod3 has a lower LL than the null model
> anova(mod.null, mod1, mod2, mod3)

   Do you get the same results from just using logLik() ?  Perhaps
anova() is scrambling things up?




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