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