<html><body><div style="color:#000; background-color:#fff; font-family:arial, helvetica, sans-serif;font-size:12pt"><div class="yiv1939687067MsoNormal">Obrigado Manoel,</div><div class="yiv1939687067MsoNormal"><br></div>
<div class="yiv1939687067MsoNormal">Consegui resolver o problema com as funções do pacote plyr
que você indicou. Obrigado pelas observações em relação às comparações múltiplas
(vou cuidar disso).</div><div class="yiv1939687067MsoNormal"><br></div>
<div class="yiv1939687067MsoNormal">Abraço.</div><div><span><br></span></div><div> </div><div><font face="arial"><span style="font-weight:bold;">Marcelo Costa<br><br></span></font><hr style="width:100%;height:2px;"><font size="2"><span style="font-family:tahoma, new york, times, serif;font-weight:bold;"></span></font><br><div style="font-family:arial, helvetica, sans-serif;text-align:center;"><font color="#111111" size="3"><b>Laboratório de Dinâmicas Ecológicas</b></font></div><div></div><div style="font-family:arial, helvetica, sans-serif;text-align:center;"><font size="3"><a rel="nofollow" target="_blank" href="http://www.pgento.ufpr.br/">Programa de Pós-graduação em Entomologia</a> - Doutorado <br>Departamento de Zoologia <br>Universidade Federal do Paraná - UFPR <br>Jardim das Américas, Curitiba, Paraná - Brasil <br>Caixa Postal 19020, CEP 81.531-980</font></div><br><hr style="width:100%;height:2px;"><font
face="arial"><br></font></div><div><br><blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; margin-top: 5px; padding-left: 5px;"> <div style="font-family: arial, helvetica, sans-serif; font-size: 12pt;"> <div style="font-family: times new roman, new york, times, serif; font-size: 12pt;"> <div dir="ltr"> <font face="Arial" size="2"> <hr size="1"> <b><span style="font-weight:bold;">De:</span></b> Manoel Galdino <mcz.fea@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> Terça-feira, 29 de Janeiro de 2013 18:10<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="yiv720541351">Você pode usar o pacote plyr. Ele permite dividir um banco em subsets e rodar uma função para cada
subset e retornar um data.frame, ou vetor, ou lista, com o resultado da função em cada subset.<div><br></div><div><br></div>
<div>Mas você tem um problema estatístico de múltiplas comparações e um teste t não é o mais recomendado.</div><div>Basicamente, se um p-valor de 0,05 quer dizer a probabilidade de se obter dados tão ou mais extremos do que o encontrado assumindo que a hipótese nula (de nenhuma diferença) é verdadeira, então 1 em cada 20 testes vai dar significativo, mesmo com a hipótese nula sendo verdadeira.</div>
<div><br></div><div>Há várias formas de corrigir isso e se você googlar por "multiple comparisons" vai achar algo.</div><div><br></div><div>De toda forma, o jeito mais simples de fazer o que você pediu é simplesmente:</div>
<div><br></div><div>reg <- lm(tamanho ~ sexo*especie)</div><div>summary(reg)</div><div><br></div><div>Você terá uma regressão com um termo de interação entre sexo e espécie e testará o efeito de sexo para cada espécie.</div>
<div><br></div><div>No exemplo que você deu, a linha abaixo faz isso:</div><div>(summary(lm(tamanho ~ sexo*especie))</div><div><br></div><div>Nesse exemplo, sexo é significativo, bem como espécies b e c (e mrelação ao baseline de sexo feminino e espécie a). Porém, não há diferença no tamanho por sexo/espécie.</div>
<div><br></div><div>abç</div><div>Manoel<br><br><div class="yiv720541351gmail_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="yiv720541351gmail_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="yiv720541351MsoNormal">Prezados membros,</div><div class="yiv720541351MsoNormal"><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="yiv720541351MsoNormal">
<br></div><div class="yiv720541351MsoNormal">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="yiv720541351MsoNormal"><br></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:arial, helvetica, sans-serif;" class="yiv720541351MsoNormal">
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="yiv720541351MsoNormal">
<br></div>
<div class="yiv720541351MsoNormal">Eu estou enviando um exemplo hipotético do formato da minha
planilha:</div><div class="yiv720541351MsoNormal"><br></div>
<div class="yiv720541351MsoNormal"><span lang="EN-US">especie<-rep(c("a","b","c"),each=6)</span></div>
<div class="yiv720541351MsoNormal"><span lang="EN-US">sexo<-rep(c("f","m"),
each = 3, len = 18)</span></div>
<div class="yiv720541351MsoNormal">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="yiv720541351MsoNormal">data<-data.frame(especie,sexo,tamanho)</div>
<div class="yiv720541351MsoNormal"> </div>
<div class="yiv720541351MsoNormal">Um abraço a todos e obrigado pela atenção.</div><span class="yiv720541351HOEnZb"><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><br clear="all"><div><br></div>-- <br>
Manoel Galdino<br><a rel="nofollow" target="_blank" href="https://sites.google.com/site/galdinomcz/">https://sites.google.com/site/galdinomcz/</a><br>
</div>
</div><br><br> </div> </div> </blockquote></div> </div></body></html>