Algo assim não seria uma solução?<div><br></div><div><div>y<-split(data,factor(data[,"especie"]))</div><div>model<-function(x) x<-summary.aov(lm(x[,3]~x[,2]))</div><div>lapply(y,model)</div><div><br></div>
<div>Abraços</div><div><br><br><div class="gmail_quote">Em 29 de janeiro de 2013 18:10, Manoel Galdino <span dir="ltr"><<a href="mailto:mcz.fea@gmail.com" target="_blank">mcz.fea@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
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 class="h5">
<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><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></div></div>_______________________________________________<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" 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.<span class="HOEnZb"><font color="#888888"><br></font></span></blockquote>
</div><span class="HOEnZb"><font color="#888888"><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>
</font></span></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>
<blockquote style="margin:0 0 0 40px;border:none;padding:0px"></blockquote><div><b style="font-size:12.727272033691406px;color:rgb(34,34,34);font-family:arial,sans-serif;background-color:rgb(255,255,255)">Luciano F. Sgarbi</b><br>
<br style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:12.727272033691406px;background-color:rgb(255,255,255)"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:12.727272033691406px;background-color:rgb(255,255,255)">Mestrando em Ecologia e Evolução</span><br style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:12.727272033691406px;background-color:rgb(255,255,255)">
<span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:12.727272033691406px;background-color:rgb(255,255,255)">Laboratório de Ecologia de Insetos (sl. 222)</span><br style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:12.727272033691406px;background-color:rgb(255,255,255)">
<div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:12.727272033691406px;background-color:rgb(255,255,255)">Departamento de Ecologia<br>Instituto de Ciências Biológicas - ICB 1<br>Universidade Federal de Goiás, campus II<br>
Goiânia-GO<br>74001-970</div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:12.727272033691406px;background-color:rgb(255,255,255)">Brazil</div></div><div><br></div><div><br></div>
</div></div>