powered by
data(Plastic)
str(Plastic) plastic.mod <- lm(cbind(tear, gloss, opacity) ~ rate*additive, data=Plastic) Anova(plastic.mod) pairs(plastic.mod)
Run the code above in your browser using DataLab