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==