<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style></head>
<body class='hmmessage'><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><div><br><br><span class="ecxApple-style-span" style="line-height:17px"><font class="ecxApple-style-span" face="'Comic Sans MS'" color="#002060">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> </div></body>
</html>