Usando aov().
http://www.leg.ufpr.br/~walmes/cursoR/cnpaf/cap17parsub-iso.R

Usando modelos mistos, lme().
http://www.leg.ufpr.br/~walmes/ensino/af722-2014-01/af722-2014-01-23parsub.html