[R-sig-ME] anova on lmer object

Doug Adams fog0 at gmx.com
Wed Feb 10 06:06:27 CET 2010


On Tue, Feb 9, 2010 at 4:32 PM, Doug Adams <fog0 at gmx.com> wrote:
> Thanks to both of you!
>

If I only pass one model into the anova function, does it assume the
nested model being tested against is the one without fixed effects?

And are the denominator degrees of freedom the ones listed in the output table:
         Df Sum Sq Mean Sq F value
division  2 103.59  51.793  1.0155

If it's 2, am I correct in thinking that it's because my fixed effect
has 3 states, the first of which would be considered the baseline
against which the others are offset -- so it would exist in the model
without fixed effects?  (So the baseline, and then the other 2
additional states for fixed effects...)

One more question:  If I want to explicitly spell out the simpler
model that doesn't have fixed effects so I can test for them, would
this be correct:
Complete Model:        lmer(outcome ~ FixedEffect + (1|RandomGrouping), data...)
Simplified Model:        lmer(outcome ~ 1 + (1|RandomGrouping))

so that "1" is like a placeholder for the first level?

Goodness, thanks everyone for your patience with my endless questions.




More information about the R-sig-mixed-models mailing list