Janeiro 2018 - R-br

previsão utilizando o lasso
by João Pedro Domingues 16 Jan '18

16 Jan '18

Boa tarde colegas, Estou tentando realizar uma previsão de vendas utilizando o algoritimo LASSO pelo pacote HDeconometrics. Alguém com experiência neste pacote poderia me ajudar por gentileza? Só preciso resolver isso para finalizar minha dissertação. Muito obrigado a todos! Data é uma matriz de series temporais com 158 observações (diarias) e 13 produtos (colunas). Eu inicialmente separo a primeira coluna como sendo o produto focal para analisar e o coloco como y, sendo a variável dependente. Já as outras colunas eu coloco como x, variáveis independentes. Então o LASSO me informa quais são relevantes. Com isso, eu separo as primeiras 148 observações para ser o training set, e as últimas 10 observações para o test set, e ver se o modelo realmente funciona. Até aqui ok. O problema: quando eu realizo o ajuste da série com o comando (lasso=ic.glmnet(x.in,y.in,crit = "bic")) ele ajusta o modelo para a série sem problemas, utilizando as variáveis que o lasso identificou como relevantes para o training set. Porém, quando eu vou executar a previsão real utilizando a última linha (previsao.lasso=predict(lasso,newdata=x.out)), ele não faz a previsão, ele simplesmente faz um ajuste igual o training set. Alguém tem alguma ideia de como resolver essa questão e conseguir prever corretamente? Grande abraço Segue abaixo os comandos utilizados: library(HDeconometrics) library(forecast) ## Inicio i = 0 y = as.matrix(Data[,i+1]) #variável dependente primeira coluna x = (Data) #cópia da base toda x[,i + 1] <- NULL #retira a variável y e fica com todas as outras variáveis x = as.matrix(x) # transforma em matriz ### separa a série em training e test set de x e y y.in=y[1:148] #training set y.out=y[-c(1:148)] #test set x.in=x[1:148,] #training set x.out=x[-c(1:148),] #test set ## ajuste do modelo e previsão LASSO lasso=ic.glmnet(x.in,y.in,crit = "bic") #ajuste do modelo com o training set previsao.lasso=predict(lasso,newdata=x.out) #previsão com o test set João Pedro Domingues

2 1

Problema ao gerar um mapa com dados estatísticos
by Edson Lira 13 Jan '18

13 Jan '18

Caros amigos, abaixo os arquivos (links) que estou usando para gerar um mapa do estado do Amazonas com o IDH Municipal de 2000 (está 2005, mas é 2000). Quando rodo a rotina para fazer a categorização do IDH (cut) me dá a mensagem: > am$idhc2000 <- cut(am$idh2000, 4) Error in seq.int(rx[1L], rx[2L], length.out = nb) : 'from' must be a finite number Além disso: Warning messages: 1: In min(x) : nenhum argumento não faltante para min; retornando Inf 2: In max(x) : nenhum argumento não faltante para max; retornando -Inf > Não consegui descobrir o que está errado. Bom dia!!! Download amazonas idh 2005.R | | | Download amazonas idh 2005.R | | | Download amazonas calha.csv | | | Download amazonas calha.csv | | | Download 13MUE250GC SIR.shx | | | Download 13MUE250GC SIR.shx | | | Download 13MUE250GC SIR.dbf | | | Download 13MUE250GC SIR.dbf | | | Download 13MUE250GC SIR.cpg | | | Download 13MUE250GC SIR.cpg | | | https://www.datafilehost.com/d/27e20ca4 # SHAPE [ ]'s Prof. Edson Lira, Me Estatístico Manaus-Amazonas

2 3

Problema na instalação do R base
by Jefferson Ferreira-Ferreira 13 Jan '18

13 Jan '18

Pessoal, estou passando por um problema estranho. O R está (aparentemente) installado na minha máquina (Ubuntu 16.04). Quando tento instalar, o sistema retorna: ele retorna "r-base-dev is already the newest version (3.4.3-1xenial0). 0 upgraded, 0 newly installed, 0 to remove and 23 not upgraded." No entanto, o RStudio não abre porque não é capaz de localizar a instalação e o comando 'which R' retorna "not found". Já desisntalei e instalei novamente, mas o mesmo continua ocorrendo. Alguma idea? Se alguém puder me ajudar, agradeço. Estou em apuros com meu trabalho sem o R. Eu tambpem postei uma pergunta aqui https://askubuntu.com/questions/995618/software-is-installed-but-which-retu… Obrigado e abraços. *Jefferson Ferreira-Ferreira, **PhD (abd)* *Geographer* *Ecosystem Dynamics Observatory <http://tscanada.wix.com/ecodyn> - EcoDyn/UNESP* *Department of **Geography * *Institute of Geosciences and Exact Sciences** (IGCE)* *São Paulo State University (UNESP)* *Rio Claro, SP - Brazil*

1 0

Instalar pacote do github
by Maurício Lordêlo 12 Jan '18

12 Jan '18

Caros, Certa vez eu consegui instalar o pacote "ggbiplot" que considero bem interessante. Não sei se pela atualização do R, ou de mudança de versão do Windows, não mais consigo instalar este pacote. Pelo que tenho registrado nos meus scripts, eu seguia os passos abaixo: library(devtools) install.packages("githubinstall") library(githubinstall) install.packages("ggbiplot-master") library(ggbiplot-master) install_github("ggbiplot","vqv") Também fiz estas duas tentativas sem sucesso install.packages("diretório de trabalho/ggbiplot-master.zip") #aqui eu fiz o download do arquivo zip e salvei no "diretório de trabalho" install.packages("https://github.com/vqv/ggbiplot.git") Alguém saberia informar como resolver este problema?

2 1

Re: [R-br] Instalar pacote do github
by sznelwar＠uol.com.br 11 Jan '18

11 Jan '18

3 2

Re: [R-br] Digest R-br, volume 85, assunto 6
by Edmar Caldas 11 Jan '18

11 Jan '18

Eu entendi assim, Vamos supor que as suas variáveis (dependentes e independentes) tem valores até 12/2017. e vc quer fazer previsão para 12/2018. as suas variáveis independentes terão que ter valores até 12/2018 senão não vai fazer a previsão.Edmar Em quarta-feira, 10 de janeiro de 2018 12:00:33 BRST, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: Enviar submissões para a lista de discussão R-br para r-br(a)listas.c3sl.ufpr.br Para se cadastrar ou descadastrar via WWW, visite o endereço https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br ou, via email, envie uma mensagem com a palavra 'help' no assunto ou corpo da mensagem para r-br-request(a)listas.c3sl.ufpr.br Você poderá entrar em contato com a pessoa que gerencia a lista pelo endereço r-br-owner(a)listas.c3sl.ufpr.br Quando responder, por favor edite sua linha Assunto assim ela será mais específica que "Re: Contents of R-br digest..." Tópicos de Hoje: 1. previsão utilizando o lasso (João Pedro Domingues) 2. Re: Distribuição de probabilidade (Wagner Bonat) 3. Re: Distribuição de probabilidade (Wagner Bonat) ---------------------------------------------------------------------- Message: 1 Date: Tue, 9 Jan 2018 20:04:49 +0000 From: João Pedro Domingues <J__P(a)hotmail.com> To: a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] previsão utilizando o lasso Message-ID: <BN6PR2001MB16668FB0A072A2FC1AD8D6B0D5100(a)BN6PR2001MB1666.namprd20.prod.outlook.com> Content-Type: text/plain; charset="iso-8859-1" Boa tarde colegas, Estou tentando realizar uma previsão de vendas utilizando o algoritimo LASSO pelo pacote HDeconometrics. Alguém com experiência neste pacote poderia me ajudar por gentileza? Só preciso resolver isso para finalizar minha dissertação. Muito obrigado a todos! Data é uma matriz de series temporais com 158 observações (diarias) e 13 produtos (colunas). Eu inicialmente separo a primeira coluna como sendo o produto focal para analisar e o coloco como y, sendo a variável dependente. Já as outras colunas eu coloco como x, variáveis independentes. Então o LASSO me informa quais são relevantes. Com isso, eu separo as primeiras 148 observações para ser o training set, e as últimas 10 observações para o test set, e ver se o modelo realmente funciona. Até aqui ok. O problema: quando eu realizo o ajuste da série com o comando (lasso=ic.glmnet(x.in,y.in,crit = "bic")) ele ajusta o modelo para a série sem problemas, utilizando as variáveis que o lasso identificou como relevantes para o training set. Porém, quando eu vou executar a previsão real utilizando a última linha (previsao.lasso=predict(lasso,newdata=x.out)), ele não faz a previsão, ele simplesmente faz um ajuste igual o training set. Alguém tem alguma ideia de como resolver essa questão e conseguir prever corretamente? Grande abraço Segue abaixo os comandos utilizados: library(HDeconometrics) library(forecast) ## Inicio i = 0 y = as.matrix(Data[,i+1]) #variável dependente primeira coluna x = (Data) #cópia da base toda x[,i + 1] <- NULL #retira a variável y e fica com todas as outras variáveis x = as.matrix(x) # transforma em matriz ### separa a série em training e test set de x e y y.in=y[1:148] #training set y.out=y[-c(1:148)] #test set x.in=x[1:148,] #training set x.out=x[-c(1:148),] #test set ## ajuste do modelo e previsão LASSO lasso=ic.glmnet(x.in,y.in,crit = "bic") #ajuste do modelo com o training set previsao.lasso=predict(lasso,newdata=x.out) #previsão com o test set João Pedro Domingues -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20180109/580bb5da/atta…> ------------------------------ Message: 2 Date: Tue, 9 Jan 2018 22:25:09 -0200 From: Wagner Bonat <wbonat(a)gmail.com> To: Wagner Wolff <wwolff(a)usp.br> Cc: a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Distribuição de probabilidade Message-ID: <CANt=4Mi4E58a3ikPDH2V9aN36b3d6D2M7HJ0V9u0W8yRHko=9w(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Olá Essas distribuições como a Weibul, Gumbel e similares são interessantes. Porém, nem sempre elas tem a tradicional parametrização de esperança e dispersão como a Tweedie e as distribuições mais clássicas (normal, gamma, normal inversa) tem. Muitas vezes a esperança ou é desconhecida ou uma função complicada de mais de um parâmetro. Eu penso que ao invés de tentar ajustar várias distribuições é mais interessante você procurar aspectos relevantes sobre o seu conjunto de dados, como por exemplo, excesso de zeros, forte assimetria, caudas pesadas, etc e verificar se a distribuição que você ajustando é capaz de descrever estes aspectos. No caso a Tweedie pode lidar com todos estes aspectos, o que significa que ela deve ajustar bem para uma grande quantidade de dados reais. O que eu quero dizer é que essa idéia de ficar comparando várias distribuições pode não ser realmente necessário. Eu exploro um pouco destas idéias em dois artigos. Eu mostro que a Tweedie ajusta muito bem, mesmo para distribuições que não fazem parte da familia, como a t e a slash. Talvez vc possa explorar idéias similares para a weibul, gumbel e similares. http://www.tandfonline.com/doi/abs/10.1080/00949655.2017.1318876?journalCod… http://journals.sagepub.com/doi/abs/10.1177/1471082X17715718 Att Wagner Em 9 de janeiro de 2018 01:11, Wagner Wolff <wwolff(a)usp.br> escreveu: > Olá Wagner > > Gostei da sua ideia sobre a familia de destribuição Tweedie. Além dessa > família você poderia indicar outras que englobem distribuições como por > exemplo, Weibull, Gumbel, Burr ... > > Ab > > 2018-01-08 18:18 GMT-02:00 Wagner Bonat via R-br <r-br(a)listas.c3sl.ufpr.br > >: > >> Quais distribuições vc tem interesse? >> A distribuição Tweedie tem como casos particulares algumas das mais >> populares distribuições como a normal, gamma e inversa normal. >> Assim, ao invés de ajustar um monte de distribuições basta ajustar uma >> que tem as outras como casos particulares. >> >> Att >> Wagner >> >> Em 8 de janeiro de 2018 10:38, Rodrigo Campos via R-br < >> r-br(a)listas.c3sl.ufpr.br> escreveu: >> >>> Bom dia Felipe, >>> Talvez o pacote "fitdistrplus" possa te ajudar! >>> Att. >>> >>> 2018-01-08 10:21 GMT-02:00 Felipe Felix Costa via R-br < >>> r-br(a)listas.c3sl.ufpr.br>: >>> >>>> olá, >>>> >>>> gostaria de saber se existe alguma forma para rodar todas as possíveis >>>> distribuições sobre um determinado conjunto de dados com variáveis >>>> continuas!? >>>> >>>> pois estou procurando e não achei nada pela internet falando dessa >>>> possibilidade. >>>> >>>> meu objetivo é verificar qual melhor a distribuição para os meus dados >>>> e, a partir dessa "melhor", estabelecer um modelo geral com todos os dados >>>> e comparar com modelos individuais de áreas isoladas. >>>> >>>> >>>> HELP >>>> >>>> Att., >>>> -- >>>> *Felipe Felix Costa* >>>> Tel. 96-98122-3077 >>>> >>>> >>>> _______________________________________________ >>>> R-br mailing list >>>> R-br(a)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 Campos >>> >>> >>> >>> _______________________________________________ >>> R-br mailing list >>> R-br(a)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. >>> >> >> >> >> -- >> Wagner Hugo Bonat >> ------------------------------------------------------------ >> ---------------------------------- >> Department of Mathematics and Computer Science (IMADA) >> University of Southern Denmark (SDU) and >> Laboratório de Estatística e Geoinformação (LEG) >> Universidade Federal do Paraná (UFPR) >> >> >> _______________________________________________ >> R-br mailing list >> R-br(a)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. >> > > > > -- > *Wagner Wolff, **PhD* > "*Luiz de Queiroz**" College of Agriculture,* > University of São Paulo > Pádua Dias avenue11 | 13418-900| Piracicaba-SP| Brazil > Phone: +55 19 982385582 <+55%2019%2098238-5582> > http://orcid.org/0000-0003-3426-308X > https://github.com/wwolff7 > http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4463141A1 > -- Wagner Hugo Bonat ---------------------------------------------------------------------------------------------- Department of Mathematics and Computer Science (IMADA) University of Southern Denmark (SDU) and Laboratório de Estatística e Geoinformação (LEG) Universidade Federal do Paraná (UFPR) -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20180109/86d7063c/atta…> ------------------------------ Message: 3 Date: Tue, 9 Jan 2018 22:28:25 -0200 From: Wagner Bonat <wbonat(a)gmail.com> To: Wagner Wolff <wwolff(a)usp.br> Cc: a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Distribuição de probabilidade Message-ID: <CANt=4Mj3Tk41SVoGUVAPTpA0h_TZCCQVyLjrYkYGWFpSHsCdHA(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Olá Essas distribuições como a Weibul, Gumbel e similares são interessantes. Porém, nem sempre elas tem a tradicional parametrização de esperança e dispersão como a Tweedie e as distribuições mais clássicas (normal, gamma, normal inversa) tem. Muitas vezes a esperança ou é desconhecida ou uma função complicada de mais de um parâmetro. Eu penso que ao invés de tentar ajustar várias distribuições é mais interessante você procurar aspectos relevantes sobre o seu conjunto de dados, como por exemplo, excesso de zeros, forte assimetria, caudas pesadas, etc e verificar se a distribuição que você ajustando é capaz de descrever estes aspectos. No caso a Tweedie pode lidar com todos estes aspectos, o que significa que ela deve ajustar bem para uma grande quantidade de dados reais. O que eu quero dizer é que essa idéia de ficar comparando várias distribuições pode não ser realmente necessário. Eu exploro um pouco destas idéias em dois artigos. Eu mostro que a Tweedie ajusta muito bem, mesmo para distribuições que não fazem parte da familia, como a t e a slash. Talvez vc possa explorar idéias similares para a weibul, gumbel e similares. http://www.tandfonline.com/doi/abs/10.1080/00949655.2017. 1318876?journalCode=gscs20 http://journals.sagepub.com/doi/abs/10.1177/1471082X17715718 Att Wagner Em 9 de janeiro de 2018 01:11, Wagner Wolff <wwolff(a)usp.br> escreveu: > Olá Wagner > > Gostei da sua ideia sobre a familia de destribuição Tweedie. Além dessa > família você poderia indicar outras que englobem distribuições como por > exemplo, Weibull, Gumbel, Burr ... > > Ab > > 2018-01-08 18:18 GMT-02:00 Wagner Bonat via R-br <r-br(a)listas.c3sl.ufpr.br > >: > >> Quais distribuições vc tem interesse? >> A distribuição Tweedie tem como casos particulares algumas das mais >> populares distribuições como a normal, gamma e inversa normal. >> Assim, ao invés de ajustar um monte de distribuições basta ajustar uma >> que tem as outras como casos particulares. >> >> Att >> Wagner >> >> Em 8 de janeiro de 2018 10:38, Rodrigo Campos via R-br < >> r-br(a)listas.c3sl.ufpr.br> escreveu: >> >>> Bom dia Felipe, >>> Talvez o pacote "fitdistrplus" possa te ajudar! >>> Att. >>> >>> 2018-01-08 10:21 GMT-02:00 Felipe Felix Costa via R-br < >>> r-br(a)listas.c3sl.ufpr.br>: >>> >>>> olá, >>>> >>>> gostaria de saber se existe alguma forma para rodar todas as possíveis >>>> distribuições sobre um determinado conjunto de dados com variáveis >>>> continuas!? >>>> >>>> pois estou procurando e não achei nada pela internet falando dessa >>>> possibilidade. >>>> >>>> meu objetivo é verificar qual melhor a distribuição para os meus dados >>>> e, a partir dessa "melhor", estabelecer um modelo geral com todos os dados >>>> e comparar com modelos individuais de áreas isoladas. >>>> >>>> >>>> HELP >>>> >>>> Att., >>>> -- >>>> *Felipe Felix Costa* >>>> Tel. 96-98122-3077 >>>> >>>> >>>> _______________________________________________ >>>> R-br mailing list >>>> R-br(a)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 Campos >>> >>> >>> >>> _______________________________________________ >>> R-br mailing list >>> R-br(a)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. >>> >> >> >> >> -- >> Wagner Hugo Bonat >> ------------------------------------------------------------ >> ---------------------------------- >> Department of Mathematics and Computer Science (IMADA) >> University of Southern Denmark (SDU) and >> Laboratório de Estatística e Geoinformação (LEG) >> Universidade Federal do Paraná (UFPR) >> >> >> _______________________________________________ >> R-br mailing list >> R-br(a)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. >> > > > > -- > *Wagner Wolff, **PhD* > "*Luiz de Queiroz**" College of Agriculture,* > University of São Paulo > Pádua Dias avenue11 | 13418-900| Piracicaba-SP| Brazil > Phone: +55 19 982385582 <+55%2019%2098238-5582> > http://orcid.org/0000-0003-3426-308X > https://github.com/wwolff7 > http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4463141A1 > -- Wagner Hugo Bonat ---------------------------------------------------------------------------------------------- Department of Mathematics and Computer Science (IMADA) University of Southern Denmark (SDU) and Laboratório de Estatística e Geoinformação (LEG) Universidade Federal do Paraná (UFPR) -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20180109/e6dfd363/atta…> ------------------------------ Subject: Legenda do Digest _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br ------------------------------ Fim da Digest R-br, volume 85, assunto 6 ****************************************

3 2

Distribuição de probabilidade
by Felipe Felix Costa 10 Jan '18

10 Jan '18

olá, gostaria de saber se existe alguma forma para rodar todas as possíveis distribuições sobre um determinado conjunto de dados com variáveis continuas!? pois estou procurando e não achei nada pela internet falando dessa possibilidade. meu objetivo é verificar qual melhor a distribuição para os meus dados e, a partir dessa "melhor", estabelecer um modelo geral com todos os dados e comparar com modelos individuais de áreas isoladas. HELP Att., -- *Felipe Felix Costa* Tel. 96-98122-3077

5 8

Re: [R-br] Salvando Gráficos num Loop
by sznelwar＠uol.com.br 06 Jan '18

06 Jan '18

1 0

Salvando Gráficos num Loop
by Yury Duarte 05 Jan '18

05 Jan '18

Bom dia colegas listeiros! Estou tendo um pequeno problema para automatizar a criação e salvamento de gráficos temáticos utilizando um loop. Quando rodo o script lendo um arquivo de cada vez (fora do loop), o gráfico é salvo sem prejuízo na pasta de destino, entretanto, quando peço para rodar o script dentro do loop, ele salva os gráficos em branco. Gostaria de entender o porque dessa peculiaridade, já que programa não me retorna nenhum aviso de falha ou problema. O código utilizado segue abaixo e também anexei dois arquivos distintos dos bancos de dados que estou trabalhando, caso alguém queira/precise rodar o script. Desde já, agradeço pela ajuda de todos! rm(list=ls()) library(hydroTSM) library(lattice) library(stringr) #----------Verificando Diretorio----------# setwd('C:\\Users\\Yury\\Desktop\\Fernandes') length(dir('C:\\Users\\Yury\\Desktop\\Fernandes')) a = grep('.txt', dir(getwd()), value = T) a = str_sub(a, end = 5) #----------Loop de todas as EM----------# for(i in 1:length(a)){ #----------Selecao Arquivo INMET----------# b = read.table(paste0('C:\\Users\\Yury\\Desktop\\Fernandes\\', a[i], '.txt'), header=TRUE) b.cor = data.frame(b$Data, b$Precipitacao, b$TempMaxima, b$TempMinima, b$Insolacao, b$Umidade.Relativa.Media) colnames(b.cor) = c("Data","RAIN","TMAX","TMIN","n","RHUM"); str(b.cor) b.cor$Data = as.Date(levels(b$Data))[b$Data] #----------Criar Vetor Completo de Datas----------# vetor.data = seq(as.Date('01/07/1961', "%d/%m/%Y"), as.Date('30/06/2015', "%d/%m/%Y"), 1) Data = data.frame(Data=vetor.data) #----------Corrigir Coluna de Datas com Falha----------# duff = merge(b.cor, Data, all = TRUE) str(duff) #----------Criando Graficos----------# pdf(paste0('C:\\Users\\Yury\\Desktop\\Fernandes\\Graphs\\', a[i], '.pdf'), width = 20, height = 14) matrixplot(dwi(duff, var.type="Days"),ColorRamp="Days", main = "Estação_b") dev.off() } Yury Duarte Engenheiro Agrônomo - ESALQ/USP

2 2

Dúvida em componentes principais - gráfico
by Ana Julia Righetto 04 Jan '18

04 Jan '18

Prezad@s, Realizei análise de componentes principais e fiquei com 3 componentes. Estou com dificuldade para fazer um gráfico similar ao biplot (quando tem apenas 2CP). Alguém poderia me ajudar? Obrigada desde já -- Ana Julia Righetto Estaticista - Universidade Estadual Paulista "Júlio de Mesquita Filho" - UNESP Mestre em Ciências (Estatística e Experimentação Agronômica) - ESALQ/USP - Piracicaba - SP Doutora em Ciências (Estatística e Experimentação Agronômica) - ESALQ/USP - Piracicaba - SP Telefone:(43) 3422-8700 E-mail: ajrighetto(a)gmail.com ajrighetto(a)usp.br Skype: anajrighetto Lattes: http://lattes.cnpq.br/5933607596600115

3 2