-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy path03_tables.do
282 lines (256 loc) · 17.6 KB
/
03_tables.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
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
use "atusall0319.dta"
local summaryx any_waiting waiting_all waiting_all_nz travel_all loinc inc2040 inc4060 inc6075 inc75100 inc100150 inc150p white black native asian pacisl multi otherrc hisp age female male nohsdip hsdip somcol coldeg married hhchild yngchild metro
local timex yr2004 yr2005 yr2006 yr2007 yr2008 yr2009 yr2010 yr2011 yr2012 yr2013 yr2014 yr2015 yr2016 yr2017 yr2018 yr2019 feb mar apr may jun jul aug sep oct nov dec mon tue wed thr fri sat
local income loinc inc2040 inc4060 inc6075 inc75100 inc100150
local demos black native asian pacisl multi otherrc hisp age female
local edu hsdip somcol coldeg unemployed
local family married hhchild yngchild
** Table 1
xtset gestfips
reg any_waiting `income' `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m1
reg any_waiting `income' metro travel_all `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m2
reg any_waiting `income' metro travel_all worktime `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m3
reg any_waiting `income' metro travel_all worktime `family' `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m4
reg any_waiting `income' metro travel_all worktime `family' `demos' `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m5
xtreg any_waiting `income' metro travel_all worktime `demos' `edu' `family' `timex' if any_time == 1, fe cluster(gestfips)
est sto m6
esttab m1 m2 m3 m4 m5 m6 using "output\table1_panelA.tex", replace cell(b(fmt(2)) se p ci) ar2(2) nostar compress nogaps obslast keep(`income') tex label addnote("Robust standard errors clustered at the state-level (in parentheses); *p < .10 **p < .05 ***p < .01.")
xtset gestfips
reg any_waiting `income' `timex' [pweight=tufnwgtp] if any_time == 1, cluster(gestfips)
est sto m1
reg any_waiting `income' metro travel_all `timex' [pweight=tufnwgtp] if any_time == 1, cluster(gestfips)
est sto m2
reg any_waiting `income' metro travel_all worktime `timex' [pweight=tufnwgtp] if any_time == 1, cluster(gestfips)
est sto m3
reg any_waiting `income' metro travel_all worktime `family' `timex' [pweight=tufnwgtp] if any_time == 1, cluster(gestfips)
est sto m4
reg any_waiting `income' metro travel_all worktime `family' `demos' `timex' [pweight=tufnwgtp] if any_time == 1, cluster(gestfips)
est sto m5
xtreg any_waiting `income' metro travel_all worktime `demos' `edu' `family' `timex' if any_time == 1, fe cluster(gestfips)
est sto m6
esttab m1 m2 m3 m4 m5 m6 using "output\table1_panelB.tex", replace cell(b(fmt(2)) se p ci) ar2(2) nostar compress nogaps obslast keep(`income') tex label addnote("Robust standard errors clustered at the state-level (in parentheses); *p < .10 **p < .05 ***p < .01.")
xtset gestfips
reg waiting_all `income' `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m1
reg waiting_all `income' metro travel_all `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m2
reg waiting_all `income' metro travel_all worktime `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m3
reg waiting_all `income' metro travel_all worktime `family' `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m4
reg waiting_all `income' metro travel_all worktime `family' `demos' `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m5
xtreg waiting_all `income' metro travel_all worktime `demos' `edu' `family' `timex', fe cluster(gestfips)
est sto m6
esttab m1 m2 m3 m4 m5 m6 using "output\table1_panelC.tex", replace cell(b(fmt(2)) se p ci) ar2(2) nostar compress nogaps obslast keep(`income') tex label addnote("Robust standard errors clustered at the state-level (in parentheses); *p < .10 **p < .05 ***p < .01.")
xtset gestfips
reg waiting_all_nz `income' `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m1
reg waiting_all_nz `income' metro travel_all `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m2
reg waiting_all_nz `income' metro travel_all worktime `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m3
reg waiting_all_nz `income' metro travel_all worktime `family' `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m4
reg waiting_all_nz `income' metro travel_all worktime `family' `demos' `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m5
xtreg waiting_all_nz `income' metro travel_all worktime `demos' `edu' `family' `timex', fe cluster(gestfips)
est sto m6
esttab m1 m2 m3 m4 m5 m6 using "output\table1_panelC.tex", replace cell(b(fmt(2)) se p ci) ar2(2) nostar compress nogaps obslast keep(`income') tex label addnote("Robust standard errors clustered at the state-level (in parentheses); *p < .10 **p < .05 ***p < .01.")
** Table SI 1
reg any_waiting i.hiloinc metro travel_all worktime `family' `demos' `educ' `timex' [pweight=tufnwgtp] if any_time == 1, cluster(gestfips)
margins i.hiloinc, post
est sto m1
reg any_waiting i.hiloinc metro travel_all worktime `family' `demos' `educ' `timex' [pweight=tufnwgtp], cluster(gestfips)
margins i.hiloinc, post
est sto m2
reg any_time i.hiloinc metro travel_all worktime `family' `demos' `educ' `timex' [pweight=tufnwgtp], cluster(gestfips)
margins i.hiloinc, post
est sto m3
esttab m1 m2 m3 using "output\tableSI1.tex", replace b(2) se(2) ar2(2) star(* 0.10 ** 0.05 *** 0.01) compress nogaps obslast tex label addnote("Robust standard errors clustered at the state-level (in parentheses); *p < .10 **p < .05 ***p < .01.")
** Table SI 2
log using "logs\summary.log", replace
estpost tabstat `summaryx' [aweight=tufnwgtp], statistics(mean sd count) columns(statistics)
est sto sumall
estpost tabstat `summaryx' [aweight=tufnwgtp] if loinc == 1, statistics(mean sd count) columns(statistics)
est sto sumlow
estpost tabstat `summaryx' [aweight=tufnwgtp] if inc150p == 1, statistics(mean sd count) columns(statistics)
est sto sumhigh
estpost tabstat `summaryx' [aweight=tufnwgtp] if loinc == 1 & any_time == 1, statistics(mean sd count) columns(statistics)
est sto sumlowany
estpost tabstat `summaryx' [aweight=tufnwgtp] if inc150p == 1 & any_time == 1, statistics(mean sd count) columns(statistics)
est sto sumhighany
estpost tabstat `summaryx' [aweight=tufnwgtp] if loinc == 1 & waiting_all > 0, statistics(mean sd count) columns(statistics)
est sto sumlownz
estpost tabstat `summaryx' [aweight=tufnwgtp] if inc150p == 1 & waiting_all > 0, statistics(mean sd count) columns(statistics)
est sto sumhighnz
esttab sumall sumlow sumhigh sumlowany sumhighany sumlownz sumhighnz using "output\t1", cell(mean(fmt(2)) sd(par)) label unstack compress obslast tex mtitles("All" "Low income" "High income") addnotes("Standard deviations in parentheses; the statistical significance of mean differences between high and low income respondents is tested using t tests. *p < .10 **p < .05 ***p < .01.") replace
foreach x of varlist `summaryx'{
reg `x' inc150p [pweight=tufnwgtp] if (inc150p == 1 | loinc == 1), r
reg `x' inc150p [pweight=tufnwgtp] if (inc150p == 1 | loinc == 1) & any_time == 1, r
reg `x' inc150p [pweight=tufnwgtp] if (inc150p == 1 | loinc == 1) & waiting_all > 0, r
}
log close
** Table SI 3
xtset gestfips
reg waiting_all_nz `income' `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m1
reg waiting_all_nz `income' metro travel_all `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m2
reg waiting_all_nz `income' metro travel_all worktime `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m3
reg waiting_all_nz `income' metro travel_all worktime `family' `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m4
reg waiting_all_nz `income' metro travel_all worktime `family' `demos' `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m5
xtreg waiting_all_nz `income' metro travel_all worktime `demos' `edu' `family' `timex', fe cluster(gestfips)
est sto m6
esttab m1 m2 m3 m4 m5 m6 using "output\tableSI3.tex", replace cell(b(fmt(2)) se p ci) ar2(2) nostar compress nogaps obslast keep(`income') tex label addnote("Robust standard errors clustered at the state-level (in parentheses); *p < .10 **p < .05 ***p < .01.")
** Table SI 4
xtset gestfips
reg waiting_all `income' `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m1
reg waiting_all `income' metro travel_all `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m2
reg waiting_all `income' metro travel_all worktime `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m3
reg waiting_all `income' metro travel_all worktime `family' `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m4
reg waiting_all `income' metro travel_all worktime `family' `demos' `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m5
xtreg waiting_all `income' metro travel_all worktime `demos' `edu' `family' `timex', fe cluster(gestfips)
est sto m6
esttab m1 m2 m3 m4 m5 m6 using "output\tableSI4.tex", replace cell(b(fmt(2)) se p ci) ar2(2) nostar compress nogaps obslast keep(`income') tex label addnote("Robust standard errors clustered at the state-level (in parentheses); *p < .10 **p < .05 ***p < .01.")
** Table SI 5
xtset gestfips
reg any_waiting `income' `timex' [pweight=tufnwgtp] if any_time == 1, cluster(gestfips)
est sto m1
reg any_waiting `income' metro travel_all `timex' [pweight=tufnwgtp] if any_time == 1, cluster(gestfips)
est sto m2
reg any_waiting `income' metro travel_all worktime `timex' [pweight=tufnwgtp] if any_time == 1, cluster(gestfips)
est sto m3
reg any_waiting `income' metro travel_all worktime `family' `timex' [pweight=tufnwgtp] if any_time == 1, cluster(gestfips)
est sto m4
reg any_waiting `income' metro travel_all worktime `family' `demos' `timex' [pweight=tufnwgtp] if any_time == 1, cluster(gestfips)
est sto m5
xtreg any_waiting `income' metro travel_all worktime `demos' `edu' `family' `timex' if any_time == 1, fe cluster(gestfips)
est sto m6
esttab m1 m2 m3 m4 m5 m6 using "output\tableSI5.tex", replace cell(b(fmt(2)) se p ci) ar2(2) nostar compress nogaps obslast keep(`income') tex label addnote("Robust standard errors clustered at the state-level (in parentheses); *p < .10 **p < .05 ***p < .01.")
** Table SI 6
xtset gestfips
reg any_waiting `income' `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m1
reg any_waiting `income' metro travel_all `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m2
reg any_waiting `income' metro travel_all worktime `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m3
reg any_waiting `income' metro travel_all worktime `family' `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m4
reg any_waiting `income' metro travel_all worktime `family' `demos' `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m5
xtreg any_waiting `income' metro travel_all worktime `demos' `edu' `family' `timex' if any_time == 1, fe cluster(gestfips)
est sto m6
esttab m1 m2 m3 m4 m5 m6 using "output\tableSI6.tex", replace cell(b(fmt(2)) se p ci) ar2(2) nostar compress nogaps obslast keep(`income') tex label addnote("Robust standard errors clustered at the state-level (in parentheses); *p < .10 **p < .05 ***p < .01.")
**Table SI 7
xtset gestfips
reg waiting_all `income' `timex' [pweight=tufnwgtp] if any_time == 1, cluster(gestfips)
est sto m1
reg waiting_all `income' metro travel_all `timex' [pweight=tufnwgtp] if any_time == 1, cluster(gestfips)
est sto m2
reg waiting_all `income' metro travel_all worktime `timex' [pweight=tufnwgtp] if any_time == 1, cluster(gestfips)
est sto m3
reg waiting_all `income' metro travel_all worktime `family' `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m4
reg waiting_all `income' metro travel_all worktime `family' `demos' `timex' [pweight=tufnwgtp] if any_time == 1, cluster(gestfips)
est sto m5
xtreg waiting_all `income' metro travel_all worktime `demos' `edu' `family' `timex' if any_time == 1, fe cluster(gestfips)
est sto m6
esttab m1 m2 m3 m4 m5 m6 using "output\tableSI7.tex", replace cell(b(fmt(2)) se p ci) ar2(2) nostar compress nogaps obslast keep(`income') tex label addnote("Robust standard errors clustered at the state-level (in parentheses); *p < .10 **p < .05 ***p < .01.")
**Table SI 8
logit any_waiting `income' `timex' [pweight=tufnwgtp] if any_time == 1, vce(cluster gestfips)
margins, dydx(`income') post
est sto m1
logit any_waiting `income' metro travel_all `timex' [pweight=tufnwgtp] if any_time == 1, vce(cluster gestfips)
margins, dydx(`income') post
est sto m2
logit any_waiting `income' metro travel_all worktime `timex' [pweight=tufnwgtp] if any_time == 1, vce(cluster gestfips)
margins, dydx(`income') post
est sto m3
logit any_waiting `income' metro travel_all worktime `family' `timex' [pweight=tufnwgtp] if any_time == 1, vce(cluster gestfips)
margins, dydx(`income') post
est sto m4
logit any_waiting `income' metro travel_all worktime `family' `demos' `timex' [pweight=tufnwgtp] if any_time == 1, vce(cluster gestfips)
margins, dydx(`income') post
est sto m5
logit any_waiting `income' metro travel_all worktime `family' `demos' `educ' `timex' [pweight=tufnwgtp] if any_time == 1, vce(cluster gestfips)
margins, dydx(`income') post
est sto m6
esttab m1 m2 m3 m4 m5 m6 using "output\tableSI8.tex", replace cell(b(fmt(2)) se p ci) pr2(2) nostar compress nogaps obslast keep(`income') tex label addnote("Robust standard errors clustered at the state-level (in parentheses); *p < .10 **p < .05 ***p < .01.")
**Table SI 9
reg any_waiting i.race_cat metro unemployed travel_all worktime `family' `timex' [pweight=tufnwgtp] if any_time == 1, cluster(gestfips)
est sto m1
reg any_waiting i.race_cat##i.hiloinc metro unemployed travel_all worktime `family' `timex' [pweight=tufnwgtp] if any_time == 1, cluster(gestfips)
est sto m2
reg waiting_all_nz i.race_cat metro unemployed travel_all worktime `family' `timex' [pweight=tufnwgtp]
est sto m3
reg waiting_all_nz i.race_cat##i.hiloinc metro unemployed travel_all worktime `family' `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto m4
esttab m1 m2 m3 m4 using "output\tableSI9.tex", replace cell(b(fmt(2)) se p ci) pr2(2) nostar compress nogaps obslast tex label addnote("Robust standard errors clustered at the state-level (in parentheses); *p < .10 **p < .05 ***p < .01.")
**Table SI 13
tobit waiting_all `income' `timex' [pweight=tufnwgtp], ll(0) vce(cluster gestfips)
margins, dydx(*) predict(ystar(0,.)) post
est sto m1
tobit waiting_all `income' metro travel_all `timex' [pweight=tufnwgtp], ll(0) vce(cluster gestfips)
margins, dydx(*) predict(ystar(0,.)) post
est sto m2
tobit waiting_all `income' metro travel_all worktime `timex' [pweight=tufnwgtp], ll(0) vce(cluster gestfips)
margins, dydx(*) predict(ystar(0,.)) post
est sto m3
tobit waiting_all `income' metro travel_all worktime `family' `timex' [pweight=tufnwgtp], ll(0) vce(cluster gestfips)
margins, dydx(*) predict(ystar(0,.)) post
est sto m4
tobit waiting_all `income' metro travel_all worktime `family' `demos' `timex' [pweight=tufnwgtp], ll(0) vce(cluster gestfips)
margins, dydx(*) predict(ystar(0,.)) post
est sto m5
tobit waiting_all `income' metro travel_all worktime `family' `demos' `educ' `timex' [pweight=tufnwgtp], ll(0) vce(cluster gestfips)
margins, dydx(*) predict(ystar(0,.)) post
est sto m6
esttab m1 m2 m3 m4 m5 m6 using "output\tableSI13.tex", replace cell(b(fmt(2)) se p ci) pr2(2) nostar compress nogaps obslast keep(`income') tex label addnote("Robust standard errors clustered at the state-level (in parentheses); *p < .10 **p < .05 ***p < .01.")
**Table SI 14
reg waiting_all `income' metro travel_all worktime `family' `demos' `educ' `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto ols
tobit waiting_all `income' metro travel_all worktime `family' `demos' `educ' `timex' [pweight=tufnwgtp], ll(0) vce(cluster gestfips)
margins, dydx(`income') predict(ystar(0,.)) post
est sto tobit1
tobit waiting_all `income' metro travel_all worktime `family' `demos' `educ' `timex' [pweight=tufnwgtp], ll(0) ul(120) vce(cluster gestfips)
margins, dydx(`income') predict(ystar(0,.)) post
est sto tobit2
poisson waiting_all `income' metro travel_all worktime `family' `demos' `educ' `timex' [pweight=tufnwgtp], vce(cluster gestfips)
margins, dydx(`income') post
est sto poisson1
cpoisson waiting_all `income' metro travel_all worktime `family' `demos' `educ' `timex' [pweight=tufnwgtp], ul vce(cluster gestfips)
margins, dydx(`income') post
est sto rcpoisson1
esttab ols tobit1 tobit2 poisson1 poisson2 using "output\tableSI14.tex", replace b(2) se(2) ar2(2) star(* 0.10 ** 0.05 *** 0.01) compress nogaps obslast tex label keep(`income') mtitles("OLS" "Tobit" "2-Way Tobit" "Poisson" "RC Poisson") addnote("Robust standard errors clustered at the state-level (in parentheses); *p < .10 **p < .05 ***p < .01.")
**Table SI15
reg waiting_all_nz `income' metro travel_all worktime `family' `demos' `educ' `timex' [pweight=tufnwgtp], cluster(gestfips)
est sto ols
tobit waiting_all_nz `income' metro travel_all worktime `family' `demos' `educ' `timex' [pweight=tufnwgtp], ll(0) vce(cluster gestfips)
margins, dydx(`income') predict(ystar(0,.)) post
est sto tobit1
tobit waiting_all_nz `income' metro travel_all worktime `family' `demos' `educ' `timex' [pweight=tufnwgtp], ll(0) ul(120) vce(cluster gestfips)
margins, dydx(`income') predict(ystar(0,.)) post
est sto tobit2
poisson waiting_all_nz `income' metro travel_all worktime `family' `demos' `educ' `timex' [pweight=tufnwgtp], vce(cluster gestfips)
margins, dydx(`income') post
est sto poisson1
cpoisson waiting_all_nz `income' metro travel_all worktime `family' `demos' `educ' `timex' [pweight=tufnwgtp], ul vce(cluster gestfips)
margins, dydx(`income') post
est sto rcpoisson1
esttab ols tobit1 tobit2 poisson1 poisson2 using "output\tableSI15.tex", replace b(2) se(2) ar2(2) star(* 0.10 ** 0.05 *** 0.01) compress nogaps obslast tex label keep(`income') mtitles("OLS" "Tobit" "2-Way Tobit" "Poisson" "RC Poisson") addnote("Robust standard errors clustered at the state-level (in parentheses); *p < .10 **p < .05 ***p < .01.")
clear all