cementStrength <- reshape(p5.21, idvar = "Mix.Rate", varying=list(2:5),
direction="long", v.names=c("TensileStrength"))
rownames(cementStrength) <- NULL
anova(lm(TensileStrength ~ Mix.Rate*time, data = cementStrength))
Run the code above in your browser using DataLab