Vanúcia, Uma outra forma de efetuar esse procedimento é o indicado na norma IEC 61400-12-1. Nela, os dados são lidos a uma frequência de 1 Hz, e depois, é calculada a média para intervalos de 10 minutos. Com isso, a probabilidade de uma média ser igual a zero é muito pequena. Veja se é possível efetuar esse procedimento nos seus dados. Saudações. Em 21 de fevereiro de 2012 23:38, Vanúcia Schumacher < vanucia-schumacher@hotmail.com> escreveu:
Vanúcia Schumacher Curso de graduação em meteorologia - UFPEL Bolsista do Programa de Educação Tutorial - PET
From: vanucia-schumacher@hotmail.com To: tura@centroin.com.br; r-help@r-project.org Date: Tue, 21 Feb 2012 22:03:09 +0000 Subject: Re: [R] HELP ERROR Weibull values must be > 0
Tura, meus dados possuem mais de 700 elementos, fica muito inviavel passa-los aqui por summarya questão e que nesses dados eu tenho valores ausentes e valores 0, li algumas sugestões em que deve-se mensurar os valores em 0 para 0.00001,sendo assim ele roda tranquilo , mas retorna um aviso de Nas, mas me retorna os parametros. Não sei se estou fazendo a coisa certa, mas como estou querendo ajustar a distribuição m dados de vento paraprodução de energia eólica, conforme varios artigos que li, essa seria a melhor distribuição. Se puder ajudar fico muito agradecida.
library(MASS)> dados<-read.table("inverno.txt",header=FALSE)> vento50<-fitdistr(dados[[1]],densfun="weibull")Mensagens de aviso perdidas:1: In dweibull(x, shape, scale, log) : NaNs produzidos2: In dweibull(x, shape, scale, log) : NaNs produzidos3: In dweibull(x, shape, scale, log) : NaNs produzidos4: In dweibull(x, shape, scale, log) : NaNs produzidos5: In dweibull(x, shape, scale, log) : NaNs produzidos6: In dweibull(x, shape, scale, log) : NaNs produzidos> vento50 shape scale 0.61473436 2.43536064 (0.01895222) (0.13202750)
Vanúcia SchumacherCurso de graduação em meteorologia - UFPELBolsista do Programa de Educação Tutorial - PET
Date: Tue, 21 Feb 2012 18:40:47 -0200 From: tura@centroin.com.br To: vanucia-schumacher@hotmail.com Subject: Re: [R] HELP ERROR Weibull values must be > 0
On 02/21/2012 01:48 PM, Vanúcia Schumacher wrote:
GUYS,
I NEED HELP WITH ERROR:
library(MASS)
dados<-read.table("mediaRGinverno.txt",header=FALSE) vento50<-fitdistr(dados[[1]],densfun="weibull") Erro em fitdistr(dados[[1]], densfun = "weibull") : Weibull values must be> 0 WHY RETURN THIS ERROR? WHAT CAN I DO?
Vanúcia
Vou lhe ajudar dentro do possível mas antes de tudo se seus dados tem valores menores ou iguais a zero eles não pode ser ajustados numa distribuição weibull pois esta só permite valores >= 0.
Que tal você me mandar o resultado do summary(dados)
Abracos Tura
[[alternative HTML version deleted]]
______________________________________________ R-help@r-project.orgmailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
_______________________________________________ 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)