use "/Users/mohammadalikadivar/Dropbox/PPPDD/Quant/Book_Data.dta", clear *bysort year region: gen region_count=_N foreach var in v2x_polyarchy v2x_liberal v2xcs_ccsi { bysort year region: egen region_`var'=total(`var') replace region_`var'=(region_`var'-`var')/(region_count-1) } label var v2x_polyarchy "Electoral Democracy" label var v2x_liberal "Liberal Element" label var v2xcs_ccsi "Civil Society" label var region_v2x_polyarchy "Electoral Democracy (regional)" label var region_v2x_liberal "Liberal Element (regional)" label var region_v2xcs_ccsi "Civil Society (regional)" stset gwf_duration, id(gwf_dem_id) failure(gwf_tta) stcox nv_regime_dur v_regime_dur past_dem GDPpc_log_1 GDPpcGrowth pop_log if gwf_dem==1 & gwf_st>1950, vce(cluster ccode) est store cox_1 stcox nv_regime_dur v_regime_dur exclpop ethfrac gwf_pres gwf_mixed past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence region_gwf_dem if gwf_dem==1 & gwf_st>1950 , vce(cluster ccode) est store cox_2 stcurve, cumhaz at1(nv_regime_dur=0) at2(nv_regime_dur=1) at3(nv_regime_dur=3) at4(nv_regime_dur=6) /// legend(label(1 "{stSerif:Years of mobilization=0}" )) /// legend(label(2 "{stSerif:Years of mobilization=1}" )) /// legend(label(3 "{stSerif:Years of mobilization=3}" )) /// legend(label(4 "{stSerif:Years of mobilization=6}" )) /// ytitle("{stSerif:Cumulative hazard of democratic breakdown}") xtitle("{stSerif:Analysis time in years}") range(0 30) title("{stSerif:Cox Proportional Hazards Regression}") /// note ("{stSerif:Based on Table 5, Model 2}") /// xlabel(0 "{stSerif:0}" 10 "{stSerif:10}" 20 "{stSerif:20}" 30 "{stSerif:30}") /// ylabel(0 "{stSerif:0}" .2 "{stSerif:.2}" .4 "{stSerif:.4}" .6 "{stSerif:.6}" , labsize(small) angle(horizontal)) graphregion(color(white)) bgcolor(white) graph export cox_curve.png, width(1600) replace margins, at (nv_regime_dur = ( 0 1 3 6 )) level(90) stcox nv_regime_dur v_regime_dur exclpop ethfrac gwf_pres gwf_mixed past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence region_gwf_dem if gwf_dem==1 & gwf_st>1950 & nv_regime_dur>0, vce(cluster ccode) est store cox_3 stcox nv_regime v_regime gwf_pres gwf_mixed past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth exclpop ethfrac /// pop_log prior_mil prior_pers prior_part post_independence region_gwf_dem if gwf_dem==1 & gwf_st>1950 , vce(cluster ccode) est store cox_4 esttab cox_* using "cox_tables.rtf", aic bic se star(* .05 ** .01) order(nv_regime_dur v_regime_dur nv_regime v_regime logoil_gas_valuePOP_1 prior_mil prior_personal prior_party post_independence gwf_presidential gwf_mixed region_gwf_dem ethfrac exclpop GDPpc_log_1 GDPpcGrowth pop_log_1 ) title("Table 4. Contentious Transition Legacy and Democratic Breakdown as of `c(current_time)' `c(current_date)'") label replace *Growth models xtset ccode year xtreg v2x_polyarchy l1.v2x_polyarchy nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_polyarchy i.year i.region if gwf_dem==1 & gwf_st>1950 , vce(cluster ccode) est store growth_1 xtreg v2x_liberal l1.v2x_liberal nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_liberal i.year i.region if gwf_dem==1 & gwf_st>1950 , vce(cluster ccode) est store growth_2 xtreg v2xcs_ccsi l1.v2xcs_ccsi nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2xcs_ccsi i.year i.region if gwf_dem==1 & gwf_st>1950 , vce(cluster ccode) est store growth_3 esttab growth_* using "growth_models.rtf", pr2 se star(+ .1 * .05 ** .01) /// order(nv_regime_dur v_regime_dur past_dem_age logoil_gas_valuePOP_1 prior_mil prior_personal prior_party /// post_independence GDPpc_log_1 GDPpcGrowth exclpop ethfrac pop_log_1 ) drop (*year *region) /// title("Table 4. Contentious Transition Legacy and Democratic Growth as of `c(current_time)' `c(current_date)'") label replace *** Duration as DV nbreg nv_regime_dur pretrans10_cl past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth pop_log /// prior_mil prior_pers prior_part post_independence region_gwf_dem exclpop ethfrac i.region if gwf_dem==1 & gwf_st>1950 & gwf_duration==1 , vce(cluster ccode) est store nb_1 esttab nb_1 using "nb_model.rtf", pr2 se star(+ .1 * .05 ** .01) /// order( pretrans10_cl past_dem_age prior_mil prior_personal prior_party /// post_independence GDPpc_log_1 GDPpcGrowth logoil_gas_valuePOP_1 exclpop ethfrac pop_log_1 ) drop (*region) /// title("Table 5. The Duration of Unarmed Mobilization as of `c(current_time)' `c(current_date)'") label replace ***Components of electoral and liberal foreach var in v2x_frassoc_thick v2x_freexp_altinf v2x_suffr v2x_elecoff { xtreg `var' l1.`var' nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_polyarchy i.year i.region if gwf_dem==1 & gwf_st>1950 , fe } foreach var in v2xcl_rol v2x_jucon v2xlg_legcon { xtreg `var' l1.`var' nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_liberal i.year i.region if gwf_dem==1 & gwf_st>1950 , fe } *** mediation effect gsem (gwf_duration <- nv_regime_dur v2x_polyarchy , family(weibull, failure(gwf_tta))) (v2x_polyarchy <- nv_regime_dur ) matrix b=e(b) gen v2x_polyarchy_1=l1.v2x_polyarchy gen v2xcs_ccsi_1=l1.v2xcs_ccsi gsem (gwf_duration <- nv_regime_dur v2x_polyarchy_1 pop_log gdpcap_log v_regime_dur /// exclpop ethfrac past_dem logoil_gas_valuePOP_1 growthrate pop_log prior_mil prior_pers prior_part post_independence region_gwf_dem, /// family(weibull, failure(gwf_tta))) (v2x_polyarchy_1 <- nv_regime_dur pop_log gdpcap_log v_regime_dur exclpop ethfrac /// past_dem logoil_gas_valuePOP_1 growthrate pop_log prior_mil prior_pers prior_part post_independence l2.region_v2x_polyarchy),from(b, skip) nlcom _b[gwf_duration:v2x_polyarchy_1]*_b[v2x_polyarchy_1:nv_regime_dur] gsem (gwf_duration <- nv_regime_dur v2x_liberal pop_log gdpcap_log v_regime_dur /// exclpop ethfrac past_dem logoil_gas_valuePOP_1 growthrate pop_log prior_mil prior_pers prior_part post_independence region_gwf_dem, /// family(weibull, failure(gwf_tta))) (v2x_liberal <- nv_regime_dur pop_log gdpcap_log v_regime_dur exclpop ethfrac /// past_dem logoil_gas_valuePOP_1 growthrate pop_log prior_mil prior_pers prior_part post_independence l2.region_v2x_liberal),from(b, skip) nlcom _b[gwf_duration:v2x_liberal]*_b[v2x_liberal:nv_regime_dur] gsem (gwf_duration <- nv_regime_dur v2xcs_ccsi_1 pop_log gdpcap_log v_regime_dur /// exclpop ethfrac past_dem logoil_gas_valuePOP_1 growthrate pop_log prior_mil prior_pers prior_part post_independence region_gwf_dem, /// family(weibull, failure(gwf_tta))) (v2xcs_ccsi_1 <- nv_regime_dur l1.v2x_polyarchy_1 pop_log gdpcap_log v_regime_dur exclpop ethfrac /// past_dem logoil_gas_valuePOP_1 growthrate pop_log prior_mil prior_pers prior_part post_independence l1.region_v2xcs_ccsi),from(b, skip) nlcom _b[gwf_duration:v2xcs_ccsi_1]*_b[v2xcs_ccsi_1:nv_regime_dur] *** Growth robustness ** log version xtreg v2x_polyarchy l1.v2x_polyarchy log_nv_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_polyarchy i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_1 xtreg v2x_liberal l1.v2x_liberal log_nv_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_liberal i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_2 xtreg v2xcs_ccsi l1.v2xcs_ccsi log_nv_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2xcs_ccsi i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_3 *** drop biggest xtreg v2x_polyarchy l1.v2x_polyarchy nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_polyarchy i.year i.region if gwf_dem==1 & gwf_st>1950 & nv_regime_dur<9 , fe est store growth_1 xtreg v2x_liberal l1.v2x_liberal nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_liberal i.year i.region if gwf_dem==1 & gwf_st>1950 & nv_regime_dur<9 , fe est store growth_2 xtreg v2xcs_ccsi l1.v2xcs_ccsi nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2xcs_ccsi i.year i.region if gwf_dem==1 & gwf_st>1950 & nv_regime_dur<9 , fe est store growth_3 *** change in sample size xtreg v2x_polyarchy l1.v2x_polyarchy nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_polyarchy i.year i.region if gwf_dem==1 & gwf_st>1960 , vce(cluster ccode) est store growth_1 xtreg v2x_liberal l1.v2x_liberal nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_liberal i.year i.region if gwf_dem==1 & gwf_st>1960 , vce(cluster ccode) est store growth_2 xtreg v2xcs_ccsi l1.v2xcs_ccsi nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2xcs_ccsi i.year i.region if gwf_dem==1 & gwf_st>1960 , vce(cluster ccode) est store growth_3 xtreg v2x_polyarchy l1.v2x_polyarchy nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_polyarchy i.year i.region if gwf_dem==1 & gwf_st>1945 , vce(cluster ccode) est store growth_1 xtreg v2x_liberal l1.v2x_liberal nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_liberal i.year i.region if gwf_dem==1 & gwf_st>1945 , vce(cluster ccode) est store growth_2 xtreg v2xcs_ccsi l1.v2xcs_ccsi nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2xcs_ccsi i.year i.region if gwf_dem==1 & gwf_st>1945 , vce(cluster ccode) xtreg v2x_polyarchy l1.v2x_polyarchy nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_polyarchy i.year i.region if gwf_dem==1 & gwf_st>1950 & gwf_st<2001, vce(cluster ccode) est store growth_1 xtreg v2x_liberal l1.v2x_liberal nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_liberal i.year i.region if gwf_dem==1 & gwf_st>1950 & gwf_st<2001, vce(cluster ccode) est store growth_2 xtreg v2xcs_ccsi l1.v2xcs_ccsi nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2xcs_ccsi i.year i.region if gwf_dem==1 & gwf_st>1950 & gwf_st<2001, vce(cluster ccode) *** cold war xtreg v2x_polyarchy l1.v2x_polyarchy nv_regime_dur v_regime_dur cold_war exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_polyarchy i.year i.region if gwf_dem==1 & gwf_st>1950 , vce(cluster ccode) est store growth_1 xtreg v2x_liberal l1.v2x_liberal nv_regime_dur v_regime_dur cold_war exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_liberal i.year i.region if gwf_dem==1 & gwf_st>1950 , vce(cluster ccode) est store growth_2 xtreg v2xcs_ccsi l1.v2xcs_ccsi nv_regime_dur v_regime_dur cold_war exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2xcs_ccsi i.year i.region if gwf_dem==1 & gwf_st>1950 , vce(cluster ccode) est store growth_3 *** imperial past xtreg v2x_polyarchy l1.v2x_polyarchy nv_regime_dur v_regime_dur pimppast prior_milper prior_milex exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_polyarchy i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_1 xtreg v2x_liberal l1.v2x_liberal nv_regime_dur v_regime_dur pimppast prior_milper prior_milex exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_liberal i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_2 xtreg v2xcs_ccsi l1.v2xcs_ccsi nv_regime_dur v_regime_dur pimppast prior_milper prior_milex exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2xcs_ccsi i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_3 *** trade Russia_Trade_percent_1 China_Trade_percent_1 EU_Trade_percent_1 US_Trade_percent_1 xtreg v2x_polyarchy l1.v2x_polyarchy nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_polyarchy i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_1 xtreg v2x_liberal l1.v2x_liberal nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_liberal i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_2 xtreg v2xcs_ccsi l1.v2xcs_ccsi nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2xcs_ccsi i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_3 *** armed conflict ongoingwarl xtreg v2x_polyarchy l1.v2x_polyarchy nv_regime_dur v_regime_dur ongoingwarl exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_polyarchy i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_1 xtreg v2x_liberal l1.v2x_liberal nv_regime_dur v_regime_dur ongoingwarl exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_liberal i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_2 xtreg v2xcs_ccsi l1.v2xcs_ccsi nv_regime_dur v_regime_dur ongoingwarl exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2xcs_ccsi i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_3 *** navco navco_nv_regime_dur navco_v_regime_dur xtreg v2x_polyarchy l1.v2x_polyarchy navco_nv_regime_dur navco_v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_polyarchy i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_1 xtreg v2x_liberal l1.v2x_liberal navco_nv_regime_dur navco_v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_liberal i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_2 xtreg v2xcs_ccsi l1.v2xcs_ccsi navco_nv_regime_dur navco_v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2xcs_ccsi i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_3 *** pre-tran vars *** aid aid_log_1 xtreg v2x_polyarchy l1.v2x_polyarchy nv_regime_dur v_regime_dur aid_log_1 exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_polyarchy i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_1 xtreg v2x_liberal l1.v2x_liberal nv_regime_dur v_regime_dur aid_log_1 exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_liberal i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_2 xtreg v2xcs_ccsi l1.v2xcs_ccsi nv_regime_dur v_regime_dur aid_log_1 exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2xcs_ccsi i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_3 *** EU_2 xtreg v2x_polyarchy l1.v2x_polyarchy alt_nv_regime_dur v_regime_dur EU_2 exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_polyarchy i.year i.region if gwf_dem==1 & gwf_st>1950 , vce(cluster ccode) est store growth_1 xtreg v2x_liberal l1.v2x_liberal alt_nv_regime_dur v_regime_dur EU_2 exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_liberal i.year i.region if gwf_dem==1 & gwf_st>1950 , vce(cluster ccode) est store growth_2 xtreg v2xcs_ccsi l1.v2xcs_ccsi alt_nv_regime_dur v_regime_dur EU_2 exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2xcs_ccsi i.year i.region if gwf_dem==1 & gwf_st>1950 , vce(cluster ccode) est store growth_3 *** alternative campaign with abeyance alt_nv_regime_dur xtreg v2x_polyarchy l1.v2x_polyarchy alt_nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_polyarchy i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_1 xtreg v2x_liberal l1.v2x_liberal alt_nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_liberal i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_2 xtreg v2xcs_ccsi l1.v2xcs_ccsi alt_nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2xcs_ccsi i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_3 *** gov rev/GDP xtreg v2x_polyarchy l1.v2x_polyarchy nv_regime_dur v_regime_dur exclpop ethfrac IMFrev past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_polyarchy i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_1 xtreg v2x_liberal l1.v2x_liberal nv_regime_dur v_regime_dur IMFrev exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_liberal i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_2 xtreg v2xcs_ccsi l1.v2xcs_ccsi nv_regime_dur v_regime_dur IMFrev exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2xcs_ccsi i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_3 *** colonies replace colruler=0 if colruler==. xtreg v2x_polyarchy l1.v2x_polyarchy nv_regime_dur v_regime_dur exclpop ethfrac i.colruler past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_polyarchy i.year i.region if gwf_dem==1 & gwf_st>1950 , vce(cluster ccode) est store growth_1 xtreg v2x_liberal l1.v2x_liberal nv_regime_dur v_regime_dur exclpop ethfrac i.colruler past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_liberal i.year i.region if gwf_dem==1 & gwf_st>1950 , vce(cluster ccode) est store growth_2 xtreg v2xcs_ccsi l1.v2xcs_ccsi nv_regime_dur v_regime_dur exclpop ethfrac i.colruler past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2xcs_ccsi i.year i.region if gwf_dem==1 & gwf_st>1950 , vce(cluster ccode) est store growth_3 *** Four lags xtreg v2x_polyarchy l1.v2x_polyarchy l2.v2x_polyarchy l3.v2x_polyarchy l4.v2x_polyarchy nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_polyarchy i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_1 xtreg v2x_liberal l1.v2x_liberal l2.v2x_liberal l3.v2x_liberal l4.v2x_liberal nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_liberal i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_2 xtreg v2xcs_ccsi l1.v2xcs_ccsi l2.v2xcs_ccsi l3.v2xcs_ccsi l4.v2xcs_ccsi nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2xcs_ccsi i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_3 *** Exclude previously democratic democracies xtreg v2x_polyarchy l1.v2x_polyarchy nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_polyarchy i.year i.region if gwf_dem==1 & gwf_st>1950 & past_dem==0, vce(cluster ccode) est store growth_1 xtreg v2x_liberal l1.v2x_liberal nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_liberal i.year i.region if gwf_dem==1 & gwf_st>1950 & past_dem==0, vce(cluster ccode) est store growth_2 xtreg v2xcs_ccsi l1.v2xcs_ccsi nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2xcs_ccsi i.year i.region if gwf_dem==1 & gwf_st>1950 & past_dem==0, vce(cluster ccode) *** Fixed effect xtreg v2x_polyarchy nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_polyarchy i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_1 xtreg v2x_liberal nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2x_liberal i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_2 xtreg v2xcs_ccsi nv_regime_dur v_regime_dur exclpop ethfrac past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth /// pop_log prior_mil prior_pers prior_part post_independence l1.region_v2xcs_ccsi i.year i.region if gwf_dem==1 & gwf_st>1950 , fe est store growth_3 *** Summary table preserve tempname outfile postfile `outfile' str20 variable mean median p5 p95 sd sd_w sd_b using summary_stats.dta, replace foreach var of varlist nv_regime v_regime gwf_pres gwf_mixed past_dem logoil_gas_valuePOP_1 GDPpc_log_1 GDPpcGrowth exclpop ethfrac /// pop_log prior_mil prior_pers prior_part post_independence region_gwf_dem v2x_polyarchy v2x_liberal v2xcs_ccsi region_v2x_polyarchy region_v2x_liberal region_v2xcs_ccsi /// { local vlabel: variable label `var' su `var' if gwf_dem==1 & gwf_st>1950, det local mean=`r(mean)' local median=`r(p50)' local p5=`r(p5)' local p95=`r(p95)' xtsum `var' if gwf_dem==1 & gwf_st>1950 post `outfile' ("`vlabel'") (`mean') (`median') (`p5') (`p95') (`r(sd)') (`r(sd_w)') (`r(sd_b)') } postclose `outfile' use summary_stats.dta, clear li restore use "/Users/alikadivar/Dropbox/Book/Quant/summary_stats.dta", clear export excel using "/Users/alikadivar/Dropbox/Book/Quant/Summary_statxlsx", firstrow(variables) *** descriptives xtset ccode year * capturing the score for electoral democracy for the first and last year gen elec_b=0 replace elec_b=l1.v2x_polyarchy if year==gwf_st+1 gen elec_e=0 replace elec_e=v2x_polyarchy if year==2010 | gwf_tta==1 bysort gwf_dem_id: egen elec_st=max(elec_b) * capturing the score for liberal democracy for the first and last year xtset ccode year gen lib_b=0 replace lib_b=l1.v2x_liberal if year==gwf_st+1 gen lib_e=0 replace lib_e=v2x_liberal if year==2010 | gwf_tta==1 bysort gwf_dem_id: egen lib_st=max(lib_b) *same for civil society gen ccsi_b=0 xtset ccode year replace ccsi_b=l1.v2xcs_ccsi if year==gwf_st+1 gen ccsi_e=0 replace ccsi_e=v2xcs_ccsi if year==2010 | gwf_tta==1 bysort gwf_dem_id: egen ccsi_st=max(ccsi_b) keep if gwf_dem==1 & gwf_st>1950 keep if year==2010 | gwf_tta==1 gen elec_dif=elec_e-elec_st gen lib_dif=lib_e-lib_st gen ccsi_dif=ccsi_e-ccsi_st gen survival=0 replace survival=1 if gwf_tta==0 drop if year<1960 drop if gwf_casename=="Tanzania 60-64" drop if gwf_casename=="Czech Republic 93-NA" drop if gwf_casename=="Ghana 56-60" su survival gwf_duration elec_dif lib_dif ccsi_dif if nv_regime_dur==0 su elec_dif lib_dif ccsi_dif if nv_regime_dur==1 su elec_dif lib_dif ccsi_dif if nv_regime_dur==2 su elec_dif lib_dif ccsi_dif if nv_regime_dur==3 su elec_dif lib_dif ccsi_dif if nv_regime_dur==4 su survival gwf_duration elec_dif lib_dif ccsi_dif if nv_regime_dur>0 & nv_regime_dur<5 su survival gwf_duration elec_dif lib_dif ccsi_dif if nv_regime_dur>4 su survival gwf_duration elec_dif lib_dif ccsi_dif if nv_regime_dur>0 & nv_regime_dur<4 su survival gwf_duration elec_dif lib_dif ccsi_dif if nv_regime_dur>3 *** Book Charts use "/Users/alikadivar/Dropbox/Country_Year_V-Dem_Full+others_STATA_v11.1/V-Dem-CY-Full+Others-v11.1.dta", clear xtset country_id year set scheme s1mono replace v2x_polyarchy=v2x_polyarchy*100 replace country_id=1000 if country_id==13 replace country_id=2000 if country_id==98 labmask country_id, values(country_name) xtline v2x_polyarchy if year>1960 & (country_name=="South Africa" | country_name=="Poland" | country_name=="Pakistan" | country_name=="Tunisia" |country_name=="Egypt") graph export polayrchy_chart.png, width(1600) replace xtline v2x_polyarchy if year>1900 & ( country_name=="Tunisia" |country_name=="Turkey" | country_name=="Israel") xtline v2x_polyarchy if year>1940 & ( country_name=="Poland" |country_name=="Hungary" |country_name=="Czech Republic")