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.