- R-br - DINF UFPR

Re: [R-br] Criação de função fortran em pacote (base)
by Junior Beleti 28 Fev '12

28 Fev '12

Olá Benilton, tentando passar um exemplo reproduzivel, o experimento em questão executa um script para a obtenção dos valores preditos utilizando os dados do pivo, por meio do pacote geoComp. Nesse script, foram relatados problemas quanto ao tempo de processamento. A função responsável por tais "problemas" foi a função volta.cokri. Dentro dessa função identifiquei que o elevado tempo de processamento se deve à chamada a função eigen. Na eigen, constatei o processamento de grandes matrizes com valores exponenciais. Como minha pesquisa tem foco no processamento paralelo, estou realizando experimentos observando formas de se decompor tais matrizes de forma a realizar uma execução paralela. Em uma das execuções, crio as matrizes m1,m2,m3 e m4, como segue m1 <- mat.cokri[[2]][1:2601,1:2601] m2 <- mat.cokri[[2]][2602:5202,1:2601] m3 <- mat.cokri[[2]][1:2601,2602:5202] m4 <- mat.cokri[[2]][2602:5202,2602:5202] sendo essas sub-matrizes da mat.cokri[[2]] que é uma matriz (5202 x 5202). O que tinha em mente seria a execução paralela de tais matrizes por meio da função mclapply. Nesse experimento, fiz o seguinte: mclapply(1:4, function(i) { g1 <- mvrnorm(n=1, mat.cokri[[1]][1:2601],m1) g2 <- mvrnorm(n=1, mat.cokri[[1]][1:2601],m2) g3 <- mvrnorm(n=1, mat.cokri[[1]][2602:5202],m3) g4 <- mvrnorm(n=1, mat.cokri[[1]][2602:5202],m4) },mc.cores=4 ) Esse seria o procedimento correto? Claro, este sendo um exemplo hipotético, visto que a execução real faz um loop com 1000 iterações. Carlos Beleti.

3 3

Re: [R-br] Criação de função fortran em pacote (base)
by Junior Beleti 28 Fev '12

28 Fev '12

Olá Benilton, eu estou confuso quanto a interpretação da função mclapply. Penso que as quatro operações (g1,g2,g3,g4) são todas executadas sequencialmente, uma após a outra, dentro de um mesmo núcleo. Ou seja, cada iteração do mclapply (que contém as 4 operações) vai para um mesmo núcleo. É isso? Se for isso, então o seguinte comando: -------------------------------------------------------------------------------- mclapply(1:1, function(i) { g1 <- mvrnorm(n=1, mat.cokri[[1]][1:2601],m1) g2 <- mvrnorm(n=1, mat.cokri[[1]][1:2601],m2) g3 <- mvrnorm(n=1, mat.cokri[[1]][2602:5202],m3) g4 <- mvrnorm(n=1, mat.cokri[[1]][2602:5202],m4) },mc.cores=4 ) -------------------------------------------------------------------------------- é apenas uma execução sequencial de g1, depois g2, depois g3 e por último g4? Agora, o seguinte comando: -------------------------------------------------------------------------------- mclapply(1:2, function(i) { g1 <- mvrnorm(n=1, mat.cokri[[1]][1:2601],m1) g2 <- mvrnorm(n=1, mat.cokri[[1]][1:2601],m2) g3 <- mvrnorm(n=1, mat.cokri[[1]][2602:5202],m3) g4 <- mvrnorm(n=1, mat.cokri[[1]][2602:5202],m4) },mc.cores=4 ) -------------------------------------------------------------------------------- executa primeiramente g1, depois g2, depois g3 e por fim g4 em um núcleo e paralelamente faz o mesmo (executa g1, depois g2, depois g3 e por fim g4) em um segundo núcleo. Ou seja, executa 2 iterações de um loop em paralelo. Está correto? Estou preparando um código reproduzível e estarei enviando em seguida. Muito obrigado por sua atenção Benilton. Carlos.

2 1

Ajuda com datas
by Daniel C Bezerra 28 Fev '12

28 Fev '12

Amigos do R, Estou avaliando a frequencia de ocorrência de uma condição entre os meses de Ago 2011 e Fev 2012. A minha dúvida é simples: se eu peço para formatar as datas na forma de mes e ano e solicito uma tabela de ocorrências o R coloca o output com base na ordem alfabética dos meses. Como fazer para colocar na ordem cronológica na tabela? table(format(data.int,"%b %Y")) Aug 2011 Dec 2011 Feb 2012 Jan 2012 Nov 2011 Oct 2011 Sep 2011 4 5 5 4 12 10 13 Alguma sugestão? Muito obrigado, D

4 7

Rosa dos Ventos
by Jodavid Ferreira 27 Fev '12

27 Fev '12

Alguém sabe como colocar Rosa dos ventos nos mapas, direto do R? -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pela comissão de Comunicação e do Ministério Jovem do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig http://jodavidferreira.com

2 1

Conversão de arquivo shape file para poligono file
by Antônio Silva 27 Fev '12

27 Fev '12

Prezados! Gostaria de fazer a conversão de um arquivo shape file para poligono file ou arquivo igual ao do pacote MAPS. Abs, Antônio Alcirley

1 0

duvida
by Edson Lira 27 Fev '12

27 Fev '12

Caros, bom domingo a todos. Estou rodando a rotina abaixo: require(xlsx) pe<-read.xlsx("/Users/edsonlira/Edinho/Edinho/Objetivo/Consulta/Jed/pe_diab.xlsx",1,encoding="macroman",noColumns=57,row.names=1) O log está me mostrando o erro: Erro em data.frame(res, ...) : row names contain missing values Qual o erro? Edson Lira Estatístico Manaus-Amazonas

3 2

inscrever na lista
by Rodolfo Souza 27 Fev '12

27 Fev '12

Olá, gostaria de me inscrever na lista. e-mail: rodolfomssouza(a)gmail.com Atenciosamente, -- Rodolfo Marcondes Silva Souza *http://lattes.cnpq.br/4486320338778858* Graduando em Agronomia > Universidade Federal Rural de Pernambuco - UFRPE > Unidade Acadêmica de Serra Talhada - UAST Fazenda Saco, S/N - Caixa Postal 063 CEP 56900-000 Serra Talhada - PE

2 2

Leitura obrigatória
by Benilton Carvalho 26 Fev '12

26 Fev '12

ow.ly/9iE4d

3 2

Problemas com encoding Linux vs Windows
by Daniel C Bezerra 26 Fev '12

26 Fev '12

Prezados amigos do R, Eu uso emacs + ess muitas vezes alterno entre Windows (XP, 32 Bits) e Linux (Ubuntu 64bit 11.04). Ambos os sistemas operacionais estão em inglês. Recentemente eu importei uma planilha de excel pelo read.xls (pacote gdata) no windows e fiz as análises sem problemas. Ao carregar o mesmo workspace (sem importar os dados novamente) no linux, vários caracteres apareceram de modo bizarro. O problema é o mesmo seja rodando o R dentro do emacs ou diretamente no terminal. Por exemplo, a palavra "não" usada em algumas variáveis aparece como "n\xe3o", embora eu consiga escrevê-la e colocá-la em títulos de gráficos, etc. Alguém tem alguma sugestão de como corrigir este problema? Desde já agradeço. Abs, D

4 6

Combinar colunas de um dataframe ou matriz
by Emmanuel Arnhold 25 Fev '12

25 Fev '12

A CMR abaixo me dá as seguintes combinações de colunas do data.frame x 1 e 2; 1 e 3; 2 e 3 Gostaria de fazer na seqüência: 1 e 1; 1 e 2; 1 e 3; 2 e 1; 2 e 2; 2 e 3; 3 e 1; 3 e 2; 3 e 3 É possível? Segue a CMR x1=c(0,0,0,1,0) x2=c(1,0,0,1,1) x3=c(1,1,1,0,0) x=data.frame(x1,x2,x3) comb<-combn(ncol(x), 2, function(y) list(x[,y])) comb

1 0