parametros weibull
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 SchumacherCurso de graduação em meteorologia - UFPELBolsista do Programa de Educação Tutorial - PET
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.
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)
participantes (3)
-
Bruno César Pasquini -
Frederico Stienmetz Alvarez -
Vanúcia Schumacher