[R-br] simulação binomial

Rodrigo Coster rcoster em gmail.com
Quinta Julho 19 10:15:14 BRT 2012


Não seria uma bernoulli (desconsiderando o fato da bernoulli(p) ser
binomial(1,p))? Seria binomial se tu tivesse o numero de dias de chuva em
uma semana, por ex.

Se teus dados estiverem como 0 e 1, um mean(x) deve bastar pra estimar o p.

2012/7/19 Natalia Martins <nsmbarreto em gmail.com>

> Prezados,
>
> Estou simulando dados diários de chuva tentei simular utilizando as
> funções fitdistr (para obtenção das estimativas dos parâmetros)
> e  rgamma, no entanto tenho bastante zeros (dias sem chuva), com isso não
> é possível a obtenção dos parâmetros.
>
> Erro apresentado no R:
>
> Erro em optim(x = c(0, 29.9, 0, 0.7, 0, 0, 8.7, 0, 0, 26.8, 0, 19.4,  :
>   valor inicial em vmmin não é finito
>
> Posto isso resolvi transformar meus dados 0=não choveu e 1=choveu e
> simular uma binomial e ver a  partir dai se choveu eu simulo o quanto
> choveu.
> No entanto para a função fitdistr não é definido a distribuição binomial.
> Sendo assim alguém saberia como eu posso obter o valor do parâmetro p?
>
> Antecipadamente meu obrigada.
> Att,
>
> Natália da Silva Martins
> Bacharel em Estatística - Universidade Estadual de Maringá/ UEM
> Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP
> Contato: (19) 8306-4743
>
>
>
> _______________________________________________
> 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/20120719/290c61d8/attachment.html>


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