-
Notifications
You must be signed in to change notification settings - Fork 0
/
CovariatesGraphalaMica.do
120 lines (90 loc) · 2.85 KB
/
CovariatesGraphalaMica.do
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
* Buenos Aires.
* Aquí creo un programa pequeño para la extracción.
capture program drop extraccion
program define extraccion
matrix define B=e(tau_cl)
scalar beta=B[1,1]
matrix define right=e(ci_r_cl)
scalar high=right[1,1]
matrix define left=e(ci_l_cl)
scalar low=left[1,1]
end
* Programa de generación de gráfico
capture program drop Graphloco
program define Graphloco
cls
* Aquí vamos a almacenar los resultados para hacer la tabla luego:
cap drop pointest
gen pointest=.
cap drop low95
gen low95=.
cap drop low90
gen low90=.
cap drop high95
gen high95=.
cap drop high90
gen high90=.
cap drop bdw
gen bdw=.
cap drop meses
gen meses=_n
replace meses=. if meses>6
gen empmasde50=(tam_emp==1 | tam_emp==2)
replace empmasde50=. if tam_emp==.
* Hacemos las regresiones para estas covariables:
* Capital Federal, Hombre, Empresa con más de 50 trabajadores, sector agrario
gen primario=(secto1==1|secto2==1|secto3==1)
replace primario=. if letra==""
gen GBA=(province1==1|province2==1|province3==1)
replace GBA=. if provi==.
local j=1
forvalues h=10(10)600 /* Aca itero sobre bdw no covariates GBA sexb2 fnac_anu empmasde50 primario*/ {
* Covariables
rdrobust GBA z, p(1) kernel(uni) level(90) h(`h')
extraccion
* IMPORTANTE:
* El cambio de signo es porque el paquete asume que los de la derecha son los tratados. En nuestro caso es al revés.
replace bdw=`h' in `j'
replace low90=-low in `j'
replace high90=-high in `j'
rdrobust GBA z, p(1) kernel(uni) h(`h')
extraccion
* IMPORTANTE:
* El cambio de signo es porque el paquete asume que los de la derecha son los tratados. En nuestro caso es al revés.
replace pointest=-beta in `j'
replace low95=-low in `j'
replace high95=-high in `j'
local ++j
}
end
* Todos
use "$out/Aumento_2003", clear
append using "$out/Aumento_Enero2004"
append using "$out/Aumento_2004"
append using "$out/Aumento_2005"
append using "$out/Aumento_2006"
append using "$out/Aumento_2007"
append using "$out/Aumento_2008"
append using "$out/Aumento_2011"
cls
* No incluye efecto fijo de persona ni de año:
Graphloco
gen zeroline = 0
local color1 maroon
local color2 0 0 255
twoway rarea low95 high95 bdw, sort lcolor(`color1'*.20) fcolor(`color1') fintensity(20) || ///
rarea low90 high90 bdw, lcolor(`color1'*.40) fcolor(`color1') fintensity(40) || ///
line pointest bdw, sort lpattern(l) lcolor(`color1') || ///
lfit zeroline bdw, lcolor(black) ///
xlabel(,labsize(medium)) ///
xscale(r(9 601)) ///
xlabel(10(50)600, labsize(vsmall) angle(45)) ///
ylabel(,glcolor(none) labsize(medium)) ///
legend(off) ///
title("") ///
ytitle(Lives in Buenos Aires , margin(r=1) size(medium)) xtitle("Bandwidth (ARS)", size(medium)) ///
graphregion(color(white)) plotregion(style(none))
graph export "$out/BuenosAiresVsBandwidth.png", replace
/* yscale(r(-0.05 .5)) ///
ylabel(-.05 0 .15 .25 .35 .5) ///
*/