10.6 Mapa de informalidad laboral

El siguiente bloque de código carga los paquetes sp, sf y tmap y realiza una serie de operaciones. En primer lugar, une (rbind) las estimaciones de los dominios observados y no observados (data_dir, data_syn) y selecciona las variables dam2, pobreza (informalidad laboral), pred_arcoseno, pred_arcoseno_EE y Cv_pred utilizando la función select(). A continuación, lee un archivo Shapefile que contiene información geoespacial del país. Luego, crea un mapa temático (tmap) utilizando la función tm_shape() y agregando capas con la función tm_polygons(). El mapa representa dos variables llamadas pobreza y pred_arcoseno, utilizando una paleta de colores llamada “YlOrRd” y establece los cortes de los intervalos de las variables con la variable brks_lp. Finalmente, la función tm_layout() establece algunos parámetros de diseño del mapa, como la relación de aspecto (asp).

library(sp)
library(sf)
library(tmap)

data_map <- rbind(data_dir, data_syn) %>% 
  select(dam2, Rd, pred_arcoseno, pred_arcoseno_EE,Cv_pred ) 


## Leer Shapefile del país
ShapeSAE <- read_sf("Recursos/Día3/Sesion3/Shape/DOM_dam2.shp") %>% 
   rename(dam2 = id_dominio) 

mapa <- tm_shape(ShapeSAE %>%
                   left_join(data_map,  by = "dam2"))

brks_lp <- quantile(data_map$pred_arcoseno)
tmap_options(check.and.fix = TRUE)
Mapa_lp <-
  mapa + tm_polygons(
    c("Rd", "pred_arcoseno"),
    breaks = brks_lp,
    title = "Mapa de tasa\nde informalidad",
    palette = "YlOrRd",
    colorNA = "white"
  ) + tm_layout(asp = 2.5)

tmap_save(
  Mapa_lp,
  "Recursos/Día3/Sesion3/0Recursos/Mapa_arcoseno.PNG",
  width = 3000,
  height = 2000,
  asp = 0
)

Mapa_lp