Dear All,
I am visiting prof. at dept. of ststistics University of Pune.
I was using lme4 package  to analyze the longitudinal data. I used the well known BtheB ("Beat the Blues) data to fit random intercept and random intercept and slope models. 
The code is:
library("lme4");
 

bdf_lmer1<-lmer(bdi~bdi.pre+time+treatment+drug+length+(1|subject),method="ML",na.action=na.omit,data=bdf_long)



summary(bdf_lmer1)



bdf_lmer2<-lmer(bdi~bdi.pre+time+treatment+drug+length+(time|subject),method="ML",na.action=na.omit,data=bdf_long)



summary(bdf_lmer2)

 

anova(bdf_lmer1, bdf_lmer2)



The output I got is somewhat different from what is given in the book " A handbook of statistical analyses using R by Brian Everitt and Torsten Hothorn". The d.f are different, AIC and BIC is also different. (The output in book is also obtained by using the same package ). Is this package recently modified?

Can the difference in the degrees of be attributed to estimation or nonestimation of variance of random effects?

Please reply at your earliest convenience.



Thanks and best regards,

Dr. Sudha Purohit.

	[[alternative HTML version deleted]]


