A função fitdisr é a que utilizo para estimação dos parâmetros para medição de ventos e é muito prática. Só lembrando, caso queira utilizar os parâmetros encontrados, salve numa variável e chame a função unlist.<div><br></div>
<div>> a = fitdistr(dados, "Weibull")</div><div>> a = unlist(a)</div><div><br></div><div>Saudações.<br><br><div class="gmail_quote">Em 18 de fevereiro de 2012 22:25, Bruno César Pasquini <span dir="ltr"><<a href="mailto:bruno.pasquini@gmail.com">bruno.pasquini@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Vanúcia,<br><br>O pacote MASS tem uma função chamada "fitdistr" que estima os parâmetros de várias distribuições por máxima verossimilhança.<br>
<br>Não me lembro de já ter ajustado uma weibull, mas é uma das distribuições possíveis de se estimar parâmetros.<br>
<br>Salvo engano, a sintaxe, tendo como base o código que você passou (e tendo instalado o pacote MASS) seria:<br><br>library(MASS) ## para carregar o pacote<br><br>fitdistr(dados, densfun="weibull") ## para fazer as estimações dos parâmetros<br>
<br>Mais detalhes você encontra pedindo um help (literalmente)...<br><br>?fitdistr<br><br>Espero que a ajude. Atenciosamente,<br><br>Bruno Pasquini.<br><br>========================================================<br><br>
<div class="gmail_quote">
Em 18 de fevereiro de 2012 23:08, Vanúcia Schumacher <span dir="ltr"><<a href="mailto:vanucia-schumacher@hotmail.com" target="_blank">vanucia-schumacher@hotmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><div class="h5">
<div><div dir="ltr">
Olá pessoal,<div><br></div><div>Gostaria de ajuda para calcular os parâmetros de weibull, pois as linhas de comando que estou utilizando são para </div><div>calcular pelo método dos momentos, porém , está me retornando valores negativos, se alguém que entenda melhor poder ajudar agradeço.</div>
<div>Talvez tenha uma maneira melhor de calcular</div><div><br></div><div><div>dados<-read.table("mediaPOAinverno.txt",header=FALSE)</div><div><br></div><div>med.weibull<-mean(dados) ## sample mean</div><div>
<br></div><div>var.weibull<-var(dados) ## sample variance</div><div><br></div><div>l.est<-med.weibull/var.weibull ## lambda estimate (corresponds to rate)</div><div><br></div><div>a.est<-((med.weibull)^2)/var.weibull ## alfa estimate</div>
</div><div><br></div><div><br></div><div>Valeu</div><div><br></div><font color="#888888"><div><br><br><span style="line-height:17px"><font color="#002060" face="'Comic Sans MS'">Vanúcia Schumacher<br style="line-height:17px">
Curso de graduação em meteorologia - UFPEL<br style="line-height:17px">Bolsista do Programa de Educação Tutorial - PET</font><br style="line-height:17px"></span></div> </font></div></div>
<br></div></div>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">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></blockquote></div><br>
<br>_______________________________________________<br>
R-br mailing list<br>
<a 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></blockquote></div><br><br clear="all"><div><br></div>-- <br>
Frederico Steinmetz Alvarez<br>Engenheiro de Produção Elétrica<div>Mestrando em Energia Eólica</div><div>Laboratório Didático Pedagógico de Fluidos - LDPFLU</div><div>Departamento de Engenharia Mecânica - DEMEC</div><div>
Universidade Federal de Pernambuco - UFPE<br>(81) 8885-1105 (Oi)<br>(48) 9143-5141 (Vivo)<br></div><br>
</div>