** replication do file for IO internal review

clear all
set more off


use io_basedata, clear

stset warday, id(warnum) f(endday)
** table 3, model 1
stcox lndeadpop pcturnsum pa10lag deminit terrain contig relcap nstates majpowboth durscale civclash if default==1, robust nohr

stset warday, id(warnum) f(enddayset)
** table 3, model 2
stcrreg lndeadpop pcturnsum pa10lag deminit terrain contig relcap nstates majpowboth durscale civclash if default==1, robust compete(enddaycq) nohr
** table 3, model 3
stcrreg lndeadpop delinten60a pcturnsum pa10lag deminit terrain contig relcap nstates majpowboth durscale civclash if default==1, robust compete(enddaycq) nohr

** table 3, model 5
stset warday, id(warnum) f(enddaycq)
stcrreg lndeadpop pcturnsum pa10lag deminit terrain contig relcap nstates majpowboth durscale civclash if default==1, robust compete(enddayset) nohr

stset warday, id(warnum) f(endday)
** table 4, model 1
stcox lndeadpop pcncsum pcnewcsum pa10lag deminit terrain contig relcap nstates majpowboth durscale civclash if default==1, robust nohr

stset warday, id(warnum) f(enddayset)
** table 4, model 2
stcrreg lndeadpop pcncsum pcnewcsum pa10lag deminit terrain contig relcap nstates majpowboth durscale civclash if default==1, robust compete(enddaycq) nohr
** table 4, model 3
stcrreg lndeadpop pcncsumc pcnewcsumc pa10lag deminit terrain contig relcap nstates majpowboth durscale civclash if default==1, robust compete(enddaycq) nohr

stset warday, id(warnum) f(enddaycq)
** table 4, model 5
stcrreg lndeadpop pcncsum pcnewcsum pa10lag deminit terrain contig relcap nstates majpowboth durscale civclash if default==1, robust compete(enddayset) nohr


use io_cntrydata, clear

stset warday, id(idvar) f(cowset2)
** table 3, model 4
stcrreg lnodeadpop bft60 pcturn pa10lag terrain contig relcap nstates majpowboth durscale civclash deminit if default==1, cluster(warnum) compete(cowalt2) nohr
** table 4, model 4
stcrreg lnodeadpop bft60 pcnc pcnewc pa10lag terrain contig relcap nstates majpowboth durscale civclash deminit if default==1, cluster(warnum) compete(cowalt2) nohr


use io_leaderdata, clear

stset obsday, id(idvar) f(turnover)
** table 5, model 1
stcox otrend60 oshare lndays irreg wondum lostdum demdum cinc tpopalt if default==1, cluster(ccode) nohr
** table 5, model 2
stcox otrend60 oshare lndays irreg wondum lostdum demdum cinc tpopalt growthlag lngdppc if default==1, cluster(ccode) nohr

stset obsday, id(idvar) f(nculp)
** table 5, model 3
stcrreg otrend60 oshare lndays irreg wondum lostdum demdum cinc tpopalt if default==1, cluster(ccode) compete(altout) nohr
** table 5, model 4
stcrreg otrend60 oshare lndays irreg wondum lostdum demdum cinc tpopalt growth lngdppc if default==1, cluster(ccode) compete(altout) nohr
** table 5, model 5
stcrreg otrend60 demtrend demdum oshare lndays irreg wondum lostdum cinc tpopalt if default==1, cluster(ccode) compete(altout) nohr

stset obsday, id(idvar) f(culpdum)
** table 5, model 6
stcrreg otrend60 oshare lndays irreg wondum lostdum demdum cinc tpopalt if default==1, cluster(ccode) compete(altout2) nohr