Para poder trabajar con los mapas, primero cargo las paqueterías: “rgdal”, “leaflet”, “dplyr”, “htmltools”, “shiny” y “RColorBrewer”.
library(rgdal)
library(leaflet)
library(RColorBrewer)
library(dplyr)
library(htmltools)
library (shiny)Después, importo el mapa desde “QGIS” al ambiente de R.
plot(mapa_oaxaca)Utilizando el reporte “ANEXO ESTADÍSTICO DE POBREZA A NIVEL MUNICIPIO 2010 Y 2015” ( ANEXO ), elegir un estado (y sus municipios), y representar el nivel de pobreza y pobreza extrema para los años 2010 y 2015.
OAXACA
Ahora, cargo el archivo: “ANEXO ESTADÍSTICO DE POBREZA A NIVEL MUNICIPIO 2010 Y 2015”, de excel a R.
library(readxl)
Indicadores_pobreza<- read_excel("Concentrado, indicadores de pobreza.xlsx",
sheet = "Oax")Para el año 2010, considerando la división (cut) por cuantiles, representar en un mapa la distribución de la pobreza de los municipios. (10%)
Indicadores_pobreza$PropPob<-as.numeric(Indicadores_pobreza$PobrezaPorcentaje2010)
tabla_oax_2010<-merge(x=mapa_oaxaca@data, y=Indicadores_pobreza, by.x="CVEGEO", by.y="Clave", sort =FALSE)
mapa_oaxaca@data$PobProp<-tabla_oax_2010$PobrezaPorcentaje2010
cuantil_oax<-cut(mapa_oaxaca@data$PobProp,5)
cortes_cuantiles_oax<- c(0, 20.1, 51.7, 83.2, Inf)
colores_<-brewer.pal(4, "YlOrRd")
colores <- colorBin( palette="YlOrRd", domain=mapa_oaxaca$PobProp, na.color="black", bins=cortes_cuantiles_oax)
textos <- paste(
"Municipio: ", mapa_oaxaca$NOMGEO, "<br/>",
"% Pobreza en 2010: ", round(mapa_oaxaca$PobProp,2)) %>% lapply(htmltools::HTML)
leaflet(data=mapa_oaxaca) %>%
addTiles() %>%
addPolygons(label = textos, fillColor = colores(mapa_oaxaca$PobProp),
fillOpacity = 0.7) %>%
addLegend(color=colores_, labels=c("0-20.1%", "20.1-51.7%", "51.7-83.2%", "83.2-99.1%"))