Agosto 2012 - R-br

Mapa de Bairros do Rio de Janeiro
by Sandro Sperandei 14 Ago '12

14 Ago '12

Srs, Estou tentando encontrar um mapa que tenha os bairros do Rio de Janeiro em um formato utilizável pelo maptools e que seja o mais atual possível. Os que poderiam existir no endereço do IBGE (http://www.ibge.gov.br/home/download/geociencias.shtm) estão vindo com a mensagem “Final inesperado do arquivo” quando o arquivo .zip é aberto. Alguém pode me ajudar? Att Sandro Sperandei

2 3

Cálculo entre todos os pares de distâncias em um objeto espacial
by ASANTOS 14 Ago '12

14 Ago '12

Boa noite pessoal, Gostaria de calcular todas as distâncias de 6 coordenadas em UTM, à partir de um objeto da classe espacial, para isto fiz: ### Cálculo da distancia entre todos os pontos em UTM ninhos<-c("N01","N02","N03","N04","N05","N06") lat<-c(0503252,0503319,0503716,0503682,0503688,0503707) long<-c(7651950,7651694,7652672,7652670,7652627,7652662) ##Cálculo que compara a primeira coordenada com a segunda e assim sucessivamente .... lat.diff<-diff(lat) long.diff<-diff(long) d=(long.diff^2+lat.diff^2)^.5 d ##Porém não era bem isto que eu queria ainda ##Transformando em um objeto da classe espacial require(sp) pontos<- SpatialPoints(cbind(lat,long)) pontos<- SpatialPoints(list(lat,long)) pontos<- SpatialPoints(data.frame(lat,long)) Agora, gostaria de baseado no objeto espacial pontos, calcular todas as combinações de distâncias entre as coordenadas fornecidas. Alguém saberia me indicar uma função em algum pacote que faz isso, Obrigado, -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304

4 4

Teste - mudança de email , desconsiderar
by marcelo claro de souza 14 Ago '12

14 Ago '12

1 0

Encontrar valores duplicados por linha em um data.frame
by Guilherme Gerhardt Mazzochini 14 Ago '12

14 Ago '12

Olá pessoal, Eu possuo um banco de dados com com quase 900 mil linhas e 23 colunas. Segue um exemplo simplificado abaixo: set.seed(321) v<-data.frame(rbind(sample(1:3,10,replace = T), sample(1:3,10,replace = T), sample(1:3,10,replace = T), sample(1:3,10,replace = T), rep(1,10))) > v X1 X2 X3 X4 X5 X6 X7 X8 X9 X10 1 3 1 2 2 3 2 1 2 3 2 2 1 1 2 3 2 3 3 1 1 1 3 1 1 2 2 2 3 2 1 3 1 4 2 3 3 2 3 1 3 3 2 1 5 1 1 1 1 1 1 1 1 1 1 Para cada linha eu quero determinar qual a coluna do primeiro valor que não é duplicado do primeiro elemento da linha (coluna 1) e salvar em um vetor o numero da coluna do primeiro elemento não duplicado. Assim, > which(v[1,]!=v[1,1])[1] [1] 2 Ou seja, a segunda coluna possui o primeiro elemento diferente do elemento na primeira coluna. Eu tentei fazer um loop para fazer isso para todas as linhas, mas devido ao grande número de linhas está demorando muito tempo. x<-numeric(nrow(v)) for(i in 1:nrow(v)){ x[i]<-which(v[i,]!=v[i,1])[1] } > x [1] 2 3 3 2 NA # o valor NA é porque na quinta linha todos os valores são iguais. Alguém tem idéia de como eu posso fazer isso de forma mais rápida?? Obrigado antecipadamente, Guilherme

2 1

Re: [R-br] Duvida no TAPPLY
by Ivan Bezerra Allaman 14 Ago '12

14 Ago '12

?rowSums, ?colSums \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0

Re: [R-br] Eliminar outliers de uma planilha.
by Ivan Bezerra Allaman 14 Ago '12

14 Ago '12

Uma alternativa seria: remove_outliers <- function(x, na.rm = TRUE, ...) { qnt <- quantile(x, probs=c(.25, .75), na.rm = na.rm, ...) H <- 2.5 * IQR(x, na.rm = na.rm) y <- x y[x < (qnt[1] - H)] <- NA y[x > (qnt[2] + H)] <- NA y } (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

5 7

microdados censo IBGE 2010
by Manoel Nascimento 14 Ago '12

14 Ago '12

Olá! Sou novo na lista e no uso do R, ainda estou me adaptando e entendendo como funciona, e já apareceram várias dúvidas. Estou precisando ler os microdados do Censo IBGE 2010. Pelo que vi, os pacotes "IBGEPesq" e "dicionariosIBGE" ainda não foram atualizados para nos fornecer os dicionários para o censo, então as diversas funções que eles fornecem ainda não servem para o censo. Por isso, fiz o seguinte (podem dar risada que é coisa de novato mesmo): 1) Fiz um script para cada conjunto de dados, que é mais ou menos o seguinte (mudando alguns dados para cada conjunto de dados -- domicílios, mortalidade, pessoas, migrações): initial.dir<-getwd() setwd("/diretorio/local/IBGE/MicrodadosCenso2010") sink("DomiciliosCenso2010.out") DomicilioAC2010 <- read.fwf("/diretorio/local/IBGE/MicrodadosCenso2010/AC/Dados/Domicilios/Amostra_Domicilios_12.txt", widths = c(2, 5, 13, 8, 3, 1, 2, 3, 2, 1, 2, 2, 1, 6, 4, 1, 2, 2, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 7, 5, 6, 4, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1)) [A mesma sintaxe é repetida para os microdados de domicílio de cada estado.] save.image("~/Documentos/cidades/IBGE/MicrodadosCenso2010/Censo2010Domicilios.RData") setwd(initial.dir) Feito deste modo, o script demora um pouco para terminar, mas funciona direitinho. Todas as tabelas são construídas de acordo com as informações que forneci. 2) Tentei fazer este script funcionar usando o "read.fwf" com a variável "col.names" para inserir os cabeçalhos de cada variável, mas ele dá um erro. (Terei que refazer mais tarde o script com esta variável para colar aqui o erro que dá.) Por isto, consegui construir as tabelas, mas elas estão sem cabeçalhos. Pergunto: há algum modo de inserir os cabeçalhos nas tabelas sem ser pela variável "col.names" do comando "read.fwf"? 3) Alguns números constantes nas tabelas são códigos, a serem substituídos por referências específicas. (Exs.: na coluna "Unidades da Federação", para a região Norte, 11=Rondônia, 12=Acre, 13=Amazonas, 14=Roraima, 15=Pará, 16=Amapá, 17=Tocantins; cada município tem um código específico; etc..) A única forma que tenho de saber se as tabelas que construí estão corretas é substituir os números pelas referências. Há algum meio simples de fazer isto? Obrigado pela atenção e paciência. Atenciosamente, Manoel Nascimento

1 0

gráfico de barras
by Vitor Aguiar 14 Ago '12

14 Ago '12

Boa tarde, gostaria de saber como fazer um simples gráfico de barras com o dataframe abaixo, onde as colunas 2 e 3 seriam barras consecutivas, com o nome na coluna 1 para esse grupo de 2 barras. Names NAgel NAcap 1 PentaE 0.14970140 0.09954859 2 D18S51 0.06891474 0.05968163 3 D21S11 0.13054984 0.11116655 4 TH01 0.26187191 0.07833215 5 D3S1358 0.06763797 0.05561891 6 FGA 0.08232084 0.08612497 Muito Obrigado, Vitor

3 4

Re: [R-br] Operação analitica de matrizes.
by Ivan Bezerra Allaman 13 Ago '12

13 Ago '12

Talvez "Ryacas " package? \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0

Re: [R-br] Operação analitica de matrizes.
by Ivan Bezerra Allaman 13 Ago '12

13 Ago '12

Uma googlada iria lhe mostrar o seguinte link: "matrix algebra r" http://personality-project.org/r/sem.appendix.1.pdf (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

3 3