Hi everybody,

I am using linear mixed model package "lmer" version 0.999999-0  (R for mac version 2.15.2) for my data analysis and 
I am facing an issue of getting different  values of (AIC, BIC and logLik) 
in the "lmer" and "anova" summaries. The same issue can be found in the provided example of sleep study:

> (fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy))
Linear mixed model fit by REML 
Formula: Reaction ~ Days + (Days | Subject) 
   Data: sleepstudy 
  AIC  BIC logLik deviance REMLdev
 1756 1775 -871.8     1752    1744
> (fm2 <- lmer(Reaction ~ Days + (1|Subject) + (0+Days|Subject), sleepstudy))
Linear mixed model fit by REML 
Formula: Reaction ~ Days + (1 | Subject) + (0 + Days | Subject) 
   Data: sleepstudy 
  AIC  BIC logLik deviance REMLdev
 1754 1770 -871.8     1752    1744
> anova(fm1, fm2)
Data: sleepstudy
Models:
fm2: Reaction ~ Days + (1 | Subject) + (0 + Days | Subject)
fm1: Reaction ~ Days + (Days | Subject)
    Df  AIC    BIC  logLik  Chisq Chi Df Pr(>Chisq)
fm2  5 1762 1778.0 -876.02                         
fm1  6 1764 1783.1 -875.99 0.0609      1     0.8051
 
Is there any explanation for these differences?
Thanks in advance for your help,
Best,
A. Naceri
 
 
	[[alternative HTML version deleted]]

