[R-sig-ME] Likelihood estimation by glm and glmer/lmer
bbolker at gmail.com
Tue Sep 5 14:24:01 CEST 2017
anova(), AIC(), etc. work fine across glm / glmer models and have since
1.0-0 (Aug 2013). More info on counting parameters (which are *numerator*
or *model* degrees of freedom, not *residual* df) is available at ...
On Tue, Sep 5, 2017 at 4:55 AM, Phillip Alday <phillip.alday at mpi.nl> wrote:
> For "true likelihood" ... you should search the list archives for
> discussions on REML vs ML estimation and D. Bates' comments on *the*
> likelihood. But, yes, if you you use REML=FALSE in lmer, you are
> estimating the likelihood.
> However, there's another problem with using AIC/BIC to compare mixed vs.
> non-mixed models, namely how to count parameters in mixed models.
> There's also been some discussion here of late with issues in counting
> parameters in mixed models. For Bayesian models using DIC and WAIC, this
> seems to be somewhat less of a problem because the effective number of
> parameters is estimated as part of the procedure (maybe Jarrod Hadfield
> or Paul Bürkner can comment/correct here), but there doesn't seem to be
> a clear answer for what the actual number of parameters in a model is.
> This is related to the degrees of freedom issue (
> All that said, there's an older post on this list that suggests that you
> can use the deviance (which is -2 log likelihood) to compare lm and lmer
> I think newer versions of lme4 even support this explicitly without any
> tl;dr: compare the likelihoods with REML=FALSE, but be careful with
> counting parameters and hence AIC, etc.
> On 09/05/2017 10:24 AM, Toni Hernandez-Matias wrote:
> > Dear all,
> > I would like to compare AIC values of a null model estimated with glm
> > function and AIC values of models that only have random effects fitted
> > glmer or lmer functions. I understand that they are comparable because
> > estimate true likelihood. Could you confirm me this?
> > Thank you very much in advance,
> > Antonio
> R-sig-mixed-models at r-project.org mailing list
[[alternative HTML version deleted]]
More information about the R-sig-mixed-models