<div dir="ltr"><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Segue exemplos que podem ser úteis.<br><br>##-----------------------------------------------------------------------------<br>## CMR.<br><br>## Fatorial 4x5 em parcela subdividida.<br>## Calagem na parcela, adubação na subparcela.<br><br>## DIC.<br>da <- expand.grid(cal=gl(4,1), adub=gl(5,1), rept=gl(4,1))<br>da$y <- rnorm(nrow(da))<br>da$parc <- with(da, interaction(cal, rept))<br>str(da)<br><br>## Quadro de anova com 2 estratos.<br><br>m0 <- aov(y~cal*adub+Error(parc), data=da)<br>summary(m0)<br><br>## Retorna mensagem de aviso mas é equivalente.<br>m0 <- aov(y~cal*adub+Error(cal:rept), data=da)<br>summary(m0)<br><br>## DBC.<br>da <- expand.grid(cal=gl(4,1), adub=gl(5,1), bloc=gl(4,1))<br>da$y <- rnorm(nrow(da))<br>da$parc <- with(da, interaction(cal, bloc))<br>str(da)<br><br>## Quadro de anova com 2 estratos.<br><br>## Com bloco de efeito fixo.<br>m0 <- aov(y~bloc+cal*adub+Error(parc), data=da)<br>m0 <- aov(y~bloc+cal*adub+Error(bloc:cal), data=da)<br>summary(m0)<br><br>## Com bloco de efeito aleatório.<br>m0 <- aov(y~cal*adub+Error(bloc+parc), data=da)<br>m0 <- aov(y~cal*adub+Error(bloc/cal), data=da)<br>summary(m0)<br><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">À disposição.<br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Walmes.<br></div>​</div>