Setembro 2013 - R-br

Re: [R-br] ajuda normatização de grid
by Eder Comunello 27 Set '13

27 Set '13

Prezado Elias, boa noite! Pra ajudar de fato seria necessário ter acesso aos dados que você está manipulando (ou uma amostra). Mesmo assim vou arriscar um palpite. Você removeu NA's com base somente na coluna 14 (gride[,14]) e por fim, normaliza de 4 a 18 (grid1[,4:18]). É provável que outras colunas além da 14 tenham valores de NA e estejam causando o erro. Para verificar o número de ocorrências de NA's em grid3: length(grid3[is.na(grid3)]) Espero que ajude, -- Éder Comunello <e <comunello.eder(a)gmail.com>comunel(a)gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]

1 0

[OFF] Formatar classe de uma distribuição
by Andre Oliveira 26 Set '13

26 Set '13

Pessoal afim de ficar bem formatado eu uso o símbolo \vdash no LaTeXpara montar as classes de uma distribuição de dados contínuos, porém estou precisando de construir um intervalo incluindo as duas extremidades para a última classe. Alguém conhece algum símbolo em LaTeX que faça isto? Eu usei \vdash \dashv mas isto é um gambiarra feia! obrigado André Oliveira Souza

1 0

Análise de Série Temporal - Modelo
by Edson Lira 26 Set '13

26 Set '13

Caros amigos, estou com uma dúvida, estou analisando uma série de dados de doadores aptos de 2000 a 2011. Observando a série e fazendo o teste das raizes, conclui que a mesma é estacionária sem haver a necessidade da diferenciação. Usando a função auto.arima do pacote forecast, ele me indica um modelo com sazonalidade, porém observando os gráficos de correlação e correlação parcial não identifico o modelo sugerido pelo pacote, alguém poderia me dar uma sugestão? No link abaixo está a saída com os gráficos. http://www.datafilehost.com/d/88912c37 [ ]'s. Edson Lira Estatístico Manaus-Amazonas

3 3

Banco de dados muito grande
by Edson Lira 26 Set '13

26 Set '13

Caros amigos, boa noite, tenho um banco de dados com 10 mb e 176 mil linhas com 11 colunas. Estou usando: tra<-read.xlsx("/Users/edsonlira/Minuza/Dados_geral.xlsx",1,header=T,encoding= "UTF-8") Me mostra o seguinte erro: Erro em .jcall("RJavaTools", "Ljava/lang/Object;", "invokeMethod", cl, : java.lang.OutOfMemoryError: Java heap space O que devo fazer? [ ]'s. Edson Lira Estatístico Manaus-Amazonas

8 15

Divida em qual analise usar para resolver um problema de pesquisa
by Alisson Lucrecio 25 Set '13

25 Set '13

Caro Colegas da r-br, Boa tarde. Estou tentando ajudar um colega de universidade, mas não sei qual a ferramenta mais apropria para resolver esse problema, assim gostaria da ajuda de vocês. Um amigo fez análise 14 elementos de terra rara em 20 fertilizante por dois métodos de abertura, assim fatorial 20 (Fertilizantes) x 2 (Métodos de abertura) com 3 repetições. Ele quer responder 2 peguntas: os fertilizantes são diferentes entre si considerando todos os elementos de terra rara? Os métodos de abertura apresentam o mesmos valor de elemento de terra rara? Se não, em quais fertilizantes, os métodos de abertura superestima ou subestima o valor quando comparado ao outro método. Obrigado Alisson Lucrécio da Costa

1 0

combinar variáveis - PNAD
by Luciane Maria Pilotto 25 Set '13

25 Set '13

Olá pessoal, gostaria de ajuda para combinar variáveis. Estou trabalhando com dados da PNAD 2008 e preciso criar uma nova variável com a combinação de outras 3 variáveis. As variáveis são: 1321 (tem direito a plano de saúde) - com respostas 1,3 e 5 1342 (dá direito à assistencia odontologica) - respostas 2 e 4 1344 (plano exclusivo para odontologia) - respostas 2 e 4 Variável a ser criada (tipo de plano): Sem plano: se 1321=5 Plano com odontologia: se 1321=1 ou 1321=3 e 1342=2 e 1344=2; 1321=1 ou 1321=3 e 1342=2 e 1344=4 Plano sem odontologia: 1321=1 ou 1321=3 e 1342=4 e 1344=4 Atenciosamente, Luciane Maria Pilotto

1 0

Combinar arrays de diferentes comprimentos em data.frame
by Thiago V. dos Santos 24 Set '13

24 Set '13

Olá pessoal, Como resultado da leitura de algumas images de satélite, tenho o seguinte array no R > dim(ncdata) [1] 18 15 16 2555 > length(ncdata) [1] 11037600 onde as dimensões representam 18 é a quantidade de longitudes, 15 é a quantidade de latitudes, 16 é a quantidade de tipos de vegetação (floresta tropical, cerrado etc), 2555 é a quantidade de datas em que as imagens foram feitas. Por outro lado, tenho os seguintes arrays: > length(nclon) [1] 18 > length(nclat) [1] 15 > length(pft) [1] 16 > length(nctime) [1] 2555 E, em suma, o que eu quero fazer é criar um data.frame com a seguinte estrutura: nclon1 nclat1 pft1 ncdate1 (primeiro valor do ncdata) nclon1 nclat1 pft1 ncddate2 (segundo valor do ncdata) . . . nclon18 nclat15 pft16 nctime2555 (último valor do ncdata) Como a quantidade de linhas é muito grande, o ideal seria até usar alguma técnica que permita paralelização usando os cores do pc. Alguém conhece uma maneira (realmente) eficiente de fazer isso? Saudações, -- Thiago V. dos Santos PhD student Land and Aé qtmospheric Science University of Minnesota http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/ind… Phone: (612) 323 9898

2 2

Re: [R-br] r baixando dados inmet
by Éder Comunello 24 Set '13

24 Set '13

Bom dia! Também tenho interesse no tópico, pois costumava baixar os arquivos um a um, diretamente no site! :D Examinando o arquivo da lista, tem uma discussão sobre o RCurl em: <http://r-br.2285057.n4.nabble.com/R-br-RCurl-td4659610.html> Seguindo a discussão de Daniel Marcelino e Raphael Saldanha, fiz umas adaptações no código postado por eles para o caso em questão: ### <begin> library(RCurl) # Links para login e dados myURL1 <- "http://www.inmet.gov.br/projetos/rede/pesquisa/inicio.php" myURL2 <- " http://www.inmet.gov.br/projetos/rede/pesquisa/gera_serie_txt.php?&mRelEsta…,,,,, " ### como é só um teste, encurtei a série pesquisada (mRelDtInicio=01/01/2012) # Dados para acesso myParams=list( mCod="seuemail", ### alterar! mSenha="suasenha", ### alterar! btnProcesso = " Acessar ") myCurl <- getCurlHandle() curlSetOpt(cookiejar="cookies.txt", useragent="Mozilla/5.0", followlocation=TRUE, curl=myCurl) login <- postForm(myURL1, .params=myParams, curl=myCurl) dados <- getURLContent(myURL2, curl=myCurl) head(strwrap(dados, 70), 30) ### <end> O código acima faz o login e baixa os dados de um único exemplo (myURL2). Tá testado e funcionando!!! Logicamente requer cadastro no BDMEP... Agora falta entrar com a limpeza e importação do código armazenado na variável dados e automatizar a geração de links para dados (basicamente mudar mRelEstacao, mRelDtInicio e .mRelDtFim em cada link). Mas acredito que já seja um começo! :D Éder Comunello ### Uma saída parcial do conteúdo armazenado em dados: > head(strwrap(dados, 70), 30) [1] "<html><meta http-equiv=\"Content-Type\" content=\"text/html;" [2] "charset=iso-8859-1\"> <head> <title>Estação - 83595</title> \r <script src=\"facebox/js/jquery.js\"" [4] "type=\"text/javascript\"></script>\r <link" [5] "href=\"facebox/styles/facebox.css\" media=\"screen\" rel=\"stylesheet\"" [6] "type=\"text/css\" >\r \r <script src=\"facebox/js/facebox.js\"" [7] "type=\"text/javascript\"></script>\r <script type=\"text/javascript\">\r" [8] "jQuery(document).ready(function($) {\r $('a[rel*=facebox]').facebox({\r" [9] "loading_image : 'loading.gif',\r close_image : 'closelabel.gif'})\r" [10] "$('area[rel*=facebox]').facebox({\r loading_image : 'loading.gif',\r" [11] "close_image : 'closelabel.gif'})\r })\r </script>\r <script" [12] "type=\"text/javascript\">\r function showFoto(src) {\r //criando um link" [13] "falso, para o LightBox captar as informações\r var a_false =" [14] "document.createElement(\"A\");\r a_false.rel = \"lightbox\";\r a_false.href" [15] "= src;\r //chamando a função do lightbox responsável pelo início das" [16] "ações\r myLightbox.start(a_false);\r //anulando o link\r a_false =" [17] "null;\r }\r </script>\r </head> <pre>-------------------- BDMEP - INMET" [18] "-------------------- Estação : AIMORES - MG (OMM: 83595) Latitude" [19] "(graus) : -19.49 Longitude (graus) : -41.07 Altitude (metros): 82.74" [20] "Estação Operante Inicio de operação: 01/06/1972 Periodo solicitado" [21] "dos dados: 01/01/2012 a 01/01/2013 Os dados listados abaixo são os" [22] "que encontram-se digitados no BDMEP -------------------- Obs.: Os" [23] "dados aparecem separados por ; (ponto e vírgula) no formato txt." [24] "Para o formato planilha XLS, <A HREF=\"instrucao.html\" TARGET=\"_top\"" [25] "rel=\"facebox\">siga as instruções</a> --------------------" [26] "Estacao;Data;Hora;TempBulboSeco;UmidadeRelativa;PressaoAtmEstacao;DirecaoVento;VelocidadeVentoNebulosidade;" [27] "83595;01/01/2012;0000;25.2;83;1004.8;14;1;8;" [28] "83595;01/01/2012;1200;26;84;1004.9;9;1;10;" [29] "83595;01/01/2012;1800;27;91;1000;9;1;10;" [30] "83595;02/01/2012;0000;25;92;1003;14;1;10;"

2 2

r baixando dados inmet
by Alisson Lucrecio 24 Set '13

24 Set '13

Caro Colegas, Boa noite. Eu estou tentando baixar uma série de dados meteorológico do inmet, assim preciso automatizar o processo, alguém sabe como o r poderia ler um arquivo dessa forma e salva-lo em .txt? O dado é proveniente do site do inmet http://www.inmet.gov.br/projetos/rede/pesquisa/gera_serie_txt.php &mRelEstacao=83595&btnProcesso=serie&mRelDtInicio=01/01/1961&mRelDtFim=01/01/2013&mAtributos=1,,,,1,1,,1,1,,,1,,,,, -------------------- BDMEP - INMET -------------------- Estação : AIMORES - MG (OMM: 83595) Latitude (graus) : -19.49 Longitude (graus) : -41.07 Altitude (metros): 82.74 Estação Operante Inicio de operação: 01/06/1972 Periodo solicitado dos dados: 01/01/1961 a 01/01/2013 Os dados listados abaixo são os que encontram-se digitados no BDMEP -------------------- Obs.: Os dados aparecem separados por ; (ponto e vírgula) no formato txt. Para o formato planilha XLS, siga as instruções -------------------- Estacao;Data;Hora;TempBulboSeco;UmidadeRelativa;PressaoAtmEstacao;DirecaoVento;VelocidadeVentoNebulosidade; 83595;01/06/1972;1200;19.8;84;1011.6;36;1;8; 83595;01/06/1972;1800;27.2;47;1007.8;9;1;9; 83595;02/06/1972;0000;23;67;1009.5;14;3;10; 83595;02/06/1972;1200;22;76;1011;32;1;4; 83595;02/06/1972;1800;29;45;1006.4;32;1;4; 83595;03/06/1972;0000;19;82;1009;14;1;10; 83595;03/06/1972;1200;22.2;78;1011.4;32;1;5; 83595;03/06/1972;1800;30.8;47;1005.6;32;1;8; 83595;04/06/1972;0000;20.2;85;1008.6;23;1;10; 83595;04/06/1972;1200;23.8;78;1010.7;36;1;8; 83595;04/06/1972;1800;30.8;47;1007.4;32;1;6; 83595;05/06/1972;0000;25;81;1008.8;14;5;10; 83595;05/06/1972;1200;22;83;1010;32;1;4; 83595;05/06/1972;1800;29;57;1007.2;14;1;4; 83595;06/06/1972;0000;24;75;1008.6;14;3;10; 83595;06/06/1972;1200;21;81;1010.2;32;1;3; 83595;06/06/1972;1800;29;54;1005.8;36;1;8; 83595;07/06/1972;0000;24.5;73;1008.4;14;1;10; 83595;07/06/1972;1200;25.4;71;1009.6;9;1;7; 83595;07/06/1972;1800;25.8;71;1007;5;3;9; 83595;08/06/1972;0000;23;76;1008.7;14;1;10; 83595;08/06/1972;1200;23.5;80;1012.3;32;1;5; 83595;08/06/1972;1800;30;52;1007.2;32;1;6; 83595;09/06/1972;0000;23;76;1009.2;14;1;10; 83595;09/06/1972;1200;24.5;77;1010.3;5;1;5; 83595;09/06/1972;1800;29;53;1006.3;5;1;6; 83595;10/06/1972;0000;24;80;1008.7;14;3;10; 83595;10/06/1972;1200;23;82;1010.8;0;0;3; 83595;10/06/1972;1800;29.8;50;1005.2;5;1;5; 83595;11/06/1972;0000;24;77;1008.4;14;1;10; 83595;11/06/1972;1200;23;80;1010;5;1;3; 83595;11/06/1972;1800;30;50;1007.1;23;1;4; 83595;12/06/1972;0000;21;77;1009;18;1;10; 83595;12/06/1972;1200;20;87;1010.4;23;1;10; Obrigado. Alisson Lucrécio da Costa

3 3

Re: [R-br] Vetor de caracter como objeto em um pacote
by Ivan Bezerra Allaman 24 Set '13

24 Set '13

Complexo!! Já tentou system(paste(....)) ou cat(..., file = Sys.getenv('APPDATA')) ?? O problema de lhe ajudar, é que não tem como gerar um CRM (eu acho) para testar aqui!! Aí fica mais difícil só na imaginção!! Abraços!!

3 2