[R-br] rodar função para cada nível de um fator
marcelo costa
marcelokosta em yahoo.com.br
Quarta Janeiro 30 10:17:20 BRST 2013
Olá Pavel,
Inicialmente a minha pergunta está relacionada com cada espécie separadamente (tem ou não tem dimorfismo). Vou utilizar essa informação para analisar a evolução do dimorfismo na família. Consegui resolver o meu problema com o pacote "plyr". Mas obrigado pelas sugestões, vou levar em consideração outras abordagens também.
abraço.
Marcelo Costa
Laboratório de Dinâmicas Ecológicas
Programa de Pós-graduação em Entomologia - Doutorado
Departamento de Zoologia - Universidade Federal do Paraná - UFPR
Jardim das Américas, Curitiba, Paraná - Brasil
Caixa Postal 19020, CEP 81.531-980
>________________________________
> De: Pavel Dodonov <pdodonov em gmail.com>
>Para: r-br em listas.c3sl.ufpr.br; marcelo costa <marcelokosta em yahoo.com.br>
>Enviadas: Quarta-feira, 30 de Janeiro de 2013 8:31
>Assunto: Re: [R-br] rodar função para cada nível de um fator
>
>
>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/8f3d9e59/attachment.html>
Mais detalhes sobre a lista de discussão R-br