reg leg1 mnp partydecentr partage partideolord lagleg1 diffsubreg1 gdppercap lninflation, cluster (id) fitstat reg leg1 two multicoal2 partydecentr partage partideolord lagleg1 diffsubreg1 gdppercap lninflation, cluster (id) fitstat gen mnpdec = mnp*partydecentr reg leg1 mnp partydecentr mnpdec partage partideolord lagleg1 diffsubreg1 gdppercap lninflation, cluster (id) fitstat lincom mnp + mnpdec lincom partydecentr + mnpdec lincom partydecentr + mnpdec* 1 lincom partydecentr + mnpdec* 2 lincom partydecentr + mnpdec* 3 lincom partydecentr + mnpdec* 4 lincom partydecentr + mnpdec* 5 reg agenda mnp partydecentr partage partideolord lagagenda diffsubregagenda gdppercap lninflation, cluster (id) fitstat reg agenda two multicoal2 partydecentr partage partideolord lagagenda diffsubregagenda gdppercap lninflation, cluster (id) fitstat reg agenda mnp partydecentr mnpdec partage partideolord lagagenda diffsubregagenda gdppercap lninflation, cluster (id) fitstat reg leg3 mnp partydecentr partage partideolord lagleg3 diffsubreg3 gdppercap lninflation, cluster (id) fitstat reg leg3 two multicoal2 partydecentr partage partideolord lagleg3 diffsubreg3 gdppercap lninflation, cluster (id) fitstat reg leg3 mnp partydecentr mnpdec partage partideolord lagleg3 diffsubreg3 gdppercap lninflation, cluster (id) fitstat *GRAFICA reg leg1 mnp partydecentr mnpdec partage partideolord lagleg1 diffsubreg1 gdppercap lninflation, cluster (id) fitstat lincom partydecentr + mnpdec*1 mat B1 = r(estimate) mat se1 = r(se) lincom partydecentr + mnpdec*2 mat B2 = r(estimate) mat se2 = r(se) lincom partydecentr + mnpdec*3 mat B3 = r(estimate) mat se3 = r(se) lincom partydecentr + mnpdec*4 mat B4 = r(estimate) mat se4 = r(se) lincom partydecentr + mnpdec*5 mat B5 = r(estimate) mat se5 = r(se) clear svmat B1 svmat B2 svmat B3 svmat B4 svmat B5 mkmat B11 B21 B31 B41 B51, mat(B) mat BT = B' clear svmat se1 svmat se2 svmat se3 svmat se4 svmat se5 mkmat se11 se21 se31 se41 se51, mat(se) mat seT = se' clear svmat BT svmat seT rename BT1 B rename seT1 se *Calculate Confidence Interval gen cil = B - 1.96*se gen ciu = B + 1.96*se *Generate Parameter ID Variable gen parmid = _n eclplot B cil ciu parmid, title(" ") xtitle("Size of Reform Coalition") xlabel(1 " 1 " 2 " 2" 3 " 3" 4 "4" 5 "5") xsca(titlegap(5))ytitle("Marginal Effect of Party Decentralization") yline(0, lw(medthick) lc(black)) legend( on order (2 "Party Decentralization" 1 "95% Confidence Interval"))