<html><body><div style="color:#000; background-color:#fff; font-family:Courier New, courier, monaco, monospace, sans-serif;font-size:10pt"><div><span><br></span></div><div>
<font size="2" face="Arial">Vanúcia,</font>
<br> </div><div><font size="2" face="Arial"><br></font></div><div><font size="2" face="Arial">Neste caso, creio que a mensagem de aviso (que não é de erro) pode ser relevada. Veja o exemplo abaixo:<br></font></div><div><font size="2" face="Arial"><br></font></div><div><font face="Courier New, courier, monaco, monospace, sans-serif"><br></font></div><div><font face="Courier New, courier, monaco, monospace, sans-serif">x<-rweibull(1e3,1,1)</font><br><font face="Courier New, courier, monaco, monospace, sans-serif"><br><br>ll<-function(th,x){<br> sh<-th[1]; sc<-th[2]<br> l<-dweibull(x,sh,sc,log=T)<br> -sum(l)<br>}<br><br># Valores iniciais -- função 'fitdistr'<br>lx <- log(x)<br>m <- mean(lx)<br>v <- var(lx)<br>shape <- 1.2/sqrt(v)<br>scale <- exp(m + 0.572/shape)<br>ini<-c(shape,scale)<br><br># Limites dos
parâmetros<br>inf<-rep(1e-3,2)<br>sup<-rep(Inf,2)<br><br><br>fit<-optim(ini,ll,method="L-BFGS-B",lower=inf,upper=sup,hessian=T,x=x)<br><br>fit$par # Estimativa<br>sqrt(diag(solve(fit$hessian))) # Erro padrão<br><br><br>library(MASS)<br>fitdistr(x,densfun="weibull")</font><br><font size="2" face="Arial"><br></font></div><div><font size="2" face="Arial"><br></font></div><div><font size="2" face="Arial"><br></font></div><div><font size="2" face="Arial">Att.,<br></font></div><div><font size="2" face="Arial">Rubem<br></font></div><div style="font-size: 10pt; font-family: "Courier New", "courier", "monaco", monospace, sans-serif"> </div> </div></body></html>