├── README.md ├── #34 Medición de la Pobreza Subjetiva.do ├── #6 #7 Construcción de indicadores sociales Cálculo de porcentajes personas aseguradas.do ├── #108 Tasa neta de asistencia edu. primaria.do ├── #41 Brecha Digital de Género.do ├── #107 Tasa neta de asistencia edu. secundaria.do ├── #106 Tasa neta de asistencia edu. inicial.do ├── #53 MINEDU-ESCALE Tasa de conclusión primaria, edades 12-13 (% del total).do ├── #39 Medición del Acceso a Programas Alimentarios.do ├── #100 Tasa de Empleo Formal e Informal de la PEA ocupada con la ENAHO en STATA 2017-2022.do ├── #62 Promedio de años de escolaridad, edades 25-34 (número de años).do ├── #96 Población en edad de trabajar según condición de actividad con la ENAHO en STATA 2017-2022.do ├── #48 Población que accede a internet en el Perú con la ENAHO en STATA.do ├── #86 Niñ@s y adolescentes que trabajan usando la ENAHO del 2019.do ├── #55 MINEDU-ESCALE Tasa de conclusión secundaria, edades 17-18 (% del total).do ├── #40 Medición de la Pobreza Relativa.do ├── #11 #12 Medición de la Pobreza Monetaria.do ├── #58 Medición de la Inclusión Financiera de la población ocupada de 18 y más años.do ├── #94 MINEDU-ESCALE Tasa Total De Matrícula Edades 3-5, 6-11 & 12-16 con la ENAHO en STATA.do ├── #63 Tenencia de Celular en el Hogar por Quintiles de Ingreso.do ├── #19 Cálculo del ingreso promedio mensual proveniente del trabajo.do ├── #29 Calculo de la población ocupada por rama de actividad y otros criterios.do ├── #30 Cálculo del ingreso promedio mensual proveniente del trabajo.do ├── #49 Medición de la tasa de analfabetismo en el Perú con la Encuesta Nacional de Hogares (ENAHO) 2019 en STATA.do ├── #14 Estimación de las Necesidades Básicas Insatisfechas.do ├── #8 #9 #10 Medición de Indicadores Sociales con la ENAHO en STATA Indicadores de Educación.do ├── #56 MINEDU-ESCALE Tasa de conclusión, ed. superior, edades 22-24 (% del total).do ├── #95 Población económicamente activa ocupada por rama de actividad con la ENAHO en STATA.do ├── #28 Construcción de Indicadores de Empleo y Desempleo.do ├── #109 #110 Gasto de bolsillo.do ├── #65 Evolución del gasto real per cápita promedio mensual con ENAHO 2016 – 2018.do ├── #52 Medición del Déficit Calórico con la ENAHO - Part1 - Norma Calorica.do ├── #21 Evolución del ingreso real per cápita promedio mensual ENAHO 2016 – 2019.do └── #52 Medición del Déficit Calórico con la ENAHO - Part2 - Calculo de calorias adquiridas.do /README.md: -------------------------------------------------------------------------------- 1 | # ENAHO-INDICADORES-STATA 2 | Sintaxis en Stata para estimar indicadores con la Encuesta Nacional de Hogares (ENAHO) del Perú. 3 | -------------------------------------------------------------------------------- /#34 Medición de la Pobreza Subjetiva.do: -------------------------------------------------------------------------------- 1 | ******************************************************************************** 2 | *POBREZA SUBJETIVA //Modulo 85 de GOBERNABILIDAD de la ENAHO 3 | /* 4 | Perú: Perfil de la pobreza por dominios geográficos, 2008-2018 5 | https://www.inei.gob.pe/media/MenuRecursivo/publicaciones_digitales/Est/Lib1699/ 6 | 5.1 Percepción del bienestar 7 | 5.2 Percepción de la situación económica de los hogares 8 | 5.3 Percepción acerca de los ingresos del hogar 9 | 5.5 Percepción acerca del nivel de vida de los Hogares 10 | */ 11 | 12 | cd "D:\ENAHO" 13 | use enaho01b-2018-2.dta, clear 14 | merge 1:1 conglome vivienda hogar using sumaria-2018.dta 15 | 16 | *Creo la variable pobreza categorizada como pobre/no pobre 17 | *sobre la base de la variable "pobreza" que ya se encuentra en el archivo 18 | gen pobre2=1 if pobreza<3 19 | replace pobre2=2 if pobreza==3 20 | *Etiquetamos los valores de la variable 21 | label define pobre2 1 "pobre" 2 "no_pobre" 22 | label value pobre2 pobre2 23 | 24 | *Percepción de bienestar 25 | tab p37 pobreza [aweight=factor07], col nofreq 26 | tab p37 pobre2 [aweight=factor07], col nofreq 27 | graph bar [aweight=factor07], over(p37) over(pobre2) /// 28 | asyvars percentages blabel(bar,format(%9.1f)) 29 | 30 | *Percepción de la situación económica 31 | tab p32 [aweight=factor07] 32 | tab p32 pobreza [aweight=factor07], col nofreq 33 | tab p32 pobre2 [aweight=factor07], col nofreq 34 | graph bar [aweight=factor07], over(p32) over(pobre2) /// 35 | asyvars percentages blabel(bar,format(%9.1f)) 36 | 37 | *Percepción acerca de los ingresos del hogar 38 | tab pobreza p38a [aweight=factor07], row nofreq 39 | tab pobre2 p38a [aweight=factor07], row nofreq 40 | graph bar [aweight=factor07], over(p38a) over(pobre2) /// 41 | asyvars percentages blabel(bar,format(%9.1f)) 42 | 43 | *Percepción acerca del nivel de vida de los Hogares 44 | tab p34 [aweight=factor07] 45 | tab p34 pobreza [aweight=factor07], col nofreq 46 | -------------------------------------------------------------------------------- /#6 #7 Construcción de indicadores sociales Cálculo de porcentajes personas aseguradas.do: -------------------------------------------------------------------------------- 1 | cd "D:\ENAHO 2017" 2 | 3 | *Abrir el archivo 4 | use enaho01a-2017-400.dta, clear 5 | 6 | *Crear variable geografica area 7 | gen area=estrato 8 | recode area (1/5=1) (6/8=2) 9 | lab def area 1 "Urbana" 2 "Rural", modify 10 | lab val area area 11 | lab var area "Area de residencia" 12 | 13 | *POBLACION CON SEGURO DE SALUD 14 | *Existen 8 variables que indican si una persona esta afiliado a un seguro 15 | sum p4191-p4198 16 | 17 | *La codificacion de las variables es 1 & 2, lo quiero cambiar a 1 & 0 18 | recode p4191-p4198 (2=0) 19 | gen seguro= p4191 + p4192 + p4193 + p4194 + p4195 + p4196 + p4197 + p4198 20 | gen asegurado=1 if seguro> 0 21 | replace asegurado=0 if seguro==0 22 | label define asegurado 0 "No" 1 "Si" 23 | label val asegurado asegurado 24 | label var asegurado "Tiene seguro de salud" 25 | 26 | *Usamos el factor de expansion de la base para obtener los resultados a nivel nacional 27 | tab asegurado [iweight=factor07] 28 | 29 | 30 | *PASANDO LOS RESULTADOS A UN DOC DE WORD 31 | 32 | *Establecemos las caracteristicas de la encuesta usando las variable 33 | *factor de expansion del modulo 300 (factora07), conglomerado y estrato 34 | svyset [pweight=factora07], psu(conglome) strata(estrato) 35 | 36 | *Empezamos a registrar los resultados en word 37 | *(recordar precisar el directorio en el que estan trabajando) 38 | asdoc, text(\b INDICADORES DE SALUD) fs(12) replace 39 | 40 | asdoc, text(\b Acceso a Seguro de Salud) fs(12) append 41 | 42 | asdoc tab asegurado [iweight=factor07] , /// 43 | title(\b Tiene seguro de salud) append fs(10) 44 | 45 | asdoc tab asegurado area [iweight=factor07] , /// 46 | title(\b Tiene seguro de salud por ambito geografico) append fs(10) 47 | 48 | asdoc tab asegurado area [iweight=factor07] , col nofreq /// 49 | title(\b Tiene seguro de salud por ambito geografico) append fs(10) 50 | 51 | asdoc tab asegurado p207 [iweight=factor07] , col nofreq /// 52 | title(\b Tiene seguro de salud por sexo) append fs(10) 53 | -------------------------------------------------------------------------------- /#108 Tasa neta de asistencia edu. primaria.do: -------------------------------------------------------------------------------- 1 | 2 | ***TASA NETA DE ASISTENCIA ESCOLAR 3 | use "D:\ENAHO 2021\enaho01a-2021-300.dta", clear 4 | 5 | /* P203 ¿Cuál es la relación de parentesco con el jefe(a) del hogar? 6 | 0. Panel 7 | 1. Jefe/Jefa 8 | 2. Esposo/Esposa 9 | 3. Hijo/Hija 10 | 4. Yerno/Nuera 11 | 5. Nieto 12 | 6. Padres/Suegros 13 | 7. Otros parientes 14 | 8. Trabajador Hogar 15 | 9. Pensionista 16 | 10. Otros no parientes 17 | 18 | P204 ¿Es miembro del hogar? 1. Si 2. No 19 | 20 | P208A ¿Qúe edad tiene en años cumplidos? (En años) 21 | 22 | P307 Actualmente, ¿Asiste a algún centro o programa de educación básica 23 | o superior bajo modalidad de educación a distancia? 1. Si 2. No 24 | 25 | P308A ¿Cuál es el grado o año de estudios en el que está matriculado? 26 | 1. Educacion inicial 27 | 2. Primaria 28 | 3. Secundaria 29 | 4. Superior no universitaria 30 | 5. Superior universitaria 31 | 6. Maestria/Doctorado 32 | 7. Básica especial */ 33 | 34 | 35 | *Var. Geograficas 36 | gen area=1 if estrato<=5 37 | replace area=2 if estrato>=6 & estrato<=8 38 | lab def area 1 "urbano" 2 "rural" 39 | lab val area area 40 | 41 | gen dpto= real(substr(ubigeo,1,2)) 42 | gen pro= substr(ubigeo,1,4) 43 | replace dpto=26 if (pro=="1501") 44 | label define dpto 1"Amazonas" 2"Ancash" 3"Apurimac" 4"Arequipa" 5"Ayacucho" /// 45 | 6"Cajamarca" 7"Callao" 8"Cusco" 9"Huancavelica" 10"Huanuco" 11"Ica" 12"Junin" 13"La_Libertad" /// 46 | 14"Lambayeque" 15"Lima 2/" 16"Loreto" 17"Madre_de_Dios" 18"Moquegua" 19"Pasco" /// 47 | 20"Piura" 21"Puno" 22"San_Martin" 23"Tacna" 24"Tumbes" 25"Ucayali" 26"MML 1/" 48 | lab val dpto dpto 49 | 50 | 51 | label define tna 0 "No" 1 "T_asistencia" 52 | 53 | *6 A 11 AÑOS 54 | gen tna_pri=. 55 | replace tna_pri=0 if codinfor!="00" & p204==1 & p203!= 8 & p203!= 9 & mes>= "04" & (p208a>=6 & p208a<=11) 56 | replace tna_pri=1 if codinfor!="00" & p204==1 & p203!= 8 & p203!= 9 & mes>= "04" & (p208a>=6 & p208a<=11)& (p307==1 & inlist(p308a,2,7)) 57 | label value tna_pri tna 58 | 59 | 60 | svyset [pweight = factora07], psu(conglome)strata(estrato) 61 | svy: mean tna_pri 62 | svy: mean tna_pri , over(area) 63 | svy: mean tna_pri , over(dpto) 64 | 65 | 66 | -------------------------------------------------------------------------------- /#41 Brecha Digital de Género.do: -------------------------------------------------------------------------------- 1 | cd "D:\ENAHO" 2 | use enaho01a-2018-300.dta, clear 3 | merge m:1 conglome vivienda hogar using sumaria-2018.dta 4 | 5 | *AREA DE RESIDENCIA 6 | gen area=1 if estrato<=5 7 | replace area=2 if estrato>=6 & estrato<=8 8 | lab def area 1 "urbano" 2 "rural" 9 | lab val area area 10 | 11 | *GRUPO DE EDAD 12 | gen gedad=1 if p208a>5 & p208a<12 13 | replace gedad=2 if p208a>11 & p208a<19 14 | replace gedad=3 if p208a>18 & p208a<25 15 | replace gedad=4 if p208a>24 & p208a<41 16 | replace gedad=5 if p208a>40 & p208a<60 17 | replace gedad=6 if p208a>59 18 | replace gedad=. if p208a==. 19 | replace gedad=. if p314a==. 20 | 21 | label define gedad 1 "6 a 11" 2 "12 a 18" 3 "19 a 24" /// 22 | 4 "25 a 40" 5 "41 a 59" 6 "60 a mas", replace 23 | label values gedad gedad 24 | 25 | gen mujer=p207==2 26 | 27 | gen gedu=1 if p301a<5 28 | replace gedu=2 if p301a>4 & p301a<7 29 | replace gedu=3 if p301a>6 & p301a<9 30 | replace gedu=4 if p301a==9 | p301a==10 | p301a==11 31 | label define gedad 1 "Primaria" 2 "Secundaria" /// 32 | 3 "Sup. No universitaria" 4 "Sup. universitaria", replace 33 | label values gedu gedu 34 | 35 | *Nacional 36 | tab p314a mujer [iw=factora07] , nofreq col 37 | 38 | *Pobreza 39 | tab p314a mujer [iw=factora07] if pobreza<3, nofreq col 40 | 41 | tab p314a mujer [iw=factora07] if pobreza==3, nofreq col 42 | 43 | *Area de residencia 44 | tab p314a mujer [iw=factora07] if area==1, nofreq col 45 | 46 | tab p314a mujer [iw=factora07] if area==2, nofreq col 47 | 48 | *Grupo de Edad 49 | tab p314a mujer [iw=factora07] if gedad==1, nofreq col 50 | 51 | tab p314a mujer [iw=factora07] if gedad==2, nofreq col 52 | 53 | tab p314a mujer [iw=factora07] if gedad==3, nofreq col 54 | 55 | tab p314a mujer [iw=factora07] if gedad==4, nofreq col 56 | 57 | tab p314a mujer [iw=factora07] if gedad==5, nofreq col 58 | 59 | tab p314a mujer [iw=factora07] if gedad==6, nofreq col 60 | 61 | *Educacion 62 | tab p314a mujer [iw=factora07] if gedu==1, nofreq col 63 | 64 | tab p314a mujer [iw=factora07] if gedu==2, nofreq col 65 | 66 | tab p314a mujer [iw=factora07] if gedu==3, nofreq col 67 | 68 | tab p314a mujer [iw=factora07] if gedu==4, nofreq col 69 | -------------------------------------------------------------------------------- /#107 Tasa neta de asistencia edu. secundaria.do: -------------------------------------------------------------------------------- 1 | 2 | ***TASA NETA DE ASISTENCIA ESCOLAR 3 | use "D:\ENAHO 2021\enaho01a-2021-300.dta", clear 4 | 5 | /* P203 ¿Cuál es la relación de parentesco con el jefe(a) del hogar? 6 | 0. Panel 7 | 1. Jefe/Jefa 8 | 2. Esposo/Esposa 9 | 3. Hijo/Hija 10 | 4. Yerno/Nuera 11 | 5. Nieto 12 | 6. Padres/Suegros 13 | 7. Otros parientes 14 | 8. Trabajador Hogar 15 | 9. Pensionista 16 | 10. Otros no parientes 17 | 18 | P204 ¿Es miembro del hogar? 1. Si 2. No 19 | 20 | P207 Sexo: 1. Hombre 2. Mujer 21 | 22 | P208A ¿Qúe edad tiene en años cumplidos? (En años) 23 | 24 | P307 Actualmente, ¿Asiste a algún centro o programa de educación básica 25 | o superior bajo modalidad de educación a distancia? 1. Si 2. No 26 | 27 | P308A ¿Cuál es el grado o año de estudios en el que está matriculado? 28 | 1. Educacion inicial 29 | 2. Primaria 30 | 3. Secundaria 31 | 4. Superior no universitaria 32 | 5. Superior universitaria 33 | 6. Maestria/Doctorado 34 | 7. Básica especial */ 35 | 36 | 37 | *Var. Geograficas 38 | gen area=1 if estrato<=5 39 | replace area=2 if estrato>=6 & estrato<=8 40 | lab def area 1 "urbano" 2 "rural" 41 | lab val area area 42 | 43 | gen dpto= real(substr(ubigeo,1,2)) 44 | gen pro= substr(ubigeo,1,4) 45 | replace dpto=26 if (pro=="1501") 46 | label define dpto 1"Amazonas" 2"Ancash" 3"Apurimac" 4"Arequipa" 5"Ayacucho" /// 47 | 6"Cajamarca" 7"Callao" 8"Cusco" 9"Huancavelica" 10"Huanuco" 11"Ica" 12"Junin" 13"La_Libertad" /// 48 | 14"Lambayeque" 15"Lima 2/" 16"Loreto" 17"Madre_de_Dios" 18"Moquegua" 19"Pasco" /// 49 | 20"Piura" 21"Puno" 22"San_Martin" 23"Tacna" 24"Tumbes" 25"Ucayali" 26"MML 1/" 50 | lab val dpto dpto 51 | 52 | 53 | label define tna 0 "No" 1 "T_asistencia" 54 | 55 | *12 A 16 AÑOS 56 | gen tna_sec=. 57 | replace tna_sec=0 if codinfor!="00" & p204==1 & p203!= 8 & p203!= 9 & mes>= "04" & (p208a>=12 & p208a<=16) 58 | replace tna_sec=1 if codinfor!="00" & p204==1 & p203!= 8 & p203!= 9 & mes>= "04" & (p208a>=12 & p208a<=16)& (p307==1 & p308a==3) 59 | label value tna_sec tna 60 | 61 | 62 | svyset [pweight = factora07], psu(conglome)strata(estrato) 63 | svy: mean tna_sec 64 | svy: mean tna_sec , over(area) 65 | svy: mean tna_sec , over(p207) 66 | svy: mean tna_sec , over(p207 area) 67 | 68 | -------------------------------------------------------------------------------- /#106 Tasa neta de asistencia edu. inicial.do: -------------------------------------------------------------------------------- 1 | 2 | ***TASA NETA DE ASISTENCIA ESCOLAR 3 | use "D:\ENAHO 2021\enaho01a-2021-300.dta", clear 4 | 5 | /* P203 ¿Cuál es la relación de parentesco con el jefe(a) del hogar? 6 | 0. Panel 7 | 1. Jefe/Jefa 8 | 2. Esposo/Esposa 9 | 3. Hijo/Hija 10 | 4. Yerno/Nuera 11 | 5. Nieto 12 | 6. Padres/Suegros 13 | 7. Otros parientes 14 | 8. Trabajador Hogar 15 | 9. Pensionista 16 | 10. Otros no parientes 17 | 18 | P204 ¿Es miembro del hogar? 1. Si 2. No 19 | 20 | P208A ¿Qúe edad tiene en años cumplidos? (En años) 21 | 22 | P307 Actualmente, ¿Asiste a algún centro o programa de educación básica 23 | o superior bajo modalidad de educación a distancia? 1. Si 2. No 24 | 25 | P308A ¿Cuál es el grado o año de estudios en el que está matriculado? 26 | 1. Educacion inicial 27 | 2. Primaria 28 | 3. Secundaria 29 | 4. Superior no universitaria 30 | 5. Superior universitaria 31 | 6. Maestria/Doctorado 32 | 7. Básica especial */ 33 | 34 | 35 | *Var. Geograficas 36 | gen area=1 if estrato<=5 37 | replace area=2 if estrato>=6 & estrato<=8 38 | lab def area 1 "urbano" 2 "rural" 39 | lab val area area 40 | 41 | gen regnat=dominio 42 | recode regnat (1/3 8 =1) (4/6=2) (7=3) 43 | label var regnat "region natural" 44 | label define regnat 1 "costa" 2 "sierra" 3 "selva" 45 | label values regnat regnat 46 | 47 | gen dpto= real(substr(ubigeo,1,2)) 48 | gen pro= substr(ubigeo,1,4) 49 | replace dpto=26 if (pro=="1501") 50 | label define dpto 1"Amazonas" 2"Ancash" 3"Apurimac" 4"Arequipa" 5"Ayacucho" /// 51 | 6"Cajamarca" 7"Callao" 8"Cusco" 9"Huancavelica" 10"Huanuco" 11"Ica" 12"Junin" 13"La_Libertad" /// 52 | 14"Lambayeque" 15"Lima 2/" 16"Loreto" 17"Madre_de_Dios" 18"Moquegua" 19"Pasco" /// 53 | 20"Piura" 21"Puno" 22"San_Martin" 23"Tacna" 24"Tumbes" 25"Ucayali" 26"MML 1/" 54 | lab val dpto dpto 55 | 56 | 57 | label define tna 0 "No" 1 "T_asistencia" 58 | 59 | *DE 3 A 5 AÑOS 60 | gen tna_ini=. 61 | replace tna_ini=0 if codinfor!="00" & p204==1 & p203!= 8 & p203!= 9 & mes>= "04" & (p208a>=3 & p208a<=5) 62 | replace tna_ini=1 if codinfor!="00" & p204==1 & p203!= 8 & p203!= 9 & mes>= "04" & (p208a>=3 & p208a<=5)& (p307==1 & p308a==1) 63 | label value tna_ini tna 64 | 65 | 66 | svyset [pweight = factora07], psu(conglome)strata(estrato) 67 | svy: mean tna_ini 68 | svy: mean tna_ini , over(area) 69 | svy: mean tna_ini , over(regnat) 70 | 71 | 72 | -------------------------------------------------------------------------------- /#53 MINEDU-ESCALE Tasa de conclusión primaria, edades 12-13 (% del total).do: -------------------------------------------------------------------------------- 1 | cd "D:\CONCLUSION" 2 | 3 | *Unir las bases del módulo salud, educación y sumaria 4 | use "enaho01a-2019-400.dta",clear 5 | merge 1:1 conglome vivienda hogar codperso using "enaho01a-2019-300.dta" 6 | keep if _merge==3 7 | drop _m 8 | merge m:1 conglome vivienda hogar using "sumaria-2019.dta", nogenerate 9 | 10 | *Estimar edad en años cumplidos al 31 de Marzo 11 | rename a*o year 12 | destring year, replace 13 | gen edad_31Marzo=year-p400a3 if p400a2<4 14 | replace edad_31Marzo=(year-p400a3-1) if p400a2>3 & p400a2<=12 15 | replace edad_31Marzo=0 if edad_31Marzo<0 16 | 17 | *Variables geográficas (departamento, urbano/rural) 18 | *Departamento (distinguir Lima Metropolitana de Lima Provincias) 19 | destring ubigeo, generate(dpto) 20 | replace dpto=dpto/10000 21 | replace dpto=round(dpto) 22 | replace dpto=26 if (dominio==8) //Lima Metropolitana 23 | replace dpto=27 if (dpto==15 & dominio!=8) //Lima Provincias 24 | 25 | label define dpto_num 1 "Amazonas" 26 | label define dpto_num 2 "Ancash", add 27 | label define dpto_num 3 "Apurimac", add 28 | label define dpto_num 4 "Arequipa", add 29 | label define dpto_num 5 "Ayacucho", add 30 | label define dpto_num 6 "Cajamarca", add 31 | label define dpto_num 7 "Callao", add 32 | label define dpto_num 8 "Cusco", add 33 | label define dpto_num 9 "Huancavelica", add 34 | label define dpto_num 10 "Huanuco", add 35 | label define dpto_num 11 "Ica", add 36 | label define dpto_num 12 "Junin", add 37 | label define dpto_num 13 "La_Libertad", add 38 | label define dpto_num 14 "Lambayeque", add 39 | *label define dpto_num 15 "Lima", add 40 | label define dpto_num 16 "Loreto", add 41 | label define dpto_num 17 "Madre_de_Dios",add 42 | label define dpto_num 18 "Moquegua", add 43 | label define dpto_num 19 "Pasco", add 44 | label define dpto_num 20 "Piura", add 45 | label define dpto_num 21 "Puno", add 46 | label define dpto_num 22 "San_Martin", add 47 | label define dpto_num 23 "Tacna", add 48 | label define dpto_num 24 "Tumbes", add 49 | label define dpto_num 25 "Ucayali", add 50 | label define dpto_num 26 "Lima_Metropolitana", add 51 | label define dpto_num 27 "Lima_Provincias", add 52 | label values dpto dpto_num 53 | 54 | *Area 55 | recode estrato (1/5=1 "Urbana")(6/8=2 "Rural"), gen(area) 56 | lab var area "Area de Residencia" 57 | 58 | *Etiquetar la variable pobreza monetaria 59 | *Usar variable "pobreza" de sumaria 60 | label define pobre_3 1 "pobre_extremo" 2 "pobre_no_extremo" 3 "no_pobre" 61 | label value pobreza pobre_3 62 | 63 | *Crear la variable lengua materna 64 | *Lengua indígena: quechua, aymara, otra lengua nativa 65 | *Lengua Castellano 66 | recode p300a (4=1 "Castellano")(1/3=2 "Indigena")(5/9=.), gen(lengua) 67 | -------------------------------------------------------------------------------- /#39 Medición del Acceso a Programas Alimentarios.do: -------------------------------------------------------------------------------- 1 | *ACCESO A PROGRAMA ALIMENTARIO 2 | *Establecer carpeta de trabajo 3 | cd "D:\ENAHO" 4 | 5 | *Bajar y unir los archivos del modulo 02, 34 & 37 6 | ******************************************************************************** 7 | 8 | use enaho01-2019-700.dta, clear 9 | gen area=estrato 10 | recode area (1/5=1) (6/8=2) 11 | label define area 1 "Urbana" 2 "Rural" 12 | 13 | label values area area 14 | 15 | *Hogares con al menos un miembro beneficiario con algún programa alimentario, 16 | *por área de residencia 17 | tab p701_09 area [iw=factor07], nofreq col 18 | 19 | 20 | ******************************************************************************** 21 | 22 | use enaho01-2019-700a.dta, clear 23 | rename p702 codperso 24 | 25 | egen c_p703=count(p703), by(conglome vivienda hogar codperso) 26 | 27 | collapse (max) c_p703 , by(conglome vivienda hogar codperso) 28 | 29 | merge 1:1 conglome vivienda hogar codperso using enaho01-2019-200.dta, nogenerate 30 | 31 | *Establecer a los residentes habituales 32 | gen resi=1 if ((p204==1 & p205==2) | (p204==2 & p206==1)) 33 | 34 | *Grupo de edad 35 | gen gedad=1 if p208a<18 36 | replace gedad=2 if p208a>17 & p208a<30 37 | replace gedad=3 if p208a>29 & p208a<60 38 | replace gedad=4 if p208a>59 39 | replace gedad=. if p208a==. 40 | 41 | label define gedad 1 "De 0 a 17 años de edad" 2 "De 18 a 29 añosde edad" /// 42 | 3 "De 30 a 59 años de edad" 4 "De 60 y más años de edad" 43 | label values gedad gedad 44 | 45 | gen b_pa=c_p703~=. 46 | label var b_pa "Beneficiarios Pro. Alimentario" 47 | 48 | *Población beneficiaria con algún programa alimentario, por grupos de edad 49 | tab b_pa gedad [iw=facpob07] if resi==1, nofreq col 50 | 51 | 52 | ******************************************************************************** 53 | use enaho01-2019-700a.dta, clear 54 | rename p702 codperso 55 | 56 | egen c_p703=count(p703), by(conglome vivienda hogar codperso) 57 | 58 | gen vd_l=p703==1 59 | gen co_p=p703==2 60 | gen de_e=p703==3 61 | gen al_e=p703==4 62 | gen wa_w=p703==5 63 | gen ot_p=p703>5 64 | 65 | collapse (max) c_p703 (sum) vd_l co_p de_e al_e wa_w ot_p /// 66 | , by(conglome vivienda hogar codperso) 67 | 68 | merge 1:1 conglome vivienda hogar codperso using enaho01-2019-200.dta, nogenerate 69 | 70 | *Distribución de la población beneficiaria de programas alimentarios 71 | *por tipo de programa que recibió 72 | gen prog=5 if c_p703~=. 73 | replace prog=1 if vd_l==1 & c_p703==1 74 | replace prog=2 if vd_l==1 & c_p703>1 75 | replace prog=3 if de_e==1 & c_p703==1 76 | replace prog=4 if al_e==1 & c_p703==1 77 | label define prog 1 "Solo Vaso de Leche" 2 "Vaso de Leche y Otros" /// 78 | 3 "Solo Desayuno Escolar" 4 "Solo Almuerzo Escolar" 5 "Otros" 79 | label value prog prog 80 | 81 | tab prog [iw=facpob07] 82 | 83 | -------------------------------------------------------------------------------- /#100 Tasa de Empleo Formal e Informal de la PEA ocupada con la ENAHO en STATA 2017-2022.do: -------------------------------------------------------------------------------- 1 | clear all 2 | cd "D:\ENAHO" 3 | /* 4 | use "enaho01a-2022-500.dta", clear 5 | append using "enaho01a-2021-500.dta" 6 | append using "enaho01a-2020-500.dta" 7 | append using "enaho01a-2019-500.dta" 8 | append using "enaho01a-2018-500.dta" 9 | append using "enaho01a-2017-500.dta" 10 | save "enaho500-2017-2022.dta", replace 11 | */ 12 | 13 | use "enaho500-2017-2022.dta", clear 14 | rename a*o year 15 | gen anio=real(year) 16 | label var year "Year" 17 | 18 | tab p500i year 19 | drop if p500i=="00" 20 | 21 | *Se establece quienes son residentes habituales 22 | gen resi=1 if ((p204==1 & p205==2) | (p204==2 & p206==1)) 23 | *P204 ¿Es miembro del hogar? 204. ¿Es miembro del hogar familiar? 1. Sí 2. No 24 | *P205 ¿Se encuentra ausente del hogar 30 días o más? 1. Sí 2. No 25 | *P206 ¿Está presente en el hogar 30 días o más? 1. Sí 2. No 26 | 27 | 28 | ******************************************************************************** 29 | *GRÁFICO N° 3.2 PERÚ: TASA DE EMPLEO FORMAL E INFORMAL Y VARIACIÓN DE LA PEA OCUPADA CON EMPLEO FORMAL E INFORMAL, 2017-2022 30 | 31 | *ocupinf= situacion de informalidad, 1: empleo informal, 2: empleo formal 32 | table ocupinf year if resi==1 [iw= fac500a], nformat(%12.0fc) 33 | 34 | gen informal=ocupinf 35 | replace informal=0 if ocupinf==2 36 | gen formal =1 if ocupinf==2 37 | replace formal =0 if ocupinf==1 38 | 39 | graph bar informal formal if resi==1 [pw= fac500a], /// 40 | over(year) per stack yscale(off) /// 41 | legend(cols(2) position(6) label(1 "Tasa de empleo informal") label(2 "Tasa de empleo formal")) /// 42 | blabel(bar, format(%4.1f) size(small) position(center) color(white)) /// 43 | bar(1,color(red)) bar(2,color(orange)) /// 44 | saving(gr3_2a.gph, replace) 45 | 46 | bys year: egen double t_informal = total(fac500a * informal) if resi==1 47 | bys year: egen float t_formal = total(fac500a * formal) if resi==1 48 | label var t_informal "Empleos Informales" 49 | label var t_formal "Empleos Formales" 50 | 51 | br t_informal t_formal 52 | format t_* %20.0fc 53 | 54 | preserve 55 | 56 | keep t_* anio 57 | duplicates drop 58 | 59 | scatter t_informal t_formal anio, /// 60 | connect(l l) lcolor(red orange) msymbol(D S) mcolor(red orange) /// 61 | mlabel(t_informal t_formal) mlabcolor(red orange) mlabposition(6 12) /// 62 | mlabsize(small small) yscale(off) xtitle("") legend(cols(2) position(6)) /// 63 | xlabel(, nogrid) ylabel(, nogrid) /// 64 | saving(gr3_2b.gph, replace) 65 | 66 | restore 67 | 68 | 69 | gr combine gr3_2a.gph gr3_2b.gph, /// 70 | title("GRÁFICO N° 3.2" "PERÚ: TASA DE EMPLEO FORMAL E INFORMAL Y VAR. DE LA PEA OCUPADA" "CON EMPLEO FORMAL E INFORMAL, 2017-2022", size(small) color(red) margin(medsmall)) /// 71 | subtitle("(Porcentaje y absolutos)", size(small)) /// 72 | xcommon commonscheme 73 | graph export gr3_2.png, replace 74 | 75 | 76 | -------------------------------------------------------------------------------- /#62 Promedio de años de escolaridad, edades 25-34 (número de años).do: -------------------------------------------------------------------------------- 1 | cd "E:\ESCOLARIDAD" 2 | 3 | use enaho01a-2019-400.dta, clear 4 | merge 1:1 conglome vivienda hogar codperso using enaho01a-2019-300.dta 5 | keep if _merge==3 6 | drop _m 7 | merge m:1 conglome vivienda hogar using sumaria-2019.dta, nogenerate 8 | 9 | *Lengua materna 10 | gen lengua=1 if p300a==4 11 | replace lengua=2 if p300a<4 12 | replace lengua=3 if p300a>5 13 | lab def lengua 1 "Castellano" 2 "Indígena" 3 "Otros", modify 14 | lab val lengua lengua 15 | 16 | *Crear variable geografica area 17 | gen area=estrato 18 | recode area (1/5=1) (6/8=2) 19 | lab def area 1 "Urbana" 2 "Rural", modify 20 | lab val area area 21 | lab var area "Area de residencia" 22 | 23 | *Tomando de la sintaxis de la tesis 24 | *POBREZA Y LOGRO EDUCATIVO EN LA REGIÓN PIURA 2015 25 | *Presentada por: Br. Econ. Charlye Casariego Acenjo 26 | *https://repositorio.unp.edu.pe/bitstream/handle/UNP/1498/ECO-CAS-ASE-2018.pdf?sequence=1&isAllowed=y 27 | generate byte X5=p301b 28 | replace X5=0 if p301a==1 | p301a==2 29 | recode X5 (1=1) (2=2) (3=3) (4=4) if p301a==3 30 | recode X5 (5=5) (6=6) if p301a==4 31 | recode X5 (1=7) (2=8) (3=9) (4=10) if p301a==5 32 | recode X5 (5=11)(6=12) if p301a==6 33 | recode X5 (1=12)(2=13)(3=14)(4=15) if p301a==7 34 | recode X5 (3=14)(4=15)(5=16) if p301a==8 35 | recode X5 (1=12)(2=13)(3=14)(4=15)(5=16) (6=17) if p301a==9 36 | recode X5 (4=15)(5=16)(6=17)(7=18) if p301a==10 37 | recode X5 (1=17)(2=18) if p301a==11 38 | g _p301c=p301c 39 | recode _p301c (0=1) 40 | replace X5=_p301c if p301b==0 & p301a!=2 41 | label value X5 X5 42 | label variable X5 "Años de estudio del individuo" 43 | 44 | *Estimar edad en años cumplidos al 31 de Marzo 45 | rename a*o year 46 | destring year, replace 47 | gen edad_31Marzo=year-p400a3 48 | replace edad_31Marzo=(year-p400a3-1) if p400a2>3 49 | replace edad_31Marzo=0 if edad_31Marzo<0 50 | 51 | gen escolt2534=X5 if (edad_31Marzo >= 25 & edad_31Marzo <= 34) 52 | 53 | sum escolt2534 [iweight=factora07] 54 | sum escolt2534 [iweight=factora07] if p207==2 55 | sum escolt2534 [iweight=factora07] if p207==1 56 | sum escolt2534 [iweight=factora07] if area==1 57 | sum escolt2534 [iweight=factora07] if area==1 & p207==2 58 | sum escolt2534 [iweight=factora07] if area==1 & p207==1 59 | sum escolt2534 [iweight=factora07] if area==2 60 | sum escolt2534 [iweight=factora07] if area==2 & p207==2 61 | sum escolt2534 [iweight=factora07] if area==2 & p207==1 62 | 63 | tab lengua 64 | sum escolt2534 [iweight=factora07] if lengua==1 65 | sum escolt2534 [iweight=factora07] if lengua==2 66 | 67 | tab pobreza 68 | sum escolt2534 [iweight=factora07] if pobreza==3 69 | sum escolt2534 [iweight=factora07] if pobreza==2 70 | sum escolt2534 [iweight=factora07] if pobreza==1 71 | -------------------------------------------------------------------------------- /#96 Población en edad de trabajar según condición de actividad con la ENAHO en STATA 2017-2022.do: -------------------------------------------------------------------------------- 1 | clear all 2 | cd "D:\ENAHO" 3 | 4 | /* 5 | use "enaho01a-2022-500.dta", clear 6 | append using "enaho01a-2021-500.dta" 7 | append using "enaho01a-2020-500.dta" 8 | append using "enaho01a-2019-500.dta" 9 | append using "enaho01a-2018-500.dta" 10 | append using "enaho01a-2017-500.dta" 11 | save "enaho500-2017-2022.dta", replace 12 | */ 13 | 14 | use "enaho500-2017-2022.dta", clear 15 | rename a*o year 16 | label var year "Year" 17 | 18 | tab p500i year 19 | drop if p500i=="00" 20 | 21 | *Se establece quienes son residentes habituales 22 | gen resi=1 if ((p204==1 & p205==2) | (p204==2 & p206==1)) 23 | *P204 ¿Es miembro del hogar? 204. ¿Es miembro del hogar familiar? 1. Sí 2. No 24 | *P205 ¿Se encuentra ausente del hogar 30 días o más? 1. Sí 2. No 25 | *P206 ¿Está presente en el hogar 30 días o más? 1. Sí 2. No 26 | 27 | 28 | 29 | ******************************************************************************** 30 | *CUADRO N° 2.1 PERÚ: DISTRIBUCIÓN DE LA POBLACIÓN EN EDAD DE TRABAJAR, 31 | *SEGÚN CONDICIÓN DE ACTIVIDAD, 2017-2022 32 | 33 | collect clear 34 | 35 | tabulate ocu500 year if resi==1 [iw= fac500a] 36 | table ocu500 year if resi==1 [iw= fac500a], nformat(%12.0fc) 37 | 38 | 39 | *collect style header edits the content of the table headers 40 | collect style header, title(hide) 41 | collect style header var, level(value) 42 | collect style header result, level(hide) 43 | collect preview 44 | 45 | putexcel set cuadros.xlsx, sheet(2_1, replace) replace 46 | putexcel A1= "CUADRO N° 2.1 PERÚ: DISTRIBUCIÓN DE LA POBLACIÓN EN EDAD DE TRABAJAR, SEGÚN CONDICIÓN DE ACTIVIDAD, 2017-2022", bold font("Calibri",11,"red") 47 | putexcel A2= "(Absoluto y porcentaje)", font("",10,"") 48 | putexcel A4 = collect 49 | putexcel A4 = "Condición de actividad" 50 | putexcel A5= "PEA Ocupada" 51 | putexcel A6= "PEA Desempleada" 52 | putexcel A7= "Desempleo oculto" 53 | putexcel A8= "Inactivo pleno" 54 | putexcel A9= "Pob. en Edad de Trabajar (PET)" 55 | 56 | putexcel A10= "Pob. Eco. Activa (PEA)" 57 | putexcel B10=formula(B5 + B6), nformat(#,###) 58 | putexcel C10=formula(C5 + C6), nformat(#,###) 59 | putexcel D10=formula(D5 + D6), nformat(#,###) 60 | putexcel E10=formula(E5 + E6), nformat(#,###) 61 | putexcel F10=formula(F5 + F6), nformat(#,###) 62 | putexcel G10=formula(G5 + G6), nformat(#,###) 63 | putexcel H10=formula(H5 + H6), nformat(#,###) 64 | 65 | putexcel A11= "Pob. Eco. Inactiva (PEI)" 66 | putexcel B11=formula(B7 + B8), nformat(#,###) 67 | putexcel C11=formula(C7 + C8), nformat(#,###) 68 | putexcel D11=formula(D7 + D8), nformat(#,###) 69 | putexcel E11=formula(E7 + E8), nformat(#,###) 70 | putexcel F11=formula(F7 + F8), nformat(#,###) 71 | putexcel G11=formula(G7 + G8), nformat(#,###) 72 | putexcel H11=formula(H7 + H8), nformat(#,###) 73 | 74 | forvalues i=5/11 { 75 | putexcel J`i'=formula(G`i'- D`i'), nformat(#,###) 76 | putexcel K`i'=formula((G`i'- D`i')*100/D`i' ), nformat(#.0) 77 | } 78 | putexcel J4= "Var. Abs. 2022/2019" 79 | putexcel K4= "Var. Por. 2022/2019" 80 | 81 | putexcel A12= "Nota: La suma de las partes puede no coincidir con el total debido al redondeo de las cifras." 82 | putexcel A13= "Fuente: INEI-Encuesta Nacional de Hogares, 2017-2022." 83 | putexcel A14= "Elaboración: MTPE-DGPE-Dirección de Investigación Socio Económico Laboral (DISEL)." 84 | 85 | putexcel A4:K4 A9:K11, bold font("",9,"") 86 | putexcel A5:K8 A12:A14, font("",9,"") 87 | putexcel A1:K14, border(all,none) 88 | -------------------------------------------------------------------------------- /#48 Población que accede a internet en el Perú con la ENAHO en STATA.do: -------------------------------------------------------------------------------- 1 | *********** 2 | *USO DE INTERNET 3 | *Población de 6 y más años de edad que hace uso de internet, según ámbito geográfico, 2019 4 | *Población de 6 y más años de edad que hace uso de internet, según nivel educativo, frecuencia de uso y ámbito geográfico, 2019 5 | 6 | *Establecer el directorio donde vamos a trabajar 7 | *Cambiar el nombre segun el folder donde van a bajar la informacion 8 | cd "D:\EDUCACION" 9 | 10 | *Bajar el modulo de educacion de la Enaho 2019 y descomprimir el folder 11 | *Guardar los pasos en un archivo log 12 | log using usointernet_enaho2019, replace 13 | 14 | *Abrir archivo 15 | use enaho01a-2019-300.dta, clear 16 | 17 | *Area de residencia 18 | gen area=estrato 19 | recode area (1/5=1) (6/8=2) 20 | lab def area 1 "Urbana" 2 "Rural", modify 21 | lab val area area 22 | lab var area "Area de residencia" 23 | 24 | *Region Natural 25 | gen regnat=1 if dominio<=3 | dominio==8 26 | replace regnat=2 if dominio>=4 & dominio<=6 27 | replace regnat=3 if dominio==7 28 | lab var regnat "Region natural" 29 | lab def regnat 1 "Costa" 2 "Sierra" 3 "Selva" 30 | lab val regnat regnat 31 | 32 | *Departamento (25) 33 | destring ubigeo, generate(dpto) 34 | replace dpto=dpto/10000 35 | replace dpto=round(dpto) 36 | label variable dpto "dpto" 37 | label define dpto 1 "Amazonas" 38 | label define dpto 2 "Ancash", add 39 | label define dpto 3 "Apurimac", add 40 | label define dpto 4 "Arequipa", add 41 | label define dpto 5 "Ayacucho", add 42 | label define dpto 6 "Cajamarca", add 43 | label define dpto 7 "Callao", add 44 | label define dpto 8 "Cusco", add 45 | label define dpto 9 "Huancavelica", add 46 | label define dpto 10 "Huanuco", add 47 | label define dpto 11 "Ica", add 48 | label define dpto 12 "Junin", add 49 | label define dpto 13 "La_Libertad", add 50 | label define dpto 14 "Lambayeque", add 51 | label define dpto 15 "Lima", add 52 | label define dpto 16 "Loreto", add 53 | label define dpto 17 "Madre_de_Dios", add 54 | label define dpto 18 "Moquegua", add 55 | label define dpto 19 "Pasco", add 56 | label define dpto 20 "Piura", add 57 | label define dpto 21 "Puno", add 58 | label define dpto 22 "San_Martin", add 59 | label define dpto 23 "Tacna", add 60 | label define dpto 24 "Tumbes", add 61 | label define dpto 25 "Ucayali", add 62 | label values dpto dpto 63 | 64 | *NIVEL EDUCATIVO 65 | gen educa=1 if p301a<5 | p301a==12 66 | replace educa=2 if p301a>4 & p301a<7 67 | replace educa=3 if p301a>6 & p301a<9 68 | replace educa=4 if p301a>8 & p301a<12 69 | label define educa 1 "Hasta primaria 1/" 2 "Secundaria" /// 70 | 3 "Superior no universitaria" 4 "Superior universitaria" 71 | label value educa educa 72 | 73 | *FRECUENCIA DE USO DE INTERNET 74 | recode p314d 1=1 2=2 3 4=3, gen(uso_internet) 75 | label define uso_internet 1 "Una vez al día" 2 "Una vez a la semana" 3 "Una vez al mes o más" 76 | label value uso_internet uso_internet 77 | 78 | *USO DE INTERNET 79 | tab p314a [iweight=factora07] if p208>5 & p204==1 80 | tab area p314a [iweight=factora07] if p208>5 & p204==1, nofreq row 81 | tab regnat p314a [iweight=factora07] if p208>5 & p204==1, nofreq row 82 | tab dpto p314a [iweight=factora07] if p208>5 & p204==1, nofreq row 83 | tab educa uso_internet [iweight=factora07] if p208>5 & p204==1, nofreq col 84 | tab educa uso_internet [iweight=factora07] if p208>5 & p204==1 & area==1, nofreq col 85 | tab educa uso_internet [iweight=factora07] if p208>5 & p204==1 & area==2, nofreq col 86 | 87 | log close 88 | -------------------------------------------------------------------------------- /#86 Niñ@s y adolescentes que trabajan usando la ENAHO del 2019.do: -------------------------------------------------------------------------------- 1 | 2 | cd "C:\ENAHO" 3 | 4 | use enaho01-2019-200, clear 5 | gen TI_5_17=0 if (p208a >= 5 & p208a <= 17) 6 | replace TI_5_17=1 if ((p208a >= 5 & p208a <= 17) & ((p210==1) | (p210==2 & (t211~=9 & t211~=11)))) 7 | lab def TI_5_17 0 "No trabaja" 1 "Si trabaja", modify 8 | lab val TI_5_17 TI_5_17 9 | 10 | gen hombre=p207==1 11 | lab def hombre 0 "Femenino" 1 "Masculino", modify 12 | lab val hombre hombre 13 | 14 | gen area=estrato 15 | recode area (1/5=1) (6/8=2) 16 | lab def area 1 "Urbana" 2 "Rural", modify 17 | lab val area area 18 | lab var area "Area de residencia" 19 | 20 | gen edad=1 if p208a>4 & p208a<14 21 | replace edad=2 if p208a>13 & p208a<18 22 | label define edad 1 "5-13" 2 "14-17" 23 | lab val edad edad 24 | 25 | *Departamento 26 | *1/ Comprende los 43 distritos que conforman la provincia de Lima. 27 | *2/ Comprende las provincias de Barranca, Cajatambo, Canta, Cañete, Huaral, Huarochirí, Huaura, Oyón y Yauyos. 28 | destring ubigeo, generate(dpto) 29 | replace dpto=dpto/10000 30 | replace dpto=round(dpto) 31 | replace dpto=151 if substr(ubigeo,1,4)=="1501" 32 | label variable dpto "Departamento" 33 | label define dpto 1 "Amazonas" 34 | label define dpto 2 "Ancash", add 35 | label define dpto 3 "Apurimac", add 36 | label define dpto 4 "Arequipa", add 37 | label define dpto 5 "Ayacucho", add 38 | label define dpto 6 "Cajamarca", add 39 | label define dpto 7 "Callao", add 40 | label define dpto 8 "Cusco", add 41 | label define dpto 9 "Huancavelica", add 42 | label define dpto 10 "Huanuco", add 43 | label define dpto 11 "Ica", add 44 | label define dpto 12 "Junin", add 45 | label define dpto 13 "La_Libertad", add 46 | label define dpto 14 "Lambayeque", add 47 | label define dpto 15 "Lima_region", add 48 | label define dpto 151 "Lima_MML", add 49 | label define dpto 16 "Loreto", add 50 | label define dpto 17 "Madre_de_Dios", add 51 | label define dpto 18 "Moquegua", add 52 | label define dpto 19 "Pasco", add 53 | label define dpto 20 "Piura", add 54 | label define dpto 21 "Puno", add 55 | label define dpto 22 "San_Martin", add 56 | label define dpto 23 "Tacna", add 57 | label define dpto 24 "Tumbes", add 58 | label define dpto 25 "Ucayali", add 59 | label values dpto dpto 60 | 61 | *Establecemos las caracteristicas de la encuesta usando las variable factor de expansion, conglomerado y estrato 62 | svyset [pweight=facpob07], psu(conglome) strata(estrato) 63 | svy: tab TI_5_17 , percent 64 | svy: tab hombre TI_5_17 , percent format(%10.3g) row 65 | svy: tab area TI_5_17 , percent format(%10.3g) row 66 | svy: tab edad TI_5_17 , percent format(%10.3g) row 67 | svy: tab dpto TI_5_17 , percent format(%10.3g) row 68 | 69 | *Graficar el porcentaje de niños, niñas y adolescentes que trabajan (% de población con edades 5-17) por departamento 70 | gen gra_TI_5_17=TI_5_17*100 71 | graph bar gra_TI_5_17 [pweight=facpob07], over(dpto, sort(1) descending label(labsize(*0.5) angle(45))) /// 72 | title("Porcentaje de niños, niñas y adolescentes que trabajan (% de población con edades 5-17)," "según Región, 2019", justification(left) size(small) color(black)) /// 73 | ytitle("Porcentaje",size(vsmall)) ylabel(,labsize(vsmall)) yline(25.97, lwidth(0.4)) /// 74 | blabel(total,size(vsmall) format(%5.3g)) /// 75 | note("Lima MML comprende los 43 distritos que conforman la provincia de Lima." "Lima region comprende las provincias de: Barranca, Cajatambo, Canta, Cañete, Huaral, Huarochirí, Huaura, Oyón y Yauyos.""Fuente: Encuesta Nacional de Hogares 2019",size(tiny)) /// 76 | saving(gra_TI_5_17, asis replace) 77 | 78 | *Cerrar archivo log 79 | log close 80 | -------------------------------------------------------------------------------- /#55 MINEDU-ESCALE Tasa de conclusión secundaria, edades 17-18 (% del total).do: -------------------------------------------------------------------------------- 1 | cd "D:\CONCLUSION" 2 | 3 | *Unir las bases del módulo salud, educación y sumaria 4 | use "enaho01a-2019-400.dta",clear 5 | merge 1:1 conglome vivienda hogar codperso using "enaho01a-2019-300.dta" 6 | keep if _merge==3 7 | drop _m 8 | merge m:1 conglome vivienda hogar using "sumaria-2019.dta", nogenerate 9 | 10 | *Estimar edad en años cumplidos al 31 de Marzo 11 | rename a*o year 12 | destring year, replace 13 | gen edad_31Marzo=year-p400a3 if p400a2<4 14 | replace edad_31Marzo=(year-p400a3-1) if p400a2>3 & p400a2<=12 15 | replace edad_31Marzo=0 if edad_31Marzo<0 16 | 17 | *Variables geográficas (departamento, urbano/rural) 18 | *Departamento (distinguir Lima Metropolitana de Lima Provincias) 19 | destring ubigeo, generate(dpto) 20 | replace dpto=dpto/10000 21 | replace dpto=round(dpto) 22 | replace dpto=26 if (dominio==8) //Lima Metropolitana 23 | replace dpto=27 if (dpto==15 & dominio!=8) //Lima Provincias 24 | label define dpto_num 1 "Amazonas" 25 | label define dpto_num 2 "Ancash", add 26 | label define dpto_num 3 "Apurimac", add 27 | label define dpto_num 4 "Arequipa", add 28 | label define dpto_num 5 "Ayacucho", add 29 | label define dpto_num 6 "Cajamarca", add 30 | label define dpto_num 7 "Callao", add 31 | label define dpto_num 8 "Cusco", add 32 | label define dpto_num 9 "Huancavelica", add 33 | label define dpto_num 10 "Huanuco", add 34 | label define dpto_num 11 "Ica", add 35 | label define dpto_num 12 "Junin", add 36 | label define dpto_num 13 "La_Libertad", add 37 | label define dpto_num 14 "Lambayeque", add 38 | *label define dpto_num 15 "Lima", add 39 | label define dpto_num 16 "Loreto", add 40 | label define dpto_num 17 "Madre_de_Dios",add 41 | label define dpto_num 18 "Moquegua", add 42 | label define dpto_num 19 "Pasco", add 43 | label define dpto_num 20 "Piura", add 44 | label define dpto_num 21 "Puno", add 45 | label define dpto_num 22 "San_Martin", add 46 | label define dpto_num 23 "Tacna", add 47 | label define dpto_num 24 "Tumbes", add 48 | label define dpto_num 25 "Ucayali", add 49 | label define dpto_num 26 "Lima_Metropolitana", add 50 | label define dpto_num 27 "Lima_Provincias", add 51 | label values dpto dpto_num 52 | 53 | *Area 54 | recode estrato (1/5=1 "Urbana")(6/8=2 "Rural"), gen(area) 55 | lab var area "Area de Residencia" 56 | 57 | *Etiquetar la variable pobreza monetaria 58 | *Usar variable "pobreza" de sumaria 59 | label define pobre_3 1 "pobre_extremo" 2 "pobre_no_extremo" 3 "no_pobre" 60 | label value pobreza pobre_3 61 | 62 | *Crear la variable lengua materna 63 | *Lengua indígena: quechua, aymara, otra lengua nativa 64 | *Lengua Castellano 65 | recode p300a (4=1 "Castellano")(1/3=2 "Indigena")(5/9=.), gen(lengua) 66 | 67 | 68 | ****************************************** 69 | *Establecer las caracteristicas de la encuesta usando las variables 70 | *factor de expansion del modulo 300 (factora07), conglomerado y estrato 71 | svyset [pweight=factora07], psu(conglome) strata(estrato) singleunit(centered) 72 | 73 | *Tasa de Conclusión Secundaria, 17 a 18 años (% del total) 74 | gen conclus1718=0 if edad_31Marzo>=17 & edad_31Marzo<=18 75 | replace conclus1718=1 if (edad_31Marzo>=17 & edad_31Marzo<=18) & (p301a>=6 & p301a<=11) 76 | label def con_l 0 "Incompleta" 1 "Completa" 77 | label val conclus1718 con_l 78 | 79 | svy: prop conclus1718 80 | svy: prop conclus1718 , over(p207) 81 | svy: prop conclus1718 , over(area) 82 | svy: prop conclus1718 , over(lengua) 83 | svy: prop conclus1718 , over(pobreza) 84 | svy: prop conclus1718 , over(dpto) 85 | -------------------------------------------------------------------------------- /#40 Medición de la Pobreza Relativa.do: -------------------------------------------------------------------------------- 1 | 2 | cd "D:\ENAHO" //Cambiar el nombre segun el folder donde van a bajar la informacion 3 | 4 | *Bajamos el modulo sumaria de la pagina del INEI 5 | copy "http://iinei.inei.gob.pe/iinei/srienaho/descarga/STATA/634-Modulo34.zip" 634-Modulo34.zip, replace 6 | *Descomprimir los archivos "Manualmente" 7 | 8 | use sumaria-2018.dta, clear 9 | 10 | *Generar la variable urbano/rural (usar la variable estrato) 11 | gen area=1 if estrato<=5 12 | replace area=2 if estrato>=6 & estrato<=8 13 | lab var area "Area de residencia" 14 | lab def area 1 "Urbano" 2 "Rural" 15 | lab val area area 16 | 17 | *Departamento 18 | destring ubigeo, generate(dpto) 19 | replace dpto=dpto/10000 20 | replace dpto=round(dpto) 21 | label variable dpto "dpto" 22 | label define dpto 1 "Amazonas" 23 | label define dpto 2 "Ancash", add 24 | label define dpto 3 "Apurimac", add 25 | label define dpto 4 "Arequipa", add 26 | label define dpto 5 "Ayacucho", add 27 | label define dpto 6 "Cajamarca", add 28 | label define dpto 7 "Callao", add 29 | label define dpto 8 "Cusco", add 30 | label define dpto 9 "Huancavelica", add 31 | label define dpto 10 "Huanuco", add 32 | label define dpto 11 "Ica", add 33 | label define dpto 12 "Junin", add 34 | label define dpto 13 "La_Libertad", add 35 | label define dpto 14 "Lambayeque", add 36 | label define dpto 15 "Lima", add 37 | label define dpto 16 "Loreto", add 38 | label define dpto 17 "Madre_de_Dios", add 39 | label define dpto 18 "Moquegua", add 40 | label define dpto 19 "Pasco", add 41 | label define dpto 20 "Piura", add 42 | label define dpto 21 "Puno", add 43 | label define dpto 22 "San_Martin", add 44 | label define dpto 23 "Tacna", add 45 | label define dpto 24 "Tumbes", add 46 | label define dpto 25 "Ucayali", add 47 | label values dpto dpto 48 | label var dpto "Departamento" 49 | 50 | *Estimamos el gasto mensual promedio de los hogares en terminos per capita 51 | gen gpcm=gashog2d/(mieperho*12) 52 | 53 | *Creamos la variable factor de expansion de la poblacion 54 | gen facpob=factor07*mieperho 55 | 56 | *Establecemos las caracteristicas de la encuesta usando las variables 57 | *factor de expansion, conglomerado y estrato 58 | svyset [pweight=facpob], psu(conglome) strata(estrato) 59 | 60 | 61 | *Para medir la pobreza relativa se debe establecer una linea de pobreza 62 | *igual para todos los hogares. Dicha linea es la media, 63 | *por ello estimamos el valor promedo del gasto mensual per capita de los hogares 64 | 65 | svy: mean gpcm //esto arroja un valor de 671.8312 para la ENAHO 2018 66 | 67 | sum gpcm [w=facpob] 68 | gen p_relativa=(gpcm=6 & estrato<=8 49 | lab var area "Area de residencia" 50 | lab def area 1 "Urbano" 2 "Rural" 51 | lab val area area 52 | 53 | *Generar la variable region regiones naturales (usar la variable dominio) 54 | gen regnat=1 if dominio<=3 | dominio==8 55 | replace regnat=2 if dominio>=4 & dominio<=6 56 | replace regnat=3 if dominio==7 57 | lab var regnat "Region natural" 58 | lab def regnat 1 "Costa" 2 "Sierra" 3 "Selva" 59 | lab val regnat regnat 60 | 61 | 62 | *Pobreza Monetaria Total 63 | *Estimamos el gasto mensual promedio de los hogares en terminos per capita 64 | gen gpcm=gashog2d/(mieperho*12) 65 | 66 | *Calculamos el factor de ponderacion a nivel de la poblacion 67 | gen facpob=factor*mieperho 68 | 69 | *Contrastamos gpcm con la linea de pobreza alimentaria (linpe) y pobreza total (linea) 70 | gen pobre3=1 if gpcm < linpe 71 | replace pobre3=2 if gpcm >= linpe & gpcm < linea 72 | replace pobre3=3 if gpcm >= linea 73 | 74 | *Etiquetamos los valores de la variable 75 | label define pobre3 1 "pobre_extremo" 2 "pobre_no_extremo" 3 "no_pobre" 76 | label value pobre3 pobre3 77 | label var pobre3 "Pobreza Monetaria" 78 | gen pobre_2=1 if gpcm < linea 79 | replace pobre_2=0 if gpcm >= linea 80 | 81 | *Etiquetamos los valores de la variable 82 | label define pobre_2 1 "pobre" 0 "no_pobre" 83 | label value pobre_2 pobre_2 84 | label var pobre_2 "Pobreza Monetaria Total" 85 | 86 | *Establecemos las caracteristicas de la encuesta usando las variable 87 | *factor de expansion, conglomerado y estrato 88 | svyset [pweight=facpob], psu(conglome) strata(estrato) 89 | 90 | 91 | *INCIDENCIA, BRECHA Y SEVERIDAD DE LA POBREZA MONETARIA 92 | *FGT(0): INCIDENCIA; FGT(1): BRECHA & FGT(2): SEVERIDAD 93 | *FGT(0) o (P0): Incidencia de la pobreza, que representa la proporción de pobres 94 | *o de pobres extremos como porcentaje del total de la población. 95 | *FGT(1) o (P1): Brecha de la pobreza, que mide la insuficiencia promedio del consumo 96 | *de los pobres respecto de la línea de pobreza, tomando en cuenta la proporción de la 97 | *población pobre en la población total. 98 | *FGT(2) o (P2): Severidad de la pobreza, que mide la desigualdad entre los pobres. 99 | 100 | *Se puede usar el comando povdeco, sepov o svy 101 | povdeco gpcm [w=facpob], varpl(linea) 102 | sepov gpcm [w=facpob], p(linea) 103 | 104 | *diferencias en el formato del output usando "prop" vs. "mean" 105 | svy: prop pobre_2 106 | svy: mean pobre_2 107 | svy: prop pobre_2, over(dpto) 108 | 109 | -------------------------------------------------------------------------------- /#58 Medición de la Inclusión Financiera de la población ocupada de 18 y más años.do: -------------------------------------------------------------------------------- 1 | *Establecer carpeta de trabajo 2 | cd "D:\ENAHO" 3 | 4 | *Bajar y unir los archivos del modulo de empleo y sumaria 5 | use enaho01a-2017-500.dta, clear 6 | merge m:1 conglome vivienda hogar using "sumaria-2017.dta" 7 | 8 | *Area urbana/rural 9 | gen area=1 if estrato<=5 10 | replace area=2 if estrato>=6 & estrato<=8 11 | lab def area 1 "urbano" 2 "rural" 12 | lab val area area 13 | 14 | *Region Natural 15 | gen regnat=1 if dominio<=3 | dominio==8 16 | replace regnat=2 if dominio>=4 & dominio<=6 17 | replace regnat=3 if dominio==7 18 | lab var regnat "Region natural" 19 | lab def regnat 1 "Costa" 2 "Sierra" 3 "Selva" 20 | lab val regnat regnat 21 | 22 | *Departamento 23 | destring ubigeo, generate(dpto) 24 | replace dpto=dpto/10000 25 | replace dpto=round(dpto) 26 | label variable dpto "Departamento" 27 | label define dpto 1 "Amazonas" 28 | label define dpto 2 "Ancash", add 29 | label define dpto 3 "Apurimac", add 30 | label define dpto 4 "Arequipa", add 31 | label define dpto 5 "Ayacucho", add 32 | label define dpto 6 "Cajamarca", add 33 | label define dpto 7 "Callao", add 34 | label define dpto 8 "Cusco", add 35 | label define dpto 9 "Huancavelica", add 36 | label define dpto 10 "Huanuco", add 37 | label define dpto 11 "Ica", add 38 | label define dpto 12 "Junin", add 39 | label define dpto 13 "La_Libertad", add 40 | label define dpto 14 "Lambayeque", add 41 | label define dpto 15 "Lima", add 42 | label define dpto 16 "Loreto", add 43 | label define dpto 17 "Madre_de_Dios", add 44 | label define dpto 18 "Moquegua", add 45 | label define dpto 19 "Pasco", add 46 | label define dpto 20 "Piura", add 47 | label define dpto 21 "Puno", add 48 | label define dpto 22 "San_Martin", add 49 | label define dpto 23 "Tacna", add 50 | label define dpto 24 "Tumbes", add 51 | label define dpto 25 "Ucayali", add 52 | label values dpto dpto 53 | 54 | *Grupos de edad 55 | gen edad=1 if p208a>17 & p208a<25 56 | replace edad=2 if p208a>24 & p208a<45 57 | replace edad=3 if p208a>44 58 | label define edad 1 "18-24" 2 "25-44" 3 "45 y mas" 59 | lab val edad edad 60 | 61 | *Nivel educativo 62 | gen educa=1 if p301a<5 | p301a==12 63 | replace educa=2 if p301a>4 & p301a<7 64 | replace educa=3 if p301a>6 & p301a<9 65 | replace educa=4 if p301a>8 & p301a<12 66 | label define educa 1 "Hasta primaria 1/" 2 "Secundaria" /// 67 | 3 "Superior no universitaria" 4 "Superior universitaria" 68 | label value educa educa 69 | 70 | *Etiquetar la variable pobreza monetaria 71 | *Usar variable "pobreza" de sumaria 72 | gen pobre=pobreza<3 73 | label define pobre 1 "pobre" 0 "no pobre" 74 | label value pobre pobre 75 | 76 | *Crear variable inclusion 77 | recode p558e6 (6=0 "no accede")(0=1 "accede"), gen(inclusion) 78 | 79 | *Establecer a los residentes habituales 80 | gen resi=1 if ((p204==1 & p205==2) | (p204==2 & p206==1)) 81 | 82 | *Perú: Inclusion financiera de la poblacion ocupada mayor a 18 ahnos 83 | tab inclusion [iw= fac500a ] if resi==1 & ocu500 == 1 84 | 85 | tab inclusion [iw= fac500a ] if resi==1 & ocu500 == 1 & p208a>17 86 | 87 | *Perú: Inclusion financiera de la poblacion ocupada mayor a 18 ahnos, por sexo 88 | tab p207 inclusion [iw= fac500a ] if resi==1 & ocu500 == 1, nofreq row 89 | 90 | *Perú: Inclusion financiera de la poblacion ocupada mayor a 18 ahnos, por edad 91 | tab edad inclusion [iw= fac500a ] if resi==1 & ocu500 == 1, nofreq row 92 | 93 | *Perú: Inclusion financiera de la poblacion ocupada mayor a 18 ahnos, por edu 94 | tab educa inclusion [iw= fac500a ] if resi==1 & ocu500 == 1, nofreq row 95 | 96 | *Perú: Inclusion financiera de la poblacion ocupada mayor a 18 ahnos, por pobreza 97 | tab pobre inclusion [iw= fac500a ] if resi==1 & ocu500 == 1, nofreq row 98 | 99 | *Perú: Inclusion financiera de la poblacion ocupada mayor a 18 ahnos, por area 100 | tab area inclusion [iw= fac500a ] if resi==1 & ocu500 == 1, nofreq row 101 | 102 | *Perú: Inclusion financiera de la poblacion ocupada mayor a 18 ahnos, por region 103 | tab regnat inclusion [iw= fac500a ] if resi==1 & ocu500 == 1, nofreq row 104 | 105 | *Perú: Inclusion financiera de la poblacion ocupada mayor a 18 ahnos, por dpto 106 | tab dpto inclusion [iw= fac500a ] if resi==1 & ocu500 == 1, nofreq row 107 | 108 | 109 | 110 | 111 | 112 | 113 | 114 | 115 | -------------------------------------------------------------------------------- /#94 MINEDU-ESCALE Tasa Total De Matrícula Edades 3-5, 6-11 & 12-16 con la ENAHO en STATA.do: -------------------------------------------------------------------------------- 1 | 2 | *MATRICULA 3 | cd "G:\ENAHO 2017" 4 | 5 | use enaho01a-2017-400.dta, clear 6 | merge 1:1 conglome vivienda hogar codperso using enaho01a-2017-300.dta 7 | keep if _merge==3 8 | drop _m 9 | 10 | /* 11 | MES Mes de procesamiento de la encuesta 12 | P203 ¿Cuál es la relación de parentesco con el jefe del hogar? 0 Panel (No Presente el año actual) 13 | 1 Jefe(a) 2 Esposo(a) 3 Hijo(a) 4 Yerno/Nuera 5 Nieto 6 Padres/Suegros 7 Otros parientes 14 | 8 Trabajador Hogar 9 Pensionista 10 Otros no parientes 15 | 16 | P204 ¿Es miembro del hogar? /1. Si 2. No /Rango 1-2 17 | 18 | P205 ¿Se encuentra ausente del hogar 30 días o más? /1. Si 2. No /Rango 1-2 19 | 20 | P206 ¿Está presente en el hogar 30 días o más? /1. Si 2. No /Rango 1-2 21 | 22 | P207 Sexo /1. Hombre 2. Mujer /Rango 1-2 23 | 24 | P300A ¿ Cuál es el idioma o lengua materna que aprendió en su niñez? 25 | 1 Quechua 2 Aymara 3 Otra lengua nativa 4 Castellano 6 Portugués 26 | 7 Otra lengua extranjera 8 Es sordomudo/a, mudo/a 9 Missing value 27 | 28 | P400A1 ¿En qué día, mes y año nació? - Día /99 Missing value /Rango : 1-31 29 | P400A2 ¿En qué día, mes y año nació? - Mes /99 Missing value /Rango : 1-12 30 | P400A3 ¿En qué día, mes y año nació? - Año /9999 Missing value /Rango : 1900 – 2021 31 | 32 | P306 Este año, ¿Está Matrículado en algún centro o programa de educación básica o superior? /1 Si 2 No 9 Missing value 33 | P307 Actualmente, ¿ Asiste a algún centro o programa de educación básica o superior? /1 Si 2 No 9 Missing value 34 | */ 35 | 36 | gen nihno=p207==1 37 | lab def nihno 0 "Niñas" 1 "Niños", modify 38 | lab val nihno nihno 39 | 40 | gen area=estrato 41 | recode area (1/5=1) (6/8=2) 42 | lab def area 1 "Urbana" 2 "Rural", modify 43 | lab val area area 44 | lab var area "Area de residencia" 45 | 46 | tab p300a, nol 47 | gen lengua=1 if p300a==4 48 | replace lengua=2 if p300a<4 49 | replace lengua=3 if p300a==6 | p300a==7 | p300a==8 50 | lab def lengua 1 "Castellano" 2 "Indígena" 3 "Otros", modify 51 | lab val lengua lengua 52 | 53 | destring ubigeo, generate(dpto) 54 | replace dpto=dpto/10000 55 | replace dpto=round(dpto) 56 | label variable dpto "Departamento" 57 | label define dpto 1 "Amazonas" 2 "Ancash" 3 "Apurimac" 4 "Arequipa" 5 "Ayacucho" /// 58 | 6 "Cajamarca" 7 "Callao" 8 "Cusco" 9 "Huancavelica" 10 "Huanuco" 11 "Ica" 12 "Junin" /// 59 | 13 "La_Libertad" 14 "Lambayeque" 15 "Lima" 16 "Loreto" 17 "Madre_de_Dios" 18 "Moquegua" /// 60 | 19 "Pasco" 20 "Piura" 21 "Puno" 22 "San_Martin" 23 "Tacna" 24 "Tumbes" 25 "Ucayali" 61 | label values dpto dpto 62 | 63 | *Se excluyen a los que participaron durante los meses de enero a marzo 64 | *-correspondiente a las vacaciones escolares y al primer mes de clases 65 | destring mes, replace 66 | gen matri_esc=1 if p306==1 & mes>=4 67 | replace matri_esc=2 if p306==2 & mes>=4 68 | label define matri_esc 1 "Matriculado" 2 "No_matriculado" 69 | label values matri_esc matri_esc 70 | label var matri_esc "Tasa total de matrícula" 71 | 72 | *Estimar edad en años cumplidos al 31 de Marzo 73 | tab a*o 74 | rename a*o year 75 | destring year, replace 76 | 77 | sum p400a1 p400a2 p400a3 78 | gen edad_31Marzo=year-p400a3 79 | replace edad_31Marzo=(year-p400a3-1) if p400a2>3 80 | replace edad_31Marzo=0 if edad_31Marzo<0 81 | replace edad_31Marzo=. if p400a2==. 82 | 83 | 84 | *Usamos el factor de expansion de la base para obtener los resultados a nivel nacional 85 | *3-5 86 | table nihno area [iweight=factora07] if (edad_31Marzo>=3 & edad_31Marzo<=5) & p203!=8 & p204==1 , stat(fvpercent matri_esc) nformat(%5.3f) 87 | table lengua [iweight=factora07] if (edad_31Marzo>=3 & edad_31Marzo<=5) & p203!=8 & p204==1 , stat(fvpercent matri_esc) nformat(%5.3f) 88 | 89 | *6-11 90 | table nihno area [iweight=factora07] if (edad_31Marzo>=6 & edad_31Marzo<=11) & p203!=8 & p204==1 , stat(fvpercent matri_esc) nformat(%5.3f) 91 | table lengua [iweight=factora07] if (edad_31Marzo>=6 & edad_31Marzo<=11) & p203!=8 & p204==1 , stat(fvpercent matri_esc) nformat(%5.3f) 92 | 93 | *12-16 94 | table nihno area [iweight=factora07] if (edad_31Marzo>=12 & edad_31Marzo<=16) & p203!=8 & p204==1 , stat(fvpercent matri_esc) nformat(%5.3f) 95 | table lengua [iweight=factora07] if (edad_31Marzo>=12 & edad_31Marzo<=16) & p203!=8 & p204==1 , stat(fvpercent matri_esc) nformat(%5.3f) 96 | 97 | -------------------------------------------------------------------------------- /#63 Tenencia de Celular en el Hogar por Quintiles de Ingreso.do: -------------------------------------------------------------------------------- 1 | /* Sintaxis extraida del documento de Trabajo Nº 21-2014 2 | INDICADORES SOCIALES Y DE DESIGUALDAD EN BASE A LA ENAHO CON STATA 3 | por Alfonso Ayala- Diciembre 22, 2014 4 | ISSN 2312-4776 5 | https://economia.unmsm.edu.pe/data/doc_trab/21-2014-OBG.pdf 6 | 7 | I. INTRODUCCIÓN 8 | II. ENTORNO DE TRABAJO EN STATA 9 | III.ANALISIS DE VARIABLES EN ENCUESTA NACIONAL DE HOGARES** 10 | IV. BREVE INTRODUCCIÓN A LA PROGRAMACIÓN EN STATA 11 | V. ANÁLISIS DE LA DISTRIBUCIÓN DEL INGRESO EN LA ENAHO */ 12 | 13 | *Establecer carpeta de trabajo 14 | cd "D:\ENAHO2013" 15 | 16 | *Unir modulo de la vivienda con archivo sumaria 17 | use enaho01-2013-100.dta, clear 18 | merge 1:1 conglome vivienda hogar using sumaria-2013.dta 19 | keep if _m==3 20 | 21 | *Calculamos el factor de ponderación a nivel de la población 22 | gen factor_personas=factor*mieperho 23 | 24 | *Cálculo de la variable region natural (pag 18) 25 | gen region_natural=. 26 | label variable region_natural "Región Natural" 27 | replace region_natural=1 if dominio<=3|dominio==8 28 | replace region_natural=2 if dominio>=4& dominio<=6 29 | replace region_natural=3 if dominio==7 30 | label define etiq_region_natural 1 "Costa" 2 "Sierra" 3 "Selva" 31 | label values region_natural etiq_region_natural 32 | 33 | *Cálculo de la variable region (pag 12/13) 34 | destring ubigeo, generate(region) 35 | replace region=region/10000 36 | replace region=round(region) 37 | 38 | label define etiq_region 1 "Amazonas" 2 "Ancash" 3 "Apurimac" 4 "Arequipa" /// 39 | 5 "Ayacucho" 6 "Cajamarca" 7 "Callao" 8 "Cusco" 9 "Huancavelica" 10 "Huanuco" /// 40 | 11 "Ica" 12 "Junin" 13 "La Libertad" 14 "Lambayeque" 15 "Lima" 16 "Loreto" /// 41 | 17 "Madre de Dios" 18 "Moquegua" 19 "Pasco" 20 "Piura" 21 "Puno" /// 42 | 22 "San Martín" 23 "Tacna" 24 "Tumbes" 25 "Ucayali" 43 | label values region etiq_region 44 | 45 | *Pobreza Monetaria Total 46 | *Estimamos el gasto mensual promedio de los hogares en terminos per capita 47 | *(ver video de estimacion de pobreza) 48 | gen gpcm=gashog2d/(mieperho*12) 49 | *Contrastamos con la linea de pobreza total (linea) y etiquetamos la variable 50 | gen pobre_2=1 if gpcm < linea 51 | replace pobre_2=0 if gpcm >= linea 52 | *etiquetamos la variable y los valores que toma la variable 53 | label var pobre_2 "Pobreza Monetaria Total" 54 | label define pobre_2 1 "pobre" 0 "no_pobre" 55 | label value pobre_2 pobre_2 56 | 57 | *Pobreza a nivel regional (pag 13/14) 58 | table region [iw=factor_personas], c(mean pobre_2) 59 | *Para calcular el indicador de pobreza nacional se agrega la opción row (pag 16) 60 | table region [iw=factor_personas], c(mean pobre_2) row 61 | 62 | *Graficar (pag 15/16) 63 | graph hbar pobre_2 [w=factor_personas], over(region, sort((mean) pobre_2)) 64 | graph hbar pobre_2 [w=factor_personas], over(region, sort((mean) pobre_2)) blabel(bar, format(%5.2f)) 65 | graph save Graph "Graph1.gph", replace 66 | graph export "Graph1.png", as(png) replace 67 | 68 | *Uso de la opción by(varlist) (pag 17) 69 | /*Mediante la opción by(varlist) Stata calcula los estadísticos según los grupos 70 | formados por la variable varlist. Por ejemplo: los estadísticos número de observaciones, 71 | la media, el valor mínimo y el valor máximo del ingreso bruto del hogar (inghog1d) 72 | por cada Región para el año 2013 se obtiene mediante: */ 73 | tabstat inghog1d [w=factor07], stat(n mean min max) by(region) 74 | 75 | *Cálculo de quintiles de ingreso (pag 23) 76 | /* xtile calcula los nq quintiles de la variable ingreso neto total (inghog2d), 77 | en este caso se pondera por el factor de expansión factor07. */ 78 | xtile quintiles_ing = inghog2d [w=factor07], nq(5) 79 | 80 | *Tenencia de bienes en la vivienda (celular) (pag 23/24) 81 | *Usamos el módulo 100 de la ENAHO (Características de la Vivienda y del Hogar) 82 | lookfor celular 83 | gen TieneCelular = p1142 84 | 85 | *Etiquetamos la variable TieneCelular y los valores de la variable: 86 | label variable TieneCelular "Tenencia de celular en el hogar" 87 | label define etiq_celular 1 "Tiene cel" 0 "No tiene" 88 | label values TieneCelular etiq_celular 89 | 90 | *Calculamos la tenencia de celular por hogar: 91 | sum TieneCelular [iw=factor07] 92 | tab TieneCelular [iw=factor07] 93 | 94 | *Calculamos la tenencia de celular por regiones naturales y regiones: 95 | *Por regiones naturales 96 | table region_natural [iw=factor07], c(mean TieneCelular) row 97 | *Otra forma de verlo 98 | tab region_natural TieneCelular [iw=factor07], nofreq row 99 | *Por regiones 100 | table region [iw=factor07], c(mean TieneCelular) row 101 | 102 | *Calculamos los quintiles de ingreso según el acceso a algún servicio o 103 | *tenencia de algún bien (pag 27) 104 | table quintiles_ing [iw=factor07], c(mean TieneCelular) 105 | 106 | -------------------------------------------------------------------------------- /#19 Cálculo del ingreso promedio mensual proveniente del trabajo.do: -------------------------------------------------------------------------------- 1 | 2 | 3 | Estimad@s, en este video comparto el cálculo del ingreso promedio mensual proveniente del trabajo con la Encuesta Nacional de Hogares (ENAHO) del 2017 en STATA. El ingreso promedio por trabajo corresponde a la PEA ocupada con ingresos mayores a cero y que provienen de su actividad principal, actividad secundaria, dependiente e independiente y puede ser monetario o no monetario. 4 | 5 | Archivo: "enaho01a-2017-500.dta" que se ubica dentro del Módulo 500 de la Encuesta Nacional de Hogares (ENAHO) 2017 6 | 7 | Variables que utilizaremos: 8 | fac500a Factor de Expansión de Empleo/Ingresos proyecciones CPV-2007 9 | 10 | ocu500 11 | 1 Ocupado 12 | 2 Desocupado Abierto 13 | 3 Desocupado Oculto 14 | 4 No PEA 15 | 16 | I524A1 Ingreso total trimestral (Imputado, deflactado, Anualizado) 17 | D529T Pago en especie dependiente (Deflactado, Anualizado) 18 | I530A Ganancia (ocupación principal independiente) (Imputado, deflactado, Anualizado) 19 | D536 Valor de los productos para su consumo (Deflactado, Anualizado) 20 | I538A1 Ingreso total (Imputado, deflactado, Anualizado) 21 | D540T Pago en especie (dependiente) (Deflactado, Anualizado) 22 | I541A Ganancia (ocupación secundaria independiente) (Imputado, deflactado, Anualizado) 23 | D543 Valor de los productos utilizados para su consumo (Deflactado, Anualizado) 24 | D544T Ingreso extraordinario (Deflactado, Anualizado) 25 | 26 | P204 ¿Es miembro del hogar? 27 | P205 ¿Se encuentra ausente del hogar 30 días o más? 28 | P206 ¿Está presente en el hogar 30 días o más? 29 | 30 | Fuente: INEI, Perú: Evolución de los Indicadores de Empleo e Ingreso por Departamento, 2007-2018. 31 | Janet Porras 32 | ******* 33 | /* 34 | INGRESO PROMEDIO MENSUAL PROVENIENTE DEL TRABAJO - 2017 35 | El ingreso promedio corresponde a la PEA ocupada con ingresos mayores a cero y 36 | que provienen de su actividad principal, actividad secundaria, 37 | dependiente e independiente y puede ser monetario o no monetario. 38 | 39 | https://www.inei.gob.pe/media/MenuRecursivo/publicaciones_digitales/Est/Lib1678/libro.pdf 40 | 41 | Archivo: 42 | ******** 43 | Usamos el modulo 500 de la ENAHO 2017 44 | 45 | Variables que utilizaremos: 46 | *************************** 47 | 48 | fac500a Factor de Expansión de Empleo/Ingresos proyecciones CPV-2007 49 | 50 | ocu500 51 | 1 Ocupado 52 | 2 Desocupado Abierto 53 | 3 Desocupado Oculto 54 | 4 No PEA 55 | 56 | 57 | I524A1 Ingreso total trimestral (Imputado, deflactado, Anualizado) 58 | D529T Pago en especie dependiente (Deflactado, Anualizado) 59 | I530A Ganancia (ocupación principal independiente) (Imputado, deflactado, Anualizado) 60 | D536 Valor de los productos para su consumo (Deflactado, Anualizado) 61 | I538A1 Ingreso total (Imputado, deflactado, Anualizado) 62 | D540T Pago en especie (dependiente) (Deflactado, Anualizado) 63 | I541A Ganancia (ocupación secundaria independiente) (Imputado, deflactado, Anualizado) 64 | D543 Valor de los productos utilizados para su consumo (Deflactado, Anualizado) 65 | D544T Ingreso extraordinario (Deflactado, Anualizado) 66 | 67 | P204 ¿Es miembro del hogar? 68 | 69 | P205 ¿Se encuentra ausente del hogar 30 días o más? 70 | 71 | P206 ¿Está presente en el hogar 30 días o más? 72 | 73 | */* 74 | 75 | 76 | use "D:\enaho01a-2017-500.dta", clear 77 | gen area=1 if estrato<=5 78 | replace area=2 if estrato>=6 & estrato<=8 79 | lab def area 1 "urbano" 2 "rural" 80 | lab val area area 81 | 82 | destring ubigeo, generate(dpto) 83 | replace dpto=dpto/10000 84 | replace dpto=round(dpto) 85 | label variable dpto "Departamento" 86 | label define dpto 1 "Amazonas" 87 | label define dpto 2 "Ancash", add 88 | label define dpto 3 "Apurimac", add 89 | label define dpto 4 "Arequipa", add 90 | label define dpto 5 "Ayacucho", add 91 | label define dpto 6 "Cajamarca", add 92 | label define dpto 7 "Callao", add 93 | label define dpto 8 "Cusco", add 94 | label define dpto 9 "Huancavelica", add 95 | label define dpto 10 "Huanuco", add 96 | label define dpto 11 "Ica", add 97 | label define dpto 12 "Junin", add 98 | label define dpto 13 "La_Libertad", add 99 | label define dpto 14 "Lambayeque", add 100 | label define dpto 15 "Lima", add 101 | label define dpto 16 "Loreto", add 102 | label define dpto 17 "Madre_de_Dios", add 103 | label define dpto 18 "Moquegua", add 104 | label define dpto 19 "Pasco", add 105 | label define dpto 20 "Piura", add 106 | label define dpto 21 "Puno", add 107 | label define dpto 22 "San_Martin", add 108 | label define dpto 23 "Tacna", add 109 | label define dpto 24 "Tumbes", add 110 | label define dpto 25 "Ucayali", add 111 | label values dpto dpto 112 | 113 | egen ingtrabw= rowtotal(i524a1 d529t i530a d536 i538a1 d540t i541a d543 d544t) 114 | gen ingtra_n= ingtrabw/12 115 | label var ingtrabw "ingreso por trabajo anual" 116 | label var ingtra_n "ingreso por trabajo mensual" 117 | 118 | *Establecer a los residentes habituales 119 | gen resi=1 if ((p204==1 & p205==2) | (p204==2 & p206==1)) 120 | 121 | keep if ocu500 == 1 & ingtra_n > 0 122 | 123 | table area [iw=fac500a] if resi==1, c(mean ingtra_n) row 124 | table dpto [iw=fac500a] if resi==1, c(mean ingtra_n) row 125 | -------------------------------------------------------------------------------- /#29 Calculo de la población ocupada por rama de actividad y otros criterios.do: -------------------------------------------------------------------------------- 1 | use "D:\ENAHO\sumaria-2018.dta", clear 2 | *Creamos la variable factor de expansion de la poblacion 3 | gen facpob=factor07*mieperho 4 | 5 | *Usando la variable pobreza generada por el INEI //con sumaria 6 | generate pobreza2=1 if pobreza<3 7 | replace pobreza2=2 if pobreza==3 8 | label define pobreza2 1 "Pobre" 2 "No pobre" 9 | label value pobreza2 pobreza2 10 | 11 | merge 1:m conglome vivienda hogar using "D:\ENAHO\enaho01a-2018-500.dta" 12 | drop if p500i=="00" 13 | 14 | *Area de residencia 15 | gen area=1 if estrato<=5 16 | replace area=2 if estrato>=6 & estrato<=8 17 | lab def area 1 "urbano" 2 "rural" 18 | lab val area area 19 | 20 | *Region natural 21 | gen region=1 if dominio>=1 & dominio<=3 22 | replace region=1 if dominio==8 23 | replace region=2 if dominio>=4 & dominio<=6 24 | replace region=3 if dominio==7 25 | label define region 1 "Costa" 2 "Sierra" 3 "Selva" 26 | lab val region region 27 | 28 | *Se establece quienes son residentes habituales 29 | gen resi=1 if ((p204==1 & p205==2) | (p204==2 & p206==1)) 30 | 31 | *Población ocupada en empleo informal 32 | tab ocu500 ocupinf [iw= fac500a] if resi==1 33 | 34 | *Población ocupada en empleo informal por Condición de pobreza 35 | tab pobreza2 ocupinf [iw=fac500a] if resi==1 & ocu500==1, nofreq row 36 | 37 | *Población ocupada por tamahno de empresa 38 | gen tamahno=1 if p512b>=1 & p512b<11 39 | replace tamahno=2 if p512b>=11 & p512b<51 40 | replace tamahno=3 if p512b>50 41 | replace tamahno=4 if p512b==. & (p512a==1 | p512a==2 ) 42 | label define tamahno 1 "De 1 a 10 trabajadores" 2 "De 11 a 50 trabajadores" /// 43 | 3 "De 51 a más trabajadores" 4 "No especificado", replace 44 | label value tamahno tamahno 45 | 46 | tab tamahno area [iw=fac500a] if resi==1 & ocu500==1 47 | tab tamahno region [iw=fac500a] if resi==1 & ocu500==1 48 | 49 | 50 | *Población ocupada, según ramas de actividad 51 | *CIIU 52 | gen ciiu_aux1 =substr("0"+string(p506r4),1,.) 53 | replace ciiu_aux1 =substr(string(p506r4),1,.) if p506r4>999 54 | gen ciiu_aux2 =substr(ciiu_aux1 ,1,2) 55 | destring ciiu_aux2, generate(ciiu_2d) 56 | gen ciiu_1d=1 if ciiu_2d<=2 57 | replace ciiu_1d=2 if ciiu_2d==3 58 | replace ciiu_1d=3 if ciiu_2d>=5 & ciiu_2d<=9 59 | replace ciiu_1d=4 if ciiu_2d>=10 & ciiu_2d<=33 60 | replace ciiu_1d=5 if ciiu_2d>=41 & ciiu_2d<=43 61 | replace ciiu_1d=6 if ciiu_2d>=45 & ciiu_2d<=47 62 | replace ciiu_1d=7 if (ciiu_2d>=49 & ciiu_2d<=53) | (ciiu_2d>=58 & ciiu_2d<=63) 63 | replace ciiu_1d=8 if ciiu_2d==84 64 | replace ciiu_1d=9 if ciiu_2d>=55 & ciiu_2d<=56 65 | replace ciiu_1d=10 if ciiu_2d==68 | (ciiu_2d>=69 & ciiu_2d<=82) 66 | replace ciiu_1d=11 if ciiu_2d==85 67 | replace ciiu_1d=12 if (ciiu_2d>=35 & ciiu_2d<=39) | (ciiu_2d>=64 & ciiu_2d<=66) | /// 68 | (ciiu_2d>=86 & ciiu_2d<=88) | (ciiu_2d>=90 & ciiu_2d<=93)| (ciiu_2d>=94 & ciiu_2d<=98) | ciiu_2d==99 69 | 70 | label var ciiu_1d "Division CIIU" 71 | la de ciiu_1d 1 "Agricultura" 2 "Pesca" 3 "Mineria" 4 "Manufactura" 5 "Construccion" /// 72 | 6 "Comercio" 7 "Transportes y Comunicaciones" 8 "Gobierno" 9 "Hoteles y Restaurantes" /// 73 | 10 "Inmobiliarias y alquileres" 11 "Ensehnanza" 12 "Otros Servicios 1/" 74 | label values ciiu_1d ciiu_1d 75 | *1/ Otros Servicios lo componen las ramas de actividad de Electricidad, Gas y Agua, 76 | *Intermediación Financiera, Actividades de Servicios Sociales y de Salud, Otras activ. 77 | *de Serv. Comunitarias, Sociales y Personales y Hogares privados con servicio doméstico. 78 | 79 | tab ciiu_1d [iw= fac500a] if resi==1 & ocu500==1, m 80 | 81 | *Población ocupada en empleo informal por Rama de Actividad 82 | gen ciiu_6c=1 if ciiu_1d<4 83 | replace ciiu_6c=2 if ciiu_1d==4 84 | replace ciiu_6c=3 if ciiu_1d==5 85 | replace ciiu_6c=4 if ciiu_1d==6 86 | replace ciiu_6c=5 if ciiu_1d==7 87 | replace ciiu_6c=6 if ciiu_1d>7 88 | 89 | label var ciiu_6c "Division CIIU-6 categorias" 90 | la de ciiu_6c 1 "Agricultura/Pesca/Mineria" 2 "Manufactura" 3 "Construccion" /// 91 | 4 "Comercio" 5 "Transportes y Comunicaciones" 6 "Otros Servicios 1/" 92 | label values ciiu_6c ciiu_6c 93 | *1/ Otros Servicios lo componen las ramas de actividad de Electricidad, Gas y Agua, 94 | *Intermediación Financiera, Actividades de Servicios Sociales y de Salud, Otras activ. 95 | *de Serv. Comunitarias, Sociales y Personales y Hogares privados con servicio doméstico. 96 | *Adicionalmente incluye Gobierno, Hoteles y Restaurantes, Inmobiliarias y alquileres y Ensehnanza 97 | 98 | tab ciiu_6c ocupinf [iw= fac500a] if resi==1 & ocu500==1, nofreq row 99 | 100 | drop _m 101 | ******************************************************************************** 102 | merge 1:1 conglome vivienda hogar codperso using "D:\ENAHO\enaho01a-2018-400.dta" 103 | keep if _m==3 104 | drop _m 105 | 106 | recode p4191 p4192 p4193 p4194 p4195 p4196 p4197 p4198 (2=0) 107 | *no incluyo a los que tienen seguro escolar 108 | gen seguro= p4191 + p4192 + p4193 + p4194 + p4195 + p4196 + p4198 109 | tab seguro 110 | replace seguro=1 if seguro>=1 & seguro<=3 111 | 112 | tab ocu500 seguro [iw=fac500a] if resi==1 & ocu500==1, row 113 | 114 | tab area seguro [iw=fac500a] if resi==1 & ocu500==1, row 115 | -------------------------------------------------------------------------------- /#30 Cálculo del ingreso promedio mensual proveniente del trabajo.do: -------------------------------------------------------------------------------- 1 | cd "D:\ENAHO" 2 | use enaho01a-2018-500.dta, clear 3 | 4 | *Area urbana/rural 5 | gen area=1 if estrato<=5 6 | replace area=2 if estrato>=6 & estrato<=8 7 | lab def area 1 "urbano" 2 "rural" 8 | lab val area area 9 | 10 | *Departamento 11 | destring ubigeo, generate(dpto) 12 | replace dpto=dpto/10000 13 | replace dpto=round(dpto) 14 | label variable dpto "Departamento" 15 | label define dpto 1 "Amazonas" 16 | label define dpto 2 "Ancash", add 17 | label define dpto 3 "Apurimac", add 18 | label define dpto 4 "Arequipa", add 19 | label define dpto 5 "Ayacucho", add 20 | label define dpto 6 "Cajamarca", add 21 | label define dpto 7 "Callao", add 22 | label define dpto 8 "Cusco", add 23 | label define dpto 9 "Huancavelica", add 24 | label define dpto 10 "Huanuco", add 25 | label define dpto 11 "Ica", add 26 | label define dpto 12 "Junin", add 27 | label define dpto 13 "La_Libertad", add 28 | label define dpto 14 "Lambayeque", add 29 | label define dpto 15 "Lima", add 30 | label define dpto 16 "Loreto", add 31 | label define dpto 17 "Madre_de_Dios", add 32 | label define dpto 18 "Moquegua", add 33 | label define dpto 19 "Pasco", add 34 | label define dpto 20 "Piura", add 35 | label define dpto 21 "Puno", add 36 | label define dpto 22 "San_Martin", add 37 | label define dpto 23 "Tacna", add 38 | label define dpto 24 "Tumbes", add 39 | label define dpto 25 "Ucayali", add 40 | label values dpto dpto 41 | 42 | *Etnia 43 | gen etnia=1 if p558c<4 | p558c==9 44 | replace etnia=2 if p558c==4 45 | replace etnia=3 if p558c==6 46 | replace etnia=4 if p558c==5 | p558c==7 47 | replace etnia=5 if p558c==8 48 | lab def etnia 1 "Indigena 1/" 2 "Negro, mulato, Afro peruano" 3 "Mestizo/a" 4 "Otro 2/" 5 "No sabe" 49 | lab val etnia etnia 50 | *1/ Incluye: Quechua, Aimara y Nativo o Indígena de la Amazonía. 51 | *2/ Incluye: Blanco y otro 52 | 53 | *Grupos de edad 54 | gen edad=1 if p208a>13 & p208a<25 55 | replace edad=2 if p208a>24 & p208a<45 56 | replace edad=3 if p208a>44 & p208a<60 57 | replace edad=4 if p208a>59 & p208a<65 58 | replace edad=5 if p208a>64 59 | label define edad 1 "14-24" 2 "25-44" 3 "45-59" 4 "60-64" 5 "65 y mas" 60 | lab val edad edad 61 | 62 | *Estado civil 63 | gen ecivil=1 if p209==1 64 | replace ecivil=2 if p209==2 65 | replace ecivil=3 if p209==3 | p209==4 | p209==5 66 | replace ecivil=4 if p209==6 67 | lab def ecivil 1 "Conviviente" 2 "Casado/a" 3 "Alguna vez unido/a 1/" 4 " Soltero/a" 68 | lab val ecivil ecivil 69 | *1/ Incluye: Separado/a, divorciado/a y viudo/a. 70 | 71 | *Establecer a los residentes habituales 72 | gen resi=1 if ((p204==1 & p205==2) | (p204==2 & p206==1)) 73 | 74 | *Población ocupada por tamahno de empresa 75 | gen tamahno=1 if p512b>=1 & p512b<11 76 | replace tamahno=2 if p512b>=11 & p512b<51 77 | replace tamahno=3 if p512b>50 78 | replace tamahno=4 if p512b==. & (p512a==1 | p512a==2 ) 79 | label define tamahno 1 "De 1 a 10 trabajadores" 2 "De 11 a 50 trabajadores" /// 80 | 3 "De 51 a más trabajadores" 4 "No especificado", replace 81 | label value tamahno tamahno 82 | 83 | /* Ingreso proveniente del trabajo 84 | i524a1 Ingreso total trimestral (Imputado, deflactado, Anualizado) 85 | d529t Pago en especie dependiente (Deflactado, Anualizado) 86 | i530a Ganancia (ocupación principal independiente) (Imputado, deflactado, Anualizado) 87 | d536 Valor de los productos para su consumo (Deflactado, Anualizado) 88 | i538a1 Ingreso total (Imputado, deflactado, Anualizado) 89 | d540t Pago en especie (dependiente) (Deflactado, Anualizado) 90 | i541a Ganancia (ocupación secundaria independiente) (Imputado, deflactado, Anualizado) 91 | d543 Valor de los productos utilizados para su consumo (Deflactado, Anualizado) 92 | d544t Ingreso extraordinario (Deflactado, Anualizado) */ 93 | 94 | egen ingtrabw= rowtotal(i524a1 d529t i530a d536 i538a1 d540t i541a d543 d544t) 95 | gen ingtra_n= ingtrabw/12 96 | label var ingtrabw "ingreso por trabajo anual" 97 | label var ingtra_n "ingreso por trabajo mensual" 98 | keep if ocu500 == 1 & ingtra_n > 0 99 | 100 | *Perú: Ingreso promedio mensual proveniente del trabajo de 101 | *la población ocupada, según etnia 102 | table etnia [iw=fac500a] if resi==1 , c(mean ingtra_n) row 103 | 104 | *Perú: Ingreso promedio mensual proveniente del trabajo, 105 | *según estado civil o conyugal 106 | table ecivil [iw=fac500a] if resi==1 , c(mean ingtra_n) row 107 | 108 | *Perú: Ingreso promedio mensual proveniente del trabajo, 109 | *según tamaño de empresa y área urbana 110 | table tamahno area [iw=fac500a] if resi==1 , c(mean ingtra_n) row 111 | 112 | *Perú: Ingreso promedio mensual proveniente del trabajo de hombres y mujeres, 113 | *según área de residencia 114 | table area p207 [iw=fac500a] if resi==1 , c(mean ingtra_n) row 115 | 116 | *Perú urbano: Ingreso promedio mensual proveniente del trabajo, 117 | *según grupos de edad 118 | table edad [iw=fac500a] if resi==1 & area==1 , c(mean ingtra_n) row 119 | 120 | *Perú: Ingreso promedio mensual proveniente del trabajo 121 | *de hombres y mujeres y Brechas, según grupos de edad 122 | table edad p207 [iw=fac500a] if resi==1 , c(mean ingtra_n) row 123 | 124 | *Perú: Brechas de género en ingreso promedio proveniente 125 | *del trabajo, según departamento 126 | table dpto p207 [iw=fac500a] if resi==1, c(mean ingtra_n) row 127 | -------------------------------------------------------------------------------- /#49 Medición de la tasa de analfabetismo en el Perú con la Encuesta Nacional de Hogares (ENAHO) 2019 en STATA.do: -------------------------------------------------------------------------------- 1 | 2 | 3 | *Cambiar el nombre segun el folder donde van a bajar la informacion 4 | cd "D:\EDUCACION" 5 | 6 | 7 | *Guardar los pasos en un archivo log 8 | log using analfabetismo_enaho2019, replace 9 | 10 | *Abrir archivo 11 | use enaho01a-2019-300.dta, clear 12 | 13 | *Grupo de edad 14 | gen gedad=1 if p208a>=15 & p208a<20 15 | replace gedad=2 if p208a>=20 & p208a<30 16 | replace gedad=3 if p208a>=30 & p208a<40 17 | replace gedad=4 if p208a>=40 & p208a<50 18 | replace gedad=5 if p208a>=50 & p208a<60 19 | replace gedad=6 if p208a>=60 20 | lab def gedad 1 "De 15 a 19 años" 2 "De 20 a 29 años" 3 "De 30 a 39 años" /// 21 | 4 "De 40 a 49 años" 5 "De 50 a 59 años" 6 "De 60 y más años" 22 | lab val gedad gedad 23 | 24 | *Area de residencia 25 | gen area=estrato 26 | recode area (1/5=1) (6/8=2) 27 | lab def area 1 "Urbana" 2 "Rural", modify 28 | lab val area area 29 | lab var area "Area de residencia" 30 | 31 | *Region Natural 32 | gen regnat=1 if dominio<=3 | dominio==8 33 | replace regnat=2 if dominio>=4 & dominio<=6 34 | replace regnat=3 if dominio==7 35 | lab var regnat "Region natural" 36 | lab def regnat 1 "Costa" 2 "Sierra" 3 "Selva" 37 | lab val regnat regnat 38 | 39 | *Dominio geografico 40 | gen dom=1 if regnat==1 & area==1 41 | replace dom=2 if regnat==1 & area==2 42 | replace dom=3 if regnat==2 & area==1 43 | replace dom=4 if regnat==2 & area==2 44 | replace dom=5 if regnat==3 & area==1 45 | replace dom=6 if regnat==3 & area==2 46 | replace dom=7 if dominio==8 47 | lab def ldom 1 "Costa urbana" 2 "Costa rural" 3 "Sierra urbana" 4 "Sierra rural" /// 48 | 5 "Selva urbana" 6 "Selva rural" 7 "Lima Metropolitana" 49 | lab val dom ldom 50 | 51 | *Departamento (25) 52 | destring ubigeo, generate(dpto) 53 | replace dpto=dpto/10000 54 | replace dpto=round(dpto) 55 | label variable dpto "dpto" 56 | label define dpto 1 "Amazonas" 57 | label define dpto 2 "Ancash", add 58 | label define dpto 3 "Apurimac", add 59 | label define dpto 4 "Arequipa", add 60 | label define dpto 5 "Ayacucho", add 61 | label define dpto 6 "Cajamarca", add 62 | label define dpto 7 "Callao", add 63 | label define dpto 8 "Cusco", add 64 | label define dpto 9 "Huancavelica", add 65 | label define dpto 10 "Huanuco", add 66 | label define dpto 11 "Ica", add 67 | label define dpto 12 "Junin", add 68 | label define dpto 13 "La_Libertad", add 69 | label define dpto 14 "Lambayeque", add 70 | label define dpto 15 "Lima", add 71 | label define dpto 16 "Loreto", add 72 | label define dpto 17 "Madre_de_Dios", add 73 | label define dpto 18 "Moquegua", add 74 | label define dpto 19 "Pasco", add 75 | label define dpto 20 "Piura", add 76 | label define dpto 21 "Puno", add 77 | label define dpto 22 "San_Martin", add 78 | label define dpto 23 "Tacna", add 79 | label define dpto 24 "Tumbes", add 80 | label define dpto 25 "Ucayali", add 81 | label values dpto dpto 82 | 83 | *Departamento (distinguiendo Provincia de Lima/1 de Region Lima/2) 84 | *1/Comprende los 43 distritos que conforman la provincia de Lima 85 | *2/Incluye las provincias de: Barranca, Cajatambo, Canta, Cañete, Huaral, Huarochiri, Huaura, Oyón y Yauyos. 86 | gen dpto_26=dpto*10 87 | replace dpto_26=151 if dom==7 88 | replace dpto_26=152 if dpto==15 & dom~=7 89 | label variable dpto_26 "Dpto (dist. Provincia de Lima de Region Lima)" 90 | label define dpto_26 10 "Amazonas" 91 | label define dpto_26 20 "Ancash", add 92 | label define dpto_26 30 "Apurimac", add 93 | label define dpto_26 40 "Arequipa", add 94 | label define dpto_26 50 "Ayacucho", add 95 | label define dpto_26 60 "Cajamarca", add 96 | label define dpto_26 70 "Callao", add 97 | label define dpto_26 80 "Cusco", add 98 | label define dpto_26 90 "Huancavelica", add 99 | label define dpto_26 100 "Huanuco", add 100 | label define dpto_26 110 "Ica", add 101 | label define dpto_26 120 "Junin", add 102 | label define dpto_26 130 "La_Libertad", add 103 | label define dpto_26 140 "Lambayeque", add 104 | label define dpto_26 151 "Provincia de Lima", add 105 | label define dpto_26 152 "Region Lima", add 106 | label define dpto_26 160 "Loreto", add 107 | label define dpto_26 170 "Madre_de_Dios", add 108 | label define dpto_26 180 "Moquegua", add 109 | label define dpto_26 190 "Pasco", add 110 | label define dpto_26 200 "Piura", add 111 | label define dpto_26 210 "Puno", add 112 | label define dpto_26 220 "San_Martin", add 113 | label define dpto_26 230 "Tacna", add 114 | label define dpto_26 240 "Tumbes", add 115 | label define dpto_26 250 "Ucayali", add 116 | label values dpto_26 dpto_26 117 | 118 | *ANALFABETISMO 119 | gen analfa=0 if p208a>=15 & p204==1 120 | replace analfa=1 if p208a>=15 & p302==2 & p204==1 121 | label define analfa 0 "Si" 1 "No" 122 | label values analfa analfa 123 | label var analfa "Sabe leer y escribir" 124 | tab analfa [iweight=factora07] 125 | tab area analfa [iweight=factora07] 126 | tab area analfa [iweight=factora07], nofreq row 127 | tab regnat analfa [iweight=factora07], nofreq row 128 | tab dpto_26 analfa [iweight=factora07], nofreq row 129 | tab gedad analfa [iweight=factora07], nofreq row 130 | tab gedad analfa [iweight=factora07] if p207==1, nofreq row 131 | tab gedad analfa [iweight=factora07] if p207==2, nofreq row 132 | 133 | *Cerramo el archivo log 134 | log close -------------------------------------------------------------------------------- /#14 Estimación de las Necesidades Básicas Insatisfechas.do: -------------------------------------------------------------------------------- 1 | cd "D:\ENAHO 2017" 2 | 3 | *Bajamos los archivos zipeados del modulo 100 y sumaria de la pagina del INEI 4 | copy "http://iinei.inei.gob.pe/iinei/srienaho/descarga/STATA/603-Modulo01.zip" 603-Modulo01.zip, replace 5 | copy "http://iinei.inei.gob.pe/iinei/srienaho/descarga/STATA/603-Modulo34.zip" 603-Modulo34.zip, replace 6 | 7 | *Descomprimir los archivos del modulo 100 y sumaria "Manualmente" 8 | 9 | use enaho01-2017-100.dta, clear 10 | /* result: resultado final de la encuesta 11 | 1: completa 12 | 2: incompleta 13 | 3: rechazo 14 | 4: ausente 15 | 5: vivienda desocupada 16 | 6: otro */ 17 | 18 | *Se trabaja solo con las encuestas completas e incompletas 19 | drop if result>2 20 | 21 | *NECESIDADES BASICAS INSATISFECHAS (ya se encuentran en el modulo 100) 22 | sum nbi* 23 | 24 | collapse (mean) nbi1 nbi2 nbi3 nbi4 nbi5, by(conglome vivienda hogar) cw 25 | 26 | *Juntamos el modulo 100 con el modulo sumaria 27 | *(ambas bases presentan informacion a nivel del hogar) 28 | merge 1:1 conglome vivienda hogar using sumaria-2017.dta, nogenerate 29 | 30 | *Creamos la variable factor de expansion de la poblacion 31 | gen facpob=factor07*mieperho 32 | 33 | *Establecemos las caracteristicas de la encuesta 34 | *usando las variable factor de expansion, conglomerado y estrato 35 | svyset [pweight=facpob], psu(conglome) strata(estrato) 36 | 37 | 38 | gen nbihog=nbi1 + nbi2 + nbi3 + nbi4 + nbi5 39 | 40 | gen NBI1_POBRE=. 41 | replace NBI1_POBRE=1 if nbihog>0 42 | replace NBI1_POBRE=0 if nbihog==0 43 | 44 | label define NBI1_POBRE 0 "ninguna NBI" 1 "al menos un NBI" 45 | label value NBI1_POBRE NBI1_POBRE 46 | label var NBI1_POBRE "Con al menos una NBI" 47 | tab NBI1_POBRE 48 | 49 | gen NBI2_POBRE=. 50 | replace NBI2_POBRE=1 if nbihog>1 51 | replace NBI2_POBRE=0 if nbihog<2 52 | 53 | label define NBI2_POBRE 0 "menos de dos NBI" 1 "al menos dos NBI" 54 | label value NBI2_POBRE NBI2_POBRE 55 | label var NBI2_POBRE "Con al menos dos NBI" 56 | tab NBI2_POBRE 57 | 58 | 59 | *VARIABLES GEOGRAFICAS (area, regnat, dpto) 60 | gen area=estrato 61 | recode area (1/5=1) (6/8=2) 62 | label define area 1 "urbano" 2 "rural" 63 | label values area area 64 | 65 | gen regnat=1 if dominio<=3 | dominio==8 66 | replace regnat=2 if dominio>=4 & dominio<=6 67 | replace regnat=3 if dominio==7 68 | 69 | lab var regnat "Region natural" 70 | lab def regnat 1 "Costa" 2 "Sierra" 3 "Selva" 71 | lab val regnat regnat 72 | 73 | destring ubigeo, generate(dpto) 74 | replace dpto=dpto/10000 75 | replace dpto=round(dpto) 76 | label variable dpto "Departamento" 77 | label define dpto 1 "Amazonas" 78 | label define dpto 2 "Ancash", add 79 | label define dpto 3 "Apurimac", add 80 | label define dpto 4 "Arequipa", add 81 | label define dpto 5 "Ayacucho", add 82 | label define dpto 6 "Cajamarca", add 83 | label define dpto 7 "Callao", add 84 | label define dpto 8 "Cusco", add 85 | label define dpto 9 "Huancavelica", add 86 | label define dpto 10 "Huanuco", add 87 | label define dpto 11 "Ica", add 88 | label define dpto 12 "Junin", add 89 | label define dpto 13 "La_Libertad", add 90 | label define dpto 14 "Lambayeque", add 91 | label define dpto 15 "Lima", add 92 | label define dpto 16 "Loreto", add 93 | label define dpto 17 "Madre_de_Dios", add 94 | label define dpto 18 "Moquegua", add 95 | label define dpto 19 "Pasco", add 96 | label define dpto 20 "Piura", add 97 | label define dpto 21 "Puno", add 98 | label define dpto 22 "San_Martin", add 99 | label define dpto 23 "Tacna", add 100 | label define dpto 24 "Tumbes", add 101 | label define dpto 25 "Ucayali", add 102 | label values dpto dpto 103 | 104 | 105 | *Cambiamos el nombre de la variable ahno y le damos 106 | *nombre a los nbi para eliminar los caracteres que no reconoce STATA 107 | rename a*o anio 108 | label var nbi1 "Poblacion en viviendas con caracteristicas fisicas inadecuadas" 109 | label var nbi2 "Poblacion en viviendas con hacinamiento" 110 | label var nbi3 "Poblacion en viviendas sin desague de ningun tipo" 111 | label var nbi4 "Poblacion en hogares con ninos (6 a 12) que no asisten a la escuela" 112 | label var nbi5 "Poblacion en hogares con alta dependencia economica" 113 | label var NBI1_POBRE "Con al menos una NBI" 114 | label var NBI2_POBRE "De 2 a 5 NBI" 115 | 116 | svy: mean nbi1 nbi2 nbi3 nbi4 nbi5 117 | outreg using nbi.doc, replace stat(b se ci) nosubstat bdec(4) varlabels /// 118 | ctitle("","Pobreza","Error Estandar","Intervalo de confianza al 95%") /// 119 | title("Indicadores de NBI - 2017") /// 120 | note(""\"Ejecutado el $S_TIME, $S_DATE"\"Fuente: ENAHO 2017") 121 | 122 | svy: mean NBI1_POBRE 123 | outreg using nbi.doc, addtable stat(b se ci) nosubstat bdec(4) varlabels /// 124 | ctitle("","Pobreza","Error Estandar","Intervalo de confianza al 95%") /// 125 | title("Porcentaje de la poblacion con al menos 1 NBI (a) - 2017") 126 | 127 | svy: mean NBI1_POBRE , over(area) 128 | outreg using nbi.doc, addtable stat(b se ci) nosubstat bdec(4) varlabels /// 129 | ctitle("Ambito","Pobreza","Error Estandar","Intervalo de confianza al 95%") /// 130 | title("Porcentaje de la poblacion con al menos 1 NBI (b) - 2017") 131 | 132 | svy: mean NBI1_POBRE , over(regnat) 133 | outreg using nbi.doc, addtable stat(b se ci) nosubstat bdec(4) varlabels /// 134 | ctitle("Ambito","Pobreza","Error Estandar","Intervalo de confianza al 95%") /// 135 | title("Porcentaje de la poblacion con al menos 1 NBI (c) - 2017") 136 | -------------------------------------------------------------------------------- /#8 #9 #10 Medición de Indicadores Sociales con la ENAHO en STATA Indicadores de Educación.do: -------------------------------------------------------------------------------- 1 | cd "D:\ENAHO 2017" 2 | 3 | *Bajamos el archivos zipeado la pagina del INEI - ENAHO 2017 4 | copy "http://iinei.inei.gob.pe/iinei/srienaho/descarga/STATA/603-Modulo03.zip" 603-Modulo03.zip, replace 5 | 6 | *Tambien se puede bajar manualmente los archivos de la pagina "microdatos" del INEI 7 | 8 | *Abrir archivo 9 | use enaho01a-2017-300.dta, clear 10 | gen area=estrato 11 | recode area (1/5=1) (6/8=2) 12 | lab def area 1 "Urbana" 2 "Rural", modify 13 | lab val area area 14 | lab var area "Area de residencia" 15 | 16 | gen regnat=dominio 17 | recode regnat (1/3 8 =1) (4/6=2) (7=3) 18 | label var regnat "region natural" 19 | label define regnat 1 "costa" 2 "sierra" 3 "selva" 20 | label values regnat regnat 21 | ************************************************************ 22 | 23 | *MATRICULA 24 | *Definimos la poblacion segun edad normativa: 25 | *de 6 a 11 anios (Primaria) & de 12 a 17 anios (Secundaria) 26 | gen edad_norma=1 if p208a>=6 & p208a<=11 27 | replace edad_norma=2 if p208a>=12 & p208a<=17 28 | label define edad_norma 1 "Primaria" 2 "Secundaria" 29 | label val edad_norma edad_norma 30 | label var edad_norma "Edad Normativa" 31 | 32 | gen matri_esc=1 if p208a>=6 & p208a<=17 & p306==1 33 | replace matri_esc=2 if p208a>=6 & p208a<=17 & p306==2 34 | label define matri_esc 1 "Matriculado" 2 "No matriculado" 35 | label values matri_esc matri_esc 36 | label var matri_esc "Matricula Escolar" 37 | *Usamos el factor de expansion de la base para obtener los resultados a nivel nacional 38 | tab matri_esc edad_norma [iweight=factora07] , col nofreq 39 | 40 | *ASISTENCIA ESCOLAR 41 | tab p307 area [iweight=factora07] if p208a>=6 & p208a<=17, col nofreq 42 | tab p307 edad_norma [iweight=factora07] 43 | tab p307 edad_norma [iweight=factora07] , col nofreq 44 | 45 | /*En este caso no se incluye a los estudiantes que cursan primaria 46 | y tienen 12 y más años. Tampoco se incluye a los estudiantes que cursan 47 | secundaria y tienen 18 y más años. Dado que el atraso escolar es medido 48 | de más de dos años, el atraso escolar empieza desde los 9 años*/ 49 | 50 | *p308c, Grado de estudios al que asiste (primaria) 51 | tab p308c p208a [iweight=factora07] if p208a>=6 & p208a<=11 52 | 53 | *p308b, Ahno de estudios al que asiste (secundaria) 54 | tab p308b p208a [iweight=factora07] if p308a==3 &p208a>=12 & p208a<=17 55 | 56 | *ANALFABETISMO 57 | gen analfa=0 if p208a>=15 & p204==1 58 | replace analfa=1 if p208a>=15 & p204==1 & p302==2 59 | tab analfa [iweight=factora07] 60 | 61 | *NIVEL DE ESCOLARIDAD 62 | tab p301a [iweight=factora07] if p208a>=18 63 | 64 | *************************************************************** 65 | *************************************************************** 66 | *PASANDO LOS RESULTADOS A UN DOC DE WORD 67 | 68 | *Establecemos las caracteristicas de la encuesta usando las variable 69 | *factor de expansion del modulo 300 (factora07), conglomerado y estrato 70 | svyset [pweight=factora07], psu(conglome) strata(estrato) 71 | 72 | *Empezamos a registrar los resultados en word 73 | *(recordar precisar el directorio en el que estan trabajando) 74 | 75 | asdoc, text(\b INDICADORES DE EDUCACION) fs(12) replace 76 | 77 | asdoc, text(\b Matricula Escolar) fs(12) append 78 | 79 | asdoc tab matri_esc edad_norma [iweight=factora07] , /// 80 | title(\b Cantidad de ninos y ninas entre 6 y 17 ahnos matriculados por nivel educativo) fs(10) append 81 | 82 | asdoc tab matri_esc edad_norma [iweight=factora07] , /// 83 | title(\b Cantidad de ninos y ninas entre 6 y 17 ahnos matriculados por nivel educativo) col nofreq fs(10) append 84 | 85 | *ASISTENCIA ESCOLAR 86 | asdoc, text(\b Asistencia Escolar) fs(12) append 87 | 88 | asdoc tab p307 area [iweight=factora07] if p208a>=6 & p208a<=17, /// 89 | title(\b Nivel de Inasistencia Escolar por Area) fs(10) append 90 | 91 | asdoc tab p307 edad_norma [iweight=factora07] , /// 92 | title(\b Nivel de Inasistencia Escolar por Edad Normativa) fs(10) append 93 | 94 | asdoc tab p307 edad_norma [iweight=factora07] , /// 95 | title(\b Nivel de Inasistencia Escolar por Edad Normativa) col nofreq fs(10) append 96 | 97 | *RETRASO ESCOLAR (PORCENTAJE DE NINOS CON MAS DE 2 AHNOS DE RETRASO ESCOLAR) 98 | asdoc, text(\b Retraso escolar) fs(12) append 99 | 100 | *p308c, Grado de estudios al que asiste (primaria) 101 | asdoc tab p308c p208a [iweight=factora07] if p208a>=6 & p208a<=11, /// 102 | title(\b Numero de estudiantes por edad y ahno de estudio que cursa - nivel primaria) fs(10) append 103 | 104 | *p308b, Ahno de estudios al que asiste (secundaria) 105 | asdoc tab p308b p208a [iweight=factora07] if p308a==3 &p208a>=12 & p208a<=17, /// 106 | title(\b Numero de estudiantes por edad y ahno de estudio que cursa - nivel secundaria) fs(10) append 107 | 108 | *ANALFABETISMO 109 | asdoc, text(\b Analfabetismo) fs(12) append 110 | 111 | asdoc tab analfa [iweight=factora07], /// 112 | title(\b Tasa de analfabetismo, personas mayores de 15 ahnos) fs(10) append 113 | 114 | *NIVEL DE ESCOLARIDAD 115 | asdoc, text(\b Nivel de Escolaridad) fs(12) append 116 | 117 | asdoc tab p301a [iweight=factora07] if p208a>=18, /// 118 | title(\b Grado de Instruccion de personas de 18 ahnos a mas) fs(10) dec(1) append 119 | asdoc tab p301a p207 [iweight=factora07] if p208a>=18, /// 120 | title(\b Grado de Instruccion de personas de 18 ahnos a mas por sexo) fs(10) dec(1) append 121 | -------------------------------------------------------------------------------- /#56 MINEDU-ESCALE Tasa de conclusión, ed. superior, edades 22-24 (% del total).do: -------------------------------------------------------------------------------- 1 | 2 | Estimad@s, en este video les muestro como estimar la Tasa de conclusión, ed. superior, edades 22-24 (% del total) con la Encuesta Nacional de Hogares (ENAHO) del 2019 en STATA. Gracias a Gustavo Espinoza Peralta por compartir el do-file. 3 | 4 | *************************************************** 5 | 6 | *Tasa de conclusión 7 | clear all 8 | 9 | cd "D:\CONCLUSION" 10 | 11 | *Unir las bases del módulo salud, educación y sumaria 12 | use "enaho01a-2019-400.dta",clear 13 | merge 1:1 conglome vivienda hogar codperso using "enaho01a-2019-300.dta" 14 | keep if _merge==3 15 | drop _m 16 | merge m:1 conglome vivienda hogar using "sumaria-2019.dta", nogenerate 17 | 18 | *Estimar edad en años cumplidos al 31 de Marzo 19 | rename a*o year 20 | destring year, replace 21 | gen edad_31Marzo=year-p400a3 if p400a2<4 22 | replace edad_31Marzo=(year-p400a3-1) if p400a2>3 & p400a2<=12 23 | replace edad_31Marzo=0 if edad_31Marzo<0 24 | 25 | *Variables geográficas (departamento, urbano/rural) 26 | *Departamento (distinguir Lima Metropolitana de Lima Provincias) 27 | destring ubigeo, generate(dpto) 28 | replace dpto=dpto/10000 29 | replace dpto=round(dpto) 30 | replace dpto=26 if (dominio==8) //Lima Metropolitana 31 | replace dpto=27 if (dpto==15 & dominio!=8) //Lima Provincias 32 | label define dpto_num 1 "Amazonas" 33 | label define dpto_num 2 "Ancash", add 34 | label define dpto_num 3 "Apurimac", add 35 | label define dpto_num 4 "Arequipa", add 36 | label define dpto_num 5 "Ayacucho", add 37 | label define dpto_num 6 "Cajamarca", add 38 | label define dpto_num 7 "Callao", add 39 | label define dpto_num 8 "Cusco", add 40 | label define dpto_num 9 "Huancavelica", add 41 | label define dpto_num 10 "Huanuco", add 42 | label define dpto_num 11 "Ica", add 43 | label define dpto_num 12 "Junin", add 44 | label define dpto_num 13 "La_Libertad", add 45 | label define dpto_num 14 "Lambayeque", add 46 | *label define dpto_num 15 "Lima", add 47 | label define dpto_num 16 "Loreto", add 48 | label define dpto_num 17 "Madre_de_Dios",add 49 | label define dpto_num 18 "Moquegua", add 50 | label define dpto_num 19 "Pasco", add 51 | label define dpto_num 20 "Piura", add 52 | label define dpto_num 21 "Puno", add 53 | label define dpto_num 22 "San_Martin", add 54 | label define dpto_num 23 "Tacna", add 55 | label define dpto_num 24 "Tumbes", add 56 | label define dpto_num 25 "Ucayali", add 57 | label define dpto_num 26 "Lima_Metropolitana", add 58 | label define dpto_num 27 "Lima_Provincias", add 59 | label values dpto dpto_num 60 | 61 | *Area 62 | recode estrato (1/5=1 "Urbana")(6/8=2 "Rural"), gen(area) 63 | lab var area "Area de Residencia" 64 | 65 | *Etiquetar la variable pobreza monetaria 66 | *Usar variable "pobreza" de sumaria 67 | label define pobre_3 1 "pobre_extremo" 2 "pobre_no_extremo" 3 "no_pobre" 68 | label value pobreza pobre_3 69 | 70 | *Crear la variable lengua materna 71 | *Lengua indígena: quechua, aymara, otra lengua nativa 72 | *Lengua Castellano 73 | recode p300a (4=1 "Castellano")(1/3=2 "Indigena")(5/9=.), gen(lengua) 74 | 75 | **************************************************** 76 | *Establecer las caracteristicas de la encuesta usando las variables 77 | *factor de expansion del modulo 300 (factora07), conglomerado y estrato 78 | svyset [pweight=factora07], psu(conglome) strata(estrato) singleunit(centered) 79 | 80 | *Tasa de conclusión, educación superior, grupo de edades 22-24 (% del total) 81 | gen conclusnyu2224=0 if edad_31Marzo>=22 & edad_31Marzo<=24 82 | replace conclusnyu2224=1 if (edad_31Marzo>=22 & edad_31Marzo<=24) & (p301a==8 | p301a==10 | p301a==11) 83 | label def con_l 0 "Incompleta" 1 "Completa" 84 | label val conclusnyu2224 con_l 85 | label var conclusnyu2224 "Tasa de conclusion, educacion superior, edades 22-24" 86 | 87 | *Estimar resultados 88 | svy: prop conclusnyu2224 89 | outreg using conclusion, stat(b se ci) bdec(4) nosubstat varlabels replace /// 90 | title("Tasa de conclusion - educacion superior - edades 22-24") /// 91 | ctitle("","","Pro","SE","Intervalo de Confianza") 92 | 93 | svy: prop conclusnyu2224 , over(p207) 94 | outreg using conclusion, stat(b se ci) bdec(4) nosubstat varlabels addtable /// 95 | title("Tasa de conclusion - educacion superior - edades 22-24, por sexo") /// 96 | ctitle("Ed. Superior","Sexo","Pro","SE","Intervalo de Confianza") 97 | 98 | svy: prop conclusnyu2224 , over(area) 99 | outreg using conclusion, stat(b se ci) bdec(4) nosubstat varlabels addtable /// 100 | title("Tasa de conclusion - educacion superior - edades 22-24, por area") /// 101 | ctitle("Ed. Superior","Area","Pro","SE","Intervalo de Confianza") 102 | 103 | svy: prop conclusnyu2224 , over(lengua) 104 | outreg using conclusion, stat(b se ci) bdec(4) nosubstat varlabels addtable /// 105 | title("Tasa de conclusion - educacion superior - edades 22-24, por lengua materna") /// 106 | ctitle("Ed. Superior","Lengua","Pro","SE","Intervalo de Confianza") 107 | 108 | svy: prop conclusnyu2224 , over(pobreza) 109 | outreg using conclusion, stat(b se ci) bdec(4) nosubstat varlabels addtable /// 110 | title("Tasa de conclusion - educacion superior - edades 22-24, por pobreza") /// 111 | ctitle("Ed. Superior","Pobreza","Pro","SE","Intervalo de Confianza") 112 | 113 | svy: prop conclusnyu2224 , over(dpto) 114 | outreg using conclusion, stat(b se ci) bdec(4) nosubstat varlabels addtable /// 115 | title("Tasa de conclusion - educacion superior - edades 22-24, por dpto") /// 116 | ctitle("Ed. Superior","Dpto","Pro","SE","Intervalo de Confianza") 117 | -------------------------------------------------------------------------------- /#95 Población económicamente activa ocupada por rama de actividad con la ENAHO en STATA.do: -------------------------------------------------------------------------------- 1 | clear all 2 | cd "D:\ENAHO" 3 | 4 | /* 5 | use "enaho01a-2022-500.dta", clear 6 | append using "enaho01a-2021-500.dta" 7 | append using "enaho01a-2020-500.dta" 8 | append using "enaho01a-2019-500.dta" 9 | append using "enaho01a-2018-500.dta" 10 | append using "enaho01a-2017-500.dta" 11 | save "enaho500-2017-2022.dta", replace 12 | */ 13 | 14 | use "enaho500-2017-2022.dta", clear 15 | rename a*o year 16 | gen anio=real(year) 17 | label var year "Year" 18 | 19 | tab p500i year 20 | drop if p500i=="00" 21 | 22 | *Se establece quienes son residentes habituales 23 | gen resi=1 if ((p204==1 & p205==2) | (p204==2 & p206==1)) 24 | *P204 ¿Es miembro del hogar? 204. ¿Es miembro del hogar familiar? 1. Sí 2. No 25 | *P205 ¿Se encuentra ausente del hogar 30 días o más? 1. Sí 2. No 26 | *P206 ¿Está presente en el hogar 30 días o más? 1. Sí 2. No 27 | 28 | ******************************************************************************** 29 | *PEA ocupada por Rama de Actividad 30 | *CIIU 31 | br p506r4 32 | 33 | gen ciiu_aux1 =substr("0"+string(p506r4),1,.) 34 | 35 | br p506r4 ciiu_aux1 36 | replace ciiu_aux1 =substr(string(p506r4),1,.) if p506r4>999 37 | br p506r4 ciiu_aux1 38 | 39 | gen ciiu_aux2 =substr(ciiu_aux1 ,1,2) 40 | br p506r4 ciiu_aux1 41 | 42 | destring ciiu_aux2, generate(ciiu_2d) 43 | br p506r4 ciiu_aux1 ciiu_2d 44 | 45 | 46 | gen ciiu_1d=1 if ciiu_2d<=2 47 | replace ciiu_1d=2 if ciiu_2d==3 48 | replace ciiu_1d=3 if ciiu_2d>=5 & ciiu_2d<=9 49 | replace ciiu_1d=4 if ciiu_2d>=10 & ciiu_2d<=33 50 | replace ciiu_1d=5 if ciiu_2d>=41 & ciiu_2d<=43 51 | replace ciiu_1d=6 if ciiu_2d>=45 & ciiu_2d<=47 52 | replace ciiu_1d=7 if (ciiu_2d>=49 & ciiu_2d<=53) | (ciiu_2d>=58 & ciiu_2d<=63) 53 | replace ciiu_1d=8 if ciiu_2d==84 54 | replace ciiu_1d=9 if ciiu_2d>=55 & ciiu_2d<=56 55 | replace ciiu_1d=10 if ciiu_2d==68 | (ciiu_2d>=69 & ciiu_2d<=82) 56 | replace ciiu_1d=11 if ciiu_2d==85 57 | replace ciiu_1d=12 if (ciiu_2d>=35 & ciiu_2d<=39) | (ciiu_2d>=64 & ciiu_2d<=66) | /// 58 | (ciiu_2d>=86 & ciiu_2d<=88) | (ciiu_2d>=90 & ciiu_2d<=93)| (ciiu_2d>=94 & ciiu_2d<=98) | ciiu_2d==99 59 | label var ciiu_1d "Division CIIU" 60 | la de ciiu_1d 1 "Agricultura" 2 "Pesca" 3 "Mineria" 4 "Manufactura" 5 "Construccion" /// 61 | 6 "Comercio" 7 "Transportes y Comunicaciones" 8 "Gobierno" 9 "Hoteles y Restaurantes" /// 62 | 10 "Inmobiliarias y alquileres" 11 "Ensehnanza" 12 "Otros Servicios 1/" 63 | label values ciiu_1d ciiu_1d 64 | br p506r4 ciiu_aux1 ciiu_2d ciiu_1d 65 | 66 | 67 | gen ciiu_5=1 if ciiu_1d<4 68 | replace ciiu_5=2 if ciiu_1d==4 69 | replace ciiu_5=3 if ciiu_1d==5 70 | replace ciiu_5=4 if ciiu_1d==6 71 | replace ciiu_5=5 if ciiu_1d>6 72 | label var ciiu_5 "Division CIIU-5 categorias" 73 | la de ciiu_5 1 "Extractiva 1/" 2 "Industria" 3 "Construccion" /// 74 | 4 "Comercio" 5 "Servicios 2/" 75 | label values ciiu_5 ciiu_5 76 | br p506r4 ciiu_aux1 ciiu_2d ciiu_1d ciiu_5 77 | 78 | 79 | table ciiu_5 year if resi==1 & ocu500==1 [iw= fac500a], nformat(%12.0fc) 80 | 81 | 82 | 83 | 84 | 85 | 86 | 87 | 88 | 89 | gen extr=ciiu_5==1 90 | gen indu=ciiu_5==2 91 | gen cons=ciiu_5==3 92 | gen come=ciiu_5==4 93 | gen serv=ciiu_5==5 94 | gen tot =ciiu_5<6 95 | 96 | bys year: egen float t_extr = total(fac500a * extr) if resi==1 & ocu500==1 97 | bys year: egen float t_indu = total(fac500a * indu) if resi==1 & ocu500==1 98 | bys year: egen long t_cons = total(fac500a * cons) if resi==1 & ocu500==1 99 | bys year: egen float t_come = total(fac500a * come) if resi==1 & ocu500==1 100 | bys year: egen float t_serv = total(fac500a * serv) if resi==1 & ocu500==1 101 | bys year: egen float t_tot = total(fac500a * tot) if resi==1 & ocu500==1 102 | label var t_extr "Extractiva 1/" 103 | label var t_indu "Industria" 104 | label var t_cons "Construccion" 105 | label var t_come "Comercio" 106 | label var t_serv "Servicios 2/" 107 | 108 | gen p_extr=t_extr*100/t_tot 109 | gen p_indu=t_indu*100/t_tot 110 | gen p_cons=t_cons*100/t_tot 111 | gen p_come=t_come*100/t_tot 112 | gen p_serv=t_serv*100/t_tot 113 | 114 | format t_* %20.0fc 115 | format p_extr p_indu p_cons p_come p_serv %4.1f 116 | 117 | keep t_* p_* anio 118 | duplicates drop 119 | twoway /// 120 | (scatter p_extr p_indu p_cons p_come p_serv anio, /// 121 | mcolor(none none none none none) /// 122 | mlabel(p_extr p_indu p_cons p_come p_serv) /// 123 | mlabsize(vsmall vsmall vsmall vsmall vsmall) /// 124 | mlabcolor(orange gold gray navy red) /// 125 | mlabposition(3 12 5 5 10) /// 126 | yaxis(2) ylabel(,nogrid) yscale(off) yscale(lpattern(blank) ) /// 127 | xtitle("") xlabel(, nogrid) ) /// 128 | (scatter t_extr t_indu t_cons t_come t_serv anio, /// 129 | connect(l l l l l) lcolor(orange gold gray navy red) /// 130 | msymbol(Dh Dh Dh Dh Dh) mcolor(orange gold gray navy red) /// 131 | mlabel(t_extr t_indu t_cons t_come t_serv) /// 132 | mlabsize(vsmall vsmall vsmall vsmall vsmall) /// 133 | mlabcolor(orange gold gray navy red) /// 134 | mlabposition(6 11 7 8 8) /// 135 | yaxis(1) ylabel(, nogrid) yscale(off) yscale(lpattern(blank)) /// 136 | xtitle("") xlabel(, nogrid) ), /// 137 | legend(on) legend(order(6 7 8 9 10) size(vsmall) cols(5) position(6)) /// 138 | yscale(off) yscale(lpattern(blank)) /// 139 | title("GRÁFICO N° 2.12" "PERÚ: PEA OCUPADA, SEGÚN RAMA DE ACTIVIDAD ECONÓMICA, 2017-2022", size(small) color(red) margin(vsmall)) /// 140 | subtitle("(Absoluto y porcentaje)", size(vsmall) margin(small)) /// 141 | note("Nota: Clasificación de ramas de actividad basada en el CIIU Rev. 4." "Fuente: INEI-ENAHO, 2017-2022.", size(vsmall)) 142 | -------------------------------------------------------------------------------- /#28 Construcción de Indicadores de Empleo y Desempleo.do: -------------------------------------------------------------------------------- 1 | cd "D:\ENAHO" 2 | copy "http://iinei.inei.gob.pe/iinei/srienaho/descarga/STATA/634-Modulo05.zip" 634-Modulo05.zip, replace 3 | copy "http://iinei.inei.gob.pe/iinei/srienaho/descarga/STATA/634-Modulo04.zip" 634-Modulo04.zip, replace 4 | copy "http://iinei.inei.gob.pe/iinei/srienaho/descarga/STATA/634-Modulo34.zip" 634-Modulo34.zip, replace 5 | 6 | *Bajar el modulo de empleo/ingreso de la ENAHO y descomprimirlo (manualmente) 7 | use enaho01a-2018-500.dta, clear 8 | 9 | /* ocu500 indicador de la PEA: 1. ocupado 2. desocupado abierto 3. desocupado oculto 4. no pea 10 | fac500 factor de expansion de empleo/ingresos 11 | p500i codigo informante del capitulo 500 12 | p204 es miembro del hogar 13 | p205 se encuentra ausente del hogar 30 dias o mas 14 | p206 esta presente en el hogar 30 dias o mas 15 | p207 sexo 16 | p208a edad 17 | p209 estado civil 18 | p301a ultimo grado de estudios 19 | p558c. Por sus antepasados y de acuerdo a sus costumbres, ¿Ud. se considera: 20 | 1. Quechua? 2. Aimara? 3. Nativo o Indígena de la Amazonía? 21 | 4. Negro/ Mulato/Zambo/Afro peruano? 5. Blanco? 6. Mestizo? 7. Otro? 8. No Sabe? 22 | 9. Perteneciente o parte de otro pueblo indígena u originario */ 23 | 24 | rename a*o year 25 | drop if p500i=="00" 26 | 27 | *Area de residencia 28 | gen area=1 if estrato<=5 29 | replace area=2 if estrato>=6 & estrato<=8 30 | lab def area 1 "urbano" 2 "rural" 31 | lab val area area 32 | 33 | *Region natural 34 | gen region=1 if dominio>=1 & dominio<=3 35 | replace region=1 if dominio==8 36 | replace region=2 if dominio>=4 & dominio<=6 37 | replace region=3 if dominio==7 38 | label define region 1 "Costa" 2 "Sierra" 3 "Selva" 39 | lab val region region 40 | 41 | *Dominio geografico 42 | gen domin02=1 if dominio>=1 & dominio<=3 & area==1 43 | replace domin02=2 if dominio>=1 & dominio<=3 & area==2 44 | replace domin02=3 if dominio>=4 & dominio<=6 & area==1 45 | replace domin02=4 if dominio>=4 & dominio<=6 & area==2 46 | replace domin02=5 if dominio==7 & area==1 47 | replace domin02=6 if dominio==7 & area==2 48 | replace domin02=7 if dominio==8 49 | label define domin02 1 "Costa_urbana" 2 "Costa_rural" 3 "Sierra_urbana" /// 50 | 4 "Sierra_rural" 5 "Selva_urbana" 6 "Selva_rural" 7 "Lima_Metropolitana" 51 | label value domin02 domin02 52 | 53 | *Departamento 54 | gen dpto= real(substr(ubigeo,1,2)) 55 | replace dpto=15 if (dpto==7) 56 | label define dpto 1"Amazonas" 2"Ancash" 3"Apurimac" 4"Arequipa" 5"Ayacucho" /// 57 | 6"Cajamarca" 8"Cusco" 9"Huancavelica" 10"Huanuco" 11"Ica" 12"Junin" /// 58 | 13"La_Libertad" 14"Lambayeque" 15"Lima" 16"Loreto" 17"Madre_de_Dios" 18"Moquegua" /// 59 | 19"Pasco" 20"Piura" 21"Puno" 22"San_Martin" 23"Tacna" 24"Tumbes" 25"Ucayali" 60 | lab val dpto dpto 61 | 62 | *Se establece quienes son residentes habituales 63 | gen resi=1 if ((p204==1 & p205==2) | (p204==2 & p206==1)) 64 | 65 | *Grupos de edad 66 | gen edad=1 if p208a>13 & p208a<25 67 | replace edad=2 if p208a>24 & p208a<60 68 | replace edad=3 if p208a>59 & p208a<65 69 | replace edad=4 if p208a>64 70 | label define edad 1 "14-24" 2 "25-59" 3 "60-64" 4 "65 y mas" 71 | lab val edad edad 72 | 73 | *Estado civil 74 | gen ecivil=1 if p209==3 | p209==4 | p209==5 75 | replace ecivil=2 if p209==6 76 | replace ecivil=3 if p209==1 | p209==2 77 | lab def ecivil 1 "Alguna vez unido/a 2/" 2 " Soltero/a" 3 "Unido/a 1/" 78 | lab val ecivil ecivil 79 | *1/ Incluye : Conviviente y casado/a 80 | *2/ Incluye: Separado/a, divorciado/a y viudo/a. 81 | 82 | *Educacion 83 | gen educ=1 if p301a<5 | p301a==12 84 | replace educ=2 if p301a==5 | p301a==6 85 | replace educ=3 if p301a==7 | p301a==8 86 | replace educ=4 if p301a==9 | p301a==10 | p301a==11 87 | lab def educ 1 "Primaria 1/" 2 "Educacion secundaria" 3 "Superior no universitaria" /// 88 | 4 "Superior universitaria" 89 | lab val educ educ 90 | *1/ Incluye sin nivel e inicial. A partir del año 2017 se incluye educación básica especial 91 | 92 | *Etnia 93 | gen etnia=1 if p558c<4 | p558c==9 94 | replace etnia=2 if p558c==4 95 | replace etnia=3 if p558c==6 96 | replace etnia=4 if p558c==5 | p558c==7 97 | replace etnia=5 if p558c==8 98 | lab def etnia 1 "Indigena 1/" 2 "Negro, mulato, Afro peruano" 3 "Mestizo/a" 4 "Otro 2/" 5 "No sabe" 99 | lab val etnia etnia 100 | *1/ Incluye: Quechua, Aimara y Nativo o Indígena de la Amazonía. 101 | *2/ Incluye: Blanco y otro 102 | 103 | *Población en Edad de Trabajar por sexo, según ámbito geográfico 104 | tab area p207 [iw= fac500a] if resi==1 105 | tab region p207 [iw= fac500a] if resi==1 106 | 107 | *Población en Edad de Trabajar, según área de residencia y grupos de edad 108 | tab edad area [iw= fac500a] if resi==1 109 | 110 | *Poblacion Economicamente Activa, según ámbito geográfico 111 | tab area ocu500 [iw= fac500a ] if ocu500<3 & resi==1 112 | tab region ocu500 [iw= fac500a ] if ocu500<3 & resi==1 113 | 114 | *Población Económicamente Activa según etnia 115 | tab etnia ocu500 [iw= fac500a ] if ocu500<3 & resi==1, nofreq col 116 | 117 | *Población Económicamente Activa, según estado civil o conyugal y sexo 118 | tab ecivil p207 [iw= fac500a ] if ocu500<3 & resi==1, nofreq col 119 | 120 | *Tasa de actividad 121 | *El cociente de la Población Económicamente Activa entre el total de Población en Edad de Trabajar 122 | gen t_act=0 if p208a>=14 123 | replace t_act=1 if ocu500==1 | ocu500==2 124 | lab def t_act 0 "" 1 "Tasa de Actividad" 125 | lab val t_act t_act 126 | 127 | tab t_act area [iw= fac500a] if resi==1, col nofreq 128 | tab t_act region [iw= fac500a] if resi==1, col nofreq 129 | 130 | *Componentes de la Población Económicamente Inactiva, según sexo 131 | gen c_pei=1 if p546==6 132 | replace c_pei=2 if p546==3 | p546==8 133 | replace c_pei=3 if p546==7 134 | replace c_pei=4 if p546==4 135 | replace c_pei=5 if p546==5 136 | lab def c_pei 1 "Vivía de su pensión o jubilación u otras rentas" 2 "Otro 1/" /// 137 | 3 "Enfermo o incapacitado" 4 "Estudiando" 5 "Quehaceres del hogar" 138 | lab val c_pei c_pei 139 | *1/ Incluye : Esperando el inicio de un trabajo dependiente, otro y no especificado 140 | 141 | tab c_pei p207 [iw= fac500a ] if ocu500>2 & resi==1, nofreq col 142 | 143 | *Población ocupada según nivel de educación alcanzado 144 | tab educ [iw= fac500a ] if ocu500==1 & resi==1 145 | 146 | *Tasa de desempleo abierto 147 | *Es la proporción de la fuerza de trabajo desocupada disponible y que busca activamente trabajo 148 | gen t_de_a=0 if ocu500==1 | ocu500==2 149 | replace t_de_a=1 if ocu500==2 150 | lab def t_de_a 0 "" 1 "Tasa de Desempleo Abierto" 151 | lab val t_de_a t_de_a 152 | 153 | *Perú urbano: Tasa de desempleo abierto por sexo, según etnia 154 | tab etnia t_de_a [iw= fac500a] if resi==1 & area==1 & p207==1, row nofreq 155 | tab etnia t_de_a [iw= fac500a] if resi==1 & area==1 & p207==2, row nofreq 156 | -------------------------------------------------------------------------------- /#109 #110 Gasto de bolsillo.do: -------------------------------------------------------------------------------- 1 | 2 | *Salud 3 | use "D:\ENAHO\enaho01a-2019-400.dta", clear 4 | 5 | /* 6 | P4025 En las últ 4 semanas, ¿Presentó Ud. Algún(a)? -No tuvo enfermedad, síntoma, recaída, accidente 7 | 8 | P414$01 En las últ 4 semanas, ¿recibió: Consulta? 9 | P414$02 En las últ 4 semanas, ¿recibió: Medicinas? 10 | P414$03 En las últ 4 semanas, ¿recibió: Análisis? 11 | P414$04 En las últ 4 semanas, ¿recibió: Rayos X, Tomografía, etc? 12 | P414$05 En las últ 4 semanas, ¿recibió: Otros exámenes (Hemodiálisis, etc.)? 13 | P414$06 En los últ 3 meses, ¿recibió: Servicio dental y conexos? 14 | P414$07 En los últ 3 meses, ¿recibió: Servicio Oftalmológico? 15 | P414$08 En los últ 3 meses, ¿recibió: Compra de Lentes? 16 | P414$09 En los últ 3 meses, ¿recibió: Vacunas? 17 | P414$10 En los últ 3 meses, ¿recibió: Control de Salud de los Niños? 18 | P414$11 En los últ 3 meses, ¿recibió : Anticonceptivos? 19 | P414$12 En los últ 3 meses, ¿recibió: Otros Gastos (ortopedia,termómetro, etc.)? 20 | P414$13 En los últ 12 meses, ¿recibió: Hospitalización? 21 | P414$14 En los últ 12 meses, ¿recibió: Intervención Quirúrgica? 22 | P414$15 En los últ 12 meses, ¿recibió: Controles por embarazo? 23 | P414$16 En los últ 12 meses, ¿recibió: Atenciones de parto? 24 | 25 | ¿El gasto realizado fue: Pagado por algún miembro de este hogar?: 26 | P4151$01 Consulta 27 | P4151$02 Medicinas/Insumos 28 | P4151$03 Análisis 29 | P4151$04 Rayos X, Tomografía, etc 30 | P4151$05 Otros Exámenes (Hemodiálisis, etc.) 31 | P4151$06 Servicio dental y conexos? 32 | P4151$07 Servicio Oftalmológico 33 | P4151$08 Compra de Lentes 34 | P4151$09 Vacunas 35 | P4151$10 Control de Salud de los Niños 36 | P4151$11 Anticonceptivos 37 | P4151$12 Otros Gastos (ortopedia, termómetro, etc.) 38 | P4151$13 Hospitalización 39 | P4151$14 Intervención Quirúrgica? 40 | P4151$15 Controles por embarazo 41 | P4151$16 Atenciones de parto 42 | 43 | ¿Cuánto fue el monto total por la compra o servicio?: 44 | P41601 Consulta 45 | P41602 Medicinas/Insumos 46 | P41603 Análisis 47 | P41604 Rayos X, Tomografía, etc 48 | P41605 Otros Exámenes (Hemodiálisis, etc.) 49 | P41606 Servicio dental y conexos? 50 | P41607 Servicio Oftalmológico 51 | P41608 Compra de Lentes 52 | P41609 Vacunas 53 | P41610 Control de Salud de los Niños 54 | P41611 Anticonceptivos 55 | P41612 Otros Gastos (ortopedia, termómetro, etc.) 56 | P41613 Hospitalización 57 | P41614 Intervención Quirúrgica 58 | P41615 Controles por embarazo 59 | P41616 Atenciones de parto 60 | 61 | i41606 (imputado, deflactado, anualizado) Servicio dental y conexos 62 | 63 | FACTOR07 Factor de Expansión anual de Población Proyecciones CPV-2007 64 | */ 65 | 66 | 67 | gen tuvo_ps=. 68 | replace tuvo_ps=1 if p4025==0 69 | label var tuvo_ps "Algún problema de salud" 70 | 71 | gen double gto01=0 72 | gen double gto02=0 73 | gen double gto03=0 74 | gen double gto04=0 75 | gen double gto05=0 76 | gen double gto06=0 77 | gen double gto07=0 78 | gen double gto08=0 79 | gen double gto09=0 80 | gen double gto10=0 81 | gen double gto11=0 82 | gen double gto12=0 83 | gen double gto13=0 84 | gen double gto14=0 85 | gen double gto15=0 86 | gen double gto16=0 87 | replace gto01=i41601 if p4151_01==1 88 | replace gto02=i41602 if p4151_02==1 89 | replace gto03=i41603 if p4151_03==1 90 | replace gto04=i41604 if p4151_04==1 91 | replace gto05=i41605 if p4151_05==1 92 | replace gto06=i41606 if p4151_06==1 93 | replace gto07=i41607 if p4151_07==1 94 | replace gto08=i41608 if p4151_08==1 95 | replace gto09=i41609 if p4151_09==1 96 | replace gto10=i41610 if p4151_10==1 97 | replace gto11=i41611 if p4151_11==1 98 | replace gto12=i41612 if p4151_12==1 99 | replace gto13=i41613 if p4151_13==1 100 | replace gto14=i41614 if p4151_14==1 101 | replace gto15=i41615 if p4151_15==1 102 | replace gto16=i41616 if p4151_16==1 103 | egen gto_usu=rowtotal(gto01 gto02 gto03 gto04 gto05 gto06 gto07 gto08 gto09 gto10 gto11 gto12 gto13 gto14 gto15 gto16) 104 | lab var gto_usu "Gasto tot en salud x usuario" 105 | 106 | *Total 107 | egen double g_salud=total(gto_usu *factor07/1000000) 108 | lab var g_salud "Total Gasto en salud, mill" 109 | sum g_salud 110 | 111 | *Tuvo algun problema de salud 112 | egen double g_ps = total(gto_usu * factor07/1000000) if tuvo_ps==1 113 | gen p_ps = g_ps*100/g_salud 114 | label var g_ps "G. salud de quienes presentan algun p.s., mill" 115 | label var p_ps "G. salud de quienes presentan algun p.s., %" 116 | sum g_ps p_ps 117 | 118 | *Por tipo de gasto (en millones) 119 | egen double g_1=total(gto01 * factor07/1000000) 120 | egen double g_2=total(gto02 * factor07/1000000) 121 | egen double g_3=total(gto03 * factor07/1000000) 122 | egen double g_4=total(gto04 * factor07/1000000) 123 | egen double g_5=total(gto05 * factor07/1000000) 124 | egen double g_6=total(gto06 * factor07/1000000) 125 | egen double g_7=total(gto07 * factor07/1000000) 126 | egen double g_8=total(gto08 * factor07/1000000) 127 | egen double g_9=total(gto09 * factor07/1000000) 128 | egen double g_10=total(gto10 * factor07/1000000) 129 | egen double g_11=total(gto11 * factor07/1000000) 130 | egen double g_12=total(gto12 * factor07/1000000) 131 | egen double g_13=total(gto13 * factor07/1000000) 132 | egen double g_14=total(gto14 * factor07/1000000) 133 | egen double g_15=total(gto15 * factor07/1000000) 134 | egen double g_16=total(gto16 * factor07/1000000) 135 | 136 | *Por tipo de gasto (Como porcentaje del gto total) 137 | gen p_1 = g_1*100/g_salud 138 | gen p_2 = g_2*100/g_salud 139 | gen p_3 = g_3*100/g_salud 140 | gen p_4 = g_4*100/g_salud 141 | gen p_5 = g_5*100/g_salud 142 | gen p_6 = g_6*100/g_salud 143 | gen p_7 = g_7*100/g_salud 144 | gen p_8 = g_8*100/g_salud 145 | gen p_9 = g_9*100/g_salud 146 | gen p_10 = g_10*100/g_salud 147 | gen p_11 = g_11*100/g_salud 148 | gen p_12 = g_12*100/g_salud 149 | gen p_13 = g_13*100/g_salud 150 | gen p_14 = g_14*100/g_salud 151 | gen p_15 = g_15*100/g_salud 152 | gen p_16 = g_16*100/g_salud 153 | 154 | label var g_1 "Consulta, mill" 155 | label var g_2 "Medicamentos, mill" 156 | label var g_3 "Análisis, mill" 157 | label var g_4 "Rayos X, Tomografía, etc" 158 | label var g_5 "Otros Exámenes (Hemodiálisis, etc.), mill" 159 | label var g_6 "Servicio dental y conexos, mill" 160 | label var g_7 "Servicio Oftalmológico, mill" 161 | label var g_8 "Compra de Lentes, mill" 162 | label var g_9 "Vacunas, mill" 163 | label var g_10 "Control de Salud de los Niños, mill" 164 | label var g_11 "Anticonceptivos, mill" 165 | label var g_12 "Otros Gastos (ortopedia, termómetro, etc.), mill" 166 | label var g_13 "Hospitalización, mill" 167 | label var g_14 "Intervención Quirúrgica, mill" 168 | label var g_15 "Controles por embarazo, mill" 169 | label var g_16 "Atenciones de parto, mill" 170 | label var p_1 "Consulta, %" 171 | label var p_2 "Medicamentos, %" 172 | label var p_3 "Análisis, %" 173 | label var p_4 "Rayos X, Tomografía, etc" 174 | label var p_5 "Otros Exámenes (Hemodiálisis, etc.), %" 175 | label var p_6 "Servicio dental y conexos, %" 176 | label var p_7 "Servicio Oftalmológico, %" 177 | label var p_8 "Compra de Lentes, %" 178 | label var p_9 "Vacunas, %" 179 | label var p_10 "Control de Salud de los Niños, %" 180 | label var p_11 "Anticonceptivos, %" 181 | label var p_12 "Otros Gastos (ortopedia, termómetro, etc.), %" 182 | label var p_13 "Hospitalización, %" 183 | label var p_14 "Intervención Quirúrgica, %" 184 | label var p_15 "Controles por embarazo, %" 185 | label var p_16 "Atenciones de parto, %" 186 | 187 | *Pasarlo a un archivo en excel 188 | keep g_* p_* 189 | duplicates drop 190 | drop if p_ps==. 191 | format p_* %5.2fc 192 | 193 | export excel using "D:\ENAHO\Gto de bolsillo 2019.xls", replace firstrow(varlabels) 194 | -------------------------------------------------------------------------------- /#65 Evolución del gasto real per cápita promedio mensual con ENAHO 2016 – 2018.do: -------------------------------------------------------------------------------- 1 | clear all 2 | cd "E:\GTO" 3 | global gto = "E:\GTO" 4 | 5 | /* Bajar los archivos sumaria de la pagina del INEI 6 | copy "http://iinei.inei.gob.pe/iinei/srienaho/descarga/STATA/546-Modulo34.zip" 546-Modulo34.zip, replace 7 | copy "http://iinei.inei.gob.pe/iinei/srienaho/descarga/STATA/603-Modulo34.zip" 603-Modulo34.zip, replace 8 | copy "http://iinei.inei.gob.pe/iinei/srienaho/descarga/STATA/634-Modulo34.zip" 634-Modulo34.zip, replace 9 | */ 10 | 11 | *Descomprimo manualmente los archivos sumaria, "deflactores_base2018_new.dta" & "despacial_ldnew.dta" 12 | *y los coloco en la misma carpeta de trabajo 13 | 14 | set mem 500m 15 | set more off 16 | use sumaria-2016, clear 17 | append using sumaria-2017 18 | append using sumaria-2018 19 | 20 | destring conglome, replace 21 | tostring conglome, replace format(%06.0f) 22 | 23 | recode gru52hd2-gashog2d (.= 0) 24 | 25 | rename a*o anio 26 | gen aniorec=real(anio) 27 | 28 | gen dpto= real(substr(ubigeo,1,2)) 29 | replace dpto=15 if (dpto==7) 30 | label define dpto 1"Amazonas" 2"Ancash" 3"Apurimac" 4"Arequipa" 5"Ayacucho" 6"Cajamarca" 8"Cusco" 9"Huancavelica" 10"Huanuco" 11"Ica" /* 31 | */12"Junin" 13"La_Libertad" 14"Lambayeque" 15"Lima" 16"Loreto" 17"Madre_de_Dios" 18"Moquegua" 19"Pasco" 20"Piura" 21"Puno" 22"San_Martin" /* 32 | */23"Tacna" 24"Tumbes" 25"Ucayali" 33 | lab val dpto dpto 34 | 35 | sort aniorec dpto 36 | merge aniorec dpto using "deflactores_base2018_new.dta" 37 | tab _m 38 | drop if _merge==2 39 | drop _m 40 | 41 | *Genero las variables de area de residencia, region natural y dominio geografico 42 | replace estrato = 1 if dominio ==8 43 | gen area = estrato <6 44 | replace area=2 if area==0 45 | label define area 2 rural 1 urbana 46 | label val area area 47 | 48 | gen domin02=1 if dominio>=1 & dominio<=3 & area==1 49 | replace domin02=2 if dominio>=1 & dominio<=3 & area==2 50 | replace domin02=3 if dominio>=4 & dominio<=6 & area==1 51 | replace domin02=4 if dominio>=4 & dominio<=6 & area==2 52 | replace domin02=5 if dominio==7 & area==1 53 | replace domin02=6 if dominio==7 & area==2 54 | replace domin02=7 if dominio==8 55 | label define domin02 1 "Costa_urbana" 2 "Costa_rural" 3 "Sierra_urbana" 4 "Sierra_rural" /* 56 | */ 5 "Selva_urbana" 6 "Selva_rural" 7 "Lima_Metropolitana" 57 | label value domin02 domin02 58 | 59 | gen region=1 if dominio>=1 & dominio<=3 60 | replace region=1 if dominio==8 61 | replace region=2 if dominio>=4 & dominio<=6 62 | replace region=3 if dominio==7 63 | label define region 1 "Costa" 2 "Sierra" 3 "Selva" 64 | label values region region 65 | 66 | gen dominioA=1 if dominio==1 & area==1 67 | replace dominioA=2 if dominio==1 & area==2 68 | replace dominioA=3 if dominio==2 & area==1 69 | replace dominioA=4 if dominio==2 & area==2 70 | replace dominioA=5 if dominio==3 & area==1 71 | replace dominioA=6 if dominio==3 & area==2 72 | replace dominioA=7 if dominio==4 & area==1 73 | replace dominioA=8 if dominio==4 & area==2 74 | replace dominioA=9 if dominio==5 & area==1 75 | replace dominioA=10 if dominio==5 & area==2 76 | replace dominioA=11 if dominio==6 & area==1 77 | replace dominioA=12 if dominio==6 & area==2 78 | replace dominioA=13 if dominio==7 & area==1 79 | replace dominioA=14 if dominio==7 & area==2 80 | replace dominioA=15 if dominio==7 & (dpto==16 | dpto==17 | dpto==25) & area==1 81 | replace dominioA=16 if dominio==7 & (dpto==16 | dpto==17 | dpto==25) & area==2 82 | replace dominioA=17 if dominio==8 & area==1 83 | replace dominioA=17 if dominio==8 & area==2 84 | label define dominioA 1 "Costa norte urbana" 2 "Costa norte rural" 3 "Costa centro urbana" 4 "Costa centro rural" /* 85 | */ 5 "Costa sur urbana" 6 "Costa sur rural" 7 "Sierra norte urbana" 8 "Sierra norte rural" 9 "Sierra centro urbana" /* 86 | */ 10 "Sierra centro rural" 11 "Sierra sur urbana" 12 "Sierra sur rural" 13 "Selva alta urbana" 14 "Selva alta rural" /* 87 | */ 15 "Selva baja urbana" 16 "Selva baja rural" 17"Lima Metropolitana" 88 | lab val dominioA dominioA 89 | 90 | drop ld 91 | 92 | sort dominioA 93 | 94 | merge dominioA using "despacial_ldnew.dta" 95 | tab _m 96 | drop _m 97 | 98 | gen factornd07=round(factor07*mieperho,1) 99 | 100 | 101 | *GASTOS REALES 102 | *******************************************************************************/ 103 | *CREANDO VARIABLES DEL GASTO DEFLACTADO A PRECIOS DE LIMA Y BASE 2018 a nivel total**/ 104 | *******************************************************************************/ 105 | *Gasto por 8 grupos de la canastas* 106 | gen gpcrg3= (gru11hd + gru12hd1 + gru12hd2 + gru13hd1 + gru13hd2 + gru13hd3 )/(12*mieperho*ld*i01) 107 | lab var gpcrg3 "Preparados dentro del hogar" 108 | 109 | gen gpcrg6 = ((g05hd + g05hd1 + g05hd2 + g05hd3 + g05hd4 + g05hd5 +g05hd6 +ig06hd)/(12*mieperho*ld*i01)) 110 | lab var gpcrg6 "Adquiridos Fuera del hogar 559" 111 | 112 | gen gpcrg8= ((sg23 + sig24)/(12*mieperho*ld*i01)) 113 | lab var gpcrg8 "Adquiridos de instituciones beneficas 602a" 114 | 115 | gen gpcrg9= ((gru14hd + gru14hd1 + gru14hd2 + gru14hd3 + gru14hd4 + gru14hd5 + sg25 + sig26)/(12*mieperho*ld*i01)) 116 | lab var gpcrg9 "Adquiridos fuera del hogar item 47 y 50 y 602" 117 | 118 | gen gpcrg10= ((gru21hd + gru22hd1 + gru22hd2 + gru23hd1 + gru23hd2 + gru23hd3 + gru24hd)/(12*mieperho*ld*i02)) 119 | lab var gpcrg10 "Vestido y calzado" 120 | 121 | gen gpcrg12= ((gru31hd + gru32hd1 + gru32hd2 + gru33hd1 + gru33hd2 + gru33hd3 + gru34hd)/(12*mieperho*ld*i03)) 122 | lab var gpcrg12 "Gasto Alquiler de vivienda y combustible" 123 | 124 | gen gpcrg14= ((gru41hd + gru42hd1 + gru42hd2 + gru43hd1 + gru43hd2 + gru43hd3 + gru44hd + sg421 + sg42d1 + sg423 + sg42d3)/(12*mieperho*ld*i04)) 125 | lab var gpcrg14 "Muebles y enseres" 126 | 127 | gen gpcrg16= ((gru51hd + gru52hd1 + gru52hd2 + gru53hd1 + gru53hd2 + gru53hd3 + gru54hd)/(12*mieperho*ld*i05)) 128 | lab var gpcrg16 "Cuidados de la salud" 129 | 130 | gen gpcrg18= ((gru61hd + gru62hd1 + gru62hd2 + gru63hd1 + gru63hd2 + gru63hd3 + gru64hd + g07hd + ig08hd + sg422 + sg42d2)/(12*mieperho*ld*i06)) 131 | lab var gpcrg18 "Transporte y comunicaciones" 132 | 133 | gen gpcrg19= ((gru71hd + gru72hd1 + gru72hd2 + gru73hd1 + gru73hd2 + gru73hd3 + gru74hd + sg42 + sg42d)/(12*mieperho*ld*i07)) 134 | lab var gpcrg19 "Esparcimiento diversión y cultura" 135 | 136 | gen gpcrg21= ((gru81hd + gru82hd1 + gru82hd2 + gru83hd1 + gru83hd2 + gru83hd3 + gru84hd)/(12*mieperho*ld*i08)) 137 | lab var gpcrg21 "Otros gastos de bienes y servicios" 138 | 139 | 140 | *RECODIFICANDO POR grupo de gastos 141 | gen gpgru2= gpcrg3 142 | lab var gpgru2 "G011.Alimentos dentro del hogar real" 143 | gen gpgru3= gpcrg6 + gpcrg8 + gpcrg9 144 | lab var gpgru3 "G012.Alimentos fuera del hogar real" 145 | gen gpgru1 = gpgru2 + gpgru3 146 | lab var gpgru1 "G01.Total en Alimentos real" 147 | 148 | gen gpgru4 = gpcrg10 149 | lab var gpgru4 "G02.Vestido y calzado real" 150 | 151 | gen gpgru5 = gpcrg12 152 | lab var gpgru5 "G03.Alquiler de Vivienda y combustible real" 153 | 154 | gen gpgru6 = gpcrg14 155 | lab var gpgru6 "G04.Muebles y enseres real" 156 | 157 | gen gpgru7= gpcrg16 158 | lab var gpgru7 "G05.Cuidados de la salud real" 159 | 160 | gen gpgru8 = gpcrg18 161 | lab var gpgru8 "G06.Transportes y comunicaciones real" 162 | 163 | gen gpgru9 = gpcrg19 164 | lab var gpgru9 "G07.Esparcimiento diversion y cultura real" 165 | 166 | gen gpgru10 = gpcrg21 167 | lab var gpgru10 "G08.otros gastos en bienes y servicios real" 168 | 169 | gen gpgru0 = gpgru1 + gpgru4 + gpgru5 + gpgru6 + gpgru7 + gpgru8 + gpgru9 + gpgru10 170 | lab var gpgru0 "Gto real promedio pc mensual" 171 | 172 | *** Salidas *** 173 | svyset [pweight = factornd07], psu(conglome)strata(estrato) 174 | 175 | *** Gasto real promedio percapita mensual *** 176 | svy:mean gpgru0, over(aniorec) 177 | 178 | svy:mean gpgru0 if area==1, over(aniorec) 179 | svy:mean gpgru0 if area==2, over(aniorec) 180 | 181 | preserve 182 | collapse (mean) gpgru0 [aw=factornd07], by(aniorec) 183 | export excel using "$gto/year", firstrow(variables) sheet("all") replace 184 | restore 185 | 186 | preserve 187 | collapse (mean) gpgru0 [aw=factornd07], by(aniorec area) 188 | export excel using "$gto/year", firstrow(variables) sheet("area") 189 | restore 190 | 191 | 192 | *** Grafico de barras 2017/2018 193 | gen temp2017=gpgru0 if aniorec==2017 194 | gen temp2018=gpgru0 if aniorec==2018 195 | lab var temp2017 "2017" 196 | lab var temp2018 "2018" 197 | 198 | graph bar temp2017 temp2018 [pweight = factornd07], over(region, label) blabel(bar, format(%9.0f)) /// 199 | legend(label(1 "2017") label(2 "2018")) /// 200 | bar(2,color(ltblue)) /// 201 | title("Perú: Gasto real promedio per cápita mensual, según regiones naturales, 2017 - 2018",size(small)) /// 202 | subtitle("(Soles constantes base=2018 a precios de Lima Metropolitana)",size(vsmall)) 203 | 204 | graph export "Graph1.png", as(png) replace 205 | 206 | -------------------------------------------------------------------------------- /#52 Medición del Déficit Calórico con la ENAHO - Part1 - Norma Calorica.do: -------------------------------------------------------------------------------- 1 | 2 | cd "D:\DCALORICO" 3 | 4 | ***PROGRAMA PARA DETERMINAR LA NORMA CALORICA 5 | *Programa STATA de estimación de los requerimientos calóricos 2010* 6 | 7 | use enaho01-2011-200.dta, clear 8 | rename a*o anio 9 | gen str13 identh = anio+ conglome+ vivienda+ hogar 10 | gen str15 identi = anio+ conglome+ vivienda+ hogar+ codperso 11 | order identh identi 12 | sort identh 13 | 14 | *rename factor facpob 15 | *facpob07 ya esta en la base 16 | svyset [pweight=facpob07], psu(conglome) strata(estrato) 17 | 18 | *variables geograficas 19 | gen area=estrato<6 if estrato!=. 20 | replace area=1 if dominio==8 & area==0 21 | 22 | gen dominio2=1 if dominio<4 & area==1 23 | replace dominio2=2 if dominio<4 & area==0 24 | replace dominio2=3 if dominio>3 & dominio<7 & area==1 25 | replace dominio2=4 if dominio>3 & dominio<7 & area==0 26 | replace dominio2=5 if dominio==7 & area==1 27 | replace dominio2=6 if dominio==7 & area==0 28 | replace dominio2=7 if dominio==8 29 | 30 | label define dominio2 1 "Costa urbana" 2 "Costa rural" 3 "Sierra urbana" /// 31 | 4 "Sierra rural" 5 "Selva urbana" 6 "Selva rural" 7 "Lima Metrop." 32 | label values dominio2 dominio2 33 | tab dominio2 34 | 35 | gen region=dominio2<3 | dominio2==7 if dominio2!=. 36 | replace region=2 if dominio2>2 & dominio2<5 & dominio2!=. 37 | replace region=3 if dominio2>4 & dominio2<7 & dominio2!=. 38 | label define region 1 "Costa" 2 "Sierra" 3 "Selva" 39 | label values region region 40 | 41 | *solo consideramos a los miembros del hogar 42 | gen miembro=p204 43 | drop if miembro!=1 44 | gen sexo=p207 45 | recode sexo 2=0 46 | drop if sexo==. 47 | *edad en meses de los ninios 48 | gen edad=p208a 49 | gen edadni=p208b if p208a==0 50 | label var edadni "edad en meses de los menores de 1 anio" 51 | 52 | *requerimientos caloricos y TMB 53 | *Fuente: Anne J. Swindale, Punam Ohri-Vachaspati (1997): 54 | *"Household food consumption indicator guide". IMPACT. 55 | 56 | * niños de 6 meses o menores 57 | * requerimientos caloricos 58 | gen reqcal=470 if sexo==1 & edadni<1 59 | replace reqcal=445 if sexo==0 & edadni<1 60 | replace reqcal=550 if sexo==1 & edadni>=1 & edadni <2 61 | replace reqcal=505 if sexo==0 & edadni>=1 & edadni <2 62 | replace reqcal=610 if sexo==1 & edadni>=2 & edadni <3 63 | replace reqcal=545 if sexo==0 & edadni>=2 & edadni <3 64 | replace reqcal=655 if sexo==1 & edadni>=3 & edadni <4 65 | replace reqcal=590 if sexo==0 & edadni>=3 & edadni <4 66 | replace reqcal=695 if sexo==1 & edadni>=4 & edadni <5 67 | replace reqcal=630 if sexo==0 & edadni>=4 & edadni <5 68 | replace reqcal=730 if sexo==1 & edadni>=5 & edadni <6 69 | replace reqcal=670 if sexo==0 & edadni>=5 & edadni <6 70 | 71 | * niños de 6 meses a 11 meses 72 | * requerimientos caloricos 73 | replace reqcal=765 if sexo==1 & edadni==6 74 | replace reqcal=720 if sexo==0 & edadni==6 75 | replace reqcal=810 if sexo==1 & edadni==7 76 | replace reqcal=750 if sexo==0 & edadni==7 77 | replace reqcal=855 if sexo==1 & edadni==8 78 | replace reqcal=800 if sexo==0 & edadni==8 79 | replace reqcal=925 if sexo==1 & edadni==9 80 | replace reqcal=865 if sexo==0 & edadni==9 81 | replace reqcal=970 if sexo==1 & edadni==10 82 | replace reqcal=905 if sexo==0 & edadni==10 83 | replace reqcal=1050 if sexo==1 & edadni==11 84 | replace reqcal=975 if sexo==0 & edadni==11 85 | 86 | * niños de 1 a 5 años 87 | * requerimientos caloricos 88 | replace reqcal=1200 if sexo==1 & edad==1 89 | replace reqcal=1140 if sexo==0 & edad==1 90 | replace reqcal=1410 if sexo==1 & edad==2 91 | replace reqcal=1310 if sexo==0 & edad==2 92 | replace reqcal=1560 if sexo==1 & edad==3 93 | replace reqcal=1440 if sexo==0 & edad==3 94 | replace reqcal=1690 if sexo==1 & edad==4 95 | replace reqcal=1540 if sexo==0 & edad==4 96 | replace reqcal=1810 if sexo==1 & edad==5 97 | replace reqcal=1630 if sexo==0 & edad==5 98 | 99 | * niños de 6 a 9 años 100 | * requerimientos caloricos 101 | replace reqcal=1900 if sexo==1 & edad==6 102 | replace reqcal=1700 if sexo==0 & edad==6 103 | replace reqcal=1990 if sexo==1 & edad==7 104 | replace reqcal=1770 if sexo==0 & edad==7 105 | replace reqcal=2070 if sexo==1 & edad==8 106 | replace reqcal=1830 if sexo==0 & edad==8 107 | replace reqcal=2150 if sexo==1 & edad==9 108 | replace reqcal=1880 if sexo==0 & edad==9 109 | 110 | *tasas de metabolismo basal por edades y sexo 111 | gen tmb=(17.5*30.6)+651 if sexo==1 & edad==10 112 | replace tmb=(12.2*31.7)+746 if sexo==0 & edad==10 113 | replace tmb=(17.5*32.4)+651 if sexo==1 & edad==11 114 | replace tmb=(12.2*35.7)+746 if sexo==0 & edad==11 115 | replace tmb=(17.5*36.5)+651 if sexo==1 & edad==12 116 | replace tmb=(12.2*40.0)+746 if sexo==0 & edad==12 117 | replace tmb=(17.5*41.4)+651 if sexo==1 & edad==13 118 | replace tmb=(12.2*41.6)+746 if sexo==0 & edad==13 119 | replace tmb=(17.5*46.9)+651 if sexo==1 & edad==14 120 | replace tmb=(12.2*47.8)+746 if sexo==0 & edad==14 121 | replace tmb=(17.5*52.3)+651 if sexo==1 & edad==15 122 | replace tmb=(12.2*48.1)+746 if sexo==0 & edad==15 123 | replace tmb=(17.5*53.1)+651 if sexo==1 & edad==16 124 | replace tmb=(12.2*49.8)+746 if sexo==0 & edad==16 125 | replace tmb=(17.5*56.3)+651 if sexo==1 & edad==17 126 | replace tmb=(12.2*50.4)+746 if sexo==0 & edad==17 127 | replace tmb=(15.3*58.2)+679 if sexo==1 & edad>17 & edad<30 128 | replace tmb=(14.7*51.0)+496 if sexo==0 & edad>17 & edad<30 129 | replace tmb=(11.6*58.2)+879 if sexo==1 & edad>29 & edad<60 130 | replace tmb=(8.7*51.0)+829 if sexo==0 & edad>29 & edad<60 131 | replace tmb=(13.5*58.2)+487 if sexo==1 & edad>59 & edad!=. 132 | replace tmb=(10.5*51.0)+596 if sexo==0 & edad>59 & edad!=. 133 | label var tmb "tasa metabolica basal" 134 | 135 | *tasas de metabolismo basal corregidos por tipo de actividad 136 | *fuente de los coeficientes de correccion por niveles de actividad: 137 | *José Maria Bengoa, Benjamin Torùn, Moisés Bahar y Nevin Scrimshaw: 138 | *Food nutrition Bulletin, vol 11, n°1. p.8 cuadro 1.The United Nations University 139 | 140 | *tmbam= actividad moderada 141 | *tmbai= actividad intensa 142 | 143 | *tmbal= actividad ligera 144 | gen tmbal=tmb*1.75 if sexo==1 & edad>9 & edad<13 145 | replace tmbal=tmb*1.64 if sexo==0 & edad>9 & edad<13 146 | replace tmbal=tmb*1.68 if sexo==1 & edad>12 & edad<15 147 | replace tmbal=tmb*1.59 if sexo==0 & edad>12 & edad<15 148 | replace tmbal=tmb*1.62 if sexo==1 & edad>14 & edad<18 149 | replace tmbal=tmb*1.55 if sexo==0 & edad>14 & edad<18 150 | replace tmbal=tmb*1.55 if sexo==1 & edad>17 & edad<66 151 | replace tmbal=tmb*1.55 if sexo==0 & edad>17 & edad<66 152 | replace tmbal=tmb*1.40 if sexo==1 & edad>65 & edad!=. 153 | replace tmbal=tmb*1.40 if sexo==0 & edad>65 & edad!=. 154 | label var tmbal "tmb con actividad ligera" 155 | 156 | *tmbam= actividad moderada 157 | gen tmbam=tmb*1.75 if sexo==1 & edad>9 & edad<13 158 | replace tmbam=tmb*1.64 if sexo==0 & edad>9 & edad<13 159 | replace tmbam=tmb*1.68 if sexo==1 & edad>12 & edad<15 160 | replace tmbam=tmb*1.59 if sexo==0 & edad>12 & edad<15 161 | replace tmbam=tmb*1.80 if sexo==1 & edad>14 & edad<18 162 | replace tmbam=tmb*1.65 if sexo==0 & edad>14 & edad<18 163 | replace tmbam=tmb*1.80 if sexo==1 & edad>17 & edad<66 164 | replace tmbam=tmb*1.65 if sexo==0 & edad>17 & edad<66 165 | replace tmbam=tmb*1.60 if sexo==1 & edad>65 & edad!=. 166 | replace tmbam=tmb*1.60 if sexo==0 & edad>65 & edad!=. 167 | label var tmbam "tmb con actividad moderada" 168 | 169 | *tmbai= actividad intensa 170 | gen tmbai=tmb*1.75 if sexo==1 & edad>9 & edad<13 171 | replace tmbai=tmb*1.64 if sexo==0 & edad>9 & edad<13 172 | replace tmbai=tmb*1.68 if sexo==1 & edad>12 & edad<15 173 | replace tmbai=tmb*1.59 if sexo==0 & edad>12 & edad<15 174 | replace tmbai=tmb*2.10 if sexo==1 & edad>14 & edad<18 175 | replace tmbai=tmb*1.80 if sexo==0 & edad>14 & edad<18 176 | replace tmbai=tmb*2.10 if sexo==1 & edad>17 & edad<66 177 | replace tmbai=tmb*1.80 if sexo==0 & edad>17 & edad<66 178 | replace tmbai=tmb*1.90 if sexo==1 & edad>65 & edad!=. 179 | replace tmbai=tmb*1.80 if sexo==0 & edad>65 & edad!=. 180 | label var tmbai "tmb con actividad intensa" 181 | 182 | *requerimientos calóricos individuales con actividad moderada 183 | gen reqcalam=reqcal 184 | replace reqcalam=tmbam if edad>9 & edad!=. 185 | label var reqcalam "r. calori. ind. act. moderada en las a. urb. y en las a. rur." 186 | note reqcalam: requerimientos caloricos individuales con actividad moderada en las areas urbanas y en las areas rurales 187 | 188 | *requerimientos calóricos individuales con actividad moderada 189 | *en las areas urbanas e intensa en las areas rurales 190 | gen reqcalad=reqcalam 191 | replace reqcalad=tmbai if area==0 & edad>9 & edad!=. 192 | label var reqcalad "r. calori. ind. act. moderada en las a. urb. e intensa en las a. rur." 193 | note reqcalad: requerimientos caloricos individuales con actividad moderada en las areas urbanas e intensa en las areas rurales 194 | 195 | *test de diferencias en los requerimientos caloricos con actividad moderada 196 | *en area urbana y en area rural 197 | svy: mean reqcalam, over(area) 198 | lincom reqcalam#1.area-reqcalam#0.area 199 | 200 | *por dominios respecto a Lima 201 | svy: mean reqcalam, over(dominio2) 202 | lincom reqcalam#7.dominio2 - reqcalam#1.dominio2 203 | lincom reqcalam#7.dominio2 - reqcalam#2.dominio2 204 | lincom reqcalam#7.dominio2 - reqcalam#3.dominio2 205 | lincom reqcalam#7.dominio2 - reqcalam#4.dominio2 206 | lincom reqcalam#7.dominio2 - reqcalam#5.dominio2 207 | lincom reqcalam#7.dominio2 - reqcalam#6.dominio2 208 | 209 | *por region natural respecto a la Costa 210 | svy: mean reqcalam, over(region) 211 | lincom reqcalam#1.region-reqcalam#2.region 212 | lincom reqcalam#1.region-reqcalam#3.region 213 | 214 | *entre sierra y selva 215 | lincom reqcalam#2.region-reqcalam#3.region 216 | 217 | collapse (sum) reqcalam reqcalad, by (identh) 218 | save normacalorica2011, replace -------------------------------------------------------------------------------- /#21 Evolución del ingreso real per cápita promedio mensual ENAHO 2016 – 2019.do: -------------------------------------------------------------------------------- 1 | cd "D:\Evolucion" //Cambiar el nombre segun el folder donde van a bajar la informacion 2 | 3 | copy "http://iinei.inei.gob.pe/iinei/srienaho/descarga/STATA/546-Modulo34.zip" 546-Modulo34.zip, replace 4 | copy "http://iinei.inei.gob.pe/iinei/srienaho/descarga/STATA/603-Modulo34.zip" 603-Modulo34.zip, replace 5 | copy "http://iinei.inei.gob.pe/iinei/srienaho/descarga/STATA/634-Modulo34.zip" 634-Modulo34.zip, replace 6 | copy "http://iinei.inei.gob.pe/iinei/srienaho/descarga/STATA/687-Modulo34.zip" 687-Modulo34.zip, replace 7 | *Descomprimo los archivos manualmente y los coloco en la carpeta 8 | 9 | 10 | *Agregamos las bases sumaria del 2016, 2017, 2018 & 2019 11 | use sumaria-2016.dta, clear 12 | append using sumaria-2017.dta 13 | append using sumaria-2018.dta 14 | append using sumaria-2019.dta 15 | 16 | destring conglome, replace 17 | tostring conglome, replace format(%06.0f) 18 | 19 | recode ingtpu01 ingtpu02 ingtpu03 ingtpu04 ingtpu05 ig03hd1 ig03hd2 ig03hd3 ig03hd4 (.= 0) 20 | 21 | *Renombro la variable ahno 22 | rename a*o anio 23 | gen aniorec=real(anio) 24 | 25 | *Genero variable departamento 26 | gen dpto= real(substr(ubigeo,1,2)) 27 | replace dpto=15 if (dpto==7) 28 | label define dpto 1"Amazonas" 2"Ancash" 3"Apurimac" 4"Arequipa" 5"Ayacucho" /* 29 | */ 6"Cajamarca" 8"Cusco" 9"Huancavelica" 10"Huanuco" 11"Ica" 12"Junin" /* 30 | */ 13"La_Libertad" 14"Lambayeque" 15"Lima" 16"Loreto" 17"Madre_de_Dios" 18"Moquegua" /* 31 | */ 19"Pasco" 20"Piura" 21"Puno" 22"San_Martin" 23"Tacna" 24"Tumbes" 25"Ucayali" 32 | lab val dpto dpto 33 | 34 | *Agrego la base que contiene los deflactores 35 | sort aniorec dpto 36 | merge aniorec dpto using deflactores_base2019_new.dta 37 | tab _m 38 | drop if _merge==2 39 | drop _m 40 | 41 | *Genero las variables de area de residencia, region natural y dominio geografico 42 | replace estrato = 1 if dominio ==8 43 | gen area = estrato <6 44 | replace area = 2 if area==0 45 | label define area 2 rural 1 urbana 46 | label val area area 47 | 48 | gen region=1 if dominio>=1 & dominio<=3 49 | replace region=1 if dominio==8 50 | replace region=2 if dominio>=4 & dominio<=6 51 | replace region=3 if dominio==7 52 | label define region 1 "Costa" 2 "Sierra" 3 "Selva" 53 | 54 | gen domin02=1 if dominio>=1 & dominio<=3 & area==1 55 | replace domin02=2 if dominio>=1 & dominio<=3 & area==2 56 | replace domin02=3 if dominio>=4 & dominio<=6 & area==1 57 | replace domin02=4 if dominio>=4 & dominio<=6 & area==2 58 | replace domin02=5 if dominio==7 & area==1 59 | replace domin02=6 if dominio==7 & area==2 60 | replace domin02=7 if dominio==8 61 | label define domin02 1 "Costa_urbana" 2 "Costa_rural" 3 "Sierra_urbana" /* 62 | */ 4 "Sierra_rural" 5 "Selva_urbana" 6 "Selva_rural" 7 "Lima_Metropolitana" 63 | label value domin02 domin02 64 | 65 | gen dominioA=1 if dominio==1 & area==1 66 | replace dominioA=2 if dominio==1 & area==2 67 | replace dominioA=3 if dominio==2 & area==1 68 | replace dominioA=4 if dominio==2 & area==2 69 | replace dominioA=5 if dominio==3 & area==1 70 | replace dominioA=6 if dominio==3 & area==2 71 | replace dominioA=7 if dominio==4 & area==1 72 | replace dominioA=8 if dominio==4 & area==2 73 | replace dominioA=9 if dominio==5 & area==1 74 | replace dominioA=10 if dominio==5 & area==2 75 | replace dominioA=11 if dominio==6 & area==1 76 | replace dominioA=12 if dominio==6 & area==2 77 | replace dominioA=13 if dominio==7 & area==1 78 | replace dominioA=14 if dominio==7 & area==2 79 | replace dominioA=15 if dominio==7 & (dpto==16 | dpto==17 | dpto==25) & area==1 80 | replace dominioA=16 if dominio==7 & (dpto==16 | dpto==17 | dpto==25) & area==2 81 | replace dominioA=17 if dominio==8 & area==1 82 | replace dominioA=17 if dominio==8 & area==2 83 | 84 | label define dominioA 1 "Costa norte urbana" 2 "Costa norte rural" 3 "Costa centro urbana" 4 "Costa centro rural" /* 85 | */ 5 "Costa sur urbana" 6 "Costa sur rural" 7 "Sierra norte urbana" 8 "Sierra norte rural" 9 "Sierra centro urbana" /* 86 | */ 10 "Sierra centro rural" 11 "Sierra sur urbana" 12 "Sierra sur rural" 13 "Selva alta urbana" 14 "Selva alta rural" /* 87 | */ 15 "Selva baja urbana" 16 "Selva baja rural" 17"Lima Metropolitana" 88 | lab val dominioA dominioA 89 | 90 | drop ld 91 | 92 | sort dominioA 93 | 94 | merge dominioA using despacial_ldnew.dta 95 | tab _m 96 | drop _m 97 | 98 | 99 | *Calculamos el factor de ponderacion 100 | gen factornd07=round(factor07*mieperho,1) 101 | 102 | svyset [pweight = factornd07], psu(conglome) 103 | 104 | *** Ingresos *** 105 | gen ipcr_2 = (ingbruhd +ingindhd)/(12*mieperho*ld*i00) 106 | gen ipcr_3 = (insedthd + ingseihd)/(12*mieperho*ld*i00) 107 | gen ipcr_4 = (pagesphd + paesechd + ingauthd + isecauhd)/(12*mieperho*ld*i00) 108 | gen ipcr_5 = (ingexthd)/(12*mieperho*ld*i00) 109 | gen ipcr_1 = (ipcr_2 + ipcr_3 + ipcr_4 + ipcr_5) 110 | 111 | gen ipcr_7 = (ingtrahd)/(12*mieperho*ld*i00) 112 | gen ipcr_8 = (ingtexhd)/(12*mieperho*ld*i00) 113 | gen ipcr_6 = (ipcr_7 + ipcr_8) 114 | 115 | gen ipcr_9 = (ingtprhd)/(12*mieperho*ld*i00) 116 | gen ipcr_10 = (ingtpuhd)/(12*mieperho*ld*i00) 117 | gen ipcr_11 = (ingtpu01)/(12*mieperho*ld*i00) 118 | gen ipcr_12 = (ingtpu03)/(12*mieperho*ld*i00) 119 | gen ipcr_13 = (ingtpu05)/(12*mieperho*ld*i00) 120 | gen ipcr_14 = (ingtpu04)/(12*mieperho*ld*i00) 121 | gen ipcr_15 = (ingtpu02)/(12*mieperho*ld*i00) 122 | gen ipcr_16 = (ingrenhd)/(12*mieperho*ld*i00) 123 | gen ipcr_17 = (ingoexhd + gru13hd3 + gru23hd3 + gru33hd3 + gru43hd3 + gru53hd3 + gru63hd3 + gru73hd3 + /* 124 | */ gru83hd3 + gru24hd +gru44hd + gru54hd + gru74hd + gru84hd + gru14hd5)/(12*mieperho*ld*i00) 125 | 126 | *ajuste por el alquiler imputado 127 | gen ipcr_18 =(ia01hd +gru34hd - ga04hd + gru64hd)/(12*mieperho*ld*i00) 128 | 129 | gen ipcr_19 = (gru13hd1 + sig24 + gru23hd1 + gru33hd1 + gru43hd1 + gru53hd1 + gru63hd1 + gru73hd1 + gru83hd1 /* 130 | */+ gru14hd3 + sig26)/(12*mieperho*ld*i00) 131 | 132 | gen ipcr_20 = (gru13hd2 + ig06hd + gru23hd2 + gru33hd2 + gru43hd2 + gru53hd2 + gru63hd2 + ig08hd + gru73hd2 + /* 133 | */ gru83hd2 + gru14hd4 + sg42d + sg42d1 + sg42d2 + sg42d3)/(12*mieperho*ld*i00) 134 | 135 | 136 | gen ipcr_0= ipcr_2 + ipcr_3 + ipcr_4 + ipcr_5+ ipcr_7 + ipcr_8 + ipcr_16 + ipcr_17 + ipcr_18 + ipcr_19 + ipcr_20 137 | 138 | label var ipcr_0 "Ingreso percapita mensual a precios de Lima monetario" 139 | label var ipcr_1 "Ingreso percapita mensual a precios de Lima monetario por trabajo" 140 | label var ipcr_2 "Ingreso percapita mensual a precios de Lima monetario por trabajo principal" 141 | label var ipcr_3 "Ingreso percapita mensual a precios de Lima monetario por trabajo secundario" 142 | label var ipcr_4 "Ingreso percapita mensual a precios de Lima pago en especie y autocon" 143 | label var ipcr_5 "Ingreso percapita mensual a precios de Lima pago extraordinario por trabajo" 144 | label var ipcr_6 "Ingreso percapita mensual a precios de Lima transferencia corriente" 145 | label var ipcr_7 "Ingreso percapita mensual a precios de Lima transferencia monetaria del pais" 146 | label var ipcr_8 "Ingreso percapita mensual a precios de Lima transferencia monetaria extranjero" 147 | label var ipcr_9 "Ingreso percapita mensual a precios de Lima transferencia monetaria privada" 148 | label var ipcr_10 "Ingreso percapita mensual a precios de Lima transferencia monetaria Publica total" 149 | label var ipcr_11 "Ingreso percapita mensual a precios de Lima transferencia monetaria Publica Juntos" 150 | label var ipcr_12 "Ingreso percapita mensual a precios de Lima transferencia monetaria Publica Pensión65" 151 | label var ipcr_13 "Ingreso percapita mensual a precios de Lima transferencia monetaria Bono Gas" 152 | label var ipcr_14 "Ingreso percapita mensual a precios de Lima transferencia monetaria Beca 18" 153 | label var ipcr_15 "Ingreso percapita mensual a precios de Lima transferencia monetaria Otros Publica" 154 | label var ipcr_16 "Ingreso percapita mensual a precios de Lima renta" 155 | label var ipcr_17 "Ingreso percapita mensual a precios de Lima extraordinario" 156 | label var ipcr_18 "Ingreso percapita mensual a precios de Lima alquiler imputado" 157 | label var ipcr_19 "Ingreso percapita mensual a precios de Lima donacion publica" 158 | label var ipcr_20 "Ingreso percapita mensual a precios de Lima donacion privada" 159 | 160 | 161 | *** Salidas *** 162 | svyset [pweight = factornd07], psu(conglome)strata(estrato) 163 | 164 | *** Ingreso real promedio per capita mensual *** 165 | svy:mean ipcr_0, over(aniorec) 166 | 167 | outreg using "D:\Evolucion\pobreza.doc", replace stat(b se ci) nosubstat bdec(1) varlabels /// 168 | ctitle("Year", "Ingreso real promedio per capita mensual","Error Estandar","Intervalo de confianza al 95%") /// 169 | title("EVOLUCION DEL INGRESO REAL PROMEDIO PER CAPITA MENSUAL") note(""\"Fuente: ENAHO 2016-2019.") 170 | 171 | svy:mean ipcr_0 if area==1, over(aniorec) 172 | outreg using "D:\Evolucion\pobreza.doc", addtable stat(b se ci) nosubstat bdec(1) /// 173 | ctitle("Year", "Ingreso real promedio per capita mensual","Error Estandar","Intervalo de confianza al 95%") /// 174 | title("EVOLUCION DEL INGRESO REAL PROMEDIO PER CAPITA MENSUAL, URBANA") note(""\"Fuente: ENAHO 2016-2019.") 175 | 176 | svy:mean ipcr_0 if area==2, over(aniorec) 177 | outreg using "D:\Evolucion\pobreza.doc", addtable stat(b se ci) nosubstat bdec(1) /// 178 | ctitle("Year", "Ingreso real promedio per capita mensual","Error Estandar","Intervalo de confianza al 95%") /// 179 | title("EVOLUCION DEL INGRESO REAL PROMEDIO PER CAPITA MENSUAL, RURAL") note(""\"Fuente: ENAHO 2016-2019.") 180 | 181 | svy:mean ipcr_0 if region==1, over(aniorec) 182 | outreg using "D:\Evolucion\pobreza.doc", addtable stat(b se ci) nosubstat bdec(1) /// 183 | ctitle("Year", "Ingreso real promedio per capita mensual","Error Estandar","Intervalo de confianza al 95%") /// 184 | title("EVOLUCION DEL INGRESO REAL PROMEDIO PER CAPITA MENSUAL, COSTA") note(""\"Fuente: ENAHO 2016-2019.") 185 | 186 | svy:mean ipcr_0 if region==2, over(aniorec) 187 | outreg using "D:\Evolucion\pobreza.doc", addtable stat(b se ci) nosubstat bdec(1) /// 188 | ctitle("Year", "Ingreso real promedio per capita mensual","Error Estandar","Intervalo de confianza al 95%") /// 189 | title("EVOLUCION DEL INGRESO REAL PROMEDIO PER CAPITA MENSUAL, SIERRA") note(""\"Fuente: ENAHO 2016-2019.") 190 | 191 | svy:mean ipcr_0 if region==3, over(aniorec) 192 | outreg using "D:\Evolucion\pobreza.doc", addtable stat(b se ci) nosubstat bdec(1) /// 193 | ctitle("Year", "Ingreso real promedio per capita mensual","Error Estandar","Intervalo de confianza al 95%") /// 194 | title("EVOLUCION DEL INGRESO REAL PROMEDIO PER CAPITA MENSUAL, SELVA") note(""\"Fuente: ENAHO 2016-2019.") 195 | 196 | -------------------------------------------------------------------------------- /#52 Medición del Déficit Calórico con la ENAHO - Part2 - Calculo de calorias adquiridas.do: -------------------------------------------------------------------------------- 1 | 2 | cd "D:\DCALORICO" 3 | 4 | ****CALCULO DE LAS CALORIAS ADQUIRIDAS POR EL HOGAR 5 | 6 | *trabajo con archivo sumaria 7 | use sumaria-2011.dta, clear 8 | compress 9 | rename a*o anio 10 | gen str13 identh=anio+conglome+vivienda+hogar 11 | gen gpcm =gashog2d/(mieperho*12) 12 | gen facpob =factor*mieperho 13 | order identh 14 | sort identh 15 | save tempo_sumaria-2011.dta, replace 16 | 17 | keep identh mieperho gpcm facpob 18 | sort identh 19 | save tempo2011, replace 20 | 21 | 22 | *Cálculo de calorias adquiridas por rubro de gasto 23 | 24 | *Usando archivo 601 25 | use enaho01-2011-601.dta, clear 26 | keep if p601b==1 27 | keep if substr(p601a,1,2)!="49" 28 | drop if substr(p601a,3,2)=="00" 29 | 30 | rename a*o anio 31 | gen str13 identh=anio+conglome+vivienda+hogar 32 | order identh 33 | sort identh 34 | merge identh using tempo2011 35 | tab _m 36 | drop if _m==2 37 | drop _m 38 | sort produc61 39 | 40 | *USANDO TABLA DE COMPOSICION DE ALIMENTOS en archivo 601 41 | merge m:1 produc61 using "tabla de composicion de alimentos.dta" 42 | tab _m 43 | drop if _m==2 44 | drop _m 45 | 46 | gen regnat=1 if (dominio>=1 & dominio<=3) 47 | replace regnat=2 if (dominio>=4 & dominio<=6) 48 | replace regnat=3 if (dominio==7) 49 | replace regnat=1 if (dominio==8) 50 | label define regnat 1 "costa" 2 "sierra" 3 "selva" 51 | label val regnat regnat 52 | 53 | gen dominio2=1 if ((dominio==1 | dominio==2 | dominio==3) & (estrato >=1 & estrato<=5)) 54 | replace dominio2=2 if ((dominio==1 | dominio==2 | dominio==3) & (estrato >=6 & estrato<=8)) 55 | replace dominio2=3 if ((dominio==4 | dominio==5 | dominio==6) & (estrato >=1 & estrato<=5)) 56 | replace dominio2=4 if ((dominio==4 | dominio==5 | dominio==6) & (estrato >=6 & estrato<=8)) 57 | replace dominio2=5 if ((dominio==7 ) & (estrato >=1 & estrato<=5)) 58 | replace dominio2=6 if ((dominio==7 ) & (estrato >=6 & estrato<=8)) 59 | replace dominio2=7 if (dominio==8 ) 60 | label define dominio2 1 "costa urbana" 2 "costa rural" 3 "sierra urbana" 4 "sierra rural" 5 "selva urbana" 6 "selva rural" 7 "lima metropolitana" 61 | label val dominio2 dominio2 62 | 63 | gen p601n=real(substr(p601a,1,2)) 64 | *para productos que no se encontraron en tabla se usan promedios de calorías por rubro 65 | replace energ= 313.6 if ((regnat == 1) & (p601n == 1) & energ==.) 66 | replace energ= 361.0 if ((regnat == 1) & (p601n == 2) & energ==.) 67 | replace energ= 359.0 if ((regnat == 1) & (p601n == 3) & energ==.) 68 | replace energ= 101.7 if ((regnat == 1) & (p601n == 4) & energ==.) 69 | replace energ= 98.4 if ((regnat == 1) & (p601n == 5) & energ==.) 70 | replace energ= 384.0 if ((regnat == 1) & (p601n == 6) & energ==.) 71 | replace energ= 141.0 if ((regnat == 1) & (p601n == 7) & energ==.) 72 | replace energ= 147.0 if ((regnat == 1) & (p601n == 8) & energ==.) 73 | replace energ= 164.6 if ((regnat == 1) & (p601n == 9) & energ==.) 74 | replace energ= 157.0 if ((regnat == 1) & (p601n == 10) & energ==.) 75 | replace energ= 366.0 if ((regnat == 1) & (p601n == 11) & energ==.) 76 | replace energ= 127.0 if ((regnat == 1) & (p601n == 12) & energ==.) 77 | replace energ= 104.0 if ((regnat == 1) & (p601n == 13) & energ==.) 78 | replace energ= 121.0 if ((regnat == 1) & (p601n == 14) & energ==.) 79 | replace energ= 363.0 if ((regnat == 1) & (p601n == 15) & energ==.) 80 | replace energ= 347.6 if ((regnat == 1) & (p601n == 16) & energ==.) 81 | replace energ= 374.0 if ((regnat == 1) & (p601n == 17) & energ==.) 82 | replace energ= 350.4 if ((regnat == 1) & (p601n == 18) & energ==.) 83 | replace energ= 355.0 if ((regnat == 1) & (p601n == 19) & energ==.) 84 | replace energ= 116.3 if ((regnat == 1) & (p601n == 20) & energ==.) 85 | replace energ= 235.3 if ((regnat == 1) & (p601n == 21) & energ==.) 86 | replace energ= 80.0 if ((regnat == 1) & (p601n == 22) & energ==.) 87 | replace energ= 889.0 if ((regnat == 1) & (p601n == 23) & energ==.) 88 | replace energ= 230.0 if ((regnat == 1) & (p601n == 24) & energ==.) 89 | replace energ= 720.0 if ((regnat == 1) & (p601n == 25) & energ==.) 90 | replace energ= 729.0 if ((regnat == 1) & (p601n == 26) & energ==.) 91 | replace energ= 67.0 if ((regnat == 1) & (p601n == 27) & energ==.) 92 | replace energ= 0.0 if ((regnat == 1) & (p601n == 28) & energ==.) 93 | replace energ= 40.0 if ((regnat == 1) & (p601n == 29) & energ==.) 94 | replace energ= 81.0 if ((regnat == 1) & (p601n == 30) & energ==.) 95 | replace energ= 339.0 if ((regnat == 1) & (p601n == 31) & energ==.) 96 | replace energ= 49.0 if ((regnat == 1) & (p601n == 32) & energ==.) 97 | replace energ= 19.0 if ((regnat == 1) & (p601n == 33) & energ==.) 98 | replace energ= 33.4 if ((regnat == 1) & (p601n == 34) & energ==.) 99 | replace energ= 129.0 if ((regnat == 1) & (p601n == 35) & energ==.) 100 | replace energ= 128.5 if ((regnat == 1) & (p601n == 36) & energ==.) 101 | replace energ= 33.0 if ((regnat == 1) & (p601n == 37) & energ==.) 102 | replace energ= 30.0 if ((regnat == 1) & (p601n == 38) & energ==.) 103 | replace energ= 34.2 if ((regnat == 1) & (p601n == 39) & energ==.) 104 | replace energ= 86.2 if ((regnat == 1) & (p601n == 40) & energ==.) 105 | replace energ= 52.0 if ((regnat == 1) & (p601n == 41) & energ==.) 106 | replace energ= 346.2 if ((regnat == 1) & (p601n == 42) & energ==.) 107 | replace energ= 355.0 if ((regnat == 1) & (p601n == 43) & energ==.) 108 | replace energ= 36.0 if ((regnat == 1) & (p601n == 44) & energ==.) 109 | replace energ= 40.0 if ((regnat == 1) & (p601n == 45) & energ==.) 110 | replace energ= 22.0 if ((regnat == 1) & (p601n == 46) & energ==.) 111 | replace energ= 149.6 if ((regnat == 1) & (p601n == 47) & energ==.) 112 | replace energ= 188.4 if ((regnat == 1) & (p601n == 48) & energ==.) 113 | replace energ= 188.4 if ((regnat == 1) & (p601n == 50) & energ==.) 114 | replace energ= 305.0 if ((regnat == 2) & (p601n == 1) & energ==.) 115 | replace energ= 361.0 if ((regnat == 2) & (p601n == 2) & energ==.) 116 | replace energ= 359.0 if ((regnat == 2) & (p601n == 3) & energ==.) 117 | replace energ= 139.1 if ((regnat == 2) & (p601n == 4) & energ==.) 118 | replace energ= 97.5 if ((regnat == 2) & (p601n == 5) & energ==.) 119 | replace energ= 384.0 if ((regnat == 2) & (p601n == 6) & energ==.) 120 | replace energ= 141.0 if ((regnat == 2) & (p601n == 7) & energ==.) 121 | replace energ= 164.4 if ((regnat == 2) & (p601n == 8) & energ==.) 122 | replace energ= 161.3 if ((regnat == 2) & (p601n == 9) & energ==.) 123 | replace energ= 157.0 if ((regnat == 2) & (p601n == 10) & energ==.) 124 | replace energ= 366.0 if ((regnat == 2) & (p601n == 11) & energ==.) 125 | replace energ= 127.0 if ((regnat == 2) & (p601n == 12) & energ==.) 126 | replace energ= 104.0 if ((regnat == 2) & (p601n == 13) & energ==.) 127 | replace energ= 121.0 if ((regnat == 2) & (p601n == 14) & energ==.) 128 | replace energ= 363.0 if ((regnat == 2) & (p601n == 15) & energ==.) 129 | replace energ= 348.9 if ((regnat == 2) & (p601n == 16) & energ==.) 130 | replace energ= 374.0 if ((regnat == 2) & (p601n == 17) & energ==.) 131 | replace energ= 343.9 if ((regnat == 2) & (p601n == 18) & energ==.) 132 | replace energ= 346.3 if ((regnat == 2) & (p601n == 19) & energ==.) 133 | replace energ= 114.0 if ((regnat == 2) & (p601n == 20) & energ==.) 134 | replace energ= 288.0 if ((regnat == 2) & (p601n == 21) & energ==.) 135 | replace energ= 80.0 if ((regnat == 2) & (p601n == 22) & energ==.) 136 | replace energ= 889.0 if ((regnat == 2) & (p601n == 23) & energ==.) 137 | replace energ= 230.0 if ((regnat == 2) & (p601n == 24) & energ==.) 138 | replace energ= 720.0 if ((regnat == 2) & (p601n == 25) & energ==.) 139 | replace energ= 729.0 if ((regnat == 2) & (p601n == 26) & energ==.) 140 | replace energ= 67.0 if ((regnat == 2) & (p601n == 27) & energ==.) 141 | replace energ= 0.0 if ((regnat == 2) & (p601n == 28) & energ==.) 142 | replace energ= 40.0 if ((regnat == 2) & (p601n == 29) & energ==.) 143 | replace energ= 212.0 if ((regnat == 2) & (p601n == 30) & energ==.) 144 | replace energ= 339.0 if ((regnat == 2) & (p601n == 31) & energ==.) 145 | replace energ= 49.0 if ((regnat == 2) & (p601n == 32) & energ==.) 146 | replace energ= 19.0 if ((regnat == 2) & (p601n == 33) & energ==.) 147 | replace energ= 33.5 if ((regnat == 2) & (p601n == 34) & energ==.) 148 | replace energ= 129.0 if ((regnat == 2) & (p601n == 35) & energ==.) 149 | replace energ= 130.3 if ((regnat == 2) & (p601n == 36) & energ==.) 150 | replace energ= 36.0 if ((regnat == 2) & (p601n == 37) & energ==.) 151 | replace energ= 30.0 if ((regnat == 2) & (p601n == 38) & energ==.) 152 | replace energ= 33.7 if ((regnat == 2) & (p601n == 39) & energ==.) 153 | replace energ= 86.9 if ((regnat == 2) & (p601n == 40) & energ==.) 154 | replace energ= 56.1 if ((regnat == 2) & (p601n == 41) & energ==.) 155 | replace energ= 390.3 if ((regnat == 2) & (p601n == 42) & energ==.) 156 | replace energ= 353.7 if ((regnat == 2) & (p601n == 43) & energ==.) 157 | replace energ= 36.0 if ((regnat == 2) & (p601n == 44) & energ==.) 158 | replace energ= 40.0 if ((regnat == 2) & (p601n == 45) & energ==.) 159 | replace energ= 22.0 if ((regnat == 2) & (p601n == 46) & energ==.) 160 | replace energ= 149.6 if ((regnat == 2) & (p601n == 47) & energ==.) 161 | replace energ= 196.9 if ((regnat == 2) & (p601n == 48) & energ==.) 162 | replace energ= 196.9 if ((regnat == 2) & (p601n == 50) & energ==.) 163 | replace energ= 329.8 if ((regnat == 3) & (p601n == 1) & energ==.) 164 | replace energ= 361.0 if ((regnat == 3) & (p601n == 2) & energ==.) 165 | replace energ= 359.0 if ((regnat == 3) & (p601n == 3) & energ==.) 166 | replace energ= 165.9 if ((regnat == 3) & (p601n == 4) & energ==.) 167 | replace energ= 97.4 if ((regnat == 3) & (p601n == 5) & energ==.) 168 | replace energ= 384.0 if ((regnat == 3) & (p601n == 6) & energ==.) 169 | replace energ= 141.0 if ((regnat == 3) & (p601n == 7) & energ==.) 170 | replace energ= 130.5 if ((regnat == 3) & (p601n == 8) & energ==.) 171 | replace energ= 168.5 if ((regnat == 3) & (p601n == 9) & energ==.) 172 | replace energ= 157.0 if ((regnat == 3) & (p601n == 10) & energ==.) 173 | replace energ= 366.0 if ((regnat == 3) & (p601n == 11) & energ==.) 174 | replace energ= 127.0 if ((regnat == 3) & (p601n == 12) & energ==.) 175 | replace energ= 104.0 if ((regnat == 3) & (p601n == 13) & energ==.) 176 | replace energ= 121.0 if ((regnat == 3) & (p601n == 14) & energ==.) 177 | replace energ= 363.0 if ((regnat == 3) & (p601n == 15) & energ==.) 178 | replace energ= 349.3 if ((regnat == 3) & (p601n == 16) & energ==.) 179 | replace energ= 374.0 if ((regnat == 3) & (p601n == 17) & energ==.) 180 | replace energ= 345.0 if ((regnat == 3) & (p601n == 18) & energ==.) 181 | replace energ= 339.9 if ((regnat == 3) & (p601n == 19) & energ==.) 182 | replace energ= 147.4 if ((regnat == 3) & (p601n == 20) & energ==.) 183 | replace energ= 235.3 if ((regnat == 3) & (p601n == 21) & energ==.) 184 | replace energ= 80.0 if ((regnat == 3) & (p601n == 22) & energ==.) 185 | replace energ= 889.0 if ((regnat == 3) & (p601n == 23) & energ==.) 186 | replace energ= 230.0 if ((regnat == 3) & (p601n == 24) & energ==.) 187 | replace energ= 720.0 if ((regnat == 3) & (p601n == 25) & energ==.) 188 | replace energ= 729.0 if ((regnat == 3) & (p601n == 26) & energ==.) 189 | replace energ= 67.0 if ((regnat == 3) & (p601n == 27) & energ==.) 190 | replace energ= 0.0 if ((regnat == 3) & (p601n == 28) & energ==.) 191 | replace energ= 40.0 if ((regnat == 3) & (p601n == 29) & energ==.) 192 | replace energ= 62.0 if ((regnat == 3) & (p601n == 30) & energ==.) 193 | replace energ= 339.0 if ((regnat == 3) & (p601n == 31) & energ==.) 194 | replace energ= 49.0 if ((regnat == 3) & (p601n == 32) & energ==.) 195 | replace energ= 19.0 if ((regnat == 3) & (p601n == 33) & energ==.) 196 | replace energ= 32.2 if ((regnat == 3) & (p601n == 34) & energ==.) 197 | replace energ= 129.0 if ((regnat == 3) & (p601n == 35) & energ==.) 198 | replace energ= 152.7 if ((regnat == 3) & (p601n == 36) & energ==.) 199 | replace energ= 34.0 if ((regnat == 3) & (p601n == 37) & energ==.) 200 | replace energ= 30.0 if ((regnat == 3) & (p601n == 38) & energ==.) 201 | replace energ= 33.3 if ((regnat == 3) & (p601n == 39) & energ==.) 202 | replace energ= 86.4 if ((regnat == 3) & (p601n == 40) & energ==.) 203 | replace energ= 73.3 if ((regnat == 3) & (p601n == 41) & energ==.) 204 | replace energ= 328.2 if ((regnat == 3) & (p601n == 42) & energ==.) 205 | replace energ= 432.6 if ((regnat == 3) & (p601n == 43) & energ==.) 206 | replace energ= 36.0 if ((regnat == 3) & (p601n == 44) & energ==.) 207 | replace energ= 40.0 if ((regnat == 3) & (p601n == 45) & energ==.) 208 | replace energ= 22.0 if ((regnat == 3) & (p601n == 46) & energ==.) 209 | replace energ= 149.6 if ((regnat == 3) & (p601n == 47) & energ==.) 210 | replace energ= 209.2 if ((regnat == 3) & (p601n == 48) & energ==.) 211 | replace energ= 209.2 if ((regnat == 3) & (p601n == 50) & energ==.) 212 | 213 | recode i601b2 i601d2 (.=0) 214 | gen cant601 = ((i601b2+ i601d2)*1000)/(mieperho*30*12) 215 | gen calor601 = (cant601*energ)/100 216 | recode i601c i601e (.=0) 217 | gen gast601 = (i601c+ i601e)/(mieperho*30*12) 218 | 219 | collapse (sum) cant601 calor601 gast601, by (identh) 220 | save tempo601_2011, replace 221 | 222 | use tempo_sumaria-2011.dta, clear 223 | merge 1:1 identh using tempo601_2011 224 | tab _m 225 | drop _m 226 | recode cant601 calor601 gast601 (.=0) 227 | save, replace 228 | 229 | 230 | **en capítulo 602 A PARTIR DE LA ENAHO 2007 231 | use enaho01-2011-602.dta, clear 232 | rename a*o anio 233 | gen str13 identh=anio+conglome+vivienda+hogar 234 | order identh 235 | sort identh 236 | merge identh using tempo2011 237 | tab _m 238 | drop if _m==2 239 | drop _m 240 | keep if p602==1 241 | 242 | *Asignacion de calorias a cada item de consumo de alimentos dentro del hogar 243 | *obtenido de instituciones beneficas con archivo 6 244 | gen ragrpcd=((p602a * p602b/(7*mieperho)) * 27.5) if p602n==1 245 | replace ragrpcd=((p602a * p602b/(7*mieperho)) * 160) if p602n==2 246 | replace ragrpcd=((p602a * p602b/(7*mieperho)) * 600) if p602n==3 247 | replace ragrpcd=((p602a * p602b/(7*mieperho)) * 600) if p602n==4 248 | replace ragrpcd=((p602a * p602b/(7*mieperho)) * 160) if p602n==5 249 | replace ragrpcd=((p602a * p602b/(7*mieperho)) * 600) if p602n==6 250 | replace ragrpcd=((p602a * p602b/(7*mieperho)) * 600) if p602n==7 251 | recode ragrpcd (.=0) 252 | 253 | gen caloria=(ragrpcd * 484.0)/100 if p602n==1 254 | replace caloria=(ragrpcd * 223.6)/100 if p602n==2 255 | replace caloria=(ragrpcd * 196.3)/100 if p602n==3 256 | replace caloria=(ragrpcd * 196.3)/100 if p602n==4 257 | replace caloria=(ragrpcd * 223.6)/100 if p602n==5 258 | replace caloria=(ragrpcd * 196.3)/100 if p602n==6 259 | replace caloria=(ragrpcd * 196.3)/100 if p602n==7 260 | recode caloria i602e1 (.=0) 261 | 262 | gen gast602=i602e1/(mieperho*30*12) 263 | 264 | gen prexrac=i602e1/(51.96*p602a*p602b) if p602==1 & (p602da==1 | p602db==2) 265 | gen precxkg=prexrac*1000/27.5 if p602n==1 266 | replace precxkg=prexrac*1000/160 if p602n==2 | p602==5 267 | replace precxkg=prexrac*1000/600 if p602n==3 | p602n==4 | p602n==6 | p602n==7 268 | 269 | *agregar (sumar sin ponderar) la variable caloria a nivel de hogar 270 | *obtención de precios por kilogramo de lo consumido de instituciones beneficas. 271 | 272 | rename (ragrpcd caloria) (cant602 calor602) 273 | collapse (sum) cant602 calor602 gast602, by (identh) 274 | save tempo602_2011, replace 275 | 276 | use tempo_sumaria-2011.dta, clear 277 | sort identh 278 | merge 1:1 identh using tempo602_2011 279 | tab _m 280 | drop _m 281 | recode cant602- gast602 (.=0) 282 | save, replace 283 | 284 | 285 | *En capítulo 602A A PARTIR DE LA ENAHO 2007 286 | use enaho01-2011-602a.dta, clear 287 | rename a*o anio 288 | gen str13 identh=anio+conglome+vivienda+hogar 289 | order identh 290 | sort identh 291 | merge identh using tempo2011 292 | tab _m 293 | drop if _m==2 294 | drop _m 295 | keep if p6021==1 296 | 297 | *Asignacion de calorias a cada item de consumo de alimentos dentro del hogar 298 | *obtenido de instituciones beneficas con archivo 6 299 | gen ragrpcd=((p602a1 * p602b1/(7*mieperho)) * 160) if p602n1==1 300 | replace ragrpcd=((p602a1 * p602b1/(7*mieperho)) * 600) if p602n1==2 301 | replace ragrpcd=((p602a1 * p602b1/(7*mieperho)) * 600) if p602n1==3 302 | replace ragrpcd=((p602a1 * p602b1/(7*mieperho)) * 600) if p602n1==4 303 | replace ragrpcd=((p602a1 * p602b1/(7*mieperho)) * 600) if p602n1==5 304 | recode ragrpcd (.=0) 305 | 306 | gen caloria=(ragrpcd * 223.6)/100 if p602n1==1 307 | replace caloria=(ragrpcd * 196.3)/100 if p602n1==2 308 | replace caloria=(ragrpcd * 196.3)/100 if p602n1==3 309 | replace caloria=(ragrpcd * 196.3)/100 if p602n1==4 310 | replace caloria=(ragrpcd * 196.3)/100 if p602n1==5 311 | recode caloria i602e3 (.=0) 312 | 313 | gen gast602A1=i602e3/(mieperho*30*12) 314 | 315 | gen prexrac=i602e3/(51.96*p602a1*p602b1) if p6021==1 & (p602d1a==1 | p602d1b==2) 316 | gen precxkg=prexrac*1000/27.5 if p602n1==1 317 | replace precxkg=prexrac*1000/600 if p602n1==2 | p602n1==3 | p602n1==4 | p602n1==5 318 | 319 | *agregar (sumar sin ponderar) la variable caloria1 a nivel de hogar 320 | *obtención de precios por kilogramo de lo consumido de instituciones beneficas. 321 | rename (ragrpcd caloria) (cant602A1 calor602A1) 322 | collapse (sum) cant602A1 calor602A1 gast602A1, by (identh) 323 | save tempo602a_2011, replace 324 | 325 | use tempo_sumaria-2011.dta, clear 326 | sort identh 327 | merge 1:1 identh using tempo602a_2011 328 | tab _m 329 | drop _m 330 | recode cant602A1- gast602A1 (.=0) 331 | save, replace 332 | 333 | 334 | *Asignacion de calorias a cada item de consumo de alimentos fuera del hogar con archivo 3 335 | *A PARTIR DEL 2007 336 | use enaho01a-2011-500.dta, clear 337 | rename a*o anio 338 | gen str13 identh=anio+conglome+vivienda+hogar 339 | gen str15 identi=anio+conglome+vivienda+hogar+codperso 340 | order identh identi 341 | 342 | gen desadg = (p559a_01 /7) * 160 if p559_01==1 343 | gen almudg = (p559a_02 /7) * 600 if p559_02==1 344 | gen cenadg = (p559a_03 /7) * 600 if p559_03==1 345 | gen otrodg1= (p559a_04 /7) * 222 if p559_04==1 346 | gen otrodg2= (p559a_05 /7) * 222 if p559_05==1 347 | gen otrodg3= (p559a_06 /7) * 222 if p559_06==1 348 | gen otrodg4= (p559a_07 /7) * 222 if p559_07==1 349 | gen otrodg5= (p559a_08 /7) * 222 if p559_08==1 350 | gen otrodg6= (p559a_09 /7) * 222 if p559_09==1 351 | gen otrodg7= (p559a_10 /7) * 222 if p559_10==1 352 | 353 | recode desadg almudg cenadg otrodg1 otrodg2 otrodg3 otrodg4 otrodg5 otrodg6 otrodg7 (.=0) 354 | gen caldesa = desadg * 223.6/100 355 | gen calalmu = almudg * 196.3/100 356 | gen calcena = cenadg * 196.3/100 357 | gen calotro1 = otrodg1 * 113.0/100 358 | gen calotro2 = otrodg2 * 113.0/100 359 | gen calotro3 = otrodg3 * 113.0/100 360 | gen calotro4 = otrodg4 * 113.0/100 361 | gen calotro5 = otrodg5 * 113.0/100 362 | gen calotro6 = otrodg6 * 113.0/100 363 | gen calotro7 = otrodg7 * 113.0/100 364 | 365 | gen calor559 = caldesa + calalmu + calcena + calotro1+ calotro2 + calotro3 + calotro4 + calotro5 + calotro6 + calotro7 366 | gen cant559 = desadg+ almudg+ cenadg+ otrodg1+ otrodg2+ otrodg3+ otrodg4+ otrodg5+ otrodg6+ otrodg7 367 | recode i559d1 i559d2 i559d3 i559d41 i559d42 i559d43 i559d44 i559d45 i559d46 i559d47 (.=0) 368 | gen gast559 = (i559d1+ i559d2+ i559d3+ i559d41+ i559d42+ i559d43+ i559d44+ i559d45 +i559d46 +i559d47)/100 369 | 370 | *solo considerar miembren=1 y parentem dif 7 y 8 371 | *agregar (sumar sin ponderar) caloria2 a nivel de hogar y luego ponerlo a percapita (con mieperho). 372 | *obtención de precios por kilogramo de lo consumido fuera del hogar (utilizando inf imp). 373 | 374 | gen pxkg_des = i559d1*1000/(51.96*p559a_01*160) if p559_01 == 1 & p559c_01 == 1 375 | gen pxkg_alm = i559d2*1000/(51.96*p559a_02*600) if p559_02 == 1 & p559c_02 == 1 376 | gen pxkg_cen = i559d3*1000/(51.96*p559a_03*600) if p559_03 == 1 & p559c_03 == 1 377 | gen pxkg_ot1 = i559d41*1000/(51.96*p559a_04*222) if p559_04 == 1 & p559c_04 == 1 378 | gen pxkg_ot2 = i559d42*1000/(51.96*p559a_05*222) if p559_05 == 1 & p559c_05 == 1 379 | gen pxkg_ot3 = i559d43*1000/(51.96*p559a_06*222) if p559_06 == 1 & p559c_06 == 1 380 | gen pxkg_ot4 = i559d44*1000/(51.96*p559a_07*222) if p559_07 == 1 & p559c_07 == 1 381 | gen pxkg_ot5 = i559d45*1000/(51.96*p559a_08*222) if p559_08 == 1 & p559c_08 == 1 382 | gen pxkg_ot6 = i559d46*1000/(51.96*p559a_09*222) if p559_09 == 1 & p559c_09 == 1 383 | gen pxkg_ot7 = i559d47*1000/(51.96*p559a_10*222) if p559_10 == 1 & p559c_10 == 1 384 | 385 | egen p559a_40 = rmean(p559a_04 - p559a_10) 386 | egen i559d40= rsum (i559d41- i559d47) 387 | gen pxkg_otr = i559d40*1000/(51.96*p559a_40*222) if ((p559_05 == 1) & (p559c_05 == 1)) 388 | 389 | keep if p204==1 & (p203!=8 | p203!=9) 390 | collapse (sum) cant559 calor559 gast559, by (identh) 391 | save tempo500_2011, replace 392 | 393 | use tempo_sumaria-2011.dta, clear 394 | sort identh 395 | merge 1:1 identh using tempo500_2011 396 | tab _m 397 | drop _m 398 | recode cant559- gast559 (.=0) 399 | save, replace 400 | 401 | ******************************************************************************** 402 | use tempo_sumaria-2011.dta, clear 403 | sort identh 404 | merge identh using normacalorica2011 405 | drop _m 406 | 407 | *EN ARCHIVO A NIVEL DE HOGARES 408 | replace calor559 = calor559 /mieperho 409 | replace cant559 = cant559 /mieperho 410 | replace gast559 = gast559 /mieperho 411 | 412 | replace reqcalad=reqcalad/mieperho 413 | replace reqcalam=reqcalam/mieperho 414 | 415 | *A partir del 2007 416 | recode calor601 calor602 calor559 calor602A1 cant601 cant602 cant559 cant602A1 /// 417 | gast601 gast602 gast559 gast602A1 (.=0) 418 | 419 | gen caloria = calor601 + calor602 + calor602A1 + calor559 420 | 421 | *Calculo de la incidencia de pobreza 422 | gen pobreC= caloria< reqcalam if caloria!=. 423 | note pobreC: gen pobreC= caloria< reqcalam if caloria!=. 424 | 425 | *Variables geograficas 426 | gen area=estrato<6 if estrato!=. 427 | replace area=1 if dominio==8 & area==0 428 | 429 | gen dominio2=1 if dominio<4 & area==1 430 | replace dominio2=2 if dominio<4 & area==0 431 | replace dominio2=3 if dominio>3 & dominio<7 & area==1 432 | replace dominio2=4 if dominio>3 & dominio<7 & area==0 433 | replace dominio2=5 if dominio==7 & area==1 434 | replace dominio2=6 if dominio==7 & area==0 435 | replace dominio2=7 if dominio==8 436 | 437 | label define dominio2 1 "Costa urbana" 2 "Costa rural" 3 "Sierra urbana" /// 438 | 4 "Sierra rural" 5 "Selva urbana" 6 "Selva rural" 7 "Lima Metrop." 439 | label values dominio2 dominio2 440 | tab dominio2 441 | 442 | gen region=dominio2<3 | dominio2==7 if dominio2!=. 443 | replace region=2 if dominio2>2 & dominio2<5 & dominio2!=. 444 | replace region=3 if dominio2>4 & dominio2<7 & dominio2!=. 445 | label define region 1 "Costa" 2 "Sierra" 3 "Selva" 446 | label values region region 447 | 448 | tab region pobreC [iw=facpob], nofreq row 449 | ***END 450 | --------------------------------------------------------------------------------