Júlia,

Essa mensagem de erro ocorre devido a algo relacionado ao subespaço das matrizes envolvidas. Porém não causa problemas na análise (pelo menos no caso balanceado). Uma forma de evitar a mensagem de erro é escrever um novo fator, veja o CMR abaixo

da <- expand.grid(A=factor(1:4), B=factor(1:3), rept=factor(1:5))
da$y <- rnorm(da$A)

m0 <- aov(y~A+B+Error(A:B), data=da)
summary(m0) # anova com mensagem de erro

da$parcela <- with(da, factor(paste(A, rept))) # escrevendo um novo fator cujos níveis correspondem à A:rept
str(da$parcela)

m1 <- aov(y~rept+A*B+Error(parcela), data=da)
summary(m1) # mesma anova sem a mensagem de erro

À disposição.
Walmes.

==========================================================================
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573
VoIP: (3361 3600) 1053 1173
e-mail: walmes@ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================