Essa mensagem de erro sempre aparece. Eu suponho que seja algo relacionado a matriz criada com cada especificação para o termo de erro. No CMR abaixo você vê o mesmo modelo declarado de duas diferentes formas. No final tem a matriz para os fatores usados no termo de R e como elas são matrizes diferentes (em conteúdo e dimensões) eu suponho que a origem do erro esteja atrelado à isso. Aí tem que abrir a fonte para ver o que de fato tá pegando. ps <- expand.grid(BL=c("I", "II", "III", "IV"), ES=c("e1", "e2"), AD=c("a1", "a2", "a3")) ps$alt <- c(58,77,38,52,44,59,30,34,85,90,73,77, 59,68,45,55,66,93,67,64,54,75,53,48) str(ps) ## modelo para experimento de parcelas subsdividas m0 <- aov(alt~BL+AD*ES+Error(BL:AD), data=ps) ## criando o fator parcela fora da aov ps$parcela <- with(ps, interaction(BL, AD)) ## mesmo modelo m1 <- aov(alt~BL+AD*ES+Error(parcela), data=ps) ## a mesma anova summary(m0) summary(m1) a <- model.matrix(~BL:AD, ps) dim(a) b <- model.matrix(~parcela, ps) dim(b) À disposição. Walmes.