I'm trying to look at multiple comparisons (across levels of my between subjects factor) for a model with one between-subject factor and one within-subject factor. I was trying to use TukeyHSD on an aovlist object, and then found out that I couldn't do that.
After some browsing, I found how to do what I wanted using
options(contrasts=c("contr.sum", "contr.poly")) lme_model=lme(dv ~ between*within, data frame, random=~1|ID, correlation=corCompSymm(form=~1|ID)) anova(lme_model, type="marginal") summary(glht(lme_model, linfct=mcp(between="Tukey")))
It seems to work (the
anova() results seem to jive with what I was getting from
aov()), but I guess the problem is that I don't really understand what
lme() is, or what it's doing. Does anyone have any documentation that I can look at?