<div><span style="font-size: medium; color: #0000ff;">Tem o conjunto de dados para rodar?</span></div>
<div> </div>
<div class="notviscode" style="display: none;">p { margin: 0; }</div>
<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>-- <br /><strong>José Wellingthon dos Santos</strong> <br />Pesquisador - Métodos Quantitativos <br />Embrapa Algodão <br />Empresa Brasileira de Pesquisa Agropecuária (Embrapa) <br />Campina Grande, PB <br /><br /><span style="color: #3333ff;">jose-wellingthon.santos@embrapa.br </span><br />Telefone: +55 (83) 3182-4401 | Fax: +55 (83) 3182-4367 <br /><span style="color: #3366ff;"><a href="http://www.embrapa.br/algodao" target="_blank">www.embrapa.br/algodao</a>| <a href="https://twitter.com/embrapa_algodao" target="_blank">https://twitter.com/embrapa_algodao</a> </span><br />Confira também: <span style="color: #3366ff;"><a href="http://www.facebook.com/agrosustentavel" target="_blank">www.facebook.com/agrosustentavel</a></span> <br /><span style="font-family: arial, helvetica, sans-serif; font-size: small;"><img style="border: 0px none;" src="http://sistemas.sede.embrapa.br/nonsec/logo.png" alt="" /></span><br /><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;"><strong>De: </strong>"Jose Claudio Faria via R-br" <r-br@listas.c3sl.ufpr.br><br /><strong>Para: </strong>"Fernando Rodrigo Bortolozo" <fernandobortolozo@gmail.com>, "a lista Brasileira oficial de discussão do programa R." <r-br@listas.c3sl.ufpr.br><br /><strong>Enviadas: </strong>Terça-feira, 6 de dezembro de 2016 16:44:16<br /><strong>Assunto: </strong>Re: [R-br] Loop script<br /><br />
<div dir="ltr">
<div>Fernando,</div>
<div> </div>
Você poderia pensar em fazer uma manova e em seguida uma aov do objeto que armazena a manova.
<div> </div>
<div>Por exemplo, no código fonte da documentação da manova você pode fazer uma aov do objeto npk2.aov</div>
<div> </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> </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 />///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\</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="../../../undefined//compose?to=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> </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> </div>
<div>Aqui vai uma pequena parte das variáveis que vou analisar, tenho 7 tratamentos, 4 repetições.  </div>
<div> </div>
<div><img style="max-width: 100%; opacity: 1;" src="../content_id?msg_id=MjQxMjQ0&cid=158d596fee3d0cad2201&folder=INBOX" alt="pasted1" /></div>
<div> </div>
<div>Aqui esta o script (Agricolae Package)</div>
<div> </div>
<div>
<div><span style="font-size: xx-small;">dados <- read.csv("~/Documents/R/Data306090DAPv2.csv", header = TRUE)</span></div>
<div><span style="font-size: xx-small;">head(dados)</span></div>
<div><span style="font-size: xx-small;">model<-aov(Cas30 ~Trat, data=dados)</span></div>
<div><span style="font-size: xx-small;">out <- HSD.test(model,"Trat", group=TRUE,console=TRUE,</span><span style="font-size: x-small;"> main="Calcio")</span></div>
<div><span style="font-size: xx-small;">par(mfrow=c(1,2))</span></div>
<div><span style="font-size: xx-small;">out$means</span></div>
<div><span style="font-size: xx-small;">par (mar = c (3,3,2,0),cex=0.9)</span></div>
<div><span style="font-size: xx-small;">bar.group(out$groups,ylim=c(0, 360),density=50,border="blue", col="blue", las=1)</span></div>
<div><span style="font-size: xx-small;">bar.err(out$means,variation="SE",horiz=FALSE, ylim = c(0, 360), col=colors()[15],space=0.3,bar=TRUE, las=1)</span></div>
<div><span style="font-size: xx-small;">out<-HSD.test(model,"Trat", group=FALSE)</span></div>
<div><span style="font-size: xx-small;">means<-out$means</span></div>
<div><span style="font-size: xx-small;">df<-df.residual(model)</span></div>
<div><span style="font-size: xx-small;">MSerror<-deviance(model)/df</span></div>
<div><span style="font-size: xx-small;">with(dados,HSD.test(Cas30,Trat,df,MSerror, group=TRUE,console=TRUE,</span><span style="font-size: x-small;"> main="Cálcio"))</span></div>
</div>
<div> </div>
<div> </div>
<div>Muito obrigado desde já.</div>
</div>
<br />_______________________________________________<br /> R-br mailing list<br /> <a href="../../../undefined//compose?to=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.</blockquote>
</div>
</div>
</div>
<br />_______________________________________________<br />R-br mailing list<br />R-br@listas.c3sl.ufpr.br<br /><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" 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" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forne�a c�digo m�nimo reproduz�vel.</div>
</div>
<div>_______________________________________________<br />R-br mailing list<br />R-br@listas.c3sl.ufpr.br<br /><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" 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" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forne�a c�digo m�nimo reproduz�vel.</div>