# NOT RUN {
data(Oslo)
table(Oslo$litho)
Oslo.mod <- lm(cbind(Cu, K, Mg, Mn, P, Zn) ~ litho, data=Oslo)
Anova(Oslo.mod)
heplot(Oslo.mod, var=c("Cu", "Mn"))
pairs(Oslo.mod)
# }
# NOT RUN {
if(require(candisc)) {
Oslo.can <- candisc(Oslo.mod)
Oslo.can
heplot(Oslo.can)
heplot3d(Oslo.can, shade=TRUE, wire=FALSE, alpha=0.5, var.col="red")
}
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab