# NOT RUN {
require(lsmeans)
nutr.aov <- aov(gain ~ (group + age + race)^2, data = nutrition)
# Summarize predictions for age group 3
nutr.lsm <- lsmeans(nutr.aov, ~ race * group,
at = list(age="3"))
lsmip(nutr.lsm, race ~ group)
# Hispanics seem exceptional; but, this doesn't test out due to very sparse data
cld(nutr.lsm, by = "group")
cld(nutr.lsm, by = "race")
# }
Run the code above in your browser using DataLab