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="gmail_quote">2013/1/29 marcelo costa <span dir="ltr"><<a href="mailto:marcelokosta@yahoo.com.br" target="_blank">marcelokosta@yahoo.com.br</a>></span><br>
<blockquote class="gmail_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="MsoNormal">Prezados membros,</div><div class="MsoNormal"><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="MsoNormal">
<br></div><div class="MsoNormal">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="MsoNormal"><br></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:arial,helvetica,sans-serif" class="MsoNormal">
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="MsoNormal">
<br></div>
<div class="MsoNormal">Eu estou enviando um exemplo hipotético do formato da minha
planilha:</div><div class="MsoNormal"><br></div>
<div class="MsoNormal"><span lang="EN-US">especie<-rep(c("a","b","c"),each=6)</span></div>
<div class="MsoNormal"><span lang="EN-US">sexo<-rep(c("f","m"),
each = 3, len = 18)</span></div>
<div class="MsoNormal">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="MsoNormal">data<-data.frame(especie,sexo,tamanho)</div>
<div class="MsoNormal"> </div>
<div class="MsoNormal">Um abraço a todos e obrigado pela atenção.</div><span class="HOEnZb"><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 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></blockquote></div><br><br clear="all"><div><br></div>-- <br>
Manoel Galdino<br><a href="https://sites.google.com/site/galdinomcz/" target="_blank">https://sites.google.com/site/galdinomcz/</a><br>
</div>