<div dir="ltr"><div>Fernando,</div><div><br></div>Você poderia pensar em fazer uma manova e em seguida uma aov do objeto que armazena a manova.<div><br></div><div>Por exemplo, no código fonte da documentação da manova você pode fazer uma aov do objeto npk2.aov</div><div><br></div><div><div>npk2 <- within(npk, foo <- rnorm(24))</div><div>(npk2.aov <- manova(cbind(yield, foo) ~ block + N*P*K, npk2))</div><div>summary(aov(npk2.aov))</div></div><div><br></div><div>HTH,</div><div class="gmail_extra"><div><div class="gmail_signature" data-smartmail="gmail_signature">///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\<br>Jose Claudio Faria<br>Estatistica<br>UESC/DCET/Brasil<br>joseclaudio.faria at <a href="http://gmail.com" target="_blank">gmail.com</a><br>Telefones:<br>55(73)3680.5545 - UESC<br>55(73)99966.9100 - VIVO<br>55(73)99100.7351 - TIM<br>55(73)98817.6159 - OI<br>55(73)98129.9942 - CLARO<br>///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\<br></div></div>
<br><div class="gmail_quote">2016-12-06 16:29 GMT-03:00 Fernando Rodrigo Bortolozo via R-br <span dir="ltr"><<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Eu estou fazendo anova e com posterior teste de médias de um experimento com muitas variáveis. Então eu tenho que repetir o script muitas vezes o mesmo script para cada variável. </div><div><br></div><div>Os colegas sabem me dizer se existe alguma maneira de fazer um Loop para um script, ou mesmo criar um elemento com todas variáveis para analisar de uma só vez?</div><div><br></div><div>Aqui vai uma pequena parte das variáveis que vou analisar, tenho 7 tratamentos, 4 repetições.  </div><div><br></div><div><img src="cid:158d596fee3d0cad2201" alt="pasted1" style="max-width:100%;opacity:1"><br></div><div><br></div><div>Aqui esta o script (Agricolae Package)</div><div><br></div><div><div><font size="1">dados <- read.csv("~/Documents/R/<wbr>Data306090DAPv2.csv", header = TRUE)</font></div><div><font size="1">head(dados)</font></div><div><font size="1">model<-aov(Cas30 ~Trat, data=dados)</font></div><div><font size="1">out <- HSD.test(model,"Trat", group=TRUE,console=TRUE,</font><span style="font-size:x-small"> main=<wbr>"Calcio")</span></div><div><font size="1">par(mfrow=c(1,2))</font></div><div><font size="1">out$means</font></div><div><font size="1">par (mar = c (3,3,2,0),cex=0.9)</font></div><div><font size="1">bar.group(out$groups,ylim=c(0, 360),density=50,border="blue", col="blue", las=1)</font></div><div><font size="1">bar.err(out$means,variation="<wbr>SE",horiz=FALSE, ylim = c(0, 360), col=colors()[15],space=0.3,<wbr>bar=TRUE, las=1)</font></div><div><font size="1">out<-HSD.test(model,"Trat", group=FALSE)</font></div><div><font size="1">means<-out$means</font></div><div><font size="1">df<-df.residual(model)</font></div><div><font size="1">MSerror<-deviance(model)/df</font></div><div><font size="1">with(dados,HSD.test(Cas30,<wbr>Trat,df,MSerror, group=TRUE,console=TRUE,</font><span style="font-size:x-small"> main=<wbr>"Cálcio"))</span></div></div><div><br></div><div><br></div><div>Muito obrigado desde já.<br></div></div>
<br>______________________________<wbr>_________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/<wbr>cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-<wbr>guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div></div>