[R-sig-ME] Very weird lmer results, compared to SAS proc mix

[hadley wickham]

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/