*==============================================================================
* FILE: 2F-Regressions
* Researcher: Julian B. Adam
* Version: July 2018
*
* Project: Works Councils reform act
*
* Data used: Establishment panel 98-04
*
* Description: This program runs regressions (base specifications (1)-(8)
* for two different dependent variables (number of vol. quits, vol. quit ratio)
*
*==============================================================================

clear
set more off
set linesize 120
cap log close
log using "$log\2F-Regressions.log", replace
use "$data\panel_1_F.dta", clear


* I. BASE SPECIFICATIONS #######################################################


bysort jahr treat_group: egen mean_kuend_an_persistent=mean(kuend_an) if betrratmean==1 & gesamt>150 & gesamt<251
label var mean_kuend_an_persistent "mean(voluntary quits)persistent WoCo"
bysort jahr treat_group: egen N_mean_kuend_an_persistent=count(mean_kuend_an) if betrrat==1 & gesamt>150 & gesamt<251
sum N_mean_kuend_an_persistent
* THIS GRAPH IS FOR 151-250 AND DOES DISPLAY CT
*number of obs. green line
count if treat_group==3& betrratmean==1
bysort jahr: count if treat_group==3 & betrratmean==1
bysort jahr: count if treat_group==4 & betrratmean==1 & gesamt<251
*number of obs. red line
count if treat_group==4 & betrrat==1 & gesamt<251
twoway (line mean_kuend_an_persistent jahr if treat_group==3 & betrrat==1, title("common trend: voluntary quits") lwidth(thick) lcolor(green) legend(label("1" "control group, 151-200"))) ///
(line mean_kuend_an_persistent jahr if treat_group==4 & betrrat==1, lwidth(thick) legend(label("2" "treatment group, 201-250"))) if jahr<2005 &gesamt>150 &gesamt<251 , xline(2001) yscale(range(1.5(0.5)4.5))
graph save $log\ct34_150250, replace asis
graph save $log\ct34_150250-1.png, replace asis
graph export $log\ct34_150250-1pers.tif, replace

*spec. (1)
xtreg kuend_an ib3.treat_group##i.post ib2001.jahr lgeschvol ///
tarif_1 west gesamt i.bran if betrrat==1 & geschart<4 & jahr<2005 ///
& treat_group>2 & treat_group<5 & betrratmean==1 & gesamt>150 & gesamt<251, fe vce(cluster idnum)

*spec. (2)
xtreg kuend_an ib3.treat_group##i.post ib2001.jahr lgeschvol lva_pc rpc_pay_1000 ///
tarif_1 west gesamt i.bran if betrrat==1 & geschart<4 & jahr<2005 ///
& treat_group>2 & treat_group<5 & betrratmean==1 & gesamt>150 & gesamt<251,fe vce(cluster idnum)

*spec. (3) manually
xtreg kuend_an treat4 post ib2001.jahr treatpost treatposttarif lgeschvol ///
tarif_1 west gesamt i.bran if betrrat==1 & geschart<4 & jahr<2005 ///
& treat_group>2 & treat_group<5 & betrratmean==1 & gesamt>150 & gesamt<251, fe vce(cluster idnum)
test (treatpost treatposttarif)

*spec. (4) manually
xtreg kuend_an treat4 post ib2001.jahr treatpost treatposttarif lgeschvol lva_pc rpc_pay_1000 ///
tarif_1 west gesamt i.bran if betrrat==1 & geschart<4 & jahr<2005 ///
& treat_group>2 & treat_group<5 & betrratmean==1 & gesamt>150 & gesamt<251, fe vce(cluster idnum)
test (treatpost treatposttarif)

*spec. (5)
xtreg kuend_an ib3.treat_group##ib2001.jahr lgeschvol ///
tarif_1 west gesamt i.bran if betrrat==1 & geschart<4 & jahr<2005 ///
& treat_group>2 & treat_group<5 & betrratmean==1 & gesamt>150 & gesamt<251,fe vce(cluster idnum)

*spec. (6)
xtreg kuend_an ib3.treat_group##ib2001.jahr lgeschvol lva_pc rpc_pay_1000 ///
tarif_1 west gesamt i.bran if betrrat==1 & geschart<4 & jahr<2005 ///
& treat_group>2 & treat_group<5 & betrratmean==1 & gesamt>150 & gesamt<251,fe vce(cluster idnum)

*spec. (7) manually
xtreg kuend_an treat4 j98 j99 j00 j02 j03 j04 j98treat j99treat j00treat j02treat j03treat ///
j04treat j98treattarif j99treattarif j00treattarif j02treattarif j03treattarif j04treattarif lgeschvol ///
tarif_1 west gesamt i.bran if betrrat==1 & geschart<4 & jahr<2005 ///
& treat_group>2 & treat_group<5 & betrratmean==1 & gesamt>150 & gesamt<251, fe vce(cluster idnum)
test (j00treat j00treattarif)
test (j02treat j02treattarif)
test (j03treat j03treattarif)
test (j04treat j04treattarif)

*spec. (8) manually
xtreg kuend_an treat4 j98 j99 j00 j02 j03 j04 j98treat j99treat j00treat j02treat j03treat ///
j04treat j98treattarif j99treattarif j00treattarif j02treattarif j03treattarif j04treattarif lgeschvol ///
lva_pc rpc_pay_1000 tarif_1 west gesamt i.bran if betrrat==1 & geschart<4 & jahr<2005 ///
& treat_group>2 & treat_group<5 & betrratmean==1 & gesamt>150 & gesamt<251, fe vce(cluster idnum)
test (j00treat j00treattarif)
test (j02treat j02treattarif)
test (j03treat j03treattarif)
test (j04treat j04treattarif)



* DV: quit ratio
*spec. (1)
xtreg kuend_ratio ib3.treat_group##i.post ib2001.jahr lgeschvol ///
tarif_1 west gesamt i.bran if betrrat==1 & geschart<4 & jahr<2005 ///
& treat_group>2 & treat_group<5 & betrratmean==1 & gesamt>150 & gesamt<251, fe vce(cluster idnum)

*spec. (2)
xtreg kuend_ratio ib3.treat_group##i.post ib2001.jahr lgeschvol lva_pc rpc_pay_1000 ///
tarif_1 west gesamt i.bran if betrrat==1 & geschart<4 & jahr<2005 ///
& treat_group>2 & treat_group<5 & betrratmean==1 & gesamt>150 & gesamt<251,fe vce(cluster idnum)

*spec. (3) manually
xtreg kuend_ratio treat4 post ib2001.jahr treatpost treatposttarif lgeschvol ///
tarif_1 west gesamt i.bran if betrrat==1 & geschart<4 & jahr<2005 ///
& treat_group>2 & treat_group<5 & betrratmean==1 & gesamt>150 & gesamt<251, fe vce(cluster idnum)
test (treatpost treatposttarif)

*spec. (4) manually
xtreg kuend_ratio treat4 post ib2001.jahr treatpost treatposttarif lgeschvol lva_pc rpc_pay_1000 ///
tarif_1 west gesamt i.bran if betrrat==1 & geschart<4 & jahr<2005 ///
& treat_group>2 & treat_group<5 & betrratmean==1 & gesamt>150 & gesamt<251, fe vce(cluster idnum)
test (treatpost treatposttarif)

*spec. (5)
xtreg kuend_ratio ib3.treat_group##ib2001.jahr lgeschvol ///
tarif_1 west gesamt i.bran if betrrat==1 & geschart<4 & jahr<2005 ///
& treat_group>2 & treat_group<5 & betrratmean==1 & gesamt>150 & gesamt<251,fe vce(cluster idnum)

*spec. (6)
xtreg kuend_ratio ib3.treat_group##ib2001.jahr lgeschvol lva_pc rpc_pay_1000 ///
tarif_1 west gesamt i.bran if betrrat==1 & geschart<4 & jahr<2005 ///
& treat_group>2 & treat_group<5 & betrratmean==1 & gesamt>150 & gesamt<251,fe vce(cluster idnum)

*spec. (7) manually
xtreg kuend_ratio treat4 j98 j99 j00 j02 j03 j04 j98treat j99treat j00treat j02treat j03treat ///
j04treat j98treattarif j99treattarif j00treattarif j02treattarif j03treattarif j04treattarif lgeschvol ///
tarif_1 west gesamt i.bran if betrrat==1 & geschart<4 & jahr<2005 ///
& treat_group>2 & treat_group<5 & betrratmean==1 & gesamt>150 & gesamt<251, fe vce(cluster idnum)
test (j00treat j00treattarif)
test (j02treat j02treattarif)
test (j03treat j03treattarif)
test (j04treat j04treattarif)

*spec. (8) manually
xtreg kuend_ratio treat4 j98 j99 j00 j02 j03 j04 j98treat j99treat j00treat j02treat j03treat ///
j04treat j98treattarif j99treattarif j00treattarif j02treattarif j03treattarif j04treattarif lgeschvol ///
lva_pc rpc_pay_1000 tarif_1 west gesamt i.bran if betrrat==1 & geschart<4 & jahr<2005 ///
& treat_group>2 & treat_group<5 & betrratmean==1 & gesamt>150 & gesamt<251, fe vce(cluster idnum)
test (j00treat j00treattarif)
test (j02treat j02treattarif)
test (j03treat j03treattarif)
test (j04treat j04treattarif)


log cl