Abril 2012 - R-br - DINF UFPR

Ajuda para evitar um loop com elegância 2
by Gustavo Marcatti 23 Abr '12

23 Abr '12

Olá pessoal, tenho que executar uma função parecida com essa no R. Mas ela normalmente demora longos 500 segundos para executar. Fiz uma versão que chama um código C, e essa é bem eficiente (demora cerca de 0,5 segundos), porém gostaria de uma função totalmente em R. Caso alguém queira ajudar, desde já agradeço. set.seed(1) vetor1 <- rep(1:4, each = 5) var1 <- abs(rnorm(4,3)) * vetor1 cond1 <-abs(rnorm(4,10)) * vetor1 var2 <- matrix(0, nrow=length(vetor1)) cond2 <- matrix(0, nrow=length(vetor1)) for(i in 1:(length(vetor1) - 1)){ if(vetor1[i] == vetor1[i+1]){ var2[i] = var1[i+1] cond2[i] = cond1[i+1] }else { var2[i] = NA cond2[i] = NA } } var2[length(vetor1)] <- NA cond2[length(vetor1)] <- NA dados <- na.exclude(data.frame(var2, cond2)) Gustavo Marcatti Eng. Florestal

2 2

R servidor Windows?
by Vinicius Brito Rocha 23 Abr '12

23 Abr '12

Senhores, é possível instalar o R em um servidor Windows para a utilização de multiplus usuários? Att. -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

2 2

bootstrap
by Maria Papa 23 Abr '12

23 Abr '12

Olá Pessoal, boa noite Estou tentando, sem sucesso, realizar uma reamostragem bootstrap para médias de 3 amostras estão organizados em uma matriz na qual, que sendo cada coluna representa uma amostra. Porém, a cada reamostragem, eu preciso que seja feita para as três amostras, ou seja, quando um dado for selecionado para formar a amostra bootstrap, deve ser retirada a linha toda. O resultado final deve ser 3 médias reamostradas, que foram obtidas em um único processo de reamostragem. Alguém sabe como fazer isso no R? Eu consegui somente quando eu realizo 3 procedimentos de reamostragem distintos. Desde já agradeço a todos célia

3 2

Dúvida de Geoestatística
by Mauro Sznelwar 21 Abr '12

21 Abr '12

Queria saber se a bibliotéca GeoR faz kriging external drift (KED), e se tem o recurso para análise multivariada, onde eu tenho mais de uma variável como variável resposta. Outra coisa, se ela é tridimensional, pois eu tenho as dita cujas em várias camadas, é uma exploração arqueológica. Caso não tenha estes recursos, quais bibliotécas as tem?

2 3

Problema com pacote parallel
by Vanessa Leotti Torman 20 Abr '12

20 Abr '12

Prezados membros da lista, Bom dia. Estou rodando uma simulação bastante complexa, pois para cada replicação tenho que rodar uma cadeia mcmc através do pacote brugs (que necessita do software OpenBUGS). Para que a execução seja mais rápida, estou usando o pacote paralell, mais especificamente a função parCapply deste pacote. Os testes com poucas replicações funcionam perfeitamente. Mas quando coloco para rodar todas (4000), mesmo após várias horas executando, o Windows dá uma mensagem que o R front end parou de funcionar. Eu percebo através do gerenciador de tarefas que apenas um dos núcleos parou de executar, então clico em fechar e acabo perdendo todos os resultados. Já tentei não usar todos os núcleos do processador, mas mesmo assim o problema acontece. Alguém teria alguma sugestão? Abraços, Vanessa

3 5

Dependencias no mclapply
by Junior Beleti 20 Abr '12

20 Abr '12

Olá a todos. Na busca de melhorias no processamento paralelo, estou com uma dúvida quanto a função mclapply do pacote multicore. Benilton já me deu várias dicas de como proceder com o mesmo, mas ainda tenho a seguinte dúvida: - Quanto a questão de dependências dentro das operações do mclapply, por exemplo, nas seguintes instruções: mclapply(1:100, function(i) { g <- mvrnorm(n=1, mat.cokri[[1]],mat.cokri[[2]]) #1 seq1 <- seq(1,nlinhas,by=2) #2 seq2 <- seq(2,nlinhas,by=2) #3 y1 <- g[seq1] #4 y2 <- g[seq2] #5 gerado <- data.frame(y1,y2) #6 compos <- agl(gerado) #7 compos1 <- cbind(compos,compos1) #8 },mc.cores=4 ) - Em minha análise, cada um dos 4 núcleos vai executar todas as instruções dentro da função, ou seja, cada núcleo vai executar as 8 instruções dentro da função. Mas a dúvida surge quanto as dependências, visto que em uma execução sequencial, todas as instruções seriam executadas sequencialmente, e a próxima iteração só iniciaria após o resultado de "compos1" ter sido atribuído. E quanto ao processo paralelo, seria respeitado tal sequencia de instruções, ou seja, a segunda iteração inicia somente após "compos1" ter sido atribuído na primeira iteração? Como o mclapply possui tal controle? Estou ciente que esta questão pode ser mais da implementação do próprio pacote, mas se alguém puder me esclarecer ficaria muito grato. Att, Carlos Beleti.

6 13

Dúvida conectar R com Interbase/firebird
by Marco Pollo 20 Abr '12

20 Abr '12

Caros, Registro os agradecimentos do email anterior, mas, alguém saberia qual seria configuração do ODBC (na plataforma Windows 32bits) para essa base? E como ficaria o comando odbcConnect, considerando minhas tentativas (infrutíferas) registradas abaixo: Configuração do meu ODBC INTERBASE/FIREBIRD : TENTATIVAS NO R: 1ª tentativa: channel <- odbcConnect(SISATER, uid = "SYSDBA", pwd = "masterkey") Erro em paste("DSN=", dsn, sep = "") : objeto 'SISATER' não encontrado 2ª tentativa: channel <- odbcConnect(SISATER.GDB, uid = "SYSDBA", pwd = "masterkey") Erro em paste("DSN=", dsn, sep = "") : objeto 'SISATER.GDB' não encontrado 3ª tentativa: channel <- odbcConnect(C:\SISATER\SISATER.GDB, uid = "SYSDBA", pwd = "masterkey") Erro: entrada inesperado em "channel <- odbcConnect(C:\" Obrigado Marco Pollo

2 1

Re: [R-br] Dependencias no mclapply
by Junior Beleti 20 Abr '12

20 Abr '12

Okay Benilton, entendido.. Obrigado, Carlos Beleti.

1 0

STARMA
by Gilbert Queiroz 20 Abr '12

20 Abr '12

Alguém pode me ajudar a encontrar alguma coisa referente a isto???? Mesmo que ainda não tenha código p/ o R. Abs. Gilbert

3 3

(sem assunto)
by Paulo Ricardo Silva Moreira 19 Abr '12

19 Abr '12

Boa tarde, Estou precisando usar a função tapply. Alguém sabe como utilizá-la? Att, Paulo Ricardo.

3 2