version 11
set matsize 1000
tsset dyadid yearid

* Table 1
quietly xi: logit  dr_at_recipient_dum  lndistance powerpc_recipient powerpc_sender    s3uni   c.s3uni#c.powerpc_recipient  c.s3uni#c.powerpc_sender  lngdppc_recipient  lngdppc_sender  westandruscoldummy dr_at_sender_dum  if  samedyad==0, 
su dr_at_recipient_dum lndistance powerpc_recipient powerpc_sender    s3uni  lngdppc_recipient  lngdppc_sender   c.s3uni#c.powerpc_recipient  c.s3uni#c.powerpc_sender  lngdppc_recipient  lngdppc_sender  westandruscoldummy dr_at_sender_dum  bilattradegdp_sender independence_sender demdem autocautoc  if e(sample)


* Table 2
xi: logit  dr_at_recipient_dum   lndistance powerpc_recipient powerpc_sender    s3uni   c.s3uni#c.powerpc_recipient  c.s3uni#c.powerpc_sender  lngdppc_recipient  lngdppc_sender  westandruscoldummy  i.year if  samedyad==0, cluster( dyadid)
estat class
margins, dydx(s3uni) at(powerpc_recipient=(-2.2 -1.2 -.2 .8 1.8 2.8 3.8 4.8)) atmeans
margins, dydx(s3uni) at(powerpc_sender=(-2.2 -1.2 -.2 .8 1.8 2.8 3.8 4.8)) atmeans
margins, dydx(*) atmeans


xi: logit  dr_at_recipient_dum   lndistance powerpc_recipient powerpc_sender    s3uni   c.s3uni#c.powerpc_recipient  c.s3uni#c.powerpc_sender  lngdppc_recipient  lngdppc_sender  westandruscoldummy  dr_at_sender_dum  i.year if  samedyad==0, cluster( dyadid)
estat class
margins, dydx(s3uni) at(powerpc_recipient=(-2.2 -1.2 -.2 .8 1.8 2.8 3.8 4.8)) atmeans
margins, dydx(s3uni) at(powerpc_sender=(-2.2 -1.2 -.2 .8 1.8 2.8 3.8 4.8)) atmeans
margins, dydx(*) atmeans

** Robustness tests
* LDV
capture drop ldr_at_recipient_dum
gen ldr_at_recipient_dum=l.dr_at_recipient_dum
xi: logit  dr_at_recipient_dum  ldr_at_recipient_dum  lndistance powerpc_recipient powerpc_sender    s3uni   c.s3uni#c.powerpc_recipient  c.s3uni#c.powerpc_sender  lngdppc_recipient  lngdppc_sender  westandruscoldummy  i.year if  samedyad==0, cluster( dyadid)
estat class
margins, dydx(s3uni) at(powerpc_recipient=(-2.2 -1.2 -.2 .8 1.8 2.8 3.8 4.8)) atmeans
margins, dydx(s3uni) at(powerpc_sender=(-2.2 -1.2 -.2 .8 1.8 2.8 3.8 4.8)) atmeans
margins, dydx(*) atmeans

* Include bilateral trade
xi: logit  dr_at_recipient_dum   lndistance powerpc_recipient powerpc_sender    s3uni   c.s3uni#c.powerpc_recipient  c.s3uni#c.powerpc_sender  lngdppc_recipient  lngdppc_sender  westandruscoldummy   bilattradegdp_sender i.year if  samedyad==0, cluster( dyadid)
estat class
margins, dydx(s3uni) at(powerpc_recipient=(-2.2 -1.2 -.2 .8 1.8 2.8 3.8 4.8)) atmeans
margins, dydx(s3uni) at(powerpc_sender=(-2.2 -1.2 -.2 .8 1.8 2.8 3.8 4.8)) atmeans
margins, dydx(*) atmeans

* Include year of independence
xi: logit  dr_at_recipient_dum   lndistance powerpc_recipient powerpc_sender    s3uni   c.s3uni#c.powerpc_recipient  c.s3uni#c.powerpc_sender  lngdppc_recipient  lngdppc_sender  westandruscoldummy   independence_sender c.independence_sender#c.independence_sender  i.year if  samedyad==0, cluster( dyadid)
estat class
margins, dydx(s3uni) at(powerpc_recipient=(-2.2 -1.2 -.2 .8 1.8 2.8 3.8 4.8)) atmeans
margins, dydx(s3uni) at(powerpc_sender=(-2.2 -1.2 -.2 .8 1.8 2.8 3.8 4.8)) atmeans
margins, dydx(independence_sender) at(independence_sender=(1816 1850 1900 1925 1950 1970 1990)) atmeans
margins, dydx(*) atmeans


* Include regime type dyad dummies
xi: logit  dr_at_recipient_dum   lndistance powerpc_recipient powerpc_sender    s3uni   c.s3uni#c.powerpc_recipient  c.s3uni#c.powerpc_sender  lngdppc_recipient  lngdppc_sender  westandruscoldummy   demdem autocautoc i.year if  samedyad==0, cluster( dyadid)
estat class
margins, dydx(s3uni) at(powerpc_recipient=(-2.2 -1.2 -.2 .8 1.8 2.8 3.8 4.8)) atmeans
margins, dydx(s3uni) at(powerpc_sender=(-2.2 -1.2 -.2 .8 1.8 2.8 3.8 4.8)) atmeans
margins, dydx(*) atmeans


* Divide into Cold War and post-Cold War period
xi: logit  dr_at_recipient_dum   lndistance powerpc_recipient powerpc_sender    s3uni   c.s3uni#c.powerpc_recipient  c.s3uni#c.powerpc_sender  lngdppc_recipient  lngdppc_sender  westandruscoldummy  i.year if  samedyad==0 & year<1990, cluster( dyadid)
estat class
margins, dydx(s3uni) at(powerpc_recipient=(-2.2 -1.2 -.2 .8 1.8 2.8 3.8 4.8)) atmeans
margins, dydx(s3uni) at(powerpc_sender=(-2.2 -1.2 -.2 .8 1.8 2.8 3.8 4.8)) atmeans
margins, dydx(*) atmeans

xi: logit  dr_at_recipient_dum   lndistance powerpc_recipient powerpc_sender    s3uni   c.s3uni#c.powerpc_recipient  c.s3uni#c.powerpc_sender  lngdppc_recipient  lngdppc_sender  westandruscoldummy  i.year if  samedyad==0 & year>=1990, cluster( dyadid)
estat class
margins, dydx(s3uni) at(powerpc_recipient=(-2.2 -1.2 -.2 .8 1.8 2.8 3.8 4.8)) atmeans
margins, dydx(s3uni) at(powerpc_sender=(-2.2 -1.2 -.2 .8 1.8 2.8 3.8 4.8)) atmeans
margins, dydx(*) atmeans