[R-sig-ME] R-sig-mixed-models Digest, Vol 53, Issue 26

Nicholas Lewin-Koh nikko at hailmail.net
Tue May 24 20:53:46 CEST 2011

Comparing non-nested models is tricky. One approach is to use the LRT,
with model A as the null
and then with model B as the null. The problem is that the distribution
of the LRT is not
known in the non-nested case. Essentially you acn simulate from each
model to generate the null distribution.
A nice review of approaches is in
Fraser et al 2010. A unified approach to model selection using the
likelihood ratio test. Methods in Ecology and Evolution 2(2) 155-162


> Message: 3
> Date: Tue, 24 May 2011 02:26:29 -0700 (PDT)
> From: Iasonas Lamprianou <lamprianou at yahoo.com>
> To: r-sig-mixed-models at r-project.org
> Subject: Re: [R-sig-ME] compare non-nested logistic models
> Message-ID: <990843.69746.qm at web120609.mail.ne1.yahoo.com>
> Content-Type: text/plain
> Dear friends, I would like to compare these two non-nested  models:
> Alt1 <- glmer(contplt2 ~ 1+A+B+(1|D),family=binomial, data=partNM3)
> Alt2 <- glmer(contplt2 ~ 1+A+C+(1|D),family=binomial, data=partNM3)
> I tried Faraway(2005) but he only deals with nested models (unless I
> missed 
> something). Can anyone drop a very quick line of help (or just a link
> which I 
> can follow?)
> Besides glmer, I tried glmmPQL and MCMCglmm and they all give the same
> results. 
> I would like to thank all the authors of the packages for giving us free 
> alternatives. Does anyone know when each of the three is mostly
> appropriate(or 
> they do exactly the same job?)
> Thanks
>  Dr. Iasonas Lamprianou
> Department of Social and Political Sciences
> University of Cyprus

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