# NOT RUN {
data(AddHealth)
# fit mlm
AH.mod <- lm(cbind(depression, anxiety) ~ grade, data=AddHealth)
Anova(AH.mod)
summary(Anova(AH.mod))
heplot(AH.mod, hypotheses="grade.L", fill=c(TRUE, FALSE))
# }
Run the code above in your browser using DataLab