<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><br><br>dados<-read.table("e:/adalgisa/magna.txt", header=TRUE)<br>dados<br>TRAT<-as.factor(dados[,1])<br>REP<-as.factor(dados[,2])<br><br>require(laercio)<br><br>colin<-3<br>for(i in colin:ncol(dados))<br>{<br> modelo<-aov(dados[,i]~TRAT+REP)<br>print(summary(modelo))<br><br>LTukey(modelo,"TRAT",conf.level=0.95)<br><br>LScottKnott(modelo,"TRAT",conf.level=0.95)<br>}<br><br>Att<br><br><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>"Jose Claudio Faria via R-br" <r-br@listas.c3sl.ufpr.br><br><b>Para: </b>"Fernando Rodrigo Bortolozo" <fernandobortolozo@gmail.com>, "a lista Brasileira oficial de discussão do programa R." <r-br@listas.c3sl.ufpr.br><br><b>Enviadas: </b>Terça-feira, 6 de dezembro de 2016 16:44:16<br><b>Assunto: </b>Re: [R-br] Loop script<br><br><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">///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\<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/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>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">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/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-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><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>