<div dir="ltr">Fernando,<br><div><br>> long_shape <- reshape2::melt(dados, id.vars = 'Trat')<br><div><div>> model <- function(df) lm(value ~ Trat, data = df)<br>> fits <- plyr::ldply(long_shape, .(variable), model)<br><br></div><div>Deste modo, "fits" é uma lista que armazena os ajustes para todas as variáveis. Minha sugestão e usar outro **ply para tratar tais ajustes e obter os testes de média e demais analíses necessárias.<br><br></div><div>Att,<br></div><div>FH<br></div><div><br><br></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">2016-12-06 13:29 GMT-06: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>