- R-br - DINF UFPR

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

Re: [R-br] Notação potências de dez (Gráficos)
by Ivan Bezerra Allaman 22 Set '13

22 Set '13

Em 2006 um sabido respondeu no r-help, sciNotation <- function(x, digits = 1) { if (length(x) > 1) { return(append(sciNotation(x[1]), sciNotation(x[-1]))) } if (!x) return(0) exponent <- floor(log10(x)) base <- round(x / 10^exponent, digits) as.expression(substitute(base %.% 10^exponent, list(base = base, exponent = exponent))) } plot(1:1000, axes = FALSE, type = "l", frame.plot = TRUE) axis(1, at = axTicks(1), label = sciNotation(axTicks(1), 1)) Deve ser o que procuras!!! (s,f,p) Allaman

1 0

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

22 Set '13

Não tive problema com as dependências! De fato agora não apareceu mesmo, ou seja, obtive o seguinte: trPaths [1] "\\Tinn-R\\tmp\\" "\\Tinn-R\\tmp\\search.txt" [3] "\\Tinn-R\\tmp\\objects.txt" "\\Tinn-R\\tmp\\file.r" [5] "\\Tinn-R\\tmp\\selection.r" "\\Tinn-R\\tmp\\block.r" [7] "\\Tinn-R\\tmp\\lines.r" "\\Tinn-R\\tmp\\reformat-input.r" [9] "\\Tinn-R\\tmp\\reformat-output.r" Abraços!

1 0

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

22 Set '13

Ok! Agora que entendi que Sys lalala é uma função. Pois bem!! Estou no windows agora e reproduzi o que você postou, e não tive NENHUM PROBLEMA, ou seja, ao rodar trPaths <- paste(paste(as.vector(Sys.getenv('APPDATA')), '\\Tinn-R\2013/09/22 - 12:20:51mp\\', sep=''), c('', 'search.txt', 'objects.txt', 'file.r', 'selection.r', 'block.r', 'lines.r', 'reformat-input.r', 'reformat-output.r'), sep='') obtive trPaths [1] "C:\\Users\\ivan\\AppData\\Roaming\\Tinn-R3/09/22 - 12:20:51mp\\" [2] "C:\\Users\\ivan\\AppData\\Roaming\\Tinn-R3/09/22 - 12:20:51mp\\search.txt" [3] "C:\\Users\\ivan\\AppData\\Roaming\\Tinn-R3/09/22 - 12:20:51mp\\objects.txt" [4] "C:\\Users\\ivan\\AppData\\Roaming\\Tinn-R3/09/22 - 12:20:51mp\\file.r" [5] "C:\\Users\\ivan\\AppData\\Roaming\\Tinn-R3/09/22 - 12:20:51mp\\selection.r" [6] "C:\\Users\\ivan\\AppData\\Roaming\\Tinn-R3/09/22 - 12:20:51mp\\block.r" [7] "C:\\Users\\ivan\\AppData\\Roaming\\Tinn-R3/09/22 - 12:20:51mp\\lines.r" [8] "C:\\Users\\ivan\\AppData\\Roaming\\Tinn-R3/09/22 - 12:20:51mp\\reformat-input.r" [9] "C:\\Users\\ivan\\AppData\\Roaming\\Tinn-R3/09/22 - 12:20:51mp\\reformat-output.r" Logo, Qual é o problema?????????????????????????????????????????

2 2