[R] lmer BIC changes between output and anova

[Darren M. Ward]

list,

i am using lmer to fit multilevel models and trying to use anova to compare the models.  however, whenever i run the anova, the AIC, BIC and loglik are different from the original model output- as below.  can someone help me out with why this is happening?  (i'm hoping the output assocaited with the anova is right!).

thank you,

darren

> unconditional<-lmer(log50 ~ 1 + (1 | Stream:Site) + (1|Stream), data)
> summary(unconditional)
Linear mixed-effects model fit by REML 
Formula: log50 ~ 1 + (1 | Stream:Site) + (1 | Stream) 
   Data: data 
    AIC    BIC logLik MLdeviance REMLdeviance
 -138.8 -132.8  72.42     -150.4       -144.8


> nosection<-lmer(log50 ~ 1 + meanlogATS + residuallogATS + (1 | Stream:Site) + (1|Stream), data)
> summary(nosection)
Linear mixed-effects model fit by REML 
Formula: log50 ~ 1 + meanlogATS + residuallogATS + (1 | Stream:Site) +      (1 | Stream) 
   Data: data 
    AIC    BIC logLik MLdeviance REMLdeviance
 -140.8 -130.7   75.4     -168.9       -150.8


> anova(unconditional, nosection)
Data: data
Models:
unconditional: log50 ~ 1 + (1 | Stream:Site) + (1 | Stream)
nosection: log50 ~ 1 + meanlogATS + residuallogATS + (1 | Stream:Site) + 
unconditional:     (1 | Stream)
              Df      AIC      BIC   logLik  Chisq Chi Df Pr(>Chisq)    
unconditional  3 -144.370 -138.294   75.185                             
nosection      5 -158.861 -148.734   84.430 18.491      2  9.657e-05 ***