Janeiro 2013 - R-br

Looping
by Fátima Lima Paula 11 Jan '13

11 Jan '13

Pessoal, tenho uma lista de nomes num banco, que pertencem a outro banco. Exemplo. lista=c("a","b","c","d","e","f","g") nome=c("a","b","a","c","d","h","a","g","e","f","g") inf=c("1","2","3","9","6","5","8","12","4","6","11") banco=cbind(nome,inf) banco=as.data.frame(banco) Quando faço, por exemplo which(banco=="a") Retorna 1 3 7 Preciso criar uma lista com os ids desses nomes no banco. Assim, a lista para o meu exemplo ficaria: lista2= 1 3 7 2 4 5 9 10 8 11 Na verdade, pensei num looping, mas não sei fazer. A minha ideia seria lista2=which(banco=="n") sendo n=cada um dos nomes da 'lista' Teria que fazer isto com todos os n da 'lista' na ordem e ir guardando na lista2 Isso é possível? Alguém poderia ajudar? Obrigada Fátima

3 5

Adicionar um intervalo (quebra - axis.break) no eixo y em gráficos no R
by Wenceslau Teixeira 11 Jan '13

11 Jan '13

*Caros colegas da lista R,* * * *Envio este post para saber se alguém pode me auxiliar ou dar algumas dicas de bibliografia sobre como colocar uma quebra (intervalo) no eixo y em gráficos.* * * *Segue o código reproduzível abaixo cujo intervalo ultrapassa o eixo y e não gostaria de aumentar os limites de Maximo no eixo y pois a maioria dos dados acaba ficando restrita num pequeno intervalo que dificulta a visualização. Eu encontrei o comando (**axis.break) mas não consegui implementar no grafico. * * * *Sds estatísticas e antecipadamente agradeço alguma orientação.* * * *Wenceslau Teixeira * *Embrapa Solos* * * *#X SÃO AS DATAS* * * *x <- c("30/04/2010","13/05/2010","28/05/2010",* * "15/06/2010","30/06/2010","15/07/2010",* * "30/07/2010","15/08/2010","30/08/2010",* * "15/09/2010","30/09/2010","15/10/2010",* * "30/10/2010","15/11/2010","30/11/2010",* * "15/12/2010","30/12/2010","15/01/2011",* * "30/01/2011","15/02/2011","28/02/2011",* * "15/03/2011","30/03/2011","15/04/2011",* * "30/04/2011","15/05/2011","30/05/2011",* * "15/06/2011","30/06/2011","15/07/2011")* * * *#Y* *y <- c(0.067,0.094,0.039,0.048,0.123,0.080,0.036,0.179,0.069,0.601,NA,NA,0.057,NA,NA,0.028,0.273,0.052,0.045,0.476,0.154,0.314,0.120,0.100,0.079,1.208,8.388,3.307,0.774,0.244) * * * *#SDXX E O DESVIO PADRAO DO ELEMENTO FE A 10 CM QUE SERA FEITO O GRAFICO* * * *SDFE10 <-c(0.018,0.028,0.013,0.047,0.132,0.056,0.011,0.177,0.008,0.000,NA,NA,0.028,NA,NA,0.000,0.387,0.030,0.000,0.622,0.037,0.157,0.054,0.030,0.031,1.907,9.203,4.629,1.041,0.189) * * * *# colocar num data.frame, eliminar NA, converter para data* *da <- data.frame(x=as.POSIXlt(x, format="%d/%m/%Y"), y=y,cex.axis=0.5, sd=SDFE10)* *#da <- na.omit(da)* *da$x* * * *#fazer o grafico * * * *plot(y~x,tck=0.01, data=da,ylab=expression(Concentração~de~Ferro~na~solução~do~solo~(mg~L^{-1})),tck=0.01,xlab="", pch=24,cex.axis=0.5, type="p", xaxt="n",ylim=c(-0.2,0.8))* *with(da, arrows(x, y-sd, x, y+sd,cex.axis=0.5, code=3, angle=90,length=0.03))* * * *axis.POSIXct(side=1, at=da$x, las=2, format="%d/%m/%Y",cex.axis=0.5)* * * *#adicionado os outros dados da prof50 cm * * * *par(new=TRUE)* * * *#Y SÃO OS DADOS DO FE50CM* * * *y1 <- c(0.097,0.029,0.026,0.053,0.056,0.014,0.015,0.122,NA,0.140,0.160,NA,0.088,NA,NA,NA,NA,NA,0.017,0.070,0.052,0.173,0.079,0.037,0.050,0.420,0.280,0.030,0.057,0.037) * * * *#SDXX E O DESVIO PADRAO FERRO A 50 CM FEITO O GRAFICO * * * *SDFE50 <- c(0.078,0.000,0.012,0.057,0.011,0.000,0.000,0.000,NA,0.000,0.000,,0.000,NA,NA,NA,NA,NA,0.000,0.027,0.029,0.062,0.019,0.006,0.000,0.114,0.227,0.000,0.000,0.000) * *# break axis* *#axis.break(2, 22, style="zigzag")* * * * * *#fazer o segundo grafico * * * *plot(y1~x, data=da,cex.axis=0.5,* * ylab=expression(Concentração~de~Ferro~na~solução~do~solo~(mg~L^{-1})),xlab="",tck=0.01, pch=19, type="p", xaxt="n",ylim=c(-0.2,0.8))* *with(da, arrows(x, y1-sd, x, y1+sd, code=3, angle=90,length = 0.03))* * * *#Adicionando uma legenda* * * *legend("topright", c("10 cm", "50 cm"), pch = c(24,19))* * * *#informação sobre O tratamento * *legend("topleft", c("Tratamento - TENS Dendê"))* * * -- Wenceslau Geraldes Teixeira Pesquisador Embrapa Embrapa Solos (www.cnps.embrapa.br) Rua Jardim Botânico, 1024 CEP 22460-000 Rio de Janeiro - RJ - Brasil Tel: 55 21 21794583 Skype: wgt007 E-mail: wenceslau.teixeira(a)embrapa.br wenceslau.teixeira(a)pq.cnpq.br CV - http://lattes.cnpq.br/4136381865367906 ________________________________________________________ II ENCONTRO BRASILEIRO DE FÍSICA DO SOLO - 2013 2nd Brazilian Soil Physics Meeting Rio de Janeiro, 06 a 10 de Maio de 2013 http://www.bspm2013.org/ -- Wenceslau Teixeira wgt007(a)gmail.com skype: wgt007

2 1

encontrar num banco a linha de um nome
by Fátima Lima Paula 11 Jan '13

11 Jan '13

Acho que estou surtando. Não consigo achar o indivíduo pelo nome em um data frame. Tenho o nome dele e quero achar a linha em que ele está. Alguém poderia me ajudar, por faor? Obrigada Fátima

3 4

data.frame para planilha de dados
by Nicolay Cunha 11 Jan '13

11 Jan '13

Prezados, Gostaria de saber se existe alguma maneira automatizada de eu reverter um data.frame para uma planilha de dados como no exemplo abaixo. Grato! Nicolay #################################################################################################################### # Data.frame # a partir disso! Ind Tipo L1 L2 L3 M1 M2 M2 C1 C2 C3 S 1 LONGO NA NA NA 20 19 18 15 14 13 30 2 MÉDIO 32 31 30 NA NA NA 16 15 13 22 3 CURTO 33 32.5 31 23 22 21.5 NA NA NA 14 #################################################################################################################### # Planilha de dados # obter isso! Ind Tipo Ant Pol S 1 LONGO M1 20 30 1 LONGO M2 19 30 1 LONGO M3 18 30 1 LONGO C1 15 30 1 LONGO C2 14 30 1 LONGO C3 13 30 2 MÉDIO L1 32 22 2 MÉDIO L2 31 22 2 MÉDIO L3 30 22 2 MÉDIO C1 16 22 2 MÉDIO C2 15 22 2 MÉDIO C3 13 22 3 CURTO L1 33 14 3 CURTO L2 32.5 14 3 CURTO L3 31 14 3 CURTO M1 23 14 3 CURTO M2 22 14 3 CURTO M3 21.5 14 -- Nicolay Leme da Cunha Biólogo, Mestre, Doutorando em Ecologia e Conservação Universidade Federal de Mato Grosso do Sul, 79070-900 Campo Grande, MS, Brasil E-mail: nicolaycunha(a)gmail.com lattes.cnpq.br/5916316648872099

2 2

Re: [R-br] Junção de dados no mclapply
by Junior Beleti 10 Jan '13

10 Jan '13

Entendido, Benilton. Apenas na instrução: X = m1[grid[idx, 2], ] Está criando um objeto X que seria um grid (malha) com os elementos de m1, certo? Mas "idx, 2" teria qual significado? Muito obrigado, Junior.

2 1

Excluir nomes sem sobrenome
by Fátima Lima Paula 10 Jan '13

10 Jan '13

Descobri que não consigo separar os nomes no bancão poque têm vários nomes que não têm sobrenome nenhum. Só que são muitos. Existe uma forma de excluir essas linhas? Aí vai o exemplo com a Maria. Não roda. Quando peço um summary no "separado", retorna que Maria tem length=1. Seria retirar todos que têm length=1. nome=c("Teresa dos Santos","Paulo da Silva","Denis Henrique de Paula","Pedro dos Santos","Pedro Maia", "Maria") idade=c("80","89","78","78","78","85") dataint=c("26","20","25","25","25","20") dados1=data.frame(nome,idade,dataint) dados1 str(dados1) nome1=dados1$nome # remover da,das,do,dos,de nome1=gsub(" d(a|e|o)s? ", " ", nome1) separado1=strsplit(as.character(nome1), split=" ") head(separado1) fa1=sapply(separado1, function(i){ n=length(i) if(n>2) return(c(i[1], paste(i[2:(n-1)],collapse=" "), i[n])) if(n==2) return(c(i[1], NA, i[n])) }) nomenovo1=t(fa1) nomenovo1=as.data.frame(nomenovo1) dim(nomenovo1) colnames(nomenovo1)=c("nome1","nome2","nome3") dados2=cbind(dados1,nomenovo1) head(dados2) Por favor, alguém pode ajudar?

3 3

Diferenças nas saídas
by Fátima Lima Paula 10 Jan '13

10 Jan '13

Pessoal, desculpem, mas não consigo rodar um script que o Walmes me passou. No exemplo, fica tudo bem, mas na real dá zebra. Acho que os objetos estão sendo considerados diferentemente. No exemplo tudo estava com factor. Passei tudo para factor. Nos exemplos os dois bancos eram data.frames. Mas quando faço a função em um a saída vem diferente da do outro. Por favor, alguém poderia me ajudar. Exemplo que apresentei nome=c("Teresa dos Santos","Paulo da Silva","Denis Henrique de Paula","Pedro dos Santos","Pedro Maia", "Maria do Carmo") idade=c("80","89","78","78","78","85") dataint=c("26","20","25","25","25","20") dados1=data.frame(nome,idade,dataint) dados1 str(dados1) nome1=dados1$nome # remover da,das,do,dos,de nome1=gsub(" d(a|e|o)s? ", " ", nome1) separado1=strsplit(as.character(nome1), split=" ") head(separado1) fa1=sapply(separado1, function(i){ n=length(i) if(n>2) return(c(i[1], paste(i[2:(n-1)],collapse=" "), i[n])) if(n==2) return(c(i[1], NA, i[n])) }) nomenovo1=t(fa1) nomenovo1=as.data.frame(nomenovo1) dim(nomenovo1) colnames(nomenovo1)=c("nome1","nome2","nome3") dados2=cbind(dados1,nomenovo1) head(dados2) Quando faço head para o fa1 no meu exemplo, retorna assim, que pelo pouco que sei de R, parece uma matriz > head(fa1) [,1] [,2] [,3] [,4] [,5] [,6] [1,] "Teresa" "Paulo" "Denis" "Pedro" "Pedro" "Maria" [2,] NA NA "Henrique" NA NA NA [3,] "Santos" "Silva" "Paula" "Santos" "Maia" "Carmo" Já quando peço o head do fa do meu bancão, retorna assim e zebra tudo a partir daí. > head(fa) [[1]] [1] "HERALDO" "TAVARES" "SILVA" [[2]] [1] "JUPITER" NA "SILVA" [[3]] [1] "JACI" NA "SANTOS" [[4]] [1] "MARCOS" "SILVA" "FILHO" [[5]] [1] "PAULO" "SILVA" "PORTO" [[6]] [1] "INACIO" "AVELINO" "BELO" Alguém tem ideia do que seja? Obrigada Fátima

3 3

Busca textual
by Clécio da Silva Ferreira 10 Jan '13

10 Jan '13

Olá. Tem uma mestranda de Linguística interessada em fazer busca textual de ocorrências de associações de uma palavra com outra na mesma sentença. Essa busca é feita em textos (livros, transcrições de falas, etc.). Existe alguma ferramenta no R para isso? Atte, xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx x Prof. Clécio da Silva Ferreira x x Depto Estatística - UFJF x xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx

4 3

Separando nomes e retirando de da dos das
by Fátima Lima Paula 09 Jan '13

09 Jan '13

Prezados, há um tempo Walmes me passou um script para separar nomes. Como vcs podem verificar o nome do meio ficou com um problema (Pedro Maia), mas eu só ia utilizar o primeiro e último nome, por isso não me importei. Mas agora preciso de tudo. Assim, eu gostaria de saber se alguém poderia me ajudar desta forma: acertando o problema que retornou para quem só tem dois nomes e retirando os "de", "das", "da", etc dos nomes. Não sei programar e nem imagino como se faz para acrescentar isso no script. Você poderia me ajudar? nome=c("Teresa dos Santos","Paulo da Silva","Denis Henrique de Paula","Pedro dos Santos","Pedro Maia", "Maria do Carmo") idade=c("80","89","78","78","78","85") dataint=c("26","20","25","25","25","20") dados=data.frame(nome,idade,dataint) dados Ulitizei seu script assim: nome=dados$nome separado=strsplit(as.character(nome),split=" ") fa=sapply(separado, function(i){ n=length(i) c(i[1],paste(i[2:(n-1)],collapse=" "), i[n]) }) nomenovo=t(fa) nomenovo=as.data.frame(nomenovo) dim(nomenovo) colnames(nomenovo)=c("nome1","nome2","nome3") dados1=cbind(dados,nomenovo) head(dados1) Obrigada Fátima

2 4

Excluindo "da" "de" ao mesmo tempo
by Fátima Lima Paula 09 Jan '13

09 Jan '13

Olá, pessoal. O Daniel Marcelino me passou esse script com a solução para extrair os "da" "de" dos nomes. Mas como sou um horror em R não consigo eliminar tudo ao mesmo tempo. Se rodo o script dele, funciona para o "de". Mas quando repito para "da" retorna um problema. Como não manjo nada, tentei usar list(c("de","da")). Não funcionou. Tentei usar ("de" | "da"). Nada. Alguém pode me ajudar, por favor. Obrigada Fátima mydata <- data.frame( name = a<-c("Ricardo de Almeida","Maria José da Cruz","Maria Tereza dos Santos","Paulo da Rocha","José Maria", "Germano das Neves"), age = b<-c(67,89,78,65,68,67), birth = c <-c(1945,1923,1934,1947,1944,1945), date = d<-c(20120907,20120703,20120301,20120805,20121004,20121209) ) mydata$name2 <-strsplit(as.character(mydata[,1]),split=c(" de "))

2 1