<html><body><div style="color:#000; background-color:#fff; font-family:lucida console, sans-serif;font-size:12pt"><div><span>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.<br></span></div><div> </div><div>Edson Lira<br>Estatístico<br>Manaus-Amazonas<br></div> <div style="font-family: lucida console, sans-serif; font-size: 12pt;"> <div style="font-family: times new roman, new york, times, serif; font-size: 12pt;"> <div dir="ltr"> <font size="2" face="Arial"> <hr size="1"> <b><span style="font-weight:bold;">De:</span></b> Rodrigo Coster <rcoster@gmail.com><br> <b><span style="font-weight: bold;">Para:</span></b> r-br@listas.c3sl.ufpr.br <br> <b><span style="font-weight: bold;">Enviadas:</span></b> Quinta-feira, 19 de Julho de 2012 9:15<br> <b><span style="font-weight: bold;">Assunto:</span></b> Re: [R-br] simulação binomial<br> </font> </div>
<br><div id="yiv976892555">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.<div><br></div><div>Se teus dados estiverem como 0 e 1, um mean(x) deve bastar pra estimar o p.<br>
<br><div class="yiv976892555gmail_quote">2012/7/19 Natalia Martins <span dir="ltr"><<a rel="nofollow" ymailto="mailto:nsmbarreto@gmail.com" target="_blank" href="mailto:nsmbarreto@gmail.com">nsmbarreto@gmail.com</a>></span><br><blockquote class="yiv976892555gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Prezados, <div><br></div><div>Estou simulando dados diários de chuva tentei simular utilizando as funções fitdistr (para obtenção das estimativas dos parâmetros)</div><div><div>e rgamma, no entanto tenho bastante zeros (dias sem chuva), com isso não é possível a obtenção dos parâmetros.</div>
<div><br></div><div>Erro apresentado no R:</div><div><br></div><div><div>Erro em optim(x = c(0, 29.9, 0, 0.7, 0, 0, 8.7, 0, 0, 26.8, 0, 19.4, : </div><div> valor inicial em vmmin não é finito</div></div><div><br></div>
<div>
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.</div><div>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?</div>
<div><br></div>Antecipadamente meu obrigada.</div><div>Att,<br><br style="font-family:Tahoma;font-size:13px;"><font size="4" face="'Times New Roman'">Natália da Silva Martins<br>Bacharel em Estatística - Universidade Estadual de Maringá/ UEM<br>
Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP<br>Contato: <a href="" rel="nofollow">(19) 8306-4743</a></font> <br><div><br></div><br>
</div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a rel="nofollow" ymailto="mailto:R-br@listas.c3sl.ufpr.br" target="_blank" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a rel="nofollow" target="_blank" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a rel="nofollow" target="_blank" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div>
</div><br>_______________________________________________<br>R-br mailing list<br><a ymailto="mailto:R-br@listas.c3sl.ufpr.br" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br><br> </div> </div> </div></body></html>