Very weird lmer results, compared to SAS proc mix
On Sat, Mar 27, 2010 at 11:25 PM, Yong Wu <wuyong88 at gmail.com> wrote:
Sorry to bother you. I am struggling in this issue for long time. Wish somebody can help me. I first used lmer to do the following analysis. fullmodel=lmer(BMI~1+exposure+(age|ID),data, REML=FALSE) ? ? ? ? ?reducemodel=lmer(BMI~1+(age|ID),data, REML=FALSE) ? ? ? ? ?anova(full,red) The "fullmodel" has AIC of 6874 and "reducemodel" has AIC of 7106, which cause "anova" analysis giving the p-value< 2.2e-16 . This result is definitely wrong
How do you know? It would be helpful if you provided the evidence you used to judge SAS correct and R incorrect. Hadley
Assistant Professor / Dobelman Family Junior Chair Department of Statistics / Rice University http://had.co.nz/