Book a Demo!
CoCalc Logo Icon
StoreFeaturesDocsShareSupportNewsAboutPoliciesSign UpSign In
robertopucp
GitHub Repository: robertopucp/1eco35_2022_2
Path: blob/main/Trabajo_grupal/WG6/Grupo_5_r.R
2714 views
1
library(haven) # leer archivos spss, stata, dbf, etc
2
library(dplyr) # limpieza de datos
3
library(stringr) # grep for regular expression
4
library(fastDummies) # crear dummy
5
library(srvyr) # libreria para declarar el dise�o muestral de una encuesta
6
library(survey)
7
8
user <- Sys.getenv("fdcc0")
9
10
setwd( paste0("C:/Users/{user}/Documentos/GitHub/1ECO35_2022_2/lab7") )
11
12
13
14
15
#Pregunta 1
16
# Establecemos las bases que usaremos enaho01_19 y enaho34_19
17
18
enaho01_2019 <- read_dta("../../../datos/2019/687-Modulo01/687-Modulo01/enaho01-2019-100.dta")
19
20
enaho34_2019 <- read_dta("../../../datos/2020/737-Modulo01/737-Modulo01/enaho01-2020-100.dta")
21
22
#Establecemos las bases que usaremos enaho01_20y enaho34_20
23
24
enaho01_2020 <- read_dta("../../../datos/2020/737-Modulo01/737-Modulo01/enaho01-2020-100.dta")
25
26
enaho34_2020 <- read_dta("../../../datos/2020/737-Modulo34/737-Modulo34/sumaria-2020.dta")
27
28
# Establecemos la base de los deflactores
29
30
deflactor_temporal <- read_dta("../../../datos/2020/737-Modulo34/737-Modulo34/ConstVarGasto-Metodologia actualizada/Gasto2020/Bases/deflactores_base2020_new.dta")
31
32
# Hacemos el merge entre ambas
33
34
enaho_merge2019 <- merge(enaho34_2019, enaho01_2019,
35
by = c("conglome", "vivienda", "hogar"),
36
all.x = T
37
)
38
39
40
enaho_merge2020 <- merge(enaho34_2020, enaho01_2020,
41
by = c("conglome", "vivienda", "hogar"),
42
all.x = T
43
)
44
45
# Ahora, vamos a realizar el append
46
enaho_append <- append(enaho_merge2019, enaho_merge2020)
47
48
#Creamose el ingreso real mensual
49
#Creamose el gasto real menusal
50
enaho_append$ingreso_mensual <- enaho_append$inghog1d / (12*enaho_append$mieperho)
51
52
enaho_append$gasto_mensual <- enaho_append$gashog2d / (12*enaho_append$mieperho)
53
54
# deflactando las variables (deflactor espacial y temporal)
55
56
# espacial
57
enaho_append$ingreso_mensual_defl <- enaho_append$ingreso_mensual * enaho_append$ld
58
59
enaho_append$gasto_mensual_defl <- enaho_append$gasto_mensual * enaho_append$ld
60
61
62
63
64
65
66
#Pregunta 2
67
#el salario por hora del trabajador dependiente
68
69
enaho01_500 <- read_dta("../../../2020/737-Modulo05/737-Modulo05/enaho01a-2020-500.dta")
70
71
# Obtenemos el salario anual del primer y segundo empleo
72
73
enaho01_500$ingreso_anual <- enaho01_500$i524e1 + enaho01_500$i538e1
74
75
# Encontramos el n�mero de hrs trabajadas en la semana
76
77
enaho01_500$horas_trab_sem <- enaho01_500$i513t + enaho01_500$i518
78
79
# Encontramos el salario x hora del trabajador
80
81
enaho01_500$salarioxhora <- enaho01_500$ingreso_anual / (enaho01_500$horas_trab_sem*52)
82
83
# reemplazamos los Na por valores 0
84
85
enaho01_500$salarioxhora[is.na(enaho01_500$salarioxhora)] = 0
86
87
88
#Pregunta 3
89
enaho02 <- read_dta("../../../enaho/2020/737-Modulo02/enaho01-2020-200.dta")
90
91
92
base1 <- enaho02%>% group_by(conglome, vivienda, hogar ) %>% summarise(edad_max = max(p208a))
93
94
enaho34 <- read_dta("../../../enaho/2020/737-Modulo34/sumaria-2020.dta")
95
96
#Hacemos el merge
97
98
num = list(enaho34)
99
merge1 = enaho02
100
101
for (i in num){
102
103
merge1 <- merge(merge1, i,
104
by = c("conglome", "vivienda", "hogar"),
105
all.x = T, suffixes = c("",".y")
106
)
107
}
108
109
names(merge1)
110
111
#Creamos la variable dummy
112
pension <- merge1 %>%mutate(g1 = ifelse(edad_max <=65,1,0))
113