Natália, porque você não atribui valores muito próximos de zero? Tipo
0,0000000001 para substituir os zeros. Acho que poderia te ajudar.
Edson Lira
Estatístico
Manaus-Amazonas
________________________________
De: Rodrigo Coster <[email protected]>
Para: [email protected]
Enviadas: Quinta-feira, 19 de Julho de 2012 9:15
Assunto: Re: [R-br] simulação binomial
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 <[email protected]>
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
>[email protected]
>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.
>
_______________________________________________
R-br mailing list
[email protected]
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.
_______________________________________________
R-br mailing list
[email protected]
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.
