7.1 Procedimiento de estimación

Este código utiliza las librerías tidyverse y magrittr para procesamiento y analizar datos.

La función readRDS() es utilizada para cargar un archivo de datos en formato RDS, que contiene las estimaciones directas y la varianza suvizada para la proporción de personas en condición de pobreza correspondientes al año 2018. Luego, se utiliza el operador %>% de la librería magrittr para encadenar la selección de las columnas de interés, que corresponden a los nombres dam2, nd, pobreza, vardir y hat_var.

library(tidyverse)
library(magrittr)

base_FH <- readRDS("Recursos/Día2/Sesion3/Data/base_FH_2018.rds") %>% 
  select(dam2, nd,  pobreza, vardir, hat_var)

Lectura de las covariables, las cuales son obtenidas previamente. Dado la diferencia entre las escalas de las variables es necesario hacer un ajuste a estas.

statelevel_predictors_df <- readRDS("Recursos/Día2/Sesion3/Data/statelevel_predictors_df_dam2.rds") %>% 
    mutate_at(.vars = c("luces_nocturnas",
                      "cubrimiento_cultivo",
                      "cubrimiento_urbano",
                      "modificacion_humana",
                      "accesibilidad_hospitales",
                      "accesibilidad_hosp_caminado"),
            function(x) as.numeric(scale(x)))

Ahora, se realiza una unión completa (full_join) entre el conjunto de datos base_FH y los predictores statelevel_predictors_df utilizando la variable dam2 como clave de unión.

Se utiliza la función tba() para imprimir las primeras 10 filas y 8 columnas del conjunto de datos resultante de la unión anterior.

La unión completa (full_join) combina los datos de ambos conjuntos, manteniendo todas las filas de ambos, y llenando con valores faltantes (NA) en caso de no encontrar coincidencias en la variable de unión (dam2 en este caso).

La función tba() imprime una tabla en formato HTML en la consola de R que muestra las primeras 10 filas y 8 columnas del conjunto de datos resultante de la unión.

base_FH <- full_join(base_FH, statelevel_predictors_df, by = "dam2" )
tba(base_FH[1:10,1:8])

dam2	nd	pobreza	vardir	hat_var	modificacion_humana	accesibilidad_hospitales	accesibilidad_hosp_caminado
00101	6796	0.2225	0.0004	0.0004	3.6127	-1.1835	-1.5653
00201	531	0.1822	0.0004	0.0039	-0.0553	0.4449	0.2100
00206	230	0.3366	0.0031	0.0041	0.5157	-0.1468	-0.1811
00301	666	0.4266	0.0043	0.0050	0.1364	0.5744	1.1660
00302	261	0.4461	0.0014	0.0029	-0.5103	0.2531	1.0880
00303	566	0.5587	0.0142	0.0075	-0.6591	0.6249	1.2229
00304	412	0.5406	0.0116	0.0042	-0.5573	1.4586	2.7337
00401	1219	0.3359	0.0010	0.0042	0.3979	-0.0833	-0.4490
00402	172	0.1496	0.0007	0.0047	-0.3661	-0.0114	-0.2863
00403	309	0.4644	0.0015	0.0031	-1.0446	0.4542	0.5702

# View(base_FH)