[R-sig-ME] Deviance tests and contrasts

David Duffy davidD at qimr.edu.au
Wed Jun 16 00:16:31 CEST 2010


On Sun, 13 Jun 2010, Adam D. I. Kramer wrote:

> d <- 
> data.frame(favcol=factor(rep(c("red","blue","yellow","purple"),each=50)),
> subj=factor(rep(1:20,each=10)), dv=rnorm(200))
>
> lmer(dv ~ favcol + (1|subj), data=d)
>
> ...in this example, with default (dummy-coded) contrasts, the "favcolpurple"
> effect tests whether purple differs from red, but I can't (easily) subtract
> the "favcolpurple" effect and recompute.
>
> ...but what I want is a more elegant way to do this for lmer models, as it
> seems to me that the "contrasts" system exists for exactly this reason.
>

Not elegant ;)

levels(d$favcol)[4] <- 'red'




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