cap log close log using "${LOGS}05_analysis", replace text set more off global NORMALIZATION "sd" global REGULPREF "regul_t0" global GERINGPREF "gering_t0" global ALLPREF "all_t1" global POT2PREF "pot2_t1" global WEIGHTS " " /*"[aweight=svb_122013]"*/ global PREFER "NachML" /*NachML Alt10_NachML */ /******************************************************************************/ /***Footnotze: Rescale Factors ***/ /******************************************************************************/ cap log close set more off use "${DATA}03_krs_panel.dta", clear keep if weiblich==9 sum bite1400_original, det gen bite1400_p90p10=bite1400_original/(r(p90)-r(p10)) gen bite1400_p95p05=bite1400_original/(r(p95)-r(p5)) foreach type in sd minmax 0max p90p10 p95p05 { sum bite1400_`type' if time==624 local meanbite_`type'=r(mean) } scalar rescale_sd_minmax=round( `meanbite_sd'/`meanbite_minmax', 0.001) scalar rescale_sd_0max= round(`meanbite_sd'/`meanbite_0max',0.001) scalar rescale_sd_p90p10= round(`meanbite_sd'/`meanbite_p90p10',0.001) scalar rescale_sd_p95p05= round(`meanbite_sd'/`meanbite_p95p05',0.001) scalar li rescale_sd_0max rescale_sd_minmax rescale_sd_p90p10 rescale_sd_p95p05 /******************************************************************************/ /***TABLE 2 and 3: Comparing Specifications for Employment and Welfare ***/ /******************************************************************************/ /* Static Panel, no timing issues */ cap log close set more off use "${DATA}03_krs_panel.dta", clear keep if weiblich==9 global TREND0 " " /* No trend*/ global TREND1 "c.trend#c.bite1400" global TREND2 "c.trend#c.bite1400 c.trend#c.trend#c.bite1400" global TREND3 "c.logtrend#c.bite1400" forval x=0/3 { /*Loop over trend specifications*/ foreach type in regul gering all pot2 sgbII sgbIII { /* Loop over 4 main outcomes */ #del ; eststo Static_`type'_t`x': xtreg S12.log_`type' NachML_mal_bite1400_${NORMALIZATION} Dt* ${TREND`x'} ${WEIGHTS}, fe vce(cluster krs); estadd scalar R2_within = e(r2_w); estadd scalar R2_between = e(r2_b); estadd scalar R2_overall = e(r2_o); #del cr } } /*Save to tables*/ local EMPLOY Static_regul_t0 Static_regul_t1 Static_regul_t2 Static_regul_t3 Static_gering_t0 Static_gering_t1 Static_gering_t2 Static_gering_t3 local REGUL Static_regul_t0 Static_regul_t1 Static_regul_t2 Static_regul_t3 local GERING Static_gering_t0 Static_gering_t1 Static_gering_t2 Static_gering_t3 local WELFARE Static_all_t0 Static_all_t1 Static_all_t2 Static_all_t3 Static_pot2_t0 Static_pot2_t1 Static_pot2_t2 Static_pot2_t3 local ALL Static_all_t0 Static_all_t1 Static_all_t2 Static_all_t3 local POT2 Static_pot2_t0 Static_pot2_t1 Static_pot2_t2 Static_pot2_t3 local UNEMP Static_sgbIII_t0 Static_sgbIII_t1 Static_sgbIII_t2 Static_sgbIII_t3 Static_sgbII_t0 Static_sgbII_t1 Static_sgbII_t2 Static_sgbII_t3 local ALGI Static_sgbIII_t0 Static_sgbIII_t1 Static_sgbIII_t2 Static_sgbIII_t3 local ALGII Static_sgbII_t0 Static_sgbII_t1 Static_sgbII_t2 Static_sgbII_t3 esttab `REGUL' using "${TABLES}Table2_Employment_Regul.tex", /// label se(3) b(3) replace mtitles(`REGUL') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) esttab `GERING' using "${TABLES}Table2_Employment_Gering.tex", /// label se(3) b(3) replace mtitles(`GERING') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) esttab `ALL' using "${TABLES}Table3_Welfare_All.tex", /// label se(3) b(3) replace mtitles(`ALL') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) esttab `POT2' using "${TABLES}Table3_Welfare_Pot2.tex", /// label se(3) b(3) replace mtitles(`POT2') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) esttab `ALGII' using "${TABLES}Table3_ALGII.tex", /// label se(3) b(3) replace mtitles(`ALGII') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) esttab `ALGI' using "${TABLES}Table3_ALGI.tex", /// label se(3) b(3) replace mtitles(`ALGI') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) /* esttab `EMPLOY' using "${TABLES}Table2_Employment.tex", /// label se(3) b(3) replace mtitles(`EMPLOY') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) esttab `WELFARE' using "${TABLES}Table3_Welfare_${NORMALIZATION}.tex", /// label se(3) b(3) replace mtitles(`WELFARE') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) esttab `UNEMP' using "${TABLES}Table3_Unemplyoment_${NORMALIZATION}.tex", /// label se(3) b(3) replace mtitles(`UNEMP') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) */ /*****************************************************/ /***TABLE 2 and 3: Alternative Spec in levels ***/ /****************************** **********************/ /* Static Panel, no timing issues */ cap log close set more off use "${DATA}03_krs_panel.dta", clear keep if weiblich==9 global TREND0 " " /* No trend*/ global TREND1 "c.trend#c.bite1400" global TREND2 "c.trend#c.bite1400 c.trend#c.trend#c.bite1400" global TREND3 "c.logtrend#c.bite1400" forval x=0/3 { /*Loop over trend specifications*/ foreach type in regul gering all pot2 sgbII sgbIII { /* Loop over 4 main outcomes */ #del ; eststo Level_`type'_t`x': xtreg log_`type' NachML_mal_bite1400_${NORMALIZATION} Dt* ${TREND`x'} ${WEIGHTS}, fe vce(cluster krs); estadd scalar R2_within = e(r2_w); estadd scalar R2_between = e(r2_b); estadd scalar R2_overall = e(r2_o); #del cr } } /*Save to tables*/ local REGUL_Level Level_regul_t0 Level_regul_t1 Level_regul_t2 Level_regul_t3 local GERING_Level Level_gering_t0 Level_gering_t1 Level_gering_t2 Level_gering_t3 local ALL_Level Level_all_t0 Level_all_t1 Level_all_t2 Level_all_t3 local POT2_Level Level_pot2_t0 Level_pot2_t1 Level_pot2_t2 Level_pot2_t3 esttab `REGUL_Level' using "${TABLES}Table2ALT_Employment_Regul_Level.tex", /// label se(3) b(3) replace mtitles(`REGUL_Level') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) esttab `GERING_Level' using "${TABLES}Table2ALT_Employment_Gering_Level.tex", /// label se(3) b(3) replace mtitles(`GERING_Level') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) esttab `ALL_Level' using "${TABLES}Table3ALT_Welfare_All_Level.tex", /// label se(3) b(3) replace mtitles(`ALL_Level') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) esttab `POT2_Level' using "${TABLES}Table3ALT_Welfare_POT2_Level.tex", /// label se(3) b(3) replace mtitles(`POT2_Level') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) /******************************************************************************/ /***TABLE 4: TIMING ISSUES ***/ /******************************************************************************/ cap log close set more off use "${DATA}03_krs_panel.dta", clear keep if weiblich==9 global TREND0 " " /* No trend*/ global TREND1 "c.trend#c.bite1400" global TREND2 "c.trend#c.bite1400 c.trend#c.trend#c.bite1400" global TREND3 "c.logtrend#c.bite1400" forval x=0/1 { foreach type in regul gering all pot2 { #del ; eststo Timing_`type'_t`x': xtreg S12.log_`type' Jul_VorML_mal_bite1400_${NORMALIZATION} Aug_VorML_mal_bite1400_${NORMALIZATION} Sep_VorML_mal_bite1400_${NORMALIZATION} Okt_VorML_mal_bite1400_${NORMALIZATION} Nov_VorML_mal_bite1400_${NORMALIZATION} Dez_VorML_mal_bite1400_${NORMALIZATION} Jan_NachML_mal_bite1400_${NORMALIZATION} Feb_NachML_mal_bite1400_${NORMALIZATION} Mar_NachML_mal_bite1400_${NORMALIZATION} Apr_NachML_mal_bite1400_${NORMALIZATION} Mai_NachML_mal_bite1400_${NORMALIZATION} Jun_NachML_mal_bite1400_${NORMALIZATION} HJ2_NachML_mal_bite1400_${NORMALIZATION} Dt* ${TREND`x'} ${WEIGHTS}, fe vce(cluster krs); #del cr estadd scalar R2_within = e(r2_w) estadd scalar R2_between = e(r2_b) estadd scalar R2_overall = e(r2_o) } } /*Save to tables*/ local TIMING Timing_${REGULPREF} Timing_${GERINGPREF} Timing_${ALLPREF} Timing_${POT2PREF} esttab `TIMING' using "${TABLES}Table4_Timing.tex", /// label se(3) b(3) replace mtitles(`TIMING') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) /******************************************************************************/ /*** TABLE 5a: PREFERRED ESTIMATES MIT NEUEM TIMING ***/ /******************************************************************************/ cap log close set more off use "${DATA}03_krs_panel.dta", clear keep if weiblich==9 global TREND0 " " /* No trend*/ global TREND1 "c.trend#c.bite1400" global TREND2 "c.trend#c.bite1400 c.trend#c.trend#c.bite1400" global TREND3 "c.logtrend#c.bite1400" forval x=0/1 { foreach type in regul gering all pot2 { #del ; eststo Prefer_`type'_t`x': xtreg S12.log_`type' Alt7_NachML_mal_bite1400_${NORMALIZATION} Dt* ${TREND`x'} ${WEIGHTS}, fe vce(cluster krs); #del cr estadd scalar R2_within = e(r2_w) estadd scalar R2_between = e(r2_b) estadd scalar R2_overall = e(r2_o) } } /*Save to tables*/ local PREFER Prefer_${REGULPREF} Prefer_${GERINGPREF} Prefer_${ALLPREF} Prefer_${POT2PREF} esttab `PREFER' using "${TABLES}Table5a_Prefered_neu.tex", /// label se(3) b(3) replace mtitles(`PREFER') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) /******************************************************************************/ /*** TABLE 5b: PREFERRED ESTIMATES MIT ALTEM TIMING ***/ /******************************************************************************/ cap log close set more off use "${DATA}03_krs_panel.dta", clear keep if weiblich==9 global TREND0 " " /* No trend*/ global TREND1 "c.trend#c.bite1400" global TREND2 "c.trend#c.bite1400 c.trend#c.trend#c.bite1400" global TREND3 "c.logtrend#c.bite1400" forval x=0/1 { foreach type in regul gering all pot2 { #del ; eststo Prefer_`type'_t`x': xtreg S12.log_`type' ${PREFER}_mal_bite1400_${NORMALIZATION} Dt* ${TREND`x'} ${WEIGHTS}, fe vce(cluster krs); #del cr estadd scalar R2_within = e(r2_w) estadd scalar R2_between = e(r2_b) estadd scalar R2_overall = e(r2_o) } } /*Save to tables*/ local PREFER Prefer_${REGULPREF} Prefer_${GERINGPREF} Prefer_${ALLPREF} Prefer_${POT2PREF} esttab `PREFER' using "${TABLES}Table5b_Prefered_alt.tex", /// label se(3) b(3) replace mtitles(`PREFER') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) /******************************************************************************/ /*** TABLE 6: AUFSTOCKER COMPOSITION ***/ /******************************************************************************/ cap log close set more off use "${DATA}03_krs_panel.dta", clear keep if weiblich==9 global TREND0 " " /* No trend*/ global TREND1 "c.trend#c.bite1400" global TREND2 "c.trend#c.bite1400 c.trend#c.trend#c.bite1400" global TREND3 "c.logtrend#c.bite1400" /* With trend differential */ forval x=0/1{ foreach type in self mini midi maxi { /* Loop over Aufstocker TYpes */ #del ; eststo Comp_`type'_t: xtreg S12.log_`type' ${PREFER}_mal_bite1400_${NORMALIZATION} Dt* ${TREND`x'} ${WEIGHTS}, fe vce(cluster krs); #del cr estadd scalar R2_within = e(r2_w) estadd scalar R2_between = e(r2_b) estadd scalar R2_overall = e(r2_o) } } /*Save to tables*/ local COMP Comp_self_t Comp_mini_t Comp_midi_t Comp_maxi_t esttab `COMP' using "${TABLES}Table6_Composition.tex", /// label se(3) b(3) replace mtitles(`COMP') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) /******************************************************************************/ /*** TABLE 7: Labor Income ***/ /******************************************************************************/ cap log close set more off use "${DATA}03_krs_panel.dta", clear keep if weiblich==9 & time==671 gen wagegrowth=(avginc2015-avginc2013)/avginc2013 eststo all_bite: reg wagegrowth bite1400_${NORMALIZATION}, robust estadd scalar R2 = e(r2) eststo all_bin: reg wagegrowth treatment2, robust estadd scalar R2 = e(r2) eststo all_level_bite: reg avginc2015 avginc2013 bite1400_${NORMALIZATION}, robust estadd scalar R2 = e(r2) eststo all_level_bin: reg avginc2015 avginc2013 treatment2, robust estadd scalar R2 = e(r2) /* gen wagechange=avginc2015-avginc2013 gen lwagechange=log(avginc2015)-log(avginc2013) gen wagegrowth=(avginc2015-avginc2013)/avginc2013 reg wagechange bite1400_${NORMALIZATION} reg wagegrowth bite1400_${NORMALIZATION} reg lwagechange bite1400_${NORMALIZATION} reg wagechange treatment2 reg wagegrowth treatment2 reg lwagechange treatment2 */ gen wagegrowth2000p=(avgincabove2k2015-avgincabove2k2013)/avgincabove2k2013 eststo p2000_bite: reg wagegrowth2000p bite1400_${NORMALIZATION}, robust estadd scalar R2 = e(r2) eststo p2000_bin: reg wagegrowth2000p treatment2, robust estadd scalar R2 = e(r2) eststo p2000_level_bite: reg avgincabove2k2015 avgincabove2k2013 bite1400_${NORMALIZATION}, robust estadd scalar R2 = e(r2) eststo p2000_level_bin: reg avgincabove2k2015 avgincabove2k2013 treatment2, robust estadd scalar R2 = e(r2) gen wagegrowth15to20=(avginc15to20in2015-avginc15to20in2013)/avginc15to20in2013 eststo f15to20_bite: reg wagegrowth15to20 bite1400_${NORMALIZATION}, robust estadd scalar R2 = e(r2) eststo f15to20_bin: reg wagegrowth15to20 treatment2, robust estadd scalar R2 = e(r2) eststo f15to20_level_bite: reg avginc15to20in2015 avginc15to20in2013 bite1400_${NORMALIZATION}, robust estadd scalar R2 = e(r2) eststo f15to20_level_bin: reg avginc15to20in2015 avginc15to20in2013 treatment2, robust estadd scalar R2 = e(r2) foreach grenze in 1400 1500 2000 { gen wagegrowth`grenze'=(avgincu`grenze'in2015-avgincu`grenze'in2013)/avgincu`grenze'in2013 eststo G`grenze'_bite: reg wagegrowth`grenze' bite1400_${NORMALIZATION} estadd scalar R2 = e(r2) eststo G`grenze'_bin: reg wagegrowth`grenze' treatment2 estadd scalar R2 = e(r2) eststo G`grenze'_level_bite: reg avgincu`grenze'in2015 avgincu`grenze'in2013 bite1400_${NORMALIZATION} estadd scalar R2 = e(r2) eststo G`grenze'_level_bin: reg avgincu`grenze'in2015 avgincu`grenze'in2013 treatment2 estadd scalar R2 = e(r2) /* gen wagechange`grenze'=avgincu`grenze'in2015-avgincu`grenze'in2013 gen lwagechange`grenze'=log(avgincu`grenze'in2015)-log(avgincu`grenze'in2013) gen wagegrowth`grenze'=(avgincu`grenze'in2015-avgincu`grenze'in2013)/avgincu`grenze'in2013 reg wagechange`grenze' bite1400_${NORMALIZATION} reg wagegrowth`grenze' bite1400_${NORMALIZATION} reg lwagechange`grenze' bite1400_${NORMALIZATION} reg wagechange`grenze' treatment2 reg wagegrowth`grenze' treatment2 reg lwagechange`grenze' treatment2 *reg wagechange`grenze' bite1500_${NORMALIZATION} *reg wagegrowth`grenze' bite1500_${NORMALIZATION} *reg wagechange`grenze' bite2000_${NORMALIZATION} *reg wagegrowth`grenze' bite2000_${NORMALIZATION} */ } /*Save to tables*/ local WAGEBITE all_bite G1400_bite G1500_bite f15to20_bite p2000_bite local WAGEBIN all_bin G1400_bin G1500_bin f15to20_bin p2000_bin local WAGEBITE2 G1400_bite G1500_bite f15to20_bite p2000_bite local WAGEBIN2 G1400_bin G1500_bin f15to20_bin p2000_bin local WAGELEVELBITE all_level_bite G1400_level_bite G1500_level_bite f15to20_level_bite p2000_level_bite esttab `WAGEBITE' using "${TABLES}Table7_Wagegrowth_bite.tex", /// label se(3) b(3) replace mtitles(`WAGEBITE') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2) esttab `WAGELEVELBITE' using "${TABLES}Table7_Wages_bite.tex", /// label se(3) b(3) replace mtitles(`WAGEBITE') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2) /* esttab `WAGEBIN' using "${TABLES}Table7_Wages_bin.tex", /// label se(3) b(3) replace mtitles(`WAGEBIN') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2) esttab `WAGEBITE2' using "${TABLES}Table7_Wages_bite2.tex", /// label se(3) b(3) replace mtitles(`WAGEBITE2') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2) esttab `WAGEBIN2' using "${TABLES}Table7_Wages_bin2.tex", /// label se(3) b(3) replace mtitles(`WAGEBIN2') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2) */ /******************************************************************************/ /*** TABLE A2: Placebo Vorjahr ***/ /******************************************************************************/ cap log close set more off use "${DATA}03_krs_panel.dta", clear keep if weiblich==9 global TREND0 " " /* No trend*/ global TREND1 "c.trend#c.bite1400" global TREND2 "c.trend#c.bite1400 c.trend#c.trend#c.bite1400" global TREND3 "c.logtrend#c.bite1400" forval y=1/6 {/* Loop over Placebo timing*/ forval x=0/1 { foreach type in all pot2 regul gering { #del ; eststo Placebo`y'_`type'_t`x': xtreg S12.log_`type' Pla`y'_NachML_mal_bite1400_${NORMALIZATION} Dt* ${TREND`x'} ${WEIGHTS}, fe vce(cluster krs); #del cr estadd scalar R2_within = e(r2_w) estadd scalar R2_between = e(r2_b) estadd scalar R2_overall = e(r2_o) } } } /*Save to tables*/ forval y=1/6 { local PLACEBO`y' Placebo`y'_${REGULPREF} Placebo`y'_${GERINGPREF} Placebo`y'_${ALLPREF} Placebo`y'_${POT2PREF} esttab `PLACEBO`y'' using "${TABLES}TableA2_Placebo`y'.tex", /// label se(3) b(3) replace mtitles(`PLACEBO`y'') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) } /* local PLACEBO12 Placebo1_${REGULPREF} Placebo2_${REGULPREF} Placebo1_${GERINGPREF} Placebo2_${GERINGPREF} /// Placebo1_${ALLPREF} Placebo2_${ALLPREF} Placebo1_${POT2PREF} Placebo2_${POT2PREF} local PLACEBO13 Placebo1_${REGULPREF} Placebo3_${REGULPREF} Placebo1_${GERINGPREF} Placebo3_${GERINGPREF} /// Placebo1_${ALLPREF} Placebo3_${ALLPREF} Placebo1_${POT2PREF} Placebo3_${POT2PREF} esttab `PLACEBO12' using "${TABLES}TableA2_Placebo1_2.tex", /// label se(3) b(3) replace mtitles(`PLACEBO12') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) esttab `PLACEBO13' using "${TABLES}TableA2_Placebo1_3.tex", /// label se(3) b(3) replace mtitles(`PLACEBO13') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) */ /******************************************************************************/ /*** TABLE A3: Different Treatment Definitions - Binary ***/ /******************************************************************************/ cap log close set more off use "${DATA}03_krs_panel.dta", clear keep if weiblich==9 foreach treat of numlist 2 3 { replace ${PREFER}_mal_treated`treat'=${PREFER}_mal_treated`treat'/100 global TREND0 " " /* No trend*/ global TREND1 "c.trend#i.treatment`treat'" global TREND2 "c.trend#i.treatment`treat' c.trend#c.trend#i.treatment`treat'" global TREND3 "c.logtrend#i.treatment`treat'" forval x=0/1 { foreach type in regul gering all pot2 { #del ; eststo Treat`treat'_`type'_t`x': xtreg S12.log_`type' ${PREFER}_mal_treated`treat' Dt* ${TREND`x'} ${WEIGHTS}, fe vce(cluster krs); #del cr estadd scalar R2_within = e(r2_w) estadd scalar R2_between = e(r2_b) estadd scalar R2_overall = e(r2_o) } /*End of Types*/ } /*End of trend `x' */ } /* End of treatment numbers */ foreach treat of numlist 2 3 { /*Save to tables*/ local TREAT`treat' Treat`treat'_${REGULPREF} Treat`treat'_${GERINGPREF} /// Treat`treat'_${ALLPREF} Treat`treat'_${POT2PREF} esttab `TREAT`treat'' using "${TABLES}TableA3_BinaryTreatment`treat'.tex", /// label se(3) b(3) replace mtitles(`TREAT`treat'') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) } /* local BOTHTREATS Treat2_${REGULPREF} Treat3_${REGULPREF} Treat2_${GERINGPREF} Treat3_${GERINGPREF} /// Treat2_${ALLPREF} Treat3_${ALLPREF} Treat2_${POT2PREF} Treat3_${POT2PREF} esttab `BOTHTREATS' using "${TABLES}TableA3_BothBinaryTreatments.tex", /// label se(3) b(3) replace mtitles(`BOTHTREATS') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) */ /******************************************************************************/ /*** TABLE A4: Different Treatment Definitions - bite thresholds ***/ /******************************************************************************/ cap log close set more off use "${DATA}03_krs_panel.dta", clear keep if weiblich==9 foreach grenze of numlist 1500 2000 { global TREND0 " " global TREND1 "c.trend#c.bite`grenze'" global TREND2 "c.trend#c.bite`grenze' c.trend#c.trend#c.bite`grenze'" global TREND3 "c.logtrend#c.bite`grenze'" forval x=0/1 { /* With trend differential */ foreach type in regul gering all pot2 { #del ; eststo Bite`grenze'_`type'_t`x': xtreg S12.log_`type' ${PREFER}_mal_bite`grenze'_${NORMALIZATION} Dt* ${TREND`x'} ${WEIGHTS}, fe vce(cluster krs); #del cr estadd scalar R2_within = e(r2_w) estadd scalar R2_between = e(r2_b) estadd scalar R2_overall = e(r2_o) } } /*Save to tables*/ local BITE`grenze' Bite`grenze'_${REGULPREF} Bite`grenze'_${GERINGPREF} /// Bite`grenze'_${ALLPREF} Bite`grenze'_${POT2PREF} esttab `BITE`grenze'' using "${TABLES}TableA4_Bite`grenze'.tex", /// label se(3) b(3) replace mtitles(`BITE`grenze'') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) } /* End of loop over bite thresholds */ /*local BOTHBITES Bite1500_${REGULPREF} Bite2000_${REGULPREF} Bite1500_${GERINGPREF} Bite2000_${GERINGPREF} /// Bite1500_${ALLPREF} Bite2000_${ALLPREF} Bite1500_${POT2PREF} Bite2000_${POT2PREF} esttab `BOTHBITES' using "${TABLES}TableA4_Bitecomparison.tex", /// label se(3) b(3) replace mtitles(`BOTHBITES') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) */ /******************************************************************************/ /*** TABLE A5: Partial Adjustemnt ***/ /******************************************************************************/ cap log close set more off use "${DATA}03_krs_panel.dta", clear keep if weiblich==9 global TREND0 " " /* No trend*/ global TREND1 "c.trend#c.bite1400" global TREND2 "c.trend#c.bite1400 c.trend#c.trend#c.bite1400" global TREND3 "c.logtrend#c.bite1400" forval x=0/1{ foreach type in regul gering all pot2 { #del ; eststo Dgmm_`type'_t`x': xtabond2 S12.log_`type' l.S12.log_`type' ${PREFER}_mal_bite1400_${NORMALIZATION} Dt* ${TREND`x'} ${WEIGHTS} , gmm(l.S12.log_`type', collapse) iv(${PREFER}_mal_bite1400_${NORMALIZATION} ${TREND`x'} Dt* ) /*${PREFER}_mal_bite1400_${NORMALIZATION} Dt* ${TREND`x'}*/ nolevel cluster(krs) robust; estadd scalar AR1inFD = e(ar1p); estadd scalar AR2inFD = e(ar2p); estadd scalar Sargan = e(sarganp); estadd scalar Hansen = e(hansenp); eststo Sgmm_`type'_t`x': xtabond2 S12.log_`type' l.S12.log_`type' ${PREFER}_mal_bite1400_${NORMALIZATION} Dt* ${TREND`x'} ${WEIGHTS} , gmm(l.S12.log_`type', collapse) iv(${PREFER}_mal_bite1400_${NORMALIZATION} ${TREND`x'} Dt*) /*${PREFER}_mal_bite1400_${NORMALIZATION} Dt* ${TREND`x'}*/ cluster(krs) robust; estadd scalar AR1inFD = e(ar1p); estadd scalar AR2inFD = e(ar2p); estadd scalar Sargan = e(sarganp); estadd scalar Hansen = e(hansenp); #del cr } local DGMM_t`x' Dgmm_regul_t`x' Dgmm_gering_t`x' Dgmm_all_t`x' Dgmm_pot2_t`x' local SGMM_t`x' Sgmm_regul_t`x' Sgmm_gering_t`x' Sgmm_all_t`x' Sgmm_pot2_t`x' /*Save to tables*/ esttab `DGMM_t`x'' using "${TABLES}TableA5_DGMM_t`x'.tex", /// label se(3) b(3) replace mtitles(`DGMM_t`x'') /// star(* 0.10 ** 0.05 *** 0.01) scalars(AR1inFD AR2inFD Sargan Hansen) esttab `SGMM_t`x'' using "${TABLES}TableA5_SGMM_t`x'.tex", /// label se(3) b(3) replace mtitles(`SGMM_t`x'') /// star(* 0.10 ** 0.05 *** 0.01) scalars(AR1inFD AR2inFD Sargan Hansen) } local DGMM Dgmm_${REGULPREF} Dgmm_${GERINGPREF} Dgmm_${ALLPREF} Dgmm_${POT2PREF} local SGMM Sgmm_${REGULPREF} Sgmm_${GERINGPREF} Sgmm_${ALLPREF} Sgmm_${POT2PREF} esttab `DGMM' using "${TABLES}TableA5_DGMM.tex", /// label se(3) b(3) replace mtitles(`DGMM') /// star(* 0.10 ** 0.05 *** 0.01) scalars(AR1inFD AR2inFD Sargan Hansen) esttab `SGMM' using "${TABLES}TableA5_SGMM.tex", /// label se(3) b(3) replace mtitles(`SGMM') /// star(* 0.10 ** 0.05 *** 0.01) scalars(AR1inFD AR2inFD Sargan Hansen) /******************************************************************************/ /*** TABLE A6: West/Ost Split sample ***/ /******************************************************************************/ cap log close set more off use "${DATA}03_krs_panel.dta", clear keep if weiblich==9 global TREND0 " " /* No trend*/ global TREND1 "c.trend#c.bite1400" global TREND2 "c.trend#c.bite1400 c.trend#c.trend#c.bite1400" global TREND3 "c.logtrend#c.bite1400" forval x=0/1 { foreach ost of numlist 0 1 { foreach type in regul gering all pot2 { #del ; eststo Ost`ost'_`type'_t`x': xtreg S12.log_`type' ${PREFER}_mal_bite1400_${NORMALIZATION} Dt* ${TREND`x'} ${WEIGHTS} if ost==`ost' , fe vce(cluster krs); #del cr estadd scalar R2_within = e(r2_w) estadd scalar R2_between = e(r2_b) estadd scalar R2_overall = e(r2_o) } } } /*Save to tables*/ local OST Ost1_${REGULPREF} Ost1_${GERINGPREF} Ost1_${ALLPREF} Ost1_${POT2PREF} local WEST Ost0_${REGULPREF} Ost0_${GERINGPREF} Ost0_${ALLPREF} Ost0_${POT2PREF} local WESTOST Ost0_${REGULPREF} Ost1_${REGULPREF} Ost0_${GERINGPREF} Ost1_${GERINGPREF} /// Ost0_${ALLPREF} Ost1_${ALLPREF} Ost0_${POT2PREF} Ost1_${POT2PREF} esttab `OST' using "${TABLES}TableA6_Ost.tex", /// label se(3) b(3) replace mtitles(`OST') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) esttab `WEST' using "${TABLES}TableA6_West.tex", /// label se(3) b(3) replace mtitles(`WEST') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) /* esttab `WESTOST' using "${TABLES}TableA6_WestOst.tex", /// label se(3) b(3) replace mtitles(`WESTOST') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) */ /******************************************************************************/ /*** TABLE A7: Männer/Frauen/Gender Cells ***/ /******************************************************************************/ cap log close set more off use "${DATA}03_krs_panel.dta", clear global G0 "weiblich==0" global G1 "weiblich==1" global G2 "weiblich!=9" global TREND0 " " /* No trend*/ global TREND1 "c.trend#c.bite1400" global TREND2 "c.trend#c.bite1400 c.trend#c.trend#c.bite1400" global TREND3 "c.logtrend#c.bite1400" forval x=0/1 { forval g=0/2 { foreach type in regul gering all pot2 { #del ; eststo G`g'_`type'_t`x': xtreg S12.log_`type' ${PREFER}_mal_bite1400_${NORMALIZATION} Dt* ${TREND`x'} ${WEIGHTS} if ${G`g'} , fe vce(cluster krs); #del cr estadd scalar R2_within = e(r2_w) estadd scalar R2_between = e(r2_b) estadd scalar R2_overall = e(r2_o) } } foreach type in regul gering all pot2 { #del ; eststo G3_`type'_t`x': xtreg S12.log_`type' ${PREFER}_mal_bite1400_${NORMALIZATION} Dt* ${TREND`x'} ${WEIGHTS} if ${G2} , fe vce(cluster genderkrsfe); #del cr estadd scalar R2_within = e(r2_w) estadd scalar R2_between = e(r2_b) estadd scalar R2_overall = e(r2_o) } } /*Save to tables*/ local MEN G0_${REGULPREF} G0_${GERINGPREF} G0_${ALLPREF} G0_${POT2PREF} local WOMEN G1_${REGULPREF} G1_${GERINGPREF} G1_${ALLPREF} G1_${POT2PREF} local MENWOMEN G0_${REGULPREF} G1_${REGULPREF} G0_${GERINGPREF} G1_${GERINGPREF} /// G0_${ALLPREF} G1_${ALLPREF} G0_${POT2PREF} G1_${POT2PREF} local GENDERCELLS G2_${REGULPREF} G2_${GERINGPREF} G2_${ALLPREF} G2_${POT2PREF} local GENDERCELLS2 G3_${REGULPREF} G3_${GERINGPREF} G3_${ALLPREF} G3_${POT2PREF} esttab `MEN' using "${TABLES}TableA7_Men.tex", /// label se(3) b(3) replace mtitles(`MEN') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) esttab `WOMEN' using "${TABLES}TableA7_Women.tex", /// label se(3) b(3) replace mtitles(`WOMEN') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) /* esttab `MENWOMEN' using "${TABLES}TableA7_MenWomen.tex", /// label se(3) b(3) replace mtitles(`MENWOMEN') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) */ esttab `GENDERCELLS' using "${TABLES}TableA7_Gendercells.tex", /// label se(3) b(3) replace mtitles(`GENDERCELLS') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) esttab `GENDERCELLS2' using "${TABLES}TableA7Alt_Gendercells2.tex", /// label se(3) b(3) replace mtitles(`GENDERCELLS2') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) /******************************************************************************/ /*** TABLE A8: AMR ***/ /******************************************************************************/ cap log close set more off use "${DATA}03_amr_panel.dta", clear keep if weiblich==9 global TREND0 " " /* No trend*/ global TREND1 "c.trend#c.bite1400" global TREND2 "c.trend#c.bite1400 c.trend#c.trend#c.bite1400" global TREND3 "c.logtrend#c.bite1400" forval x=0/1{ foreach type in regul gering all pot2 { #del ; eststo Amr_`type'_t`x': xtreg S12.log_`type' ${PREFER}_mal_bite1400_${NORMALIZATION} Dt* ${TREND`x'} ${WEIGHTS}, fe vce(cluster amr); #del cr estadd scalar R2_within = e(r2_w) estadd scalar R2_between = e(r2_b) estadd scalar R2_overall = e(r2_o) } } /*Save to tables*/ local AMR Amr_${REGULPREF} Amr_${GERINGPREF} Amr_${ALLPREF} Amr_${POT2PREF} esttab `AMR' using "${TABLES}TableA8_AMR.tex", /// label se(3) b(3) replace mtitles(`AMR') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) /******************************************************************************/ /*** TABLE A9: SPATIAL CORRELATION ***/ /******************************************************************************/ /******************************************************************************/ /***TABLE 0: SPATIAL ***/ /******************************************************************************/ /* cap log close set more off forval x=1/3 { use "${DATA}03_krs_panel.dta", clear merge m:1 krs using "${DATA}02_krs_W`x'.dta" mkmat W* if weiblich==0 &time==624, mat(W) keep if weiblich==9 global TREND0 " " /* No trend*/ global TREND1 "c.trend#c.bite1400" global TREND2 "c.trend#c.bite1400 c.trend#c.trend#c.bite1400" global TREND3 "c.logtrend#c.bite1400" #del; xtreg log_regul NachML_mal_bite1400_${NORMALIZATION} Dt*, fe vce(cluster krs) ; xsmle log_regul NachML_mal_bite1400_${NORMALIZATION} , fe type(both) vce(cluster krs) wmat(W) model(sar) iter(15); xsmle log_regul NachML_mal_bite1400_${NORMALIZATION} , fe type(both) vce(cluster krs) wmat(W) model(sdm) iter(15); xsmle log_regul NachML_mal_bite1400_${NORMALIZATION} , fe type(both) vce(cluster krs) emat(W) model(sem) iter(15); xsmle log_regul NachML_mal_bite1400_${NORMALIZATION} , fe type(both) vce(cluster krs) wmat(W) emat(W) model(sac) iter(15); #del cr } /******************************************************************************/ /***TABLE 0: TIMING ***/ /******************************************************************************/ cap log close set more off use "${DATA}03_krs_panel.dta", clear keep if weiblich==9 global TREND0 " " /* No trend*/ global TREND1 "c.trend#c.bite1400" global TREND2 "c.trend#c.bite1400 c.trend#c.trend#c.bite1400" global TREND3 "c.logtrend#c.bite1400" forval x=0/3 { /*Loop over trend specifications*/ local TIMING`x' foreach type in regul gering all pot2 { #del ; eststo Timing_`type'_t`x': xtreg log_`type' /*HJ1_VorML_mal_bite1400_${NORMALIZATION} Jan_VorML_mal_bite1400_${NORMALIZATION} Feb_VorML_mal_bite1400_${NORMALIZATION} Mar_VorML_mal_bite1400_${NORMALIZATION} Apr_VorML_mal_bite1400_${NORMALIZATION} Mai_VorML_mal_bite1400_${NORMALIZATION} Jun_VorML_mal_bite1400_${NORMALIZATION} */ Jul_VorML_mal_bite1400_${NORMALIZATION} Aug_VorML_mal_bite1400_${NORMALIZATION} Sep_VorML_mal_bite1400_${NORMALIZATION} Okt_VorML_mal_bite1400_${NORMALIZATION} Nov_VorML_mal_bite1400_${NORMALIZATION} Dez_VorML_mal_bite1400_${NORMALIZATION} /*HJ2_VorML_mal_bite1400_${NORMALIZATION}*/ Jan_NachML_mal_bite1400_${NORMALIZATION} Feb_NachML_mal_bite1400_${NORMALIZATION} Mar_NachML_mal_bite1400_${NORMALIZATION} Apr_NachML_mal_bite1400_${NORMALIZATION} Mai_NachML_mal_bite1400_${NORMALIZATION} Jun_NachML_mal_bite1400_${NORMALIZATION} HJ2_NachML_mal_bite1400_${NORMALIZATION} Dt* ${TREND`x'} ${WEIGHTS}, fe vce(cluster krs); #del cr estadd scalar R2_within = e(r2_w) estadd scalar R2_between = e(r2_b) estadd scalar R2_overall = e(r2_o) local TIMING`x' `TIMING`x'' Timing_`type'_t`x' } esttab `TIMING`x'' using "${TABLES}Table0_Timing`x'.tex", /// label se(3) b(3) replace mtitles(`TIMING`x'') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) } local TIMINGPREF Timing_${REGULPREF} Timing_${GERINGPREF} Timing_${ALLPREF} Timing_${POT2PREF} esttab `TIMINGPREF' using "${TABLES}Table0_Timingpref.tex", /// label se(3) b(3) replace mtitles(`TIMINGPREF') /// star(* 0.10 ** 0.05 *** 0.01) scalars(R2_within R2_between R2_overall) /******************************************************************************/ /*** TABLE 0: Partial Adjustemnt & Timing ***/ /******************************************************************************/ cap log close set more off use "${DATA}03_krs_panel.dta", clear keep if weiblich==9 global TREND0 " " /* No trend*/ global TREND1 "c.trend#c.bite1400" global TREND2 "c.trend#c.bite1400 c.trend#c.trend#c.bite1400" global TREND3 "c.logtrend#c.bite1400" forval x=0/1{ foreach type in regul gering all pot2 { #del ; eststo Dgmm_`type'_t`x': xtabond2 log_`type' l.log_`type' Jul_VorML_mal_bite1400_${NORMALIZATION} Aug_VorML_mal_bite1400_${NORMALIZATION} Sep_VorML_mal_bite1400_${NORMALIZATION} Okt_VorML_mal_bite1400_${NORMALIZATION} Nov_VorML_mal_bite1400_${NORMALIZATION} Dez_VorML_mal_bite1400_${NORMALIZATION} /*HJ2_VorML_mal_bite1400_${NORMALIZATION}*/ Jan_NachML_mal_bite1400_${NORMALIZATION} Feb_NachML_mal_bite1400_${NORMALIZATION} Mar_NachML_mal_bite1400_${NORMALIZATION} Apr_NachML_mal_bite1400_${NORMALIZATION} Mai_NachML_mal_bite1400_${NORMALIZATION} Jun_NachML_mal_bite1400_${NORMALIZATION} HJ2_NachML_mal_bite1400_${NORMALIZATION} Dt* ${TREND`x'} ${WEIGHTS} , gmm(l.log_`type', collapse) iv(Jul_VorML_mal_bite1400_${NORMALIZATION} Aug_VorML_mal_bite1400_${NORMALIZATION} Sep_VorML_mal_bite1400_${NORMALIZATION} Okt_VorML_mal_bite1400_${NORMALIZATION} Nov_VorML_mal_bite1400_${NORMALIZATION} Dez_VorML_mal_bite1400_${NORMALIZATION} /*HJ2_VorML_mal_bite1400_${NORMALIZATION}*/ Jan_NachML_mal_bite1400_${NORMALIZATION} Feb_NachML_mal_bite1400_${NORMALIZATION} Mar_NachML_mal_bite1400_${NORMALIZATION} Apr_NachML_mal_bite1400_${NORMALIZATION} Mai_NachML_mal_bite1400_${NORMALIZATION} Jun_NachML_mal_bite1400_${NORMALIZATION} HJ2_NachML_mal_bite1400_${NORMALIZATION} Dt* ${TREND`x'}) nolevel cluster(krs) robust; estadd scalar AR1inFD = e(ar1p); estadd scalar AR2inFD = e(ar2p); estadd scalar Sargan = e(sarganp); estadd scalar Hansen = e(hansenp); eststo Sgmm_`type'_t`x': xtabond2 log_`type' l.log_`type' Jul_VorML_mal_bite1400_${NORMALIZATION} Aug_VorML_mal_bite1400_${NORMALIZATION} Sep_VorML_mal_bite1400_${NORMALIZATION} Okt_VorML_mal_bite1400_${NORMALIZATION} Nov_VorML_mal_bite1400_${NORMALIZATION} Dez_VorML_mal_bite1400_${NORMALIZATION} /*HJ2_VorML_mal_bite1400_${NORMALIZATION}*/ Jan_NachML_mal_bite1400_${NORMALIZATION} Feb_NachML_mal_bite1400_${NORMALIZATION} Mar_NachML_mal_bite1400_${NORMALIZATION} Apr_NachML_mal_bite1400_${NORMALIZATION} Mai_NachML_mal_bite1400_${NORMALIZATION} Jun_NachML_mal_bite1400_${NORMALIZATION} HJ2_NachML_mal_bite1400_${NORMALIZATION} Dt* ${TREND`x'} ${WEIGHTS} , gmm(l.log_`type', collapse) iv(Jul_VorML_mal_bite1400_${NORMALIZATION} Aug_VorML_mal_bite1400_${NORMALIZATION} Sep_VorML_mal_bite1400_${NORMALIZATION} Okt_VorML_mal_bite1400_${NORMALIZATION} Nov_VorML_mal_bite1400_${NORMALIZATION} Dez_VorML_mal_bite1400_${NORMALIZATION} /*HJ2_VorML_mal_bite1400_${NORMALIZATION}*/ Jan_NachML_mal_bite1400_${NORMALIZATION} Feb_NachML_mal_bite1400_${NORMALIZATION} Mar_NachML_mal_bite1400_${NORMALIZATION} Apr_NachML_mal_bite1400_${NORMALIZATION} Mai_NachML_mal_bite1400_${NORMALIZATION} Jun_NachML_mal_bite1400_${NORMALIZATION} HJ2_NachML_mal_bite1400_${NORMALIZATION} Dt* ${TREND`x'}) cluster(krs) robust; estadd scalar AR1inFD = e(ar1p); estadd scalar AR2inFD = e(ar2p); estadd scalar Sargan = e(sarganp); estadd scalar Hansen = e(hansenp); #del cr } local DGMM_t`x' Dgmm_regul_t`x' Dgmm_gering_t`x' Dgmm_all_t`x' Dgmm_pot2_t`x' local SGMM_t`x' Sgmm_regul_t`x' Sgmm_gering_t`x' Sgmm_all_t`x' Sgmm_pot2_t`x' /*Save to tables*/ esttab `DGMM_t`x'' using "${TABLES}Table0_DGMM_t`x'.tex", /// label se(3) b(3) replace mtitles(`DGMM_t`x'') /// star(* 0.10 ** 0.05 *** 0.01) scalars(AR1inFD AR2inFD Sargan Hansen) esttab `SGMM_t`x'' using "${TABLES}Table0_SGMM_t`x'.tex", /// label se(3) b(3) replace mtitles(`SGMM_t`x'') /// star(* 0.10 ** 0.05 *** 0.01) scalars(AR1inFD AR2inFD Sargan Hansen) }