Marcelo, Seu código não atende à diretriz deste grupo em Guia de Postagem no R-br <http://www.leg.ufpr.br/r-br-guia>, em particular, a solicitação de fornecer o CMR (*código mínimo reproduzível*). As chances de você obter respostas são maiores quando essa recomendação é seguida. HTH On Sun, Mar 29, 2020 at 9:05 AM Marcelo Magalhães por (R-br) < r-br@listas.c3sl.ufpr.br> wrote:
Prezados e Prezadas,
Estou trabalhando com o IVS (Índice de Vulnerabilidade Social), mais precisamente com as UDHs (Unidades de Habitação) para o Município de Jaboatão dos Guararapes - PE. A base Completa + Shapefiles estão no site do Ipea (http://ivs.ipea.gov.br/index.php/pt/biblioteca). Gostaria de fazer um mapa do Índice de Vulnerabilidade Social do Município para saber quais as localidades mais vulneráveis.
Tentei unir pelo código do município e pela variável UDH_Atlas no shapefile e udh e travei.
Ao que parece:
1. UDHs são colocados como "strings" pelo R, ao colocar a função as.numeric(), há uma mudança nas linhas.
2. UDHs são coordenadas únicas (não se separam em Long. e Lat.), ou temos que separá-las em Long e Lat?
Agradeço desde já a ajuda,
Marcelo Nery Magalhães
Script:
#Bibliotecas
library(ggplot2)
library(rgdal)
#Shape
Jab<-readOGR("C:/Users/usuario/Documents/IVS", "UDHs_RM_Todas")
head(Jab@data)
#Transformando string em numérico
Jab$UDH_ATLAS<-as.numeric(Jab$UDH_ATLAS)
#Base de dados/tabela
ivs<-read.csv("C:/Users/usuario/Documents/IVS/jab.csv", sep = ";", head=T, dec = ",")
head(ivs)
ivs<-na.omit(ivs)
names(Jab@data)
names(ivs)
#Transformando os nomes das variáveis em um só para cada objeto
names(ivs)[12]<-c("Cod_mun_6")
names(Jab)[1]<-c("udh") # tentei criar uma chave com essa variável, já que com o código do município não estava funcionando
names(Jab)
names(ivs)
#Ordenando
ivs<-ivs[order(ivs$Cod_mun_6),]
Jab<-Jab@data[order(Jab$Cod_mun_6),]
head(ivs)
head(Jab)
#Combinando através da Chave do Código do Município
j<-merge(Jab,ivs,by.y="Cod_mun_6") # Houve problema ao criar “j”
head(j)
dim(Jab)
dim(ivs)
dim(j)
fortify(Jab)
j.f<-fortify(Jab,region="udh")
j$ivs
#Categorizando
j$ivscat<-cut(j$ivs, breaks = c(0,0.200,0.300,0.400,0.500,1), labels=c("0 a 0,200 - Muito Baixa", "0,200 a 0,300 - Baixa", "0,300 a 0,400 - Média", "0,400 a 0,500 - Alta", "0,500 a 1 - Muito Alta"), include.lowest=T)
names(j)
names(j.f)
#Problema: udh não possui latitude e lobgitude (como fazer, então?)
jabo<-merge(j, j.f, by.y="udh")
_______________________________________________ R-br mailing list R-br@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.