Novembro 2011 - R-br

[NEWS] Mudanças nas versões do R e melhora de performance
by Fernando Mayer 16 Fev '12

16 Fev '12

Não sei se todos já estão sabendo, mas li isso hoje e achei importante repassar para a lista: Sobre as mudanças de releases do R: R 2.14 to be released on October 31; R 2.13 patch on September 13 http://blog.revolutionanalytics.com/2011/09/r-214-to-be-released-on-october… Brian Ripley on The R Development Process http://blog.revolutionanalytics.com/2011/08/brian-ripley-on-the-r-developme… Sobre a melhora de performance: With byte compiler, R 2.14 will be even faster http://blog.revolutionanalytics.com/2011/08/with-byte-compiler-r-214-will-b… []s, --- Fernando Mayer Dept. Ecologia e Zoologia - ECZ/CCB Universidade Federal de Santa Catarina - UFSC URL: http://sites.google.com/site/fernandomayer e-mail: fernandomayer [@] gmail.com

4 4

Criando um processo de pontos sobre uma área regular
by ASANTOS 10 Jan '12

10 Jan '12

Boa tarde pessoal, Estou tentando criar um processo de pontos aleatório sobre uma área regular representada por árvores regularmente espaçadas, sendo que o processo de pontos só pode ocorrer onde não existam árvores e estou quebrando a cabeça com isso, bom primeiro criei a área com as árvores: require(plotrix) ##Cria o limite da área limx<-c(2,2,20.5,20.5) limy<-c(2,20.5,20.5,2) plot(limx,limy) lim=cbind(limx,limy) lim<-as.matrix(lim) polygon(lim) ##Aqui criei a posição das árvores e com área igual ao seu diâmetro sq=seq(2.5,20,2.5) for (k in sq) { for(j in sq){ draw.circle(k,j,radius=0.3,col="red") }} Bom agora eu queria jogar sobre essa área o processo de pontos: ##Criando o processo de pontos x <- runif(n=500,min=0, max=20) y <- runif(n=500,min=0, max=20) xy<-cbind(x,y) points(xy) ## Porém, os pontos também caem nas mesmas posições ocupadas pelas árvores (obviamente), onde gostaria de saber se alguém conheceria alguma solução para eu fazer com que não haja sobreposição entre os raios que representam os diâmetros das árvores e os pontos criados? Na verdade, preciso definir uma maneira de fazer com que as coordenadas contidas nos diâmetros das árvores não sejam validas no momento da geração do processo de pontos, atualmente estou tentando soluções com o pacote spatstat mais sem sucesso, pois não consigo fazer com que a área criada seja um objeto ppp valido. Obrigado, -- Alexandre DOS SANTOS Engenheiro Florestal, Msc. Laboratório de Entomologia Florestal Departamento de Entomologia Universidade Federal de Lavras Caixa Postal 3037 37200-000 - Lavras - Minas Gerais - Brasil Tel: +55 35 92230304

6 9

[Dúvida] Combinações de strings de um vetor.
by Pedro Rafael 18 Dez '11

18 Dez '11

Pessoal existe algum objeto no R para realizar combinações 2 a 2, 3 a 3. etc? Exemplo: Considere o vetor: vetor <- c("a","b","c","d") Queria obter todas as combinações possíveis com os elementos do vetor, ou seja, uma saída semelhante a: [1] "ab" "ac" "ad" "bc" "bd", "dc" Até poderia pensar em programar isto do zero usando o R mas será que já existe um objeto para fazer isto pronto e eficiente? Obrigado à todos desde já. -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

4 11

Ajuda
by Gilbert Queiroz 05 Dez '11

05 Dez '11

Pessoal, tenho a seguinte parte de um script que estou desenvolvendo: n<-c(87, 108, 80, 25) err.m<-matrix(NA, choose(nx,2),1) err.m nx for(i in 1:nx) { for(j in 2:nx) { if(i<j){ err<-sqrt((410.35/(n[i]+0.5))+(410.35/(n[j]+0.5))) print(err) } } } Informações úteis: nx = 4 (já foi definido em outra parte do script) err.m é uma matriz-coluna (6 x 1) O quero fazer é colocar os resultados de err dentro da matriz err.m Como faço isso??? O comando print me mostra que err está calculando o que quero. Mas me interessa colocar os resultados em err.m Abs.

5 6

lquadmath
by FHRB Toledo 30 Nov '11

30 Nov '11

Senhores, Uso Debian com R 2.14, desde o upgrade instalar (e/ou atualizar) pacotes com implementações acusa o seguinte "erro /usr/bin/ld: cannot find -lquadmath", vi alguma coisa na internet desse problema com R/Debian e as bibliotecas do gcc! Referência: http://groups.google.com/group/rattle-users/browse_thread/thread/f4d004a094… . Como proceder, pois seguindo o recomendado pela referência nada acontece! att, FH

2 4

Perguntas dentro de function
by Guilherme Biz 30 Nov '11

30 Nov '11

Bom dia, Gostaria de saber como fazer uma pergunta dentro de uma function e assim continuar esta mesma função. Por exemplo: Dentro da função é perguntada o peso da pessoa e com a resposta dada continuar a função. Note que a continuação desta função vai depender da resposta. Att, Guilherme Biz

4 4

setwidth: configuração automática de options("width") em terminal
by Jakson Alves de Aquino 30 Nov '11

30 Nov '11

Olá! Problema: Quando o R está sendo usado em emulador de terminal, a opção "width" não é atualizada automaticamente quando o terminal é redimensionado. Solução: Capturar o sinal SIGWINCH (que avisa ao aplicativo que o terminal foi redimensionado) e modificar o valor de options("width"). Estou lançando a versão para teste do pacote setwidth com a implementação da solução mencionada: http://www.lepem.ufc.br/jaa/setwidth_0.9-0.tar.gz Para testar, basta instalar o pacote R CMD INSTALL setwidth_0.9-0.tar.gz e, no R, digitar library(setwidth) Ao ser carregado o pacote, fica habilitada a captura do sinal SIGWINCH. Não há nenhuma função em R visível para o usuário. Agradeço relatos de bug e comentários. -- Jakson Alves de Aquino Universidade Federal do Ceará Departamento de Ciências Sociais www.lepem.ufc.br/aquino.php

5 11

Uso do ifelse para variáveis qualitativas
by Edson Lira 30 Nov '11

30 Nov '11

Como faço para usar o ifelse em variáveis qualitativas. Na rotina abaixo estou tendo um problema. tip_doa$bairro_resid<-ifelse(tip_doa$BAIRRORESIDENCIA=="NA", tip_doa$BAIRRO_NAO_TABELADO,tip_doa$BAIRRORESIDENCIA) Não indica erro, mas vejam table. > table(tip_doa$bairro_resid) 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 372 557 179 1956 531 108 537 252 60 417 57 861 2 328 512 3870 173 3 123 395 332 2046 995 141 240 580 306 1285 41 161 242 1072 261 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 883 581 254 365 114 54 117 245 256 154 561 479 146 3 2 173 681 92 246 828 140 403 1166 592 109 328 183 12 271 1032 117 490 34 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 193 354 616 364 177 713 443 1204 366 1 309 35 513 216 117 20 1 292 765 > Se eu usar NA ao invés de "NA" me dá a saída: > table(tip_doa$bairro_resid) character(0) [].s Edson Lira Estatístico Manaus-Amazonas

4 3

excluir variável específica de uma tabela de dados
by Rodrigo Plei 30 Nov '11

30 Nov '11

Prezados, Eu tenho uma tabela com 32 variáveis dispostas em colunas, todas com o mesmo número de linhas. Gostaria de saber como posso fazer uma segunda tabela com essa tabela original, mas excluindo somente a primeira variável. Agradeço desde logo pela ajuda. 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

2 2

Dúvida em simulação de processos não homogêneo
by Victor Eduardo 29 Nov '11

29 Nov '11

Pessoal, estou com dúvida no script abaixo, pois as vezes fica me dando como resposta NULL e não estou conseguindo identificar o erro. O meu objetivo é Implementar uma funçao que, para um dado T , retorna os instantes das ocorrências de um processo de Poisson nao-homogeneo com funcao de intensidade lambda(t) = 3/(t + 1) lambda_t<-function(t){ if(t<0){ print("tempo negativo") break() } lambdat<-3/(t+1) return(lambdat) } pois_naohomog<-function(){ lambda<-3 i<-1 T<-0 saida<-NULL T<-T+exponencial(lambda) u<-runif(1) while(u<lambda_t(T)/lambda){ saida[i]<-T i<-i+1 T<-T+exponencial(lambda) } return(saida) }

2 1