- R-br - DINF UFPR

manipulação de colunas
by Rodrigo Plei 02 Mar '12

02 Mar '12

Prezados, Como faço para colocar em uma única colunas os dados que estão em diversas colunas com o mesmo número de linhas? Por exemplo, tenho um arquivo com uma tablela de dados assim col1 col 2 col 3 col 4 154 122 258 155 128 137 259 144 111 128 233 143 Precio que fique assim: coluna 154 128 111 122 137 128 258 259 233 155 144 143 Agradeço desde logo. Rodrigo -- =8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8--> Rodrigo Silvestre Martins, PhD Bolsista FAPESP Pós-Doutorado Laboratório de Ecossistemas Pesqueiros (LabPesq) Universidade de São Paulo, Instituto Oceanográfico Praça do Oceanográfico, 191. Cidade Universitaria (sala 107-A/B) Butantã - São Paulo/SP, Brasil 05508-900 Tel: +55 11 3091 6549 Email: rodrigo.plei(a)gmail.com ; ocersm(a)lycos.com; rsmartins(a)usp.br CV Lattes: http://lattes.cnpq.br/5350064124902777

3 4

[Dúvida] Atribuindo projeção no mapa.
by Pedro Rafael 02 Mar '12

02 Mar '12

plot(mapa, col = colors[findInterval(dados$ano_2008, brks,all.inside=TRUE)], axes = TRUE) Senhores, gostaria de inserir o sistema de referência espacial (4291) na hora de adicionar o vetor de shapfile ao R. PS: O .shp não possui sistema de referência atribuido. -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

2 1

Duvida sobre como modificar dados em um data.frame
by Augusto Ribas 02 Mar '12

02 Mar '12

Ola pessoal. Estou com uma duvida simples aqui. Eu estava analisando dados aqui e vi um grafico estranho, dai vi que os dados estavam digitados errados. Mas o erro é simples. Existem 2 fatores, e os niveis de um fator estão trocados, mas somente dentro de um nivel do segundo fator, o resto esta tudo correto. Basicamente existe um vetor e eu tenho que trocar A por B e vive versa. Ai eu consegui arrumar imendando um monte de comando mas ficou uma coisa enorme e especifica para algo que parece tão simples. ai pensei se alguém tem uma solução mais economica e geral que a que eu fiz? ### #Exemplo de dados ### resposta<-c(c(rnorm(10,5),rnorm(10,8)),c(rnorm(10,8),rnorm(10,5))) fator1<-c(rep(c("A","B"),each=10),rep(c("A","B"),each=10)) fator2<-rep(c("1","2"),each=20) dados<-data.frame(resposta,fator1,fator2) ### #Grafico Errado ### library(lattice) bwplot(resposta~fator1|fator2,data=dados) ### #Minha Solução ### dados[which(dados$fator2%in%"2"),"fator1"]<-sapply(dados[which(dados$fator2%in%"2"),"fator1"],function(x){ if(x=="B"){ replace(x, x == "B", "A") } else { replace(x, x == "A", "B") } }) ### #Grafico Certo ### bwplot(resposta~fator1|fator2,data=dados) -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

4 4

R + GPU (Era: Dá pra Rodar o R nos Tablets ANDROID??)
by Fernando Mayer 02 Mar '12

02 Mar '12

Jogando lenha na fogueira (: gputools promete otimizações incríveis usando a tecnologia CULA da nVidia http://brainarray.mbni.med.umich.edu/Brainarray/Rgpgpu/ Ainda não testei. Alguém? []s, --- Fernando Mayer Universidade Federal de Santa Catarina - UFSC Departamento de Ecologia e Zoologia - ECZ/CCB URL: http://sites.google.com/site/fernandomayer e-mail: fernandomayer [@] gmail.com 2012/3/1 Benilton Carvalho <beniltoncarvalho(a)gmail.com>: > programar para usar a GPU para fazer processamentos usualmente feitos em CPU. b > > 2012/3/1 Rodrigo Sant'Ana <rodrigo.gringo(a)gmail.com>: >> Concordo, em termos de custo beneficio final, ainda vale mais a pena >> manter um Notebook, ou até mesmo >> netbook, e concordo tbm, e foi exatamente isso que quis dizer no email >> anterior, que para "produção" não vale a pena ainda. >> Mas penso que como ferramenta secundária/auxiliar ao seu PC de >> produção possa trazer algum beneficio por agora, só isso. >> >> Mas você está certo. Já havia entendido o que você quis dizer, só >> estava ressaltando que esta >> é uma área de inovação para atuação do R e até mesmo como uma área de >> desenvolvimento em ascensão. >> >> Quando você diz em programar para GPU, você diz para optimizarmos >> funções para utilizar melhor as funções gráficas da placa, ou >> para utilizar o processador da placa de video como um processador fisico? >> >> >> >> >> _______________ >> Rodrigo Sant'Ana >> Oceanógrafo >> >> >> >> Em 1 de março de 2012 18:36, Benilton Carvalho >> <beniltoncarvalho(a)gmail.com> escreveu: >>> Rodrigo, >>> >>> o que eu quero dizer eh que, se vc for usar um dispositivo movel com >>> teclado externo e um processador que nao seja compativel com i386 ou >>> x86_64, entao (hoje!) eh melhor vc usar um laptop ou, se vc estiver >>> mesmo preocupado com tamanho/peso, use um netbook que custa 1/3 do >>> valor e vai fazer exatamente o q tablet faria em termos de R... >>> >>> Se voce estiver nos EUA, o Asus Transformer Prime vai te custar, pelo >>> menos, USD 499 (tablet) + USD 149 (dock) = USD 648.... E, por esse >>> valor, vc compra algum Toshiba com processador i5 e 6GB RAM.... Vai >>> ter ate' quem prefira comprar um Dell com processador i3 com os 499 e >>> um netbook por alguns 200 dolares... Ai' vai do custo-beneficio de >>> cada um... >>> >>> Se for para eu produzir algo de util sem um teclado externo, eu nao >>> consigo ficar mais de 10min usando um tablet. Se for para usar um >>> tablet com teclado externo, prefiro um (net/note)book. (preferencia >>> pessoal: mais aplicativos disponiveis, melhor suporte de >>> software/hardware e por ai' vai) >>> >>> O mesmo aplica-se (novamente, *pelo menos por enquanto*) ao caso do >>> Ubuntu com seu desktop para Android: se vc nao tiver um teclado >>> externo e um monitor, continua "inutil" para nossas aplicacoes. >>> >>> Concordo que, cada vez mais, teremos `a disposicao equipamentos mais >>> portateis. Entretanto, para "producao" ainda ficaremos dependentes de >>> teclado e bateria por mais alguns (poucos, espero) anos... >>> >>> Enquanto isso, todo o exercicio de rodar R (ou equivalentes de outros >>> setores, etc) sob Android ou iOS vale (como ja' dito pelo Leonard) >>> pela prova do conceito e inovacoes que ainda virao disso... Mas, ainda >>> nao para "producao"... not yet.. >>> >>> Comentario provocativo: os comentarios acima tambem valem para os >>> "beneficios" da programacao para GPU vis-a-vis prog. para CPU... >>> >>> b >>> >>> >>> 2012/3/1 Rodrigo Sant'Ana <rodrigo.gringo(a)gmail.com>: >>>> Pô Benilton concordo com você em alguns pontos. Mas acho que podemos >>>> rodar algumas rotinas em tablets sim. Esse Asus vem com um dock de >>>> teclado bastante interessante, e talvez, para algumas matérias básicas >>>> e scripts simples, como mesmo calcular uma média, desvio padrão, >>>> tabela de frequencia, histogramas, teste t, e outros tipos de >>>> inferências mais comuns e simples porém didáticas e necessárias para >>>> uma sala de aula, talvez não no nível de aprofundamento de um curso de >>>> graduação mais direcionado, mas em matérias inicias, ou em mini-cursos >>>> com bases de dados sucintas utilizadas para exemplificação, seja >>>> interessante. >>>> >>>> Entendo que com grandes banco de dados não seja o caminha, mas >>>> utilizar um tablet para acesso remoto via ssh, onde o trabalho pesado >>>> será executado por um computador mais potente nos dê mobilidade, é um >>>> tecnologia interessante. Acho que a explanação do Leonardo é válida, >>>> talvez teremos celulares como computadores pessoais, realizando >>>> trabalho em conjunto com computadores mais potentes.. >>>> >>>> Acho bastante válido, mas entendo sua colocação Benilton, com certeza, >>>> no dia a dia de analises mais densas, não seja tão interessante. >>>> >>>> Abração >>>> _______________ >>>> Rodrigo Sant'Ana >>>> Oceanógrafo >>>> >>>> >>>> >>>> Em 1 de março de 2012 15:38, Leandro Marino >>>> <leandromarino(a)leandromarino.com.br> escreveu: >>>>> Para rodar o R via ssh você terá ter ter uma conexão ativa a internet o que >>>>> muitas vezes não é bom. >>>>> >>>>> Alguns softwares, como por exemplo, o QuantumGis já estão com versões para >>>>> Android. Acho que é um caminho sim. >>>>> >>>>> Cabe salientar que a Canonical (desenvolvedora do Ubuntu) já anunciou versão >>>>> do Ubuntu para Celulares e Tablets. Inclusive é possível ter um telefone da >>>>> motorola que tenha webtop como Atrix rodando o ubuntu. >>>>> >>>>> Com o desenvolvimento do Windows 8 suportando processadores ARM, e linux, >>>>> acredito que seja uma futura convergencia. Acredito que em breve os >>>>> celulares serão nossos computadores pessoais... :) >>>>> >>>>> >>>>> >>>>> Atenciosamente, >>>>> Leandro Marino >>>>> http://www.leandromarino.com.br (Fotógrafo) >>>>> http://est.leandromarino.com.br/Blog (Estatístico) >>>>> Cel.: + 55 21 9845-7707 >>>>> Cel.: + 55 21 8777-7907 >>>>> >>>>> >>>>> >>>>> Em 1 de março de 2012 15:22, Rodrigo Sant'Ana <rodrigo.gringo(a)gmail.com> >>>>> escreveu: >>>>> >>>>>> Depende Leonard, >>>>>> >>>>>> existem tablets hoje em dia bastante robustos, um exemplo é o Asus >>>>>> Transformer Prime (sem querer fazer propaganda. Que vem com um >>>>>> processador Quad-core da nVidia, num é a solução, nem se compara a um >>>>>> PC. Mas ter o R rodando nele, para rotinas rápidas, talvez seja >>>>>> funcional para quem trabalha viajando muito! >>>>>> >>>>>> Só quis colocar a ideia, pois lembrei destas novas gerações de >>>>>> tablets, e até mesmo, de uma necessidade que eu compartilho, pois acho >>>>>> interessante ter esta mobilidade no uso do R, em um tablet. Para aulas >>>>>> em geral, para demonstrações de códigos em palestras, para uso em >>>>>> qualquer ocasião. >>>>>> >>>>>> Outra coisa, achei bem interessante essa parte do acesso ao R via ssh, >>>>>> talvez seja um tópico interessante para iniciar a discussão aqui na >>>>>> lista, com demonstrações e tudo mais! Pois eu desconheço completamente >>>>>> e acho super interessante e importante. >>>>>> >>>>>> Abraço a todos! >>>>>> >>>>>> _______________ >>>>>> Rodrigo Sant'Ana >>>>>> Oceanógrafo >>>>>> >>>>>> >>>>>> >>>>>> Em 29 de fevereiro de 2012 13:06, Benilton Carvalho >>>>>> <beniltoncarvalho(a)gmail.com> escreveu: >>>>>> > Concordo contigo, Leonard... mas como o OP nao perguntou minha >>>>>> > opiniao, ne'? rsrsrs >>>>>> > _______________________________________________ >>>>>> > 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. >>>>>> _______________________________________________ >>>>>> 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. >>>>> >>>>> >>>>> >>>>> _______________________________________________ >>>>> 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. >>>> _______________________________________________ >>>> 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. >>> _______________________________________________ >>> 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. >> _______________________________________________ >> 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. > _______________________________________________ > 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.

3 3

Agrupamento de Tocher
by Melissa Oda 01 Mar '12

01 Mar '12

Alguém conhece um função no R para o agrupamento de Tocher? Obrigada Melissa

1 0

[Dúvida] Uso de linguagem funcional para estatística.
by Pedro Rafael 29 Fev '12

29 Fev '12

Alguem aqui domina alguma linguagem de programação funcional? Existe alguma vantagem de quem é da estatística ou matemática também aprender profundamente alguma linguagem funcional como é o caso de Haskell? Atualmente tudo que preciso consigo em R e C mas ultimamente venho lendo sobre linguagem funcional e cálculo lambda e me parece bem interessante para quem trabalha matemática.... -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

7 22

ebooks com 50% desconto
by Benilton Carvalho 29 Fev '12

29 Fev '12

Celebrando a Strata Conference deste ano, a O'Reilly esta' oferecendo 50% de desconto nos livros relacionados a dados (e analise). Divirtam-se: http://oreilly.com/data/index.html b

2 1

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