clear all set more off cd "E:\Outbreak_factors_data\bmjopen-2020-041516.R1 - Research data" use "bmjopen-2020-041516.R1.dta", clear drop if cityc==420100 // Wuhan global outcome "cured_xday" ******************************************************************************** *** Define variables: Connection with Wuhan ******************************************************************************** local connect_begin_date "110" local connect_end_date "124" egen from_wh=rmean(from`connect_begin_date'-from`connect_end_date') egen to_wh =rmean(to`connect_begin_date' -to`connect_end_date') replace from_wh=from_wh*100 la var from_wh "Connection, Wuhan to local" ******************************************************************************** *** Global ******************************************************************************** global outreg2_options "lab keep($control_key) sortvar($control_key) nocons bd(3) sd(3) rd(3) ad(3) adds(Dependent variable mean, r(mean))" global FE "i.procode" local end_date "128 204 211 218 225 303 310 317 324 331" global regsave_option "ci level(90) nose autoid" ******************************************************************************** *** Regressions: Main ******************************************************************************** global control_key "pct_65ab" global control_oth "share_wh15 pop_size emp_rate pct_college municipality" global control "$control_key $control_oth" // sum $control g i=1 foreach date in `end_date' { // qui reg $outcome`date' $control qui reg $outcome`date' $control closed_xday`date' qui sum $outcome`date' if e(sample), mean // if i==1 outreg2 using $outcome.xls, replace $outreg2_options // if i!=1 outreg2 using $outcome.xls, append $outreg2_options if i==1 regsave $control_key using f2&3_data_$outcome, replace $regsave_option if i!=1 regsave $control_key using f2&3_data_$outcome, append $regsave_option qui replace i=i+1 } ******************************************************************************** global control_key "from_wh" global control_oth "pct_65ab pop_size emp_rate pct_college municipality" global control "$control_key $control_oth" // sum $control foreach date in `end_date' { // qui reg $outcome`date' $control qui reg $outcome`date' $control closed_xday`date' qui sum $outcome`date' if e(sample), mean // outreg2 using $outcome.xls, append $outreg2_options regsave $control_key using f2&3_data_$outcome, append $regsave_option qui replace i=i+1 } ******************************************************************************** global control_key "num_bed" global control_oth "pct_65ab share_wh15 pop_size emp_rate pct_college municipality" global control "$control_key $control_oth" // sum $control foreach date in `end_date' { // qui reg $outcome`date' $control qui reg $outcome`date' $control closed_xday`date' qui sum $outcome`date' if e(sample), mean // outreg2 using $outcome.xls, append $outreg2_options regsave $control_key using f2&3_data_$outcome, append $regsave_option } ******************************************************************************** global control_key "gov_index" global control_oth "pct_65ab share_wh15 pop_size emp_rate pct_college municipality" global control "$control_key $control_oth" // sum $control foreach date in `end_date' { // qui reg $outcome`date' $control qui reg $outcome`date' $control closed_xday`date' qui sum $outcome`date' if e(sample), mean // outreg2 using $outcome.xls, append $outreg2_options regsave $control_key using f2&3_data_$outcome, append $regsave_option } ******************************************************************************** *** Regressions: Sub index ******************************************************************************** global control_key "gov_eff" global control "$control_key $control_oth" // sum $control foreach date in `end_date' { // qui reg $outcome`date' $control qui reg $outcome`date' $control closed_xday`date' qui sum $outcome`date' if e(sample), mean // outreg2 using $outcome.xls, append $outreg2_options regsave $control_key using f2&3_data_$outcome, append $regsave_option } global control_key "gov_law" global control "$control_key $control_oth" // sum $control foreach date in `end_date' { // qui reg $outcome`date' $control qui reg $outcome`date' $control closed_xday`date' qui sum $outcome`date' if e(sample), mean // outreg2 using $outcome.xls, append $outreg2_options regsave $control_key using f2&3_data_$outcome, append $regsave_option } global control_key "gov_trans" global control "$control_key $control_oth" // sum $control foreach date in `end_date' { // qui reg $outcome`date' $control qui reg $outcome`date' $control closed_xday`date' qui sum $outcome`date' if e(sample), mean // outreg2 using $outcome.xls, append $outreg2_options regsave $control_key using f2&3_data_$outcome, append $regsave_option } global control_key "gov_size" global control "$control_key $control_oth" // sum $control foreach date in `end_date' { // qui reg $outcome`date' $control qui reg $outcome`date' $control closed_xday`date' qui sum $outcome`date' if e(sample), mean // outreg2 using $outcome.xls, append $outreg2_options regsave $control_key using f2&3_data_$outcome, append $regsave_option }