*** Parte 2 *** Este programa debe utilizarse con el sofware Stata *** Se pega a la base de hogares las bases de datos generadas anteriormente . *** Si lo que se desea es solo obtener la información de personas y hogares pobres, sin las características de la vivienda, *** se recomienda hacer una base derivada de la base original de hogares que contenga solamente: folio, tamaño del hogar, estrato, y factor . ************************************************************************************************************************************************************************************ clear set mem 400m set more off use "C:\ENIGH\Dataorg\ENIGH2002\Stata\vivi02.dta", clear gen str12 id_hogar=substr(folio,1,12) drop folio rename id_hogar folio keep folio estrato tam_hog factor sort folio save "C:\ENIGH\Comite\Pobreza 2002\Data_Stata\hogares_4var02.dta", replace merge folio using "C:\ENIGH\Comite\Pobreza 2002\Data_Stata\ingmon02def.dta" tab _merge drop _merge sort folio merge folio using "C:\ENIGH\Comite\Pobreza 2002\Data_Stata\gasmon02def.dta" tab _merge drop _merge sort folio merge folio using "C:\ENIGH\Comite\Pobreza 2002\Data_Stata\auto_nomon02def.dta" tab _merge drop _merge sort folio merge folio using "C:\ENIGH\Comite\Pobreza 2002\Data_Stata\esp_nomon02def.dta" tab _merge drop _merge sort folio merge folio using "C:\ENIGH\Comite\Pobreza 2002\Data_Stata\reg_alqnomon02def.dta" tab _merge drop _merge save "C:\ENIGH\Comite\Pobreza 2002\Data_Stata\hog_inggas02", replace summ gasm alimentom alco_tabm vestcalzm viviendam limpiezam enseresm cristalm saludm transpubm transform comunicam edurecrm educabm cuidpersm accpersm otrosm regdadosm gasnomona alimenta alco_taba vestcalza viviendaa limpiezaa enseresa cristala saluda transpuba transfora comunicaa edurecra educaba cuidpersa accpersa otrosa regdadosa gasnomone alimente alco_tabe vestcalze viviendae enserese cristale transpube transfore comunicae edurecre educabe cuidperse accperse otrose regdadose gasnomonr alco_tabr vestcalzr viviendar enseresr cristalr saludr transpubr transforr comunicar edurecrr educabr cuidpersr accpersr otrosr regdadosr est_alq vivrecib vivprest ingmonn perceps otrings regrecibs destring estrato, replace gen rururb=1 if estrato>2 & estrato~=. replace rururb=0 if estrato<=2 *** El valor "0" de la variable rururb representa al estrato urbano, y el valor "1" al estrato rural. egen double gasmon=rsum(alimentom alco_tabm vestcalzm viviendam limpiezam enseresm cristalm saludm transpubm transform comunicam edurecrm cuidpersm accpersm otrosm) egen double autocons=rsum( alimenta alco_taba vestcalza viviendaa limpiezaa enseresa cristala saluda transpuba transfora comunicaa edurecra cuidpersa accpersa otrosa) egen double pago_esp=rsum(alimente alco_tabe vestcalze viviendae limpiezae enserese cristale salude transpube transfore comunicae edurecre cuidperse accperse otrose) egen double reg_esp=rsum(alimentr alco_tabr vestcalzr viviendar limpiezar enseresr cristalr saludr transpubr transforr comunicar edurecrr cuidpersr accpersr otrosr) egen double nomon=rsum(autocons pago_esp reg_esp est_alq) egen double regdados=rsum(regdadosm regdadosa regdadose) gen double regdadosn= -1 * regdados gen double reg_espn = -1 * reg_esp *** ict : ingreso corriente total *** gct : gasto corriente total *** intt: ingreso neto total *** gnt: gasto neto total egen double ict=rsum(ingmonn nomon) egen double gct=rsum(gasmon nomon) egen double intt=rsum(ingmonn nomon regdadosn reg_espn) egen double gnt=rsum(gasmon nomon regdadosn reg_espn) *** Información per cápita gen double ictpc= ict/tam_hog gen double gctpc= gct/tam_hog gen double intpc= intt/tam_hog gen double gntpc= gnt/tam_hog summ gasmon autocons pago_esp reg_esp nomon regdados ict gct intt gnt intpc gntpc ictpc gctpc *** Generacion del factor de expansion para personas sort rururb egen double gcalim=rsum( alimentom alimenta alimente alimentr) egen double gnalim=rsum( alimentom alimenta alimente) gen double gcalimpc= gcalim/tam_hog gen double gnalimpc= gnalim/tam_hog ************************************************************************************************************* *** POBREZA ALIMENTARIA gen pobalim=1 if intpc<=672.25 & rururb==0 replace pobalim=0 if intpc>672.25 & rururb==0 & intpc~=. replace pobalim=1 if intpc<=494.77 & rururb==1 replace pobalim=0 if intpc>494.77 & rururb==1 & intpc~=. *** Numero y porcentaje de hogares en pobreza alimentaria gen double factorp=factor*tam_hog tab pobalim [w=factor] tab pobalim [w=factor] if rururb==0 tab pobalim [w=factor] if rururb==1 *** Numero y porcentaje de personas en pobreza alimentaria tab pobalim [w=factorp] tab pobalim [w=factorp] if rururb==0 tab pobalim [w=factorp] if rururb==1 ************************************************************************************************************* sort folio drop in 1 sort rururb intpc gen fac_acum0=sum(factor) if rururb==0 gen fac_acum1=sum(factor) if rururb==1 mvencode fac_acum0 fac_acum1, mv(0) egen fac_acum=rsum(fac_acum0 fac_acum1) drop fac_acum0 fac_acum1 *** Construcción del decil movil *** Generando el factor de expansión acumulado a nivel hogar,se busca el hogar pivote. Este hogar será aquel que tenga un ingreso neto per capita igual *** o inmediatamente superior a la línea de pobreza alimentaria correspondiente a su estrato. *** Se toma el factor acumulado que se tiene para dicho hogar y a partir de este se genera el decil movil. *** La idea principal del decil movil es dejar en el centro, dentro de un decil, al hogar de referencia. gen double decmovil=1 if (fac_acum>(1340669 - (15662728/20)) & fac_acum<(1340669 + (15662728/20)) & rururb==0) replace decmovil=1 if (fac_acum>(2562346-(8987441/20)) & fac_acum<(2562346+(8987441/20)) & rururb==1) *** Para obtener la pobreza de capacidades y patrimonial se restringe el estudio a los gastos de interés . *** Para pobreza de capacidades los rubros de gasto que deben tomarse en cuenta son: salud, educacion y alimentacion. egen double gasmoncap=rsum(alimentom saludm educabm) egen double autoconcap=rsum(alimenta saluda educaba) egen double pagoespcap=rsum(alimente salude educabe) egen double nomoncap=rsum(autoconcap pagoespcap) egen double gntcap=rsum(gasmoncap nomoncap) *** Para pobreza de patrimonio los rubros de gasto que deben tomarse en cuenta son: salud, educacion, alimentacion, vestido-calzado, vivienda, transporte publico. egen double gasmonpat=rsum(alimentom vestcalzm viviendam saludm transpubm educabm) egen double autoconpat=rsum(alimenta vestcalza viviendaa saluda transpuba educaba) egen double pagoespat=rsum(alimente vestcalze viviendae salude transpube educabe) egen double nomonpat=rsum(autoconpat pagoespat est_alq) egen double gntpat=rsum(gasmonpat nomonpat) sort rururb decmovil by rururb: means gnalim gntcap gntpat [w=factor] if decmovil==1 *** Con las medias obtenidas anteriormente se divide gntcap y gntpat entre gnalim para obtener el inverso del coeficente de Engel *** Esta operación puede realizarse en una hoja de cálculo para obtener los valores por los cuales deben ser multiplicadas las líneas alimentarias y ser derivadas *** las líneas de capacidades y de patrimonio . ********************************************************************************************************* *** POBREZA DE CAPACIDADES gen pobcap=1 if (intpc<= (672.25*1.179) & rururb==0) replace pobcap=0 if (intpc>(672.25*1.179) & rururb==0) & intpc~=. replace pobcap=1 if (intpc<=(494.77*1.187) & rururb==1) replace pobcap=0 if (intpc>(494.77*1.187) & rururb==1) & intpc~=. *** Número y porcentaje de hogares en pobreza de capacidades tab pobcap [w=factor] tab pobcap [w=factor] if rururb==0 tab pobcap [w=factor] if rururb==1 *** Número y porcentaje de personas en pobreza de capacidades tab pobcap [w=factorp] tab pobcap [w=factorp] if rururb==0 tab pobcap [w=factorp] if rururb==1 *** POBREZA DE PATRIMONIO gen pobpatr=1 if (intpc<=(672.25*2.034) & rururb==0) replace pobpatr=0 if (intpc>(672.25*2.034) & rururb==0) & intpc~=. replace pobpatr=1 if (intpc<=(494.77*1.913) & rururb==1) replace pobpatr=0 if (intpc>(494.77*1.913) & rururb==1) & intpc~=. *** Número y porcentaje de hogares en pobreza de patrimonio tab pobpatr [w=factor] tab pobpatr [w=factor] if rururb==0 tab pobpatr [w=factor] if rururb==1 **** Número y porcentaje de personas en pobreza de patrimonio tab pobpatr [w=factorp] tab pobpatr [w=factorp] if rururb==0 tab pobpatr [w=factorp] if rururb==1 save, replace