Distribuição densidade de probabilidade de Weibull
Olá Pessoal, Sou usuário recente do R e como tal ainda tenho muitas dúvidas básicas na utilização do programa. Tenho uma série de dados de velocidade de vento na qual quero ajustar os parâmetros c e K da distribuição de Wuibull, no entanto quando uso as funções abaixo, da um erro que x deve ser um vetor numérico, no entanto se eu colocar um valor numérico ele gera uma série de números pseudoaleatórios, O que eu quero abrir os meus dados dentro da função. Se alguém puder me ajudar eu agradeço. Vento50m <- dweibull(Vento, shape=2, scale = 1) fitdistr(Vento50m, densfun="weibull") ATT. -- *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)
Não entedi! Você quer gerar dados de uma Weibull ou você quer estimar os parâmetros a partir de uma amostra que você tem? Esclareça o que você quer e leia com atenção a documentação das funções. À 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 ==========================================================================
Eu quero estimar os parâmetros a partir de uma série de dados que eu já tenho. No entanto os comando para a função da distribuição de weibull no R não aceita que eu insira o nome do arquivo na função, apenas um vetor numérico, por exemplo: fitdistr(n=100, densfun="weibull") ou Vento50m <- dweibull(n=100, shape=2, scale = 1) nestas duas formas o programa gera números aleatórios. Obrigado Em 25 de novembro de 2011 13:12, Walmes Zeviani <walmeszeviani@gmail.com>escreveu:
Não entedi! Você quer gerar dados de uma Weibull ou você quer estimar os parâmetros a partir de uma amostra que você tem? Esclareça o que você quer e leia com atenção a documentação das funções.
À 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)
Mas se vc quer estimar os parâmtros, pq o shape=2 e scale=1? Tente meudado<- read.table("caminhodo arquivo) # veja o help do read.table fitdistr(meudado, densfun="weibul") Deve funcionar... Em 25 de novembro de 2011 13:22, Rodrigo Muniz <uenf.rodrigo@gmail.com>escreveu:
Eu quero estimar os parâmetros a partir de uma série de dados que eu já tenho. No entanto os comando para a função da distribuição de weibull no R não aceita que eu insira o nome do arquivo na função, apenas um vetor numérico, por exemplo:
fitdistr(n=100, densfun="weibull") ou Vento50m <- dweibull(n=100, shape=2, scale = 1)
nestas duas formas o programa gera números aleatórios. Obrigado
Em 25 de novembro de 2011 13:12, Walmes Zeviani <walmeszeviani@gmail.com>escreveu:
Não entedi! Você quer gerar dados de uma Weibull ou você quer estimar os parâmetros a partir de uma amostra que você tem? Esclareça o que você quer e leia com atenção a documentação das funções.
À 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)
_______________________________________________ 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.
-- Fernando A.B. Colugnati
Esse procedimento de read.table eu já fix, até consigo plotar gráficos com os meus dados, no entanto a função fitdistr(meudado, densfun="weibull") dá um erro quando insere o nome do arquivo de dados, o erro diz que o meu dado não é um vetor numérico. Se eu substituir o meu meu dado por um valor numérico, a função roda, porem com uma série de números aleatório. OBG Em 25 de novembro de 2011 13:37, Fernando Colugnati <fernando@ipti.org.br>escreveu:
Mas se vc quer estimar os parâmtros, pq o shape=2 e scale=1? Tente
meudado<- read.table("caminhodo arquivo) # veja o help do read.table
fitdistr(meudado, densfun="weibul")
Deve funcionar...
Em 25 de novembro de 2011 13:22, Rodrigo Muniz <uenf.rodrigo@gmail.com>escreveu:
Eu quero estimar os parâmetros a partir de uma série de dados que eu já
tenho. No entanto os comando para a função da distribuição de weibull no R não aceita que eu insira o nome do arquivo na função, apenas um vetor numérico, por exemplo:
fitdistr(n=100, densfun="weibull") ou Vento50m <- dweibull(n=100, shape=2, scale = 1)
nestas duas formas o programa gera números aleatórios. Obrigado
Em 25 de novembro de 2011 13:12, Walmes Zeviani <walmeszeviani@gmail.com>escreveu:
Não entedi! Você quer gerar dados de uma Weibull ou você quer estimar os parâmetros a partir de uma amostra que você tem? Esclareça o que você quer e leia com atenção a documentação das funções.
À 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)
_______________________________________________ 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.
-- Fernando A.B. Colugnati
_______________________________________________ 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)
experimente descobrir quais sao as classes de dados q vc tem em 'meudado' str(meudado) ja' te dara' uma ideia... se vc vir algo como 'data.frame' e observar que o primeiro elemento eh o conjunto de dados de interesse, vc deve querer usar algo como: fitdistr(meudado[[1]], densfun="weibull") (acima funcionara' se is.numeric(meudado[[1]]) devolver TRUE) b 2011/11/25 Rodrigo Muniz <uenf.rodrigo@gmail.com>:
Esse procedimento de read.table eu já fix, até consigo plotar gráficos com os meus dados, no entanto a função fitdistr(meudado, densfun="weibull") dá um erro quando insere o nome do arquivo de dados, o erro diz que o meu dado não é um vetor numérico. Se eu substituir o meu meu dado por um valor numérico, a função roda, porem com uma série de números aleatório. OBG
Em 25 de novembro de 2011 13:37, Fernando Colugnati <fernando@ipti.org.br> escreveu:
Mas se vc quer estimar os parâmtros, pq o shape=2 e scale=1? Tente meudado<- read.table("caminhodo arquivo) # veja o help do read.table fitdistr(meudado, densfun="weibul") Deve funcionar...
Em 25 de novembro de 2011 13:22, Rodrigo Muniz <uenf.rodrigo@gmail.com> escreveu:
Eu quero estimar os parâmetros a partir de uma série de dados que eu já tenho. No entanto os comando para a função da distribuição de weibull no R não aceita que eu insira o nome do arquivo na função, apenas um vetor numérico, por exemplo: fitdistr(n=100, densfun="weibull") ou Vento50m <- dweibull(n=100, shape=2, scale = 1) nestas duas formas o programa gera números aleatórios. Obrigado Em 25 de novembro de 2011 13:12, Walmes Zeviani <walmeszeviani@gmail.com> escreveu:
Não entedi! Você quer gerar dados de uma Weibull ou você quer estimar os parâmetros a partir de uma amostra que você tem? Esclareça o que você quer e leia com atenção a documentação das funções.
À 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)
_______________________________________________ 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.
-- Fernando A.B. Colugnati
_______________________________________________ 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)
_______________________________________________ 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.
Um vetor númerico é uma coluna do seu data.frame para o qual a função class() retorna "numeric" ou "integer". É o vetor que você tem que passar para a fitdistr(). Se você tivesse fornecido um CMR, com uma amostra da estrutura dos seus dados, tudo teria-se resolvido na segunda mensagem. Para um entendimento sobre os tipos de objetos do R, consulte os textos disponíveis no CRAN e dê uma olhada nesse script ( http://www.leg.ufpr.br/~walmes/cursoR/cnpaf/cap01objetos-iso.R) para saber quais os tipos de objetos e como acessá-los. À 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 ==========================================================================
participantes (4)
-
Benilton Carvalho -
Fernando Colugnati -
Rodrigo Muniz -
Walmes Zeviani