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 (Pessoal) Tel (19) 3375-1196 (Residêncial) Tel (19) 3447-8537 (Trabalho)
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@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í 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@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@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@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)
A mensagem de erro diz "Weibull values must be > 0". O teu summary() mostra que o valor mínimo é 0. Só o X50m que não é assim, justamente é o que você conseguiu ajustar. O vento realmento foi zero? Ou foi um valor tão pequeno que o aparelho de medição censurou para zero? Não seria o caso de usar outra distribuição que tenha suporte no [0,inf), ou somar uma quantidade pequena aos valores zeros (0+0.000001), ou usar uma mistura de distribuições (uma para X=0, outra para X>0), ou usar a Weibull com censura? Leia o guia de postagem! À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
Muito obrigado, vou tentar seguir suas orientações. Em 4 de dezembro de 2011 22:39, Walmes Zeviani <walmeszeviani@gmail.com>escreveu:
A mensagem de erro diz "Weibull values must be > 0". O teu summary() mostra que o valor mínimo é 0. Só o X50m que não é assim, justamente é o que você conseguiu ajustar. O vento realmento foi zero? Ou foi um valor tão pequeno que o aparelho de medição censurou para zero? Não seria o caso de usar outra distribuição que tenha suporte no [0,inf), ou somar uma quantidade pequena aos valores zeros (0+0.000001), ou usar uma mistura de distribuições (uma para X=0, outra para X>0), ou usar a Weibull com censura? Leia o guia de postagem!
À disposição. Walmes.
========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
_______________________________________________ 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.
-- *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)
participantes (3)
-
Cesar Rabak -
Rodrigo Muniz -
Walmes Zeviani