<div dir="ltr"><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Euriann,</div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Basicamente a automação disso pode ser feita assim:</div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div class="gmail_default"><span style="font-family:monospace,monospace"># Conjunto de ensaios em DBC em 6 locais (gerado por simulação).<br>da <- expand.grid(loc = gl(6, 1),<br>                  rpt = gl(3, 1),<br>                  blc = gl(4, 1),<br>                  trt = gl(5, 1))<br>da$prod <- rnorm(nrow(da))<br><br># Análise conjunta.<br>m0 <- aov(terms(prod ~ loc/rpt/blc + trt, keep.order = TRUE),<br>          data = da)<br>anova(m0)<br><br># Análises por local.<br>fits <- lapply(split(da, f = da$loc),<br>               FUN = aov,<br>               formula = terms(prod ~ rpt/blc + trt, keep.order = TRUE))<br>lapply(fits, anova)</span></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">À disposição.</div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Walmes.<br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div>​</div>