[R-sig-ME] Deviance tests and contrasts

[Kay Cecil Cichini]

hello adam, david and others,

i recently was trying to figure out the same problem. i think that 
anova(l2,l0) is not appropiate because the models ar not nested.

i ended up testing anova(l2,l1) for the interaction and anova(l1,l0) for 
the main effect of sex. any thoughts on this?

yours,
kay


David Duffy schrieb:
> On Tue, 8 Jun 2010, Adam D. I. Kramer wrote:
> 
>>
>> l2 <- lmer(score ~ sex * condition + (1 + condition | subject), data=d)
>>
>> ...the main effect for sex, which has only one contrast?
> 
> Being very simple minded, I would have run something like:
> l2 <- lmer(score ~ sex * condition + (1 + condition | subject), data=d)
> l1 <- lmer(score ~ sex + condition + (1 + condition | subject), data=d)
> l0 <- lmer(score ~ condition + (1 + condition | subject), data=d,
>            subset=complete.cases(sex))
> anova(l2, l1, l0)
> 
> 
>