[R-br] parametros weibull

[Bruno César Pasquini]

Vanúcia,

O pacote MASS tem uma função chamada "fitdistr" que estima os parâmetros de
várias distribuições por máxima verossimilhança.

Não me lembro de já ter ajustado uma weibull, mas é uma das distribuições
possíveis de se estimar parâmetros.

Salvo engano, a sintaxe, tendo como base o código que você passou (e tendo
instalado o pacote MASS) seria:

library(MASS) ## para carregar o pacote

fitdistr(dados, densfun="weibull") ## para fazer as estimações dos
parâmetros

Mais detalhes você encontra pedindo um help (literalmente)...

?fitdistr

Espero que a ajude. Atenciosamente,

Bruno Pasquini.

========================================================

Em 18 de fevereiro de 2012 23:08, Vanúcia Schumacher <
vanucia-schumacher em hotmail.com> escreveu:

>  Olá pessoal,
>
> Gostaria de ajuda para calcular os parâmetros de weibull, pois as linhas
> de comando que estou utilizando são para
> calcular pelo método dos momentos, porém , está me retornando valores
> negativos, se alguém que entenda melhor poder ajudar agradeço.
> Talvez tenha uma maneira melhor de calcular
>
> dados<-read.table("mediaPOAinverno.txt",header=FALSE)
>
> med.weibull<-mean(dados) ## sample mean
>
> var.weibull<-var(dados) ## sample variance
>
> l.est<-med.weibull/var.weibull ## lambda estimate (corresponds to rate)
>
> a.est<-((med.weibull)^2)/var.weibull ## alfa estimate
>
>
> Valeu
>
>
>
> Vanúcia Schumacher
> Curso de graduação em meteorologia - UFPEL
> Bolsista do Programa de Educação Tutorial - PET
>
> _______________________________________________
> R-br mailing list
> R-br em 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120218/85af3041/attachment.html>