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
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@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@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.
Estima com a glm() ou programe sua própria função de verossimilhança. É simples. À 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@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
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 <rcoster@gmail.com> Para: r-br@listas.c3sl.ufpr.br 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 <nsmbarreto@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@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.
_______________________________________________ R-br mailing list R-br@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.
Dados de chuva geralmente são correlacionados no tempo, sera q uma bernoulli ou binomial satisfaz? Em 19 de julho de 2012 12:56, Edson Lira <edinhoestat@yahoo.com.br>escreveu:
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 <rcoster@gmail.com> *Para:* r-br@listas.c3sl.ufpr.br *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 <nsmbarreto@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@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.
_______________________________________________ R-br mailing list R-br@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.
_______________________________________________ R-br mailing list R-br@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.
-- Luis Iván Ortiz Valencia Doutorando Saúde Pública - Epidemiologia, IESC, UFRJ Estatístico Msc. Spatial Analyst Msc.
participantes (5)
-
Edson Lira -
Luis Iván Ortiz Valencia -
Natalia Martins -
Rodrigo Coster -
Walmes Zeviani