[R-br] Distribuição de Weibull

[Cesar Rabak]

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) 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
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.