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. >
_______________________________________________ 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.