Boa noite pessoal.<div>Estou com o seguinte script:</div><div><div>require(MASS)</div><div>getwd()</div><div>setwd("C:/Users/Rodrigo/Documents/Pos_graduacao/DoutoradoUSP/Geoestatistica")</div><div>Vento <- read.table("Weibull_full.txt", head=T)</div>
<div><br></div><div>WVento50m <- fitdistr(Vento[[1]], densfun="weibull")</div><div>WVento30m <- fitdistr(Vento[[2]], densfun="weibull")</div><div>WVento15m <- fitdistr(Vento[[3]], densfun="weibull")</div>
<div>WVento10m <- fitdistr(Vento[[4]], densfun="weibull")</div><div>WVento3m <- fitdistr(Vento[[5]], densfun="weibull")</div><div><br></div><div>O arquivo Weibull_full é um data frame com 5 coluna sendo a variável velocidade do vento em cinco diferentes </div>
<div>níveis de altura. Estou tentanto calcular os parâmetros de forma e escala da distribuição de Weibull para cada nível </div><div>de altura. No entanto, só consigo calcular para o nível de altura de 50 m (WVento50m), quando mando rodar os outros </div>
<div>níveis dá o seguinte erro:</div><div> <span class="Apple-style-span" style="font-family: 'Lucida Console'; font-size: 13px; line-height: 17px; white-space: pre-wrap; ">Error in fitdistr(Vento[[2]], densfun = "weibull") : </span></div>
<span class="Apple-style-span" style="border-collapse: separate; font-family: Arial; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; font-size: medium; "><span class="Apple-style-span" style="font-family: 'Lucida Console'; font-size: 13px; line-height: 17px; text-align: -webkit-left; white-space: pre-wrap; "><pre tabindex="0" class="GD40030CLR" style="font-family: 'Lucida Console'; font-size: 10pt !important; outline-style: none; outline-width: initial; outline-color: initial; border-top-style: none; border-right-style: none; border-bottom-style: none; border-left-style: none; border-width: initial; border-color: initial; white-space: pre-wrap !important; margin-top: 0px; margin-bottom: 0px; line-height: 1.3; ">
<span class="GD40030CGR"> Weibull values must be > 0</span></pre></span></span><div><br></div><div>Já conferi o arquivo dos dado para ver se há algo errado com os dados, mais não há. Tentei fazer a análise separada</div>
<div>mais deu o mesmo erro. Um observação que eu faço é que embora eu não consiga rodar a função "fitdistr", eu consigo plotar</div><div>os gráficos de histograma normalmente para todos os níveis de altura. </div>
<div>Se alguém puder me ajudar eu agradeço.</div><div>Abraço a todos.</div>-- <br><b>Rodrigo A. Muniz</b><div>Eng. Agrônomo. Ms Produção Vegetal (UENF)</div><div>Doutorando em Engenharia de Sistemas Agrícolas (ESALQ/USP)</div>
<div>Cel <a href="tel:%2819%29%208300-4333" value="+551983004333" target="_blank">(19) 8300-4333</a> (Pessoal)</div><div>Tel <a href="tel:%2819%29%203375-1196" value="+551933751196" target="_blank">(19) 3375-1196</a> (Residêncial)</div>
<div>Tel <a href="tel:%2819%29%203447-8537" value="+551934478537" target="_blank">(19) 3447-8537</a> (Trabalho)</div><br>
</div>