[R-br] rodar função para cada nível de um fator

[Pavel Dodonov]

Salve...
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.
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.
Só uma sugestão. ;)
Abs...
- Pavel
-------
Pavel Dodonov
Biologist, PhD student in Ecology and Natural Resources, São Carlos Federal
University (UFSCar), SP, Brazil
http://ufscar-br.academia.edu/pdodonov
*"The highest function of ecology is understanding consequences." - Frank
Herbert, Dune, 1965*


2013/1/29 marcelo costa <marcelokosta em yahoo.com.br>

>  Prezados membros,
>
> Peço a ajuda de vocês para o seguinte problema:  Eu tenho uma planilha
> com dois fatores (“espécie” e “sexo”) e uma  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.  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 é!).
>
> Eu queria fazer um teste T para cada nível do fator "especie". Alguém tem
> alguma ideia de como eu poderia fazer isso sem precisar fazer um subset
> para cada espécie?
>
> 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.
>
> Eu estou enviando um exemplo hipotético do formato da minha planilha:
>
> especie<-rep(c("a","b","c"),each=6)
> sexo<-rep(c("f","m"), each = 3, len = 18)
>
> 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)
> data<-data.frame(especie,sexo,tamanho)
>
> Um abraço a todos e obrigado pela atenção.
>
> Marcelo Costa
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130130/70126b91/attachment.html>