12 Fev
2015
12 Fev
'15
15:23
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.