[R-sig-ME] Deviance tests and contrasts

David Duffy davidD at qimr.edu.au
Wed Jun 9 12:30:28 CEST 2010


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)



-- 
| David Duffy (MBBS PhD)                                         ,-_|\
| email: davidD at qimr.edu.au  ph: INT+61+7+3362-0217 fax: -0101  /     *
| Epidemiology Unit, Queensland Institute of Medical Research   \_,-._/
| 300 Herston Rd, Brisbane, Queensland 4029, Australia  GPG 4D0B994A v




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