data(beets)
## maybe str(beets) ; plot(beets) ...
beets$bh <- with(beets, interaction(block, harvest))
summary(aov(yield~block+sow+harvest+Error(bh), beets))
summary(aov(sugpct~block+sow+harvest+Error(bh), beets))
Run the code above in your browser using DataLab