I found that you can do the same thing with 'aov' as well. Sorry for any confusion. :)
> model.aov <- aov(Score.changes ~ Therapy + Block, data=table1) > summary(model.aov) Df Sum Sq Mean Sq F value Pr(>F) Therapy 2 260.93 130.47 15.259 0.001861 ** Block 4 438.00 109.50 12.807 0.001484 ** Residuals 8 68.40 8.55 --- Signif. codes: 0 ?**?0.001 ?*?0.01 ??0.05 ??0.1 ??1 [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.