[R-br] Anova, fatorial triplo, interacao dupla

[walmes .]

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.
​
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150212/a9bd171c/attachment.html>