24 Out
2017
24 Out
'17
13:33
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.