Euriann, Basicamente a automação disso pode ser feita assim:
# Conjunto de ensaios em DBC em 6 locais (gerado por simulação). da <- expand.grid(loc = gl(6, 1), rpt = gl(3, 1), blc = gl(4, 1), trt = gl(5, 1)) da$prod <- rnorm(nrow(da)) # Análise conjunta. m0 <- aov(terms(prod ~ loc/rpt/blc + trt, keep.order = TRUE), data = da) anova(m0) # Análises por local. fits <- lapply(split(da, f = da$loc), FUN = aov, formula = terms(prod ~ rpt/blc + trt, keep.order = TRUE)) lapply(fits, anova) À disposição. Walmes.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.