Boa tarde!

Pessoal, tenho um experimento em DIC, com 6 repetições, em esquema de parcelas subdivididas, onde tenho nas parcelas um fatorial (2 variedades X 2 vinhaças) e nas subparcelas aplicou-se aluminio (4 níveis: 60, 45, 35 e 20%). 

Minha dúvida: como defino no R um modelo inicial para analisar isso?

Seria algo talvez parecido com isso:

aov(y ~ variedade*vinhaça*aluminio + Error(parcelas/subparcelas))

?