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">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 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>_______________________________________________<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>