[R-br] Analise dados proporção e contagem

Walmes Zeviani walmeszeviani em gmail.com
Sexta Dezembro 14 18:32:23 BRST 2012


Razão sexual é uma variável dicotômica (M=0, F=1), então você pode usar a
dsitribuição binomial e verificar diferenças com relação ao parâmetro p.
Use a função glm() com família binomial, o n ou size é a prole de cada mãe
e x o número de fêmeas (ou machos, como quiser). Já o tamanho da prole,
apesar de ser uma contagem, não atende os pressupostos estabeleciodos pela
distribuição Poisson que é a independência entre eventos. Facilmente se
verifica isso porque na barriga da mãe existe uma limitação de espaço, a
mãe tem uma limitação de recursos biológicos para produção dos óvulos que
são fecundados e gerados por ela, etc, pessoas experientes na biologia da
reprodução podem apontar outras razões. Tenho observado o mesmo com dados
vegetais que no caso são de números de frutos por planta, sementes por
fruto, folhas por ramo, etc. O fato é que se for usar distribuição Poisson
você possivelmente observar subdispersão, então é mais adequado que use uma
quasipoisson na glm(). Se você não é familiar com modelos lineares
generalizados procure fontes de informação. Sobre código, algo útil
disponível em

http://www.leg.ufpr.br/~walmes/cursoR/cnpaf/cap23propor-iso.R
http://www.leg.ufpr.br/~walmes/cursoR/cnpaf/cap24cont-iso.R

À disposição.
Walmes.

==========================================================================
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573
VoIP: (3361 3600) 1053 1173
e-mail: walmes em ufpr.br
skype: walmeszeviani
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121214/1fd58dd6/attachment.html>


Mais detalhes sobre a lista de discussão R-br