Ayudantia 3

Esta es una vista previa del archivo. Inicie sesión para ver el archivo original

* AYUDANTIA 3
* b) 	Descarguemos los modulos sumarias y modulo vivienda de la base de datos ENAHO 2006, 2008 y 2018. 
* 		Haga merge en cada uno de los años y luego haga append con las bases anuales. 
*		Finalmente haga merge a nivel distrital con la base de datos: terremotos.dta. 
*		¿Por que los merge se han hecho de manera distinta?
* Definamos ruta de trabajo
cd "G:\Mi unidad\Cursos y Talleres dictados\UC - Ayudante - Microeconometria\2019-II\Ayudantia 3\Bases/"
* EMPECEMOS UNIENDO AÑO POR AÑO
* 2006
use "enaho01-2006-100", clear
merge 1:1 conglome vivienda hogar using sumaria-2006 
keep if _m==3
drop _m
gen year=2006
save data_2006 , replace
* 2008
use enaho01-2008-100 , clear
merge 1:1 conglome vivienda hogar using sumaria-2008
keep if _m==3
drop _m
gen year=2008
save data_2008 , replace
* 2018
use enaho01-2018-100 , clear
merge 1:1 conglome vivienda hogar using sumaria-2018
keep if _m==3
drop _m
gen year=2018
save data_2018 , replace
* LUEGO UNIMOS LOS AÑOS
use data_2006 , replace
append using data_2008, force
append using data_2018, force
* UNIENDO CON LA BASE DE TRATAMIENTO
merge m:1 ubigeo using terremotos.dta
keep if _merge==3
drop _merge
* c) Haga una Tabla de Estadísticos Descriptivos (N, media, desv. estándar, mínimo y máximo) por años con las siguientes variables: 
* Tasa de pobreza, consumo per cápita (en niveles y en logaritmos) y porcentaje de viviendas con techos de material malo. 
* Ahora repita el análisis pero comparando a los "tratados" con el resto del país. ¿Que estamos observando en los resultados?
* Tasa de pobreza
gen pobre=1 if pobreza==1 | pobreza==2
replace pobre=0 if pobreza==3
* Consumo per cápita (en niveles y en logaritmos)
gen consumo_pc_USD=(gashog2d/mieperho)/3.3
gen log_consumo_pc=log(consumo_pc_USD)
* Porcentaje de viviendas con techos de material malo
tab p103a
gen techo_mala_calidad=0 if p103a==1 | p103a==2 | p103a==3 | p103a==4 | p103a==5
replace techo_mala_calidad=1 if techo_mala_calidad==. & p103a!=.
tabstat pobre log_consumo_pc techo_mala_calidad , by(year) stats(N mean sd min max) col(stats)
tabstat pobre log_consumo_pc techo_mala_calidad if year==2006, by(T) stats(N mean sd min max)
* d) Una la base de datos con controles.dta y ahora replique lo anterior pero quedandose unicamente con los tratados y con los controles. 
* ¿Que diferencia observa con lo anterior? Sea observador.
drop T
 merge m:m ubigeo using controles.dta
keep if _merge==3 
drop _merge
tabstat pobre log_consumo_pc techo_mala_calidad if year==2006, by(T) stats(N mean sd min max)
 
*e) Replique c) y d) pero ahora haciendo un test de medias solo para 2006. 
* Presente los resultados en una Tabla. En la literatura esto se suele llamar "Tabla de Balance". 
* ¿Porque nos importa esta Tabla?
 
 ssc install balancetable
 
 balancetable T pobre log_consumo_pc consumo_pc_USD techo_mala_calidad if year==2006 using "myfile.xls", ctitles("Control group" "Treatment group" "Difference") vce(robust) pval replace 
 
 * f) Usando solo los años 2006 y 2008 cree una variable dummy terremoto si es tratado, cree una variable dummy "post" si es el año 2008. 
 * Haga una Tabla cruzada para esas dos dummies, mostrando el promedio del porcentaje de viviendas de mala calidad en cada una. 
 * ¿Que nos esta indicando este resultado?
gen post=1 if year==2008
replace post=0 if year==2006
table post T , c(mean pobre)
table post T , c(mean log_consumo_pc)
table post T , c(mean techo_mala_calidad)
 
* g) Estime el siguiente modelo:
gen T_post=T*post
 reg pobre T_post T post, r
outreg2 using "reg1.xls" , replace
 reg log_consumo_pc T_post T post, r
outreg2 using "reg1.xls" , append
 reg techo_mala_calidad T_post T post, r
outreg2 using "reg1.xls" , append
 
 
* h) reghdfe
* USANDO AREG
 areg pobre T_post post , vce(robust) absorb(ubigeo)
outreg2 using "reg2.xls" , replace
 areg log_consumo_pc T_post post , vce(robust) absorb(ubigeo)
outreg2 using "reg2.xls" , append
 areg techo_mala_calidad T_post post , vce(robust) absorb(ubigeo)
outreg2 using "reg2.xls" , append
 
 * USANDO reghdfe
 reghdfe log_consumo_pc T_post post , vce(robust) absorb(FE1=ubigeo)