Autor : Diego Cubas Chavez


1. Cargado de Bases de Datos

  • Carga de bases de Datos desde los años 2011 a 2022

Los archivos en formato .dta descargados de INEI fueron convertidos a formato .csv y subidos a un repositorio de Github Con el fin de lograr ejecutar el código sin necesidad de configurar los archivos en la carpeta local. Solo se necesita conexión a internet de cualquier usuario. La data es integra luego del preprocesamiento de datos, se contrastó con los resultados obtenidos en Stata.

  • Código Utilizado para convertir los archivos de formato .dta a formato .csv

data2011to2022 <- dir(“ruta de la carpeta”,full.names=T) %>% map_df(read_dta,encoding=‘latin1’)

data2011to2022 %>% select(aÑo,conglome,estrato,factor07, mieperho,gashog2d,linea,linpe,pobreza) %>% mutate(aÑo=replace(aÑo,is.na(aÑo),2013))

write.csv(data2011to2022,“sumaria-2011a2022.csv”,row.names = FALSE)

base_datos1<-read.csv("https://raw.githubusercontent.com/YoDaN8/YoDaN8/main/sumaria-2011a2016.csv" , encoding = 'latin1')
base_datos2<-read.csv("https://raw.githubusercontent.com/YoDaN8/YoDaN8/main/sumaria-2017a2022.csv" , encoding = 'latin1')
base_datos3<-read.csv("https://raw.githubusercontent.com/YoDaN8/YoDaN8/main/sumaria-2021a2022.csv" , encoding = 'latin1')

La data se partió en 2 , desde 2011 a 2016 y desde 2017 a 2022 , dado que hay límite de 25MB al subir un archivo.


2. Metodología

1. Calculando la tasa de probreza monetaria y extrema

basetotal <- rbind(base_datos1,base_datos2)
matriz_poblacion = matrix(nrow = 2 , ncol = 12)
i=2011
for (a in 2011:2022){
  # Filtramos la base_datos1 por año para obtener la base de datos sumaria-"año"
  bd<-basetotal %>% filter(aÑo==i)
  # Creamos el factor de expansión
  bd$facpob<-bd$factor07*bd$mieperho
  # Creamos la variable diseño_i, que indica la complejidad de la encuesta 
  diseño_bd <- svydesign(ids=~1, weights = ~bd$facpob, data = bd)
  # Creamos una tabla de la variable pobreza teniendo en cuenta el diseño de la encuesta
  tab_cruzada <- svytable(~pobreza, diseño_bd)
  # Transformamos de tabla a matriz
  matriz_t<-as.matrix(tab_cruzada)
  # Calculando la población en situación de pobreza extrema, pobreza y total 
  Pobre_ext<-matriz_t[1,1]
  Pobre_mon<-Pobre_ext+matriz_t[2,1]
  Total<-Pobre_mon+matriz_t[3,1]
  # Calculando la tasa de probreza monetaria y extrema
  Tasa_pobreza_monetaria<-round(Pobre_mon/Total*100,digits = 2)
  Tasa_pobreza_extrema<-round(Pobre_ext/Total*100,digits=2)
  # Creamos una matriz con las dos tasas de pobreza y cambiamos nombre de filas y columna
  matriz_pob <- matrix(data = c(Tasa_pobreza_monetaria, Tasa_pobreza_extrema),2,1)
  row.names(matriz_pob)<-c("Tasa_pobreza_monetaria","Tasa_pobreza_extrema")
  colnames(matriz_pob)<-a
  matriz_poblacion[,i-2010] <- matriz_pob
  i<-i+1
}
colnames(matriz_poblacion) <- c(paste('Año',2011:2022))
categorias <- data.frame(Tasas = c("Tasa_pobreza_monetaria","Tasa_pobreza_extrema" ))
matriz_final <- cbind(categorias,matriz_poblacion)

2. Prueba estadística (Welch)

Diferencia entre las medias de la pobreza monetaria y media de la pobreza monetaria extrema

#Creando variables
base_datos3$Pob_ext<-ifelse(base_datos3$pobreza==1,1,0)
base_datos3$Pob_mon<-ifelse(base_datos3$pobreza==3,0,1)
#Prueba estadística de la diferenci de la tasa de pobreza total
test <- t.test(c(filter(base_datos3,aÑo=="2021")$Pob_mon),c(filter(base_datos3,aÑo=="2022")$Pob_mon))
test_format <- tidy(test)

3. Resultados

1. Tasa de pobreza Monetaria y Tasa de pobreza extrema en Perú , desde el año 2011 al año 2022

Tasas Año 2011 Año 2012 Año 2013 Año 2014 Año 2015 Año 2016 Año 2017 Año 2018 Año 2019 Año 2020 Año 2021 Año 2022
Tasa de pobreza monetaria 27.82 25.81 23.91 22.73 21.77 20.74 21.7 20.49 20.19 30.13 25.87 27.52
Tasa de pobreza extrema 6.34 6.01 4.73 4.28 4.07 3.77 3.79 2.81 2.85 5.15 4.12 5.01

2. Prueba t de dos muestras de Welch

Estimado 1 Estimado 2 Estimado 3 Estadístico p-value Parametros Lim. Inf Lim. Sup
-0.02 0.19 0.21 c(t = -6.63) 0 c(df = 68347.08) -0.03 -0.01

Para ambas pruebas estadísticas se obtiene un p-valor igual a 0.000, que es menor al 0.05, por tanto, se rechaza la hipótesis nula de diferencia no significativa y se acepta la hipótesis alternativa de diferenci significativa. Podemos decir que tanto la tasa de pobreza y pobreza extrema en el 2021 son significativamente diferentes para los valores en el año 2022, bajo un nivel de significancia del 5%,lo mismo que un nivel de confianza del 95%.

LS0tDQp0aXRsZTogIkFuw6FsaXNpcyBkZSBwb2JyZXphIG1vbmV0YXJpYSB5IHBvYnJlemEgZXh0cmVtYSBlbiBQZXLDuiAsIGRlc2RlIGxvcyBhw7FvcyAyMDExIGEgMjAyMiINCm91dHB1dDogb2lsYWJzOjpsYWJfcmVwb3J0DQotLS0NCg0KKipBdXRvciA6KiogW0RpZWdvIEN1YmFzIENoYXZlel0oaHR0cHM6Ly93d3cubGlua2VkaW4uY29tL2luL2RpZWdvY3ViYXNjaGF2ZXovKQ0KDQpgYGB7ciAsIGluY2x1ZGU9RkFMU0V9DQojQ29ycmVyIGVzdGUgY29kaWdvIGZ1ZXJhIGRlIFJNRCAsIHBhcmEgcG9kZXIgdXNhciBsYSBsaWJyZXJpYSBvaWxhYnMgZGUgR2l0SHViDQojaW5zdGFsbC5wYWNrYWdlcygicmVtb3RlcyIpDQojcmVtb3Rlczo6aW5zdGFsbF9naXRodWIoImFuZHJld3BicmF5L29pbGFicyIpDQpgYGANCg0KYGBge3Igc2V0dXAsIGluY2x1ZGU9RkFMU0V9DQprbml0cjo6b3B0c19jaHVuayRzZXQoZWNobyA9IFRSVUUpDQoNCmxpYnJhcnkoImRwbHlyIikNCmxpYnJhcnkoIm9pbGFicyIpDQpsaWJyYXJ5KCdicm9vbScpDQpsaWJyYXJ5KCdzdXJ2ZXknKQ0KbGlicmFyeSgnaHRtbFRhYmxlJykNCmBgYA0KDQotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0NCg0KIyMgKioxLiBDYXJnYWRvIGRlIEJhc2VzIGRlIERhdG9zKiogDQoNCi0gICAqKkNhcmdhIGRlIGJhc2VzIGRlIERhdG9zIGRlc2RlIGxvcyBhw7FvcyAyMDExIGEgMjAyMioqDQoNCkxvcyBhcmNoaXZvcyBlbiBmb3JtYXRvICouZHRhKiBkZXNjYXJnYWRvcyBkZSBbSU5FSV0oaHR0cHM6Ly9wcm95ZWN0b3MuaW5laS5nb2IucGUvbWljcm9kYXRvcy8pIGZ1ZXJvbiBjb252ZXJ0aWRvcyBhIGZvcm1hdG8gKi5jc3YqIHkgc3ViaWRvcyBhIHVuIHJlcG9zaXRvcmlvIGRlIFtHaXRodWJdKGh0dHBzOi8vZ2l0aHViLmNvbS9Zb0RhTjgvWW9EYU44KSBDb24gZWwgZmluIGRlIGxvZ3JhciBlamVjdXRhciBlbCBjw7NkaWdvIHNpbiBuZWNlc2lkYWQgZGUgY29uZmlndXJhciBsb3MgYXJjaGl2b3MgZW4gbGEgY2FycGV0YSBsb2NhbC4gU29sbyBzZSBuZWNlc2l0YSBjb25leGnDs24gYSBpbnRlcm5ldCBkZSBjdWFscXVpZXIgdXN1YXJpby4gTGEgZGF0YSBlcyBpbnRlZ3JhIGx1ZWdvIGRlbCBwcmVwcm9jZXNhbWllbnRvIGRlIGRhdG9zLCBzZSBjb250cmFzdMOzIGNvbiBsb3MgcmVzdWx0YWRvcyBvYnRlbmlkb3MgZW4gU3RhdGEuDQoNCi0gICAqKkPDs2RpZ28gVXRpbGl6YWRvIHBhcmEgY29udmVydGlyIGxvcyBhcmNoaXZvcyBkZSBmb3JtYXRvIC5kdGEgYSBmb3JtYXRvIC5jc3YqKg0KDQoqKipkYXRhMjAxMXRvMjAyMiBcPC0gZGlyKCJydXRhIGRlIGxhIGNhcnBldGEiLGZ1bGwubmFtZXM9VCkgJVw+JSBtYXBfZGYocmVhZF9kdGEsZW5jb2Rpbmc9J2xhdGluMScpKioqDQoNCioqKmRhdGEyMDExdG8yMDIyICVcPiUgc2VsZWN0KGHDkW8sY29uZ2xvbWUsZXN0cmF0byxmYWN0b3IwNywgbWllcGVyaG8sZ2FzaG9nMmQsbGluZWEsbGlucGUscG9icmV6YSkgJVw+JSBtdXRhdGUoYcORbz1yZXBsYWNlKGHDkW8saXMubmEoYcORbyksMjAxMykpKioqDQoNCioqKndyaXRlLmNzdihkYXRhMjAxMXRvMjAyMiwic3VtYXJpYS0yMDExYTIwMjIuY3N2Iixyb3cubmFtZXMgPSBGQUxTRSkqKioNCg0KYGBge3J9DQpiYXNlX2RhdG9zMTwtcmVhZC5jc3YoImh0dHBzOi8vcmF3LmdpdGh1YnVzZXJjb250ZW50LmNvbS9Zb0RhTjgvWW9EYU44L21haW4vc3VtYXJpYS0yMDExYTIwMTYuY3N2IiAsIGVuY29kaW5nID0gJ2xhdGluMScpDQpiYXNlX2RhdG9zMjwtcmVhZC5jc3YoImh0dHBzOi8vcmF3LmdpdGh1YnVzZXJjb250ZW50LmNvbS9Zb0RhTjgvWW9EYU44L21haW4vc3VtYXJpYS0yMDE3YTIwMjIuY3N2IiAsIGVuY29kaW5nID0gJ2xhdGluMScpDQpiYXNlX2RhdG9zMzwtcmVhZC5jc3YoImh0dHBzOi8vcmF3LmdpdGh1YnVzZXJjb250ZW50LmNvbS9Zb0RhTjgvWW9EYU44L21haW4vc3VtYXJpYS0yMDIxYTIwMjIuY3N2IiAsIGVuY29kaW5nID0gJ2xhdGluMScpDQoNCmBgYA0KDQpMYSBkYXRhIHNlIHBhcnRpw7MgZW4gMiAsIGRlc2RlIDIwMTEgYSAyMDE2IHkgZGVzZGUgMjAxNyBhIDIwMjIgLCBkYWRvIHF1ZSBoYXkgbMOtbWl0ZSBkZSAyNU1CIGFsIHN1YmlyIHVuIGFyY2hpdm8uDQoNCi0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQ0KDQojIyAqKjIuIE1ldG9kb2xvZ8OtYSoqDQoNCiMjIyAqKjEuIENhbGN1bGFuZG8gbGEgdGFzYSBkZSBwcm9icmV6YSBtb25ldGFyaWEgeSBleHRyZW1hKioNCg0KYGBge3J9DQpiYXNldG90YWwgPC0gcmJpbmQoYmFzZV9kYXRvczEsYmFzZV9kYXRvczIpDQptYXRyaXpfcG9ibGFjaW9uID0gbWF0cml4KG5yb3cgPSAyICwgbmNvbCA9IDEyKQ0KaT0yMDExDQpmb3IgKGEgaW4gMjAxMToyMDIyKXsNCiAgIyBGaWx0cmFtb3MgbGEgYmFzZV9kYXRvczEgcG9yIGHDsW8gcGFyYSBvYnRlbmVyIGxhIGJhc2UgZGUgZGF0b3Mgc3VtYXJpYS0iYcOxbyINCiAgYmQ8LWJhc2V0b3RhbCAlPiUgZmlsdGVyKGHDkW89PWkpDQogICMgQ3JlYW1vcyBlbCBmYWN0b3IgZGUgZXhwYW5zacOzbg0KICBiZCRmYWNwb2I8LWJkJGZhY3RvcjA3KmJkJG1pZXBlcmhvDQogICMgQ3JlYW1vcyBsYSB2YXJpYWJsZSBkaXNlw7FvX2ksIHF1ZSBpbmRpY2EgbGEgY29tcGxlamlkYWQgZGUgbGEgZW5jdWVzdGEgDQogIGRpc2XDsW9fYmQgPC0gc3Z5ZGVzaWduKGlkcz1+MSwgd2VpZ2h0cyA9IH5iZCRmYWNwb2IsIGRhdGEgPSBiZCkNCiAgIyBDcmVhbW9zIHVuYSB0YWJsYSBkZSBsYSB2YXJpYWJsZSBwb2JyZXphIHRlbmllbmRvIGVuIGN1ZW50YSBlbCBkaXNlw7FvIGRlIGxhIGVuY3Vlc3RhDQogIHRhYl9jcnV6YWRhIDwtIHN2eXRhYmxlKH5wb2JyZXphLCBkaXNlw7FvX2JkKQ0KICAjIFRyYW5zZm9ybWFtb3MgZGUgdGFibGEgYSBtYXRyaXoNCiAgbWF0cml6X3Q8LWFzLm1hdHJpeCh0YWJfY3J1emFkYSkNCiAgIyBDYWxjdWxhbmRvIGxhIHBvYmxhY2nDs24gZW4gc2l0dWFjacOzbiBkZSBwb2JyZXphIGV4dHJlbWEsIHBvYnJlemEgeSB0b3RhbCANCiAgUG9icmVfZXh0PC1tYXRyaXpfdFsxLDFdDQogIFBvYnJlX21vbjwtUG9icmVfZXh0K21hdHJpel90WzIsMV0NCiAgVG90YWw8LVBvYnJlX21vbittYXRyaXpfdFszLDFdDQogICMgQ2FsY3VsYW5kbyBsYSB0YXNhIGRlIHByb2JyZXphIG1vbmV0YXJpYSB5IGV4dHJlbWENCiAgVGFzYV9wb2JyZXphX21vbmV0YXJpYTwtcm91bmQoUG9icmVfbW9uL1RvdGFsKjEwMCxkaWdpdHMgPSAyKQ0KICBUYXNhX3BvYnJlemFfZXh0cmVtYTwtcm91bmQoUG9icmVfZXh0L1RvdGFsKjEwMCxkaWdpdHM9MikNCiAgIyBDcmVhbW9zIHVuYSBtYXRyaXogY29uIGxhcyBkb3MgdGFzYXMgZGUgcG9icmV6YSB5IGNhbWJpYW1vcyBub21icmUgZGUgZmlsYXMgeSBjb2x1bW5hDQogIG1hdHJpel9wb2IgPC0gbWF0cml4KGRhdGEgPSBjKFRhc2FfcG9icmV6YV9tb25ldGFyaWEsIFRhc2FfcG9icmV6YV9leHRyZW1hKSwyLDEpDQogIHJvdy5uYW1lcyhtYXRyaXpfcG9iKTwtYygiVGFzYV9wb2JyZXphX21vbmV0YXJpYSIsIlRhc2FfcG9icmV6YV9leHRyZW1hIikNCiAgY29sbmFtZXMobWF0cml6X3BvYik8LWENCiAgbWF0cml6X3BvYmxhY2lvblssaS0yMDEwXSA8LSBtYXRyaXpfcG9iDQogIGk8LWkrMQ0KfQ0KY29sbmFtZXMobWF0cml6X3BvYmxhY2lvbikgPC0gYyhwYXN0ZSgnQcOxbycsMjAxMToyMDIyKSkNCmNhdGVnb3JpYXMgPC0gZGF0YS5mcmFtZShUYXNhcyA9IGMoIlRhc2FfcG9icmV6YV9tb25ldGFyaWEiLCJUYXNhX3BvYnJlemFfZXh0cmVtYSIgKSkNCm1hdHJpel9maW5hbCA8LSBjYmluZChjYXRlZ29yaWFzLG1hdHJpel9wb2JsYWNpb24pDQpgYGANCg0KLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tDQoNCiMjIyAqKjIuIFBydWViYSBlc3RhZMOtc3RpY2EgKFdlbGNoKSoqDQoNCioqRGlmZXJlbmNpYSBlbnRyZSBsYXMgbWVkaWFzIGRlIGxhIHBvYnJlemEgbW9uZXRhcmlhIHkgbWVkaWEgZGUgbGEgcG9icmV6YSBtb25ldGFyaWEgZXh0cmVtYSoqDQoNCmBgYHtyfQ0KI0NyZWFuZG8gdmFyaWFibGVzDQpiYXNlX2RhdG9zMyRQb2JfZXh0PC1pZmVsc2UoYmFzZV9kYXRvczMkcG9icmV6YT09MSwxLDApDQpiYXNlX2RhdG9zMyRQb2JfbW9uPC1pZmVsc2UoYmFzZV9kYXRvczMkcG9icmV6YT09MywwLDEpDQojUHJ1ZWJhIGVzdGFkw61zdGljYSBkZSBsYSBkaWZlcmVuY2kgZGUgbGEgdGFzYSBkZSBwb2JyZXphIHRvdGFsDQp0ZXN0IDwtIHQudGVzdChjKGZpbHRlcihiYXNlX2RhdG9zMyxhw5FvPT0iMjAyMSIpJFBvYl9tb24pLGMoZmlsdGVyKGJhc2VfZGF0b3MzLGHDkW89PSIyMDIyIikkUG9iX21vbikpDQp0ZXN0X2Zvcm1hdCA8LSB0aWR5KHRlc3QpDQpgYGANCg0KLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tDQoNCiMjICoqMy4gUmVzdWx0YWRvcyoqDQoNCiMjIyAqKjEuIFRhc2EgZGUgcG9icmV6YSBNb25ldGFyaWEgeSBUYXNhIGRlIHBvYnJlemEgZXh0cmVtYSBlbiBQZXLDuiAsIGRlc2RlIGVsIGHDsW8gMjAxMSBhbCBhw7FvIDIwMjIqKg0KDQpgYGB7ciAsIGVjaG89RkFMU0V9DQpjb2xuYW1lcyhtYXRyaXpfcG9ibGFjaW9uKSA8LSBjKHBhc3RlKCdBw7FvJywyMDExOjIwMjIpKQ0KY2F0ZWdvcmlhcyA8LSBkYXRhLmZyYW1lKFRhc2FzID0gYygiVGFzYSBkZSBwb2JyZXphIG1vbmV0YXJpYSIsIlRhc2EgZGUgcG9icmV6YSBleHRyZW1hIiApKQ0KbWF0cml6X2ZpbmFsIDwtIGNiaW5kKGNhdGVnb3JpYXMsbWF0cml6X3BvYmxhY2lvbikNCmh0bWxUYWJsZTo6aHRtbFRhYmxlKG1hdHJpel9maW5hbCAsIHJuYW1lcyA9IEYgKQ0KYGBgDQoNCiMjIyAqKjIuIFBydWViYSB0IGRlIGRvcyBtdWVzdHJhcyBkZSBXZWxjaCoqDQoNCmBgYHtyICwgZWNobz1GfQ0KdGVzdF9mb3JtYXQgPC0gcm91bmQodGVzdF9mb3JtYXRbLDE6OF0sMikNCmNvbG5hbWVzKHRlc3RfZm9ybWF0KSA8LSBjKCJFc3RpbWFkbyAxIiAsICJFc3RpbWFkbyAyIiAsIkVzdGltYWRvIDMiLCAiRXN0YWTDrXN0aWNvIiAsInAtdmFsdWUiICwgICJQYXJhbWV0cm9zIiAsIkxpbS4gSW5mIiAsIkxpbS4gU3VwIikNCmh0bWxUYWJsZTo6aHRtbFRhYmxlKHRlc3RfZm9ybWF0ICwgcm5hbWVzID0gRikgDQpgYGANCg0KKipQYXJhIGFtYmFzIHBydWViYXMgZXN0YWTDrXN0aWNhcyBzZSBvYnRpZW5lIHVuIHAtdmFsb3IgaWd1YWwgYSAwLjAwMCwgcXVlIGVzIG1lbm9yIGFsIDAuMDUsIHBvciB0YW50bywgc2UgcmVjaGF6YSBsYSBoaXDDs3Rlc2lzIG51bGEgZGUgZGlmZXJlbmNpYSBubyBzaWduaWZpY2F0aXZhIHkgc2UgYWNlcHRhIGxhIGhpcMOzdGVzaXMgYWx0ZXJuYXRpdmEgZGUgZGlmZXJlbmNpIHNpZ25pZmljYXRpdmEuIFBvZGVtb3MgZGVjaXIgcXVlIHRhbnRvIGxhIHRhc2EgZGUgcG9icmV6YSB5IHBvYnJlemEgZXh0cmVtYSBlbiBlbCAyMDIxIHNvbiBzaWduaWZpY2F0aXZhbWVudGUgZGlmZXJlbnRlcyBwYXJhIGxvcyB2YWxvcmVzIGVuIGVsIGHDsW8gMjAyMiwgYmFqbyB1biBuaXZlbCBkZSBzaWduaWZpY2FuY2lhIGRlbCA1JSxsbyBtaXNtbyBxdWUgdW4gbml2ZWwgZGUgY29uZmlhbnphIGRlbCA5NSUuKioNCg==