Inciso 8

Se importa el mapa y los datos:

mapa_rep<-readOGR("C:\\Users\\richa\\OneDrive - CIDE\\1er Semestre-DESKTOP-RPQUJJA\\Manejo de bases de datos\\Actividad 4", layer="republica",encoding = "UTF-8")
## OGR data source with driver: ESRI Shapefile 
## Source: "C:\Users\richa\OneDrive - CIDE\1er Semestre-DESKTOP-RPQUJJA\Manejo de bases de datos\Actividad 4", layer: "republica"
## with 32 features
## It has 3 fields
datos12<-  read_excel("C:\\Users\\richa\\OneDrive - CIDE\\1er Semestre-DESKTOP-RPQUJJA\\Manejo de bases de datos\\Actividad 4\\Coneval.xlsx", 
                      sheet= 8, 
                      range= anchored("J4", dim = c(32,1)), 
                      col_names = c("Educacion"), na = "*")

mapa_rep@data$Educacion <- as.numeric(datos12$Educacion)

Se generan los cuantiles:

cut(mapa_rep@data$Educacion,4)
##  [1] (15.3,21]   (15.3,21]   (9.49,15.3] (15.3,21]   (9.49,15.3] (15.3,21]  
##  [7] (26.8,32.5] (15.3,21]   (9.49,15.3] (15.3,21]   (21,26.8]   (21,26.8]  
## [13] (15.3,21]   (15.3,21]   (9.49,15.3] (26.8,32.5] (15.3,21]   (15.3,21]  
## [19] (9.49,15.3] (26.8,32.5] (21,26.8]   (15.3,21]   (15.3,21]   (15.3,21]  
## [25] (15.3,21]   (15.3,21]   (15.3,21]   (15.3,21]   (15.3,21]   (26.8,32.5]
## [31] (21,26.8]   (21,26.8]  
## Levels: (9.49,15.3] (15.3,21] (21,26.8] (26.8,32.5]

Se genera el mapa con los cuantiles y una paleta de colores:

my_colors_rep_edu <- brewer.pal(5, "PuBuGn")
my_colors_rep_edu <- colorRampPalette(my_colors_rep_edu)(4)
cuantil_rep_edu <- cut(mapa_rep@data$Educacion,4)
my_colors_rep_edu <- my_colors_rep_edu[as.numeric(cuantil_rep_edu)]
plot(mapa_rep, col=my_colors_rep_edu, bg = "white")

Se agrega la leyenda para identificar los cuantiles:

my_colors_rep_edu <- brewer.pal(5, "PuBuGn")
my_colors_rep_edu <- colorRampPalette(my_colors_rep_edu)(4)
cuantil_rep_edu <- cut(mapa_rep@data$Educacion,4)
my_colors_rep_edu <- my_colors_rep_edu[as.numeric(cuantil_rep_edu)]
plot(mapa_rep, col=my_colors_rep_edu, bg = "white")

texto_rep_edu <- c("9.49%-15.30%", "15.30%-21.00%", "21.00%-26.80%", "26.80%-32.50%")
legend("topright",texto_rep_edu, cex = .5, bty = "o", 
       col = c("#BDC9E1", "#67A9CF", "#1C9099", "#016C59"), 
       pch = 19)

Para hacerlo interactivo:

textos_rep_edu <- paste(
  "Estado : ",mapa_rep$NOMGEO ,"<br/>",
  "% Pob. con rezago educativo: ", round(mapa_rep$Educacion,2) ) %>% lapply(htmltools::HTML)

cortes_rep_edu<- c(9.49,15.3,21.0,26.8,32.5,Inf)
colores_rep_edu<- colorBin( palette="PuBuGn", domain=mapa_rep$Educacion , na.color="transparent", bins=cortes_rep_edu) 

leaflet(data=mapa_rep) %>% 
  addTiles() %>% 
  addPolygons(label = textos_rep_edu,fillColor = colores_rep_edu(mapa_rep$Educacion),
              fillOpacity = 0.9)