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]]

