*PITF GLOBAL MODEL REPLICATION
*
*Variables used in the PITF global model:
*
*SFTPTV2A='Regime type'
*LOGIM='Infant Mortality'
*MACCAT='Bordering States with Major Civil or Ethnic Conflict (4+ vs. 0-3)'
*DISP4CAT='State-Led Discrimination (yes vs. no)'
*
*Variables used in Fearon & Laitin model:
*
*LOG_GDPC='Per capita income, logged'
*LOGTPOP='Population, logged'
*LOGMTN ='% Mountainous, logged'
*FEANCTIG ='Noncontiguous state (Yes vs. No)'
*FLOIL ='Oil export revenues from fossil fuels of total export (>=1/3  vs. <1/3 )'
*SFTPDUR2 ='Instability dummy (three-step change in Polity IV scale, Yes vs. No)'
*ETHHERF ='Ethnic fractionalization'
*RELHREL ='Religious fractionalization'
*ANOCRACY ='Anocracy dummy (-5<=POLX<=5 vs. POLX<-5)'
*DEMOCRACY ='Democracy dummy (POLX>5 vs. POLX<-5)'
*
*Indicators to be used for subseting:
*
*MISS_IND ='subset to 0 to restrict global model to same sample size as F/L'
*GLB_IND = 'Y to run model for All Problems'
*CWAR_IND = 'Y to run model for Civil War Onsets'
*REG_IND = 'Y to run model for Adverse Regime Change Onsets'
*

quietly tabulate sftptv2a, gen(sftptv2a)
label variable sftptv2a1 "Full autocracy"
label variable sftptv2a2 "Full democracy"
label variable sftptv2a3 "Partial autocracy"
label variable sftptv2a4 "Partial democracy with factionalism"
label variable sftptv2a5 "Partial democracy without factionalism"
label variable sftptv2a6 "Transition (-88)"
quietly tabulate byregn2, gen(byregn2_)
label variable byregn2_1 "AF"
label variable byregn2_2 "AS"
label variable byregn2_3 "EU"
label variable byregn2_4 "LA"
label variable byregn2_5 "ME"
label variable maccat "Armed conflict in 4+ bordering states"
label variable sftpcons "Instability onset indicator (DV)"
label variable logim "Infant mortality"
label variable disp4cat "State-led discrimination"
egen stratida = concat(year byregn2 group) if sample==1
egen stratidb = concat(year byregn2 group) if sample==2
egen stratidc = concat(year byregn2 group) if sample==3
label variable stratida "Stratum ID for Sample A"
label variable stratidb "Stratum ID for Sample B"
label variable stratidc "Stratum ID for Sample C"

**TABLE 1

**Full Problem Set

clogit sftpcons sftptv2a3 sftptv2a4 sftptv2a5 sftptv2a2 sftptv2a6 logim maccat disp4cat if sample==1 & glb_ind=="Y", group(stratida) or
clogit sftpcons sftptv2a3 sftptv2a4 sftptv2a5 sftptv2a2 sftptv2a6 logim maccat disp4cat if sample==2 & glb_ind=="Y", group(stratidb) or
clogit sftpcons sftptv2a3 sftptv2a4 sftptv2a5 sftptv2a2 sftptv2a6 logim maccat disp4cat if sample==3 & glb_ind=="Y", group(stratidc) or

**Civil War Onsets

clogit sftpcons sftptv2a3 sftptv2a4 sftptv2a5 sftptv2a2 sftptv2a6 logim maccat disp4cat if sample==1 & cwar_ind=="Y", group(stratida) or
clogit sftpcons sftptv2a3 sftptv2a4 sftptv2a5 sftptv2a2 sftptv2a6 logim maccat disp4cat if sample==2 & cwar_ind=="Y", group(stratidb) or
clogit sftpcons sftptv2a3 sftptv2a4 sftptv2a5 sftptv2a2 sftptv2a6 logim maccat disp4cat if sample==3 & cwar_ind=="Y", group(stratidc) or

**Adverse Regime Change Onsets

clogit sftpcons sftptv2a3 sftptv2a4 sftptv2a5 sftptv2a2 sftptv2a6 logim maccat disp4cat if sample==1 & reg_ind=="Y", group(stratida) or
clogit sftpcons sftptv2a3 sftptv2a4 sftptv2a5 sftptv2a2 sftptv2a6 logim maccat disp4cat if sample==2 & reg_ind=="Y", group(stratidb) or
clogit sftpcons sftptv2a3 sftptv2a4 sftptv2a5 sftptv2a2 sftptv2a6 logim maccat disp4cat if sample==3 & reg_ind=="Y", group(stratidc) or

**TABLE 3

**All Problems

****Fearon & Laitin Model

clogit sftpcons log_gdpc logtpop logmtn feanctig floil sftpdur2 ethherf relhrel anocracy democracy if sample==1 & glb_ind=="Y" & miss_ind==0, group(stratida) or
clogit sftpcons log_gdpc logtpop logmtn feanctig floil sftpdur2 ethherf relhrel anocracy democracy if sample==2 & glb_ind=="Y" & miss_ind==0, group(stratidb) or
clogit sftpcons log_gdpc logtpop logmtn feanctig floil sftpdur2 ethherf relhrel anocracy democracy if sample==3 & glb_ind=="Y" & miss_ind==0, group(stratidc) or

****PITF Model

clogit sftpcons sftptv2a3 sftptv2a4 sftptv2a5 sftptv2a2 sftptv2a6 logim maccat disp4cat if sample==1 & glb_ind=="Y" & miss_ind==0, group(stratida) or
clogit sftpcons sftptv2a3 sftptv2a4 sftptv2a5 sftptv2a2 sftptv2a6 logim maccat disp4cat if sample==2 & glb_ind=="Y" & miss_ind==0, group(stratidb) or
clogit sftpcons sftptv2a3 sftptv2a4 sftptv2a5 sftptv2a2 sftptv2a6 logim maccat disp4cat if sample==3 & glb_ind=="Y" & miss_ind==0, group(stratidc) or

**Civil War Onsets

****Fearon & Laitin Model

clogit sftpcons log_gdpc logtpop logmtn feanctig floil sftpdur2 ethherf relhrel anocracy democracy if sample==1 & cwar_ind=="Y" & miss_ind==0, group(stratida) or
clogit sftpcons log_gdpc logtpop logmtn feanctig floil sftpdur2 ethherf relhrel anocracy democracy if sample==2 & cwar_ind=="Y" & miss_ind==0, group(stratidb) or
clogit sftpcons log_gdpc logtpop logmtn feanctig floil sftpdur2 ethherf relhrel anocracy democracy if sample==3 & cwar_ind=="Y" & miss_ind==0, group(stratidc) or

****PITF Model

clogit sftpcons sftptv2a3 sftptv2a4 sftptv2a5 sftptv2a2 sftptv2a6 logim maccat disp4cat if sample==1 & cwar_ind=="Y" & miss_ind==0, group(stratida) or
clogit sftpcons sftptv2a3 sftptv2a4 sftptv2a5 sftptv2a2 sftptv2a6 logim maccat disp4cat if sample==2 & cwar_ind=="Y" & miss_ind==0, group(stratidb) or
clogit sftpcons sftptv2a3 sftptv2a4 sftptv2a5 sftptv2a2 sftptv2a6 logim maccat disp4cat if sample==3 & cwar_ind=="Y" & miss_ind==0, group(stratidc) or

**Adverse Regime Change Onsets

****Fearon & Laitin Model

clogit sftpcons log_gdpc logtpop logmtn feanctig floil sftpdur2 ethherf relhrel anocracy democracy if sample==1 & reg_ind=="Y" & miss_ind==0, group(stratida) or
clogit sftpcons log_gdpc logtpop logmtn feanctig floil sftpdur2 ethherf relhrel anocracy democracy if sample==2 & reg_ind=="Y" & miss_ind==0, group(stratidb) or
clogit sftpcons log_gdpc logtpop logmtn feanctig floil sftpdur2 ethherf relhrel anocracy democracy if sample==3 & reg_ind=="Y" & miss_ind==0, group(stratidc) or

****PITF Model

clogit sftpcons sftptv2a3 sftptv2a4 sftptv2a5 sftptv2a2 sftptv2a6 logim maccat disp4cat if sample==1 & reg_ind=="Y" & miss_ind==0, group(stratida) or
clogit sftpcons sftptv2a3 sftptv2a4 sftptv2a5 sftptv2a2 sftptv2a6 logim maccat disp4cat if sample==2 & reg_ind=="Y" & miss_ind==0, group(stratidb) or
clogit sftpcons sftptv2a3 sftptv2a4 sftptv2a5 sftptv2a2 sftptv2a6 logim maccat disp4cat if sample==3 & reg_ind=="Y" & miss_ind==0, group(stratidc) or