op <- options(contrasts = c("contr.helmert", "contr.poly"))
# AOV
( npk.aov <- aov(yield ~ block + N*P*K, npk) )
summary(npk.aov)
pprint(summary(npk.aov),
format = "summary.aov")
# MANOVA
npk2 <- within(npk, foo <- rnorm(24))
( npk2.aov <- manova(cbind(yield, foo) ~ block + N*P*K, npk2) )
summary.aov(npk2.aov)
pprint(summary.aov(npk2.aov),
format = "summary.aovlist")
options(op) # reset to previous
Run the code above in your browser using DataLab