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.
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
