<html><head><style type='text/css'>p { margin: 0; }</style></head><body><div style='font-family: arial,helvetica,sans-serif; font-size: 12pt; color: #000000'>Caro Fernado,<br>O loop é muito simples. Segue um script em que analisei 11 variáveis de uma só vez e ainda apliquei dois testes de comparações de medias. As duas primeiras colunas sao respectivamente <br>Trat e Repetiçoes, as demais são as variavies.<br>Att<br><div><span></span>-- <br><b>José Wellingthon dos Santos</b> <br>Pesquisador - Métodos Quantitativos <br>Embrapa Algodão <br>Empresa Brasileira de Pesquisa Agropecuária (Embrapa) <br>Campina Grande, PB <br><br><font color="#3333ff">jose-wellingthon.santos@embrapa.br </font><br>Telefone: +55 (83) 3182-4401 | Fax: +55 (83) 3182-4367 <br><font color="#3366ff">www.embrapa.br/algodao<u> </u>| https://twitter.com/embrapa_algodao </font><br>Confira também: <font color="#3366ff">www.facebook.com/agrosustentavel</font> <br><font size="2" face="arial,helvetica,sans-serif"><img style="border: 0px none;" src="http://sistemas.sede.embrapa.br/nonsec/logo.png"></font><br><span></span><br></div><hr id="zwchr"><div style="color:#000;font-weight:normal;font-style:normal;text-decoration:none;font-family:Helvetica,Arial,sans-serif;font-size:12pt;"><b>De: </b>"Fernando Rodrigo Bortolozo via R-br" <r-br@listas.c3sl.ufpr.br><br><b>Para: </b>r-br@listas.c3sl.ufpr.br<br><b>Enviadas: </b>Terça-feira, 6 de dezembro de 2016 16:29:24<br><b>Assunto: </b>[R-br] Loop script<br><br><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" class="" 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/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="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="SE",horiz=FALSE, ylim = c(0, 360), col=colors()[15],space=0.3,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,Trat,df,MSerror, group=TRUE,console=TRUE,</font><span style="font-size:x-small"> main="Cálcio"))</span></div></div><div><br></div><div><br></div><div>Muito obrigado desde já.<br></div></div>
<br>_______________________________________________<br>R-br mailing list<br>R-br@listas.c3sl.ufpr.br<br>https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br<br>Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo m�nimo reproduz�vel.</div><br></div></body></html>