Agosto 2016 - R-br

Geração de series sintéticas de precipitação no R
by Alisson Lucrécio 22 Ago '16

22 Ago '16

Caro Colegas, Bom dia. Algum colega do tem um script para geração de séries sintéticas de precipitação no R. Baseado nessa metodologia. Segue o link. https://www.google.com.br/url?sa=t&rct=j&q=&esrc=s&source=web&cd=1&cad=rja&… Muito obrigado. -- Alisson Lucrecio da Costa

1 0

previsão dinâmica
by Edimeire Alexandra Pinto 19 Ago '16

19 Ago '16

OI gente. Preciso de uma ajudinha. Sei que parece bobo, mas nãoconsigo estender a janela de previsões sequenciais com largura de , com largura fixa de tamanho 8.Criar um modelo com sample de 1:70 e prever 71 até 78,depois reestimar a equação com sample de 1:72 e previsão de 73:80, reestimar omodelo de 1:73 e previsão de 74:83 eassim por adiante até fechar o range em 84, última observação.Vejamos um data frame já existente no R como nome de Canadaque são séries que vão do primeiro trimestre de 1980 ao último trimestre de2004:require("dynlm") library("vars") data("Canada") summary(Canada) class(Canada)colnames(Canada) suponha que eu definaque minha amostra para estimar o modelo com sequencia em j que vai até 10,assim teríamos um modelo estimado de 01/1980 a 3/1996, depois outros modelos estimado de 01/1980a 4/1996, depoismais outro modelo estimado de 01/1980 a 1/1997, assim por adiante até 4/1998, ou seja, paraem 4/1998. No R seria mais ou menos assim:eq<-list() for(j in1:10){ eq(i)<-dynlm(e~L(e,1)+prod+rw+u,data=window(Canada,start=c(1980,1),end=c(1996,2),freq=4,extend=j)) o problema é que a previsão tem de ser sequencial tanto nostart quanto no end e com tamanho de fixo de 8, e também com uma janela que vaiindo sequencialmente em j de 1 a 18. Por exemplo, a previsão teria de serassim: a primeira de 04/1996 a 4/1998 (observe que o primeiro modelo termina em3/1996, conforme falei antes), depois 01/1997 a 2/1999 (observe que o primeiromodelo termina em 1/1997, conforme falei antes), depois faria outra previsão de2/1997 a 4/1999 e assim por adiante até chegar na última observação que é04/2000. Note que aumentamos sempre de 8+j em que j vai de 1 a 18. No R, não consigo desenvolver, pois seria algo mais ou menosassim: ## sei que está errado, mas é só para ter ideiafor( i in1:18){ predict(eq(i), newdata=window(Canada,start=c(1996,4+i),end=c(1998,4+i+8),freq=4,extend=i) ## 8 é o tamanho fixo da janela de previsão Se tiver ficado confuso, desculpe-me, mas quis apenasmostrar que quero um previsão dinâmica e não estática, tendo e, vista que aequação que sempre estimada à medida que o tamanho da amostra vai variando.Qualquer dica é válida!Obrigada, gente.

3 3

Re: [R-br] Digest R-br, volume 68, assunto 17
by ana paula coelho madeira 19 Ago '16

19 Ago '16

Prezado Alexandro, Os comandos em questão tentam apenas a normalidade. Tenho interesse em ver ficar a aderência de uma distribuição qualquer. Agradeço > From: r-br-request(a)listas.c3sl.ufpr.br > Subject: Digest R-br, volume 68, assunto 17 > To: r-br(a)listas.c3sl.ufpr.br > Date: Wed, 17 Aug 2016 12:00:02 -0300 > > 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 dinâmica (Edimeire Alexandra Pinto) > 2. Re: Anderson Darling (Alexandro (Yahoo)) > 3. OFF TOPIC - Banco de dados (Andre Oliveira) > 4. Re: OFF TOPIC - Banco de dados (luc.souza(a)usp.br) > 5. RES: OFF TOPIC - Banco de dados (Leonard Mendonça de Assis) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Tue, 16 Aug 2016 21:53:38 +0000 (UTC) > From: Edimeire Alexandra Pinto <economatistica(a)yahoo.com.br> > To: a lista Brasileira oficial de discussão do programa R. > <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] previsão dinâmica > Message-ID: > <228175548.15741202.1471384418129.JavaMail.yahoo(a)mail.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > OI gente. > > Preciso de uma ajudinha. Sei que parece bobo, mas nãoconsigo estender a janela de previsões sequenciais com largura de , com largura fixa de tamanho 8.Criar um modelo com sample de 1:70 e prever 71 até 78,depois reestimar a equação com sample de 1:72 e previsão de 73:80, reestimar omodelo de 1:73 e previsão de 74:83 eassim por adiante até fechar o range em 84, última observação.Vejamos um data frame já existente no R como nome de Canadaque são séries que vão do primeiro trimestre de 1980 ao último trimestre de2004:require("dynlm") > > library("vars") > > data("Canada") > > summary(Canada) > > class(Canada)colnames(Canada) > > > > suponha que eu definaque minha amostra para estimar o modelo com sequencia em j que vai até 10,assim teríamos um modelo estimado de 01/1980 a 3/1996, depois outros modelos estimado de 01/1980a 4/1996, depoismais outro modelo estimado de 01/1980 a 1/1997, assim por adiante até 4/1998, ou seja, paraem 4/1998. No R seria mais ou menos assim:eq<-list() > > for(j in1:10){ > > eq(i)<-dynlm(e~L(e,1)+prod+rw+u,data=window(Canada,start=c(1980,1),end=c(1996,2),freq=4,extend=j)) > > o problema é que a previsão tem de ser sequencial tanto nostart quanto no end e com tamanho de fixo de 8, e também com uma janela que vaiindo sequencialmente em j de 1 a 18. Por exemplo, a previsão teria de serassim: a primeira de 04/1996 a 4/1998 (observe que o primeiro modelo termina em3/1996, conforme falei antes), depois 01/1997 a 2/1999 (observe que o primeiromodelo termina em 1/1997, conforme falei antes), depois faria outra previsão de2/1997 a 4/1999 e assim por adiante até chegar na última observação que é04/2000. Note que aumentamos sempre de 8+j em que j vai de 1 a 18. No R, não consigo desenvolver, pois seria algo mais ou menosassim: > > ## sei que está errado, mas é só para ter ideiafor( i in1:18){ > > predict(eq(i), newdata=window(Canada,start=c(1996,4+i),end=c(1998,4+i+8),freq=4,extend=i) ## 8 é o tamanho fixo da janela de previsão > > Se tiver ficado confuso, desculpe-me, mas quis apenasmostrar que quero um previsão dinâmica e não estática, tendo e, vista que aequação que sempre estimada à medida que o tamanho da amostra vai variando.Qualquer dica é válida!Obrigada, gente. > > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160816/39b43423/atta…> > > ------------------------------ > > Message: 2 > Date: Tue, 16 Aug 2016 22:04:57 +0000 (UTC) > From: "Alexandro (Yahoo)" <vl.alexandro(a)yahoo.com.br> > To: ana paula coelho madeira <apcmadeira(a)hotmail.com>, a lista > Brasileira oficial de discussão do programa R. > <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Anderson Darling > Message-ID: > <1102947158.17317609.1471385097903.JavaMail.yahoo(a)mail.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > #Segue um exemplo!require(ADGofTest)#?ad.test > set.seed( 123 )x <- runif( 100 ) > ad.test( x )$p.valuead.test( x, pnorm, 0, 1 )$p.value #Espero ter ajudado! > Alexandro Vieira Lopes > > > > De: ana paula coelho madeira via R-br <r-br(a)listas.c3sl.ufpr.br> > Para: lista R <r-br(a)listas.c3sl.ufpr.br> > Enviadas: Segunda-feira, 15 de Agosto de 2016 15:28 > Assunto: [R-br] Anderson Darling > > Boa tarde. > Tenho um conjunto de dados nos quais estou trabalhando com diferentes fdp's. Gostaria de aplicar diferentes testes de aderência, entre eles o Anderson Darling. No entanto, consegui apenas o ks.test. No pacote Nortest encontrei o ad.test para testar a normalidade, que não é o caso. Alguém tem algum material para me indicar? > Desde já agradeço. > Ana Paula > _______________________________________________ > 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� c�igo m�imo reproduz�el. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160816/268f052c/atta…> > > ------------------------------ > > Message: 3 > Date: Wed, 17 Aug 2016 03:21:29 +0000 (UTC) > From: Andre Oliveira <andreolsouza(a)yahoo.com.br> > To: a lista Brasileira oficial de discussão do programa R. > <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] OFF TOPIC - Banco de dados > Message-ID: > <811409596.15775329.1471404089455.JavaMail.yahoo(a)mail.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > Boa noite,tenho um projeto com aluno de ensino médio sobre produzir interface em shiny para visualizar dados. É uma feira onde o tema é alimentação e pensei em explorar banco de dados de produção de milho, soja ou outro produtos por hectare ao longo de anos no Brasil, via shiny. Alguém saberia dizer onde posso encontrar banco de dados disponíveis com varias culturas e vários anos? Este aluno é do técnico em informática e vem trabalhando com estatística onde há uma integração entre as áreas. > obrigado > André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160817/1a282242/atta…> > > ------------------------------ > > Message: 4 > Date: Wed, 17 Aug 2016 06:47:50 -0300 (BRT) > From: luc.souza(a)usp.br > To: Andre Oliveira <andreolsouza(a)yahoo.com.br>, "a lista Brasileira > oficial de discussão do programa R." <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] OFF TOPIC - Banco de dados > Message-ID: <3E69A4AD-23EB-49AC-B77D-2914AD7D210C(a)usp.br> > Content-Type: text/plain; charset="utf-8" > > Na conab tem uma série histórica de vários produtos e no ibge Tb. > > Enviado do meu iPhone > > > Em 17 de ago de 2016, às 00:21, Andre Oliveira via R-br <r-br(a)listas.c3sl.ufpr.br> escreveu: > > > > Boa noite, > > tenho um projeto com aluno de ensino médio sobre produzir interface em shiny para visualizar dados. É uma feira onde o tema é alimentação e pensei em explorar banco de dados de produção de milho, soja ou outro produtos por hectare ao longo de anos no Brasil, via shiny. Alguém saberia dizer onde posso encontrar banco de dados disponíveis com varias culturas e vários anos? Este aluno é do técnico em informática e vem trabalhando com estatística onde há uma integração entre as áreas. > > obrigado > > > > André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES > > _______________________________________________ > > 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. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160817/6a558ec2/atta…> > > ------------------------------ > > Message: 5 > Date: Wed, 17 Aug 2016 09:37:57 -0300 > From: Leonard Mendonça de Assis <assis.leonard(a)gmail.com> > To: 'a lista Brasileira oficial de discussão do programa R.' > <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] RES: OFF TOPIC - Banco de dados > Message-ID: <00ad01d1f884$30927090$91b751b0$(a)gmail.com> > Content-Type: text/plain; charset="utf-8" > > André, eu sei que o CEPEA/USP disponibiliza séries temporais com a cotação de diversos produtos agrícolas. Área plantada eu desconheço algum órgão que possa ter isso. Uma tentativa seria ministério da Agricultura ou censo agropecuário do ibge > > > > De: R-br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de Andre Oliveira via R-br > Enviada em: quarta-feira, 17 de agosto de 2016 00:21 > Para: a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> > Assunto: [R-br] OFF TOPIC - Banco de dados > > > > Boa noite, > > tenho um projeto com aluno de ensino médio sobre produzir interface em shiny para visualizar dados. É uma feira onde o tema é alimentação e pensei em explorar banco de dados de produção de milho, soja ou outro produtos por hectare ao longo de anos no Brasil, via shiny. Alguém saberia dizer onde posso encontrar banco de dados disponíveis com varias culturas e vários anos? Este aluno é do técnico em informática e vem trabalhando com estatística onde há uma integração entre as áreas. > > > > obrigado > > > > > > André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160817/bcb2784f/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 68, assunto 17 > *****************************************

2 1

Re: [R-br] O R não abre no linux mint
by Sergio Pinto 18 Ago '16

18 Ago '16

Prezado Marcelo Obrigado pela ajuda. Infelizmente o problema persiste. Continua aparecendo a mesma saída no terminal. Att Sergio Medeiros Pinto Doutorando em Ciências Médicas UERJ Prof. de Educação Fisica Em 17/08/2016 21:50, "Marcelo Laia" <marcelolaia(a)gmail.com> escreveu: http://stackoverflow.com/questions/20967256/error-on- package-and-r-paths-and-installations On 17/08/16 at 08:53pm, Sergio Pinto via R-br wrote: > Prezados > tenho um computador 32 bits com linux mint e instalei o pacote r-base. > O ícone do R foi criado mas não abre o programa. > No terminal aparece a seguinte mensagem: > > sergio@sergio-Infoway ~ $ R > /usr/local/bin/R: line 236: /usr/local/lib/R/etc/ldpaths: No such file or > directory > /usr/local/bin/R: line 259: /usr/local/lib/R/bin/exec/R: No such file or > directory > > > O meu sistema é > Linux sergio-Infoway 3.13.0-39-generic #66-Ubuntu SMP Tue Oct 28 13:31:23 > UTC 2014 i686 i686 i686 GNU/Linux > > Poderiam me ajudar? > > Desde já agradeço a atenção > > *Sergio Medeiros Pinto* > Doutorando em Ciências Médicas UERJ > Prof. de Educação Física > _______________________________________________ > 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. -- Marcelo

2 1

Girar rótulo da lattice
by Andre Oliveira 18 Ago '16

18 Ago '16

Olá,como posso girar os rótulos da lattice? par(las=1) não deu certo! require(lattice) par(las=1) xyplot(GP~factor(Ano)|S,pch=19, ylab="Ganho de peso (g)",xlab="", main="Ganho de peso por sexo") André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

1 0

O R não abre no linux mint
by Sergio Pinto 18 Ago '16

18 Ago '16

Prezados tenho um computador 32 bits com linux mint e instalei o pacote r-base. O ícone do R foi criado mas não abre o programa. No terminal aparece a seguinte mensagem: sergio@sergio-Infoway ~ $ R /usr/local/bin/R: line 236: /usr/local/lib/R/etc/ldpaths: No such file or directory /usr/local/bin/R: line 259: /usr/local/lib/R/bin/exec/R: No such file or directory O meu sistema é Linux sergio-Infoway 3.13.0-39-generic #66-Ubuntu SMP Tue Oct 28 13:31:23 UTC 2014 i686 i686 i686 GNU/Linux Poderiam me ajudar? Desde já agradeço a atenção *Sergio Medeiros Pinto* Doutorando em Ciências Médicas UERJ Prof. de Educação Física

2 1

OFF TOPIC - Banco de dados
by Andre Oliveira 17 Ago '16

17 Ago '16

Boa noite,tenho um projeto com aluno de ensino médio sobre produzir interface em shiny para visualizar dados. É uma feira onde o tema é alimentação e pensei em explorar banco de dados de produção de milho, soja ou outro produtos por hectare ao longo de anos no Brasil, via shiny. Alguém saberia dizer onde posso encontrar banco de dados disponíveis com varias culturas e vários anos? Este aluno é do técnico em informática e vem trabalhando com estatística onde há uma integração entre as áreas. obrigado André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

5 4

Anderson Darling
by ana paula coelho madeira 16 Ago '16

16 Ago '16

Boa tarde. Tenho um conjunto de dados nos quais estou trabalhando com diferentes fdp's. Gostaria de aplicar diferentes testes de aderência, entre eles o Anderson Darling. No entanto, consegui apenas o ks.test. No pacote Nortest encontrei o ad.test para testar a normalidade, que não é o caso. Alguém tem algum material para me indicar? Desde já agradeço. Ana Paula

2 1

Pacote Forecast Error
by Analia Macedo 15 Ago '16

15 Ago '16

Olá pessoal, um amigo me passou o pacote forecast por email para que eu fizesse o download mas quando tento instalar o pacote recebo essa mensagem > install.packages("~/TCC/forecast_7.1.tar.gz", repos = NULL, type = "source") Installing package into ‘C:/Users/anali_000/Documents/R/win-library/3.2’ (as ‘lib’ is unspecified) * installing *source* package 'forecast' ... ** package 'forecast' successfully unpacked and MD5 sums checked ** libs Warning: running command 'make -f "Makevars.win" -f "C:/PROGRA~1/R/R-32~1.0/etc/x64/Makeconf" -f "C:/PROGRA~1/R/R-32~1.0/share/make/winshlib.mk" SHLIB_LDFLAGS='$(SHLIB_CXXLDFLAGS)' SHLIB_LD='$(SHLIB_CXXLD)' SHLIB="forecast.dll" WIN=64 TCLBIN=64 OBJECTS="calcBATS.o calcTBATS.o etsTargetFunction.o etsTargetFunctionWrapper.o etscalc.o etspolyroot.o makeBATSMatrices.o makeTBATSMatrices.o updateMatrices.o updateTBATSMatrices.o"' had status 127 ERROR: compilation failed for package 'forecast' * removing 'C:/Users/anali_000/Documents/R/win-library/3.2/forecast' Warning in install.packages : running command '"C:/PROGRA~1/R/R-32~1.0/bin/x64/R" CMD INSTALL -l "C:\Users\anali_000\Documents\R\win-library\3.2" "C:/Users/anali_000/Documents/TCC/forecast_7.1.tar.gz"' had status 1 Warning in install.packages : installation of package ‘C:/Users/anali_000/Documents/TCC/forecast_7.1.tar.gz’ had non-zero exit status Será que algum de vocês poderia me ajudar? Em anexo to mandando o zip do arquivo. Meu amigo disse que no dele funcionou sem problemas a instalação do pacote.. Preciso disso pro meu TCC :/ -- *Att. * *Anália Araujo Macedo*

3 3

Recortar .shp para plotar mapa
by Andre Oliveira 15 Ago '16

15 Ago '16

Boa noite,como posso plotar este mapa de forma que não fique tão pequeno? Tem uma borda desnecessária! download.file("ftp://geoftp.ibge.gov.br/organizacao_do_territorio/malhas_territoriais/malhas_municipais/municipio_2014/ES/es_municipios.zip",destfile="es_municipios.zip") unzip("estados_2010.zip") mapaUF = readShapePoly("32MUE250GC_SIR.shp") plot(mapaUF,axes=T) André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

2 1