Path: blob/main/Trabajo_grupal/WG6/Grupo_9_r_Pregunta1.R
2714 views
### Pregunta 1 #### Grupo 9 ####12pacman::p_load(haven,dplyr, stringr)345"1.0 Set Directorio"67user <- Sys.getenv("USERNAME") # username89setwd( paste0("C:/Users/",user,"/Documents/GitHub/1ECO35_2022_2/Grupo9") )101112"2.0 Cargar dataset de ENAHO"1314enaho01 <- read_dta("../../enaho/2020/737-Modulo01/737-Modulo01/enaho01-2020-100.dta")15sumaria <- read_dta("../../enaho/2020/737-Modulo34/737-Modulo34/sumaria-2020.dta")161718"3.0 Juntar bases"1920enaho_2020 <- merge(enaho01, sumaria,21by = c("conglome", "vivienda", "hogar"),22all.x = T, suffixes = c("",".y")23)24"4.0 Filtrar variables a usar"2526enaho_2020 <- enaho_2020[ ,c("conglome", "vivienda", "hogar" ,"ubigeo",27"aÑo" , "mieperho", "inghog1d",28"gashog2d", "ld") ]293031###### Repetir procedimiento con data 2019 ######32"2.0"33enaho01 <- read_dta("../../enaho/2019/687-Modulo01/687-Modulo01/enaho01-2019-100.dta")34sumaria <- read_dta("../../enaho/2019/687-Modulo34/687-Modulo34/sumaria-2019.dta")353637"3.0"38enaho_2019 <- merge(enaho01, sumaria,39by = c("conglome", "vivienda", "hogar",40"ubigeo", "dominio", "estrato"),41all.x = T, suffixes = c("",".y")42)4344"4.0"45enaho_2019 <- enaho_2019[ ,c("conglome", "vivienda", "hogar" ,"ubigeo",46"aÑo" , "mieperho", "inghog1d",47"gashog2d", "ld") ]4849###################################################5051"5.0 Append de base 2019 y 2020"52enaho_append <- bind_rows(enaho_2020, enaho_2019)5354#Año:55enaho_append <- enaho_append %>% rename("aniorec"="aÑo")56#Departamento:57enaho_append['dpto'] <- paste(as.numeric(substr(enaho_append$ubigeo, 1, 2)))58#Nota, se coloca as.numeric para convertirlo en numero y59#hacer merge con base de deflactor606162"6.0 Deflactar las variables"63#Deflactor espacial:64# ld de base6566#Deflactor temporal:67deflactor <- read_dta("../../enaho/2020/737-Modulo34/737-Modulo34/ConstVarGasto-Metodologia actualizada/Gasto2020/Bases/deflactores_base2020_new.dta")68deflactor <- deflactor[ , c("dpto", "aniorec", "i00")]6970#Unir deflactor a la base principal:71enaho <- merge(enaho_append, deflactor,72by = c( "dpto", "aniorec"),73all.x=T)747576"7.0 Dividir las variables de ingreso y gasto por mieperho, 12, ld e i00."7778enaho <- enaho %>%79mutate(80ingreso_month_pc = enaho$inghog1d/(enaho$mieperho* enaho$ld * enaho$i00*12),81gasto_month_pc = enaho$gashog2d/(enaho$mieperho* enaho$ld * enaho$i00*12)82)8384#Ingreso per cápita = ingreso_month_pc8586#Gasto per cápita = gasto_month_pc8788"8.0 Comprobación son los mismos valores que en Python"89enaho <- enaho[order(enaho$aniorec, enaho$conglome, enaho$vivienda, enaho$hogar),]90View(enaho)919293