* Elite appointment under Kim regime **

*** descriptive statistics ***
estpost tabstat appt party_back gov_back mil_back other_back age tenure ln_gdp gdp_growth, listwise statistics(count mean sd max min) columns(statistics)
esttab, cells("count mean sd max min") nomtitle nonumber

*** appt models ***
qui:nbreg appt mil_back if kim_era==1, cluster(id) nolog
est sto m1
qui:nbreg appt party_back if kim_era==1, cluster(id) nolog
est sto m2
qui:nbreg appt gov_back if kim_era==1, cluster(id) nolog
est sto m3
qui:nbreg appt other_back if kim_era==1, cluster(id) nolog
est sto m4
qui:nbreg appt mil_back party_back if kim_era==1, cluster(id) nolog
est sto m5
qui:nbreg appt mil_back party_back other_back if kim_era==1, cluster(id) nolog
est sto m6
qui:nbreg appt mil_back party_back other_back age tenure ln_gdp gdp_growth if kim_era==1, cluster(id) nolog
est sto m7
esttab m1 m2 m3 m4 m5 m6 m7, se nonumber title("Elite appointment in the Kim Jong Un era") mtitle("1" "2" "3" "4" "5" "6" "7") star(* 0.05 ** 0.01) varwidth(10) modelwidth(10) nobaselevel
est clear

*** after 2016 ***
qui:nbreg appt mil_back##after2016 party_back##after2016 other_back##after2016 if kim_era==1, cluster(id) nolog
est sto m8
qui:nbreg appt mil_back##after2016 party_back##after2016 other_back##after2016 age tenure ln_gdp gdp_growth if kim_era==1, cluster(id) nolog
est sto m9
esttab m8 m9, se nonumber title("Elite appointment in the Kim Jong Un era") mtitle("1" "2" ) star(* 0.05 ** 0.01) varwidth(20) modelwidth(20) nobaselevel
est clear

*** Comparison of results between Poisson and Negative Binomial models ***

nbreg appt mil_back party_back other_back age tenure ln_gdp gdp_growth if kim_era==1, cluster(id) nolog
est sto m1
poisson appt mil_back party_back other_back age tenure ln_gdp gdp_growth if kim_era==1, cluster(id) nolog
est sto m2
nbreg appt mil_back##after2016 party_back##after2016 other_back##after2016 age tenure ln_gdp gdp_growth if kim_era==1, cluster(id) nolog
est sto m3
poisson appt mil_back##after2016 party_back##after2016 other_back##after2016 age tenure ln_gdp gdp_growth if kim_era==1, cluster(id) nolog
est sto m4
esttab m1 m2 m3 m4, se nonumber title("model comparison") mtitle("poisson" "negative binomial" "poisson" "negative binomial" ) star(* 0.05 ** 0.01) varwidth(20) modelwidth(20) nobaselevel
est clear