<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">Considere o código abaixo.</div><div class="gmail_default"><span style="font-family:monospace,monospace"><br></span></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><span style="font-family:monospace,monospace">library(doBy)<br><br># 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># Médias amostrais.<br>aggregate(prod ~ trt, data = da, FUN = mean)<br><br># Médias ajustadas.<br>LSmeans(m0, effect = "trt")<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)<br><br># Médias amostrais por local.<br>aggregate(prod ~ trt + loc, data = da, FUN = mean)<br><br># Médias ajustadas.<br>lapply(fits, FUN = LSmeans, effect = "trt")<br><br># Para o delineamento usado nesse exemplo, as médias ajustadas coincidem<br># com as médias amostrais.</span><br><br></div>À<div style="font-family:trebuchet ms,sans-serif;display:inline" class="gmail_default">disposição.</div><div style="font-family:trebuchet ms,sans-serif;display:inline" class="gmail_default"><br></div><div style="font-family:trebuchet ms,sans-serif;display:inline" class="gmail_default">Walmes.</div><br><br></div>