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.

> a = fitdistr(dados, "Weibull")
> a = unlist(a)

Saudações.

Em 18 de fevereiro de 2012 22:25, Bruno César Pasquini <bruno.pasquini@gmail.com> escreveu:
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@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@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.


_______________________________________________
R-br mailing list
R-br@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.



--
Frederico Steinmetz Alvarez
Engenheiro de Produção Elétrica
Mestrando em Energia Eólica
Laboratório Didático Pedagógico de Fluidos - LDPFLU
Departamento de Engenharia Mecânica - DEMEC
Universidade Federal de Pernambuco - UFPE
(81) 8885-1105 (Oi)
(48) 9143-5141 (Vivo)