Dezembro 2011 - R-br

Python para R
by Daniel Marcelino 04 Jun '12

04 Jun '12

Caros, alguem da lista conhece programação python? Estou há dias tentando desenvolver uma rotina em R para explorar dados eleitorais, mas não tenho conseguido ir muito longe. Eu tenho um modelo do que mais ou menos eu preciso, mas a rotina está em python. Alguem encara me ajudar nisso? Daniel ----------------------------------------------------- My contacts: Phone 1+514 343 6111 Ext 3799 Skype: d.marcelino 3200 Jean Brillant, Office C5075 Montreal, QC; H3T 1N8 Canada -----------------------------------------------------

3 2

[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

Leitura de arquivos
by Emerson Cotta Bodevan 04 Jan '12

04 Jan '12

Prezados, tenho um arquivo .txt com uma sequência enorme de números. Por ex.: 0120304405560085673947684658495604932570932847509836536793586593420574092398564230562439065423095624359842350234567230945 7340127459042759834265093248560147237561092436509864230543609586723096756759867243670934675.... . . . ...87350274523475982734596239467358927609650972345792034759237459623767235967526 chego a ter milhares de números! Veja que não há espaço entre os dígitos da sequência! Gostaria de ler esse arquivo para um vetor mas, ao ler, quero que o R entenda que cada dígito deva ser armazenado numa posição diferente no vetor criado. Por ex. (no arquivo acima): 0 1 2 0 3 0 4 4 0 5 5 6 0 0 8... cada um pertencendo a uma posição no vetor criado. Agradeço qualquer ajuda, -- *Emerson*

4 7

Dúvida: Erro em shapiro test
by Flávia Machado 31 Dez '11

31 Dez '11

Olá! Estou tendo problemas com o teste shapiro wilk. Inseri o arquivo txt (com pontos ao invés de vírgulas) dos meus dados com o comando "dados<-read.table(file.choose(), header=T). Os dados foram carregados e ao digitar dados todos os objetos apareceram normalmente, mas quando fui fazer o teste de shapiro apareceu " Erro em sort(x[complete.cases(x)]) : objeto 'FRUTO' não encontrado". Para o teste de shapiro, digitei: "shapiro.test(FRUTO)" Estou achando bem estranho porque a pouco tempo fiz o mesmo teste, da mesma forma, com outra série de dados e deu certo. Hoje usei estes mesmos dados apra testar e deu o mesmo erro. Vocês saberiam me dizer como solucionar esse problema? Obs.: sou bem iniciante no R, acho que deu para perceber. Atenciosamente, -- Flávia Machado

2 3

modelo multinomial logistico
by Luiz André Zardo 30 Dez '11

30 Dez '11

Prezados, estou usando a função multinom para ajustar um modelo logistico. Porém não sei como calcular os valores da estatistica da razão de verossimilhança, do pseudo-R2 e da estatistica de Wald. Alguém saberia me indicar como conseguir esses calculos??? ou indicar algum outro pacote que possa disponibilisar essas estatisticas?? Obrigado Luiz Zardo

2 1

aggregate
by Edson Lira 30 Dez '11

30 Dez '11

Caros amigos do R, na rotina abaixo gostaria de calcular além da média, o desvio padrão e o CV=(sd/mean). Como fazer estes cálculos? ## example with character variables and NAs testDF <- data.frame(v1 = c(1,3,5,7,8,3,5,NA,4,5,7,9), v2 = c(11,33,55,77,88,33,55,NA,44,55,77,99) ) by1 <- c("red","blue",1,2,NA,"big",1,2,"red",1,NA,12) by2 <- c("wet","dry",99,95,NA,"damp",95,99,"red",99,NA,NA) aggregate(x = testDF, by = list(by1, by2), FUN = "mean") # and if you want to treat NAs as a group fby1 <- factor(by1, exclude = "") fby2 <- factor(by2, exclude = "") aggregate(x = testDF, by = list(fby1, fby2), FUN = "mean") FELIZ 2012 a todos os amigos e familiares da lista. Edson Lira Estatístico Manaus-Amazonas

2 1

Retirar elemento de uma lista
by Rodrigo Coster 29 Dez '11

29 Dez '11

Povo, como faço para retirar apenas um elemento de uma lista? No caso, eu tenho essa lista: list(fruta='melancia',x=matrix(rpois(100,5),10),y=matrix(rpois(200,7),10)) Eu quero uma lista só com x e y (no caso a minha lista tem um numero variado de 'termos', não apenas 3). É o equivalente a índice negativo em vetor, mas não funciona em list :( []'s

3 3

Encontrar o maior valor por linha entre colunas
by Daniel Marcelino 29 Dez '11

29 Dez '11

Olá a todos. Eu tenho um pequeno problema que é justamente encontrar qual é o maior valor entre as colunar em cada linha. Como eu poderia obter uma nova coluna ou mesmo um objeto com o nome da coluna para cada valor encontrado. Alguém tem alguma sugestão? Daniel

7 10

[Dúvida] Contar NA em Subset
by Pedro Rafael 28 Dez '11

28 Dez '11

Pessoal como fazer um comando subset contar os elementos NA? Considere o data.frame abaixo: dados = data.frame(x = c(1,2,3,7,12), y = c(0,10,2,3,NA)) Com este data.frame chamado dados gostaria de filtrar com o subset os elementos de y que são diferente de 10. Para o meu exemplo NA também é diferente de 10 e gostaria de fosse contabilizado. Com o mando abaixo o resultado segue logo após. > subset(dados, y != 10) x y 1 1 0 3 3 2 4 7 3 Veja que o NA não foi contabilizado. -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

4 6