Você pode obter o quadro de anova com SQ fatiadas para solo/nc e com com a linha do resíduo do modelo completo, o seu avi1. Para isso você não pode abandonar nenhum termo. Veja o exemplo.

da <- expand.grid(bloc=gl(4,1),
                  A=gl(2,1), B=gl(3,1), C=gl(4,1))
da$y <- rnorm(nrow(da))

m0 <- aov(y~bloc+A*B*C, data=da)
anova(m0)

## Fatiar B/C.
m1 <- aov(y~bloc+A*B+B/C+A:B:C, data=da)
anova(m1)

m1$assign
names(coef(m1))[m1$assign==5]

summary(m1, split=list("B:C"=list(
                           "B1"=c(1,4,7),
                           "B2"=c(2,5,8),
                           "B3"=c(3,6,9))))

À disposição.

Walmes.