I believe the recent discussion about AIC or p-values has missed a
crucial practical issue.
The AIC statistic reported by a default call to lmer() has NOTHING to do
with the choice of fixed effects. lmer() uses reml to define the fit of
a model. REML (residual ML or restricted ML) is a log
Hi Philip,
Thanks very much for this, i was completely unaware. I have read various
papers using lmer to calculate the AIC statistic and none have mentioned this?
I have just run through a random section of my models with this correction,
however the AIC / BIC values are the same with the
Chris Mcowen chrismco...@... writes:
Hi Philip,
Thanks very much for this, i was completely unaware. I have read various
papers using lmer to calculate the
AIC statistic and none have mentioned this?
I have just run through a random section of my models with this correction,
however
They are described as
“nearly” interchangeable because the ‘REML’ argument only applies
to calls to ‘lmer’ and the ‘nAGQ’ argument only applies to calls
to ‘glmer’
I am using lmer?
Thanks
Chris
On 5 Aug 2010, at 14:16, Manuel Morales wrote:
REML does not apply for glmer fits:
-sig-eco] AIC / BIC vs P-Values in lmer
I have just tried it with REML=FALSE and once again there is no difference in
the AIC/BIC values between the two models? I have given two examples this time
but have tried it with 10 models with no difference.
Thanks,
Chris
1
MODEL WITH REML=FALSE
on behalf of Ben Bolker
Sent: Thu 05/08/2010 3:13 PM
To: r-sig-ecology@r-project.org
Subject: [R-sig-eco] AIC / BIC vs P-Values in lmer
Thanks (forehead slap -- I knew that but it escaped me -- Manuel
Morales also pointed this out, off-list).
Isn't the difference between
(1|order/family
-project.org on behalf of Ben Bolker
Sent: Thu 05/08/2010 3:13 PM
To: r-sig-ecology@r-project.org
Subject: [R-sig-eco] AIC / BIC vs P-Values in lmer
Thanks (forehead slap -- I knew that but it escaped me -- Manuel
Morales also pointed this out, off-list).
Isn't the difference between
(1|order