<html><body><div style="color:#000; background-color:#fff; font-family:lucida console, sans-serif;font-size:12pt"><div><span>Uma alternativa, talvez ajude.</span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: 'lucida console', sans-serif; background-color: transparent; font-style: normal; "><span><br></span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: 'lucida console', sans-serif; background-color: transparent; font-style: normal; "><span>Veja a rotina abaixo:</span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: 'lucida console', sans-serif; background-color: transparent; font-style: normal; "><span><br></span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: 'lucida console', sans-serif; background-color: transparent; font-style: normal; "><span><div style="margin: 0.0px 0.0px 0.0px 0.0px; font: 16.0px Arial"><span style="color: #0d2099">require(</span>agricolae<span
style="color: #0d2099">)</span></div>
<div style="margin: 0.0px 0.0px 0.0px 0.0px; font: 16.0px Arial"><span style="color: #0d2099">?</span>kruskal</div>
<div style="margin: 0.0px 0.0px 0.0px 0.0px; font: 16.0px Arial; min-height: 18.0px"><br></div>
<div style="margin: 0.0px 0.0px 0.0px 0.0px; font: 13.0px Courier">comparison<span style="color: #0d2099"><-kruskal(</span>tamanho<span style="color: #0d2099">,</span>especie<span style="color: #0d2099">,</span>group<span style="color: #0d2099">=</span><span style="color: #c3892a">TRUE</span><span style="color: #0d2099">,</span> main<span style="color: #0d2099">=</span><span style="color: #af130b">"dada"</span><span style="color: #0d2099">)</span></div>
<div style="margin: 0.0px 0.0px 0.0px 0.0px; font: 13.0px Courier">comparison<span style="color: #0d2099"><-kruskal(</span>tamanho<span style="color: #0d2099">,</span>sexo<span style="color: #0d2099">,</span>group<span style="color: #0d2099">=</span><span style="color: #c3892a">TRUE</span><span style="color: #0d2099">,</span> main<span style="color: #0d2099">=</span><span style="color: #af130b">"data"</span><span style="color: #0d2099">)</span></div></span></div><div></div><div> </div><div>Edson Lira<br>Estatístico<br>Manaus-Amazonas<br></div> <div style="font-size: 12pt; font-family: 'lucida console', sans-serif; "> <div style="font-size: 12pt; font-family: 'times new roman', 'new york', times, serif; "> <div dir="ltr"> <font size="2" face="Arial"> <hr size="1"> <b><span style="font-weight:bold;">De:</span></b> Pavel Dodonov <pdodonov@gmail.com><br> <b><span style="font-weight: bold;">Para:</span></b> r-br@listas.c3sl.ufpr.br; marcelo
costa <marcelokosta@yahoo.com.br> <br> <b><span style="font-weight: bold;">Enviadas:</span></b> Quarta-feira, 30 de Janeiro de 2013 6:31<br> <b><span style="font-weight: bold;">Assunto:</span></b> Re: [R-br] rodar função para cada nível de um fator<br> </font> </div> <br><div id="yiv2055976776">Salve...<div>Uma dúvida... sua pergunta está relacionada a cada espécie separadamente, ou você quer saber se em média os machos são maiores/menores que as fêmeas independentemente da espécie? Se for a segunda opção, você poderia fazer un modelo misto (pacote nlme, função lme). Ele vai comprar o tamanho de machos e fêmeas levando em conta a variação entre as espécies.</div>
<div>Uma outra sugestão (pra reduzir o problema das comparações múltiplas) seria, ao invés de comparar espécie por espécie, dividir as espécies em grupos funcionais ou famílias, o que for mais interessante/fizer mais sentido, e rodar um modelo misto para cada grupo.</div>
<div>Só uma sugestão. ;)</div><div>Abs...</div><div>- Pavel<br clear="all"><div><div>-------</div>Pavel Dodonov<br><font size="1">Biologist, PhD student in Ecology and Natural Resources, São Carlos Federal University (UFSCar), SP, Brazil</font><div>
<font size="1"><a rel="nofollow" target="_blank" href="http://ufscar-br.academia.edu/pdodonov">http://ufscar-br.academia.edu/pdodonov</a><br></font><div><i><b><font size="1">"The highest function of ecology is understanding consequences." - Frank Herbert, Dune, 1965</font></b></i></div>
</div></div>
<br><br><div class="yiv2055976776gmail_quote">2013/1/29 marcelo costa <span dir="ltr"><<a rel="nofollow" ymailto="mailto:marcelokosta@yahoo.com.br" target="_blank" href="mailto:marcelokosta@yahoo.com.br">marcelokosta@yahoo.com.br</a>></span><br><blockquote class="yiv2055976776gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div><div style="font-size: 12pt; font-family: arial, helvetica, sans-serif; ">
<div class="yiv2055976776MsoNormal">Prezados membros,</div><div class="yiv2055976776MsoNormal"><br></div>Peço a ajuda de vocês para o seguinte problema:<span> </span>Eu tenho uma planilha
com dois fatores (“espécie” e “sexo”) e uma <span> </span>variável contínua (“tamanho”). Quero analisar se
existe diferença no tamanho de machos e fêmeas para cada espécie. Assim,
preciso fazer um teste T para cada espécie. <span> </span>O problema é que eu tenho 288 espécies (uma
planilha de 6638 linhas) e, portanto, fazer um subset para cada espécie não me parece uma solução muito boa (e provavelmente não é!). <div style="font-style: normal; font-size: 16px; background-color: transparent; font-family: arial, helvetica, sans-serif; " class="yiv2055976776MsoNormal">
<br></div><div class="yiv2055976776MsoNormal">Eu queria fazer um teste T para cada nível do fator "especie".<span> </span>Alguém tem alguma
ideia de como eu poderia fazer isso sem precisar fazer um subset para cada
espécie? <br></div><div style="font-style: normal; font-size: 16px; background-color: transparent; font-family: arial, helvetica, sans-serif; " class="yiv2055976776MsoNormal"><br></div><div style="font-style: normal; font-size: 16px; background-color: transparent; font-family: arial, helvetica, sans-serif; " class="yiv2055976776MsoNormal">
A princípio eu achei que isso seria fácil (e provavelmente seja), no entanto o meu pouco conhecimento sobre o R não está me ajudando. <br></div><div style="font-style: normal; font-size: 16px; background-color: transparent; font-family: arial, helvetica, sans-serif; " class="yiv2055976776MsoNormal">
<br></div>
<div class="yiv2055976776MsoNormal">Eu estou enviando um exemplo hipotético do formato da minha
planilha:</div><div class="yiv2055976776MsoNormal"><br></div>
<div class="yiv2055976776MsoNormal"><span lang="EN-US">especie<-rep(c("a","b","c"),each=6)</span></div>
<div class="yiv2055976776MsoNormal"><span lang="EN-US">sexo<-rep(c("f","m"),
each = 3, len = 18)</span></div>
<div class="yiv2055976776MsoNormal">tamanho<-c(1.5,1.6,1.8,1.3,1.1,1.2,2.5,2.7,2.6,2.1,2.2,2.1,3.5,3.6,3.7,3.1,3.3,3.2)</div>
<div class="yiv2055976776MsoNormal">data<-data.frame(especie,sexo,tamanho)</div>
<div class="yiv2055976776MsoNormal"> </div>
<div class="yiv2055976776MsoNormal">Um abraço a todos e obrigado pela atenção.</div><span class="yiv2055976776HOEnZb"><font color="#888888">
<div><span></span></div><div> </div><div><font face="arial"><span style="font-weight:bold;">Marcelo Costa<br></span></font><br><font face="arial"><br></font></div></font></span></div></div><br>_______________________________________________<br>
R-br mailing list<br>
<a rel="nofollow" ymailto="mailto:R-br@listas.c3sl.ufpr.br" target="_blank" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a rel="nofollow" target="_blank" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a rel="nofollow" target="_blank" href="http://www.leg.ufpr.br/r-br-guia">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><a ymailto="mailto:R-br@listas.c3sl.ufpr.br" 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" 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.<br><br> </div> </div> </div></body></html>