[R-br] Distribuição de Weibull

[Rodrigo Muniz]

Aí segue o resultado do summary.

> summary(Vento)      X50m             X30m             X15m             X10m             X3m
 Min.   : 0.166   Min.   : 0.000   Min.   : 0.000   Min.   : 0.000
Min.   :0.000
 1st Qu.: 3.225   1st Qu.: 2.370   1st Qu.: 1.567   1st Qu.: 1.154
1st Qu.:0.576
 Median : 4.686   Median : 3.570   Median : 2.494   Median : 1.918
Median :1.194
 Mean   : 4.899   Mean   : 3.887   Mean   : 2.807   Mean   : 2.205
Mean   :1.411
 3rd Qu.: 6.240   3rd Qu.: 5.062   3rd Qu.: 3.715   3rd Qu.: 2.977
3rd Qu.:2.052
 Max.   :17.870   Max.   :13.890   Max.   :12.210   Max.   :11.520
Max.   :9.680



Em 4 de dezembro de 2011 21:32, Cesar Rabak <cesar.rabak em gmail.com>escreveu:

> Rodrigo,
>
> Você não seguiu as recomendações de postar um CMR:
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-**guia<http://www.leg.ufpr.br/r-br-guia>)
> e forneça código mínimo reproduzível.
>
> Para que a gente tenha uma idéia mínima do que pode estar acontecendo,
> seria necessário entender o que há no seu dataframe "Vento".
>
> poste aqui o resultado de:
>
> > summary(Vento)
>
> E a gente segue daí.
>
> --
> Cesar Rabak
> GNU/Linux User 52247.
> Get counted: http://counter.li.org/
>
>
> Em 4/12/2011 20:14, Rodrigo Muniz escreveu:
>
>> Boa noite pessoal.
>> Estou com o seguinte script:
>> require(MASS)
>> getwd()
>> setwd("C:/Users/Rodrigo/**Documents/Pos_graduacao/**
>> DoutoradoUSP/Geoestatistica")
>> Vento <- read.table("Weibull_full.txt", head=T)
>>
>> WVento50m <- fitdistr(Vento[[1]], densfun="weibull")
>> WVento30m <- fitdistr(Vento[[2]], densfun="weibull")
>> WVento15m <- fitdistr(Vento[[3]], densfun="weibull")
>> WVento10m <- fitdistr(Vento[[4]], densfun="weibull")
>> WVento3m <- fitdistr(Vento[[5]], densfun="weibull")
>>
>> O arquivo Weibull_full é um data frame com 5 coluna sendo a variável
>> velocidade do vento em cinco diferentes
>> níveis de altura. Estou tentanto calcular os parâmetros de forma e
>> escala da distribuição de Weibull para cada nível
>> de altura. No entanto, só consigo calcular para o nível de altura de 50
>> m (WVento50m), quando mando rodar os outros
>> níveis dá o seguinte erro:
>> Error in fitdistr(Vento[[2]], densfun = "weibull") :
>>
>>   Weibull values must be>  0
>>
>>
>> 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
>> 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
>> os gráficos de histograma normalmente para todos os níveis de altura.
>> Se alguém puder me ajudar eu agradeço.
>> Abraço a todos.
>> --
>> *Rodrigo A. Muniz*
>>
>> Eng. Agrônomo. Ms Produção Vegetal (UENF)
>> Doutorando em Engenharia de Sistemas Agrícolas (ESALQ/USP)
>> Cel (19) 8300-4333 <tel:%2819%29%208300-4333> (Pessoal)
>> Tel (19) 3375-1196 <tel:%2819%29%203375-1196> (Residêncial)
>> Tel (19) 3447-8537 <tel:%2819%29%203447-8537> (Trabalho)
>>
>>
>>
>> ______________________________**_________________
>> R-br mailing list
>> R-br em listas.c3sl.ufpr.br
>> https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-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<http://www.leg.ufpr.br/r-br-guia>)
>> e forneça código mínimo reproduzível.
>>
>
>
> ______________________________**_________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-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<http://www.leg.ufpr.br/r-br-guia>)
> e forneça código mínimo reproduzível.
>



-- 
*Rodrigo A. Muniz*
Eng. Agrônomo. Ms Produção Vegetal (UENF)
Doutorando em Engenharia de Sistemas Agrícolas (ESALQ/USP)
Cel (19) 8300-4333 (Pessoal)
Tel (19) 3375-1196 (Residêncial)
Tel (19) 3447-8537 (Trabalho)
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111204/1824d961/attachment.html>