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))
#
Walmes, eu consegui fazer pelo comando lapply a anova individual dos 5 locais.
Só que agora não estou conseguindo colocasr os termos da fonte de variação em
ordem:
# ajuste do modelo de efeitos fixos para bloco dentro de repeticao
#Anova Individual
mi <- lapply(split(da, f=da$local), aov,