************************************************************************************************
*********** "To Link or Not to Link? Agenda Change in International Bargaining" ****************
**************************British Journal of Political Science**********************************
******************************author: Heather Elko McKibben*************************************
************************************************************************************************


use McKibben_BJPS_Linking-paper.dta, replace

***********Create "No-Agreement Costs" Variable**********
******(Interaction of of EVI and GDP per capita**********

**Create an index using "EVI" and gdp per capita 
generate evi_inv=evi*-1+4.8125
generate interim=evi_inv*ln_gdp_capita
generate no_agree_costs=interim*-1+36.80302429199219

**Create interaction terms
generate interaction1=no_agree_costs*diff_valued_agenda
generate interaction2=lnco2*diff_valued_agenda

****************************************************
********************Model 1*************************
****************Issue Addition**********************
****************************************************

*Model 1A:
*RUN ISSUE ADDITION MODEL, CALCULATE MARGINAL EFFECTS OF INTERACTION AT MEAN OF "DIFFERENTLY-VALUED ISSUES" VARIABLE
gllamm il_add diff_valued_agenda no_agree_costs interaction1 lnco2 high_level contact_group lngdp, i(i_bargain bali) link(logit) 

preserve

set seed 339487731

drawnorm SN_b1-SN_b11, n(10000) means(e(b)) cov(e(V)) clear

postutil clear
postfile mypost prob_hat0 lo0 hi0 prob_hat1 lo1 hi1 diff_hat diff_lo diff_hi using "McKibben_BJPS_Linking_interim", replace

noisily display "start"

local a= 0
while `a' <= 36.80302429199219 { 

{

*CALCULATION AT MEAN OF DIFF-VALUED AGENDA
scalar h_diff_valued_agenda =  -.0320825
scalar h_lnco2 =  9.331776
scalar h_high_level = 0
scalar h_contact_group = 0
scalar h_lngdp =  23.85371
scalar h_ln_gdp_capita =  8.342014
scalar h_bali = 0
scalar h_constant = 1 
scalar h_no_agree_costs = 23.94389       


    generate x_betahat0  = SN_b1*h_diff_valued_agenda + SN_b2*`a' + SN_b3*h_diff_valued_agenda*`a' + SN_b4*h_lnco2 + SN_b5*h_high_level + SN_b6*h_contact_group + SN_b7*h_lngdp + SN_b8*h_constant
                            
    generate x_betahat1  = SN_b1*(h_diff_valued_agenda+1) + SN_b2*`a' + SN_b3*(h_diff_valued_agenda+1)*`a' + SN_b4*h_lnco2 + SN_b5*h_high_level + SN_b6*h_contact_group + SN_b7*h_lngdp + SN_b8*h_constant
    
    gen prob0 =normal(x_betahat0)
    gen prob1=normal(x_betahat1)
    gen diff=prob1-prob0
    
    egen probhat0 =mean(prob0)
    egen probhat1=mean(prob1)
    egen diffhat=mean(diff)
    
    tempname prob_hat0 lo0 hi0 prob_hat1 lo1 hi1 diff_hat diff_lo diff_hi 

    _pctile prob0, p(2.5,97.5) 
    scalar `lo0' = r(r1)
    scalar `hi0' = r(r2)  
    
    _pctile prob1, p(2.5,97.5)
    scalar `lo1'= r(r1)
    scalar `hi1'= r(r2) 
    
    _pctile diff, p(2.5,97.5)
    scalar `diff_lo'= r(r1)
    scalar `diff_hi'= r(r2) 
   
    scalar `prob_hat0'=probhat0
    scalar `prob_hat1'=probhat1
    scalar `diff_hat'=diffhat
    
    post mypost (`prob_hat0') (`lo0') (`hi0') (`prob_hat1') (`lo1') (`hi1') (`diff_hat') (`diff_lo') (`diff_hi')
   
    }      
    
    drop x_betahat0 x_betahat1 prob0 prob1 diff probhat0 probhat1 diffhat
    
    local a=(`a')+ 1
    
    display "." _c
    
} 

display ""
postclose mypost

restore
merge using McKibben_BJPS_Linking_interim


gen yline=0
gen MV = (0 + (_n-1))
replace MV=. if (0 + (_n-1)) > 36.80302


line diff_hat MV, clwidth(thick) ///
         || line diff_lo MV, lpattern(dash) lcolor(gray) clwidth(thin) ///
         || line diff_hi MV, lpattern(dash) lcolor(gray) clwidth(thin)  ///
	     || line yline MV, lpattern(solid) lcolor(black) clwidth(thin) ///
         || ,   ///
           xlabel(0 36.80302, nogrid labsize(2)) ///
           ylabel(, axis(1) nogrid labsize(2)) ///
		   legend(off) ///
           xtitle("Costliness of No Agreement", size(2.5)) ///
           ytitle("Marginal Effect of an Increase in Differently-Valued Nature of the Agenda", size(2.5)) ///
		   title("Bargaining Over More Differently Valued Issues") ///
		   subtitle("Marginal Effect at Different No-Agreement Values") ///
		   note("*Note: 95% Confidence Intervals Indicated", size(2)) ///
           xsca(titlegap(4)) ysca(titlegap(4)) ///
           scheme(s2mono) graphregion(fcolor(white) ilcolor(white) lcolor(white))

		   


*Model 1B:
*RUN ISSUE ADDITION MODEL, CALCULATE MARGINAL EFFECTS OF INTERACTION AT MIN OF "DIFFERENTLY-VALUED ISSUES" VARIABLE
use McKibben_BJPS_Linking-paper.dta, replace

**Create an index using "EVI" and gdp per capita 
generate evi_inv=evi*-1+4.8125
generate interim=evi_inv*ln_gdp_capita
generate no_agree_costs=interim*-1+36.80302429199219

**Create interaction terms
generate interaction1=no_agree_costs*diff_valued_agenda
generate interaction2=lnco2*diff_valued_agenda

gllamm il_add diff_valued_agenda no_agree_costs interaction1 lnco2 high_level contact_group lngdp, i(i_bargain bali) link(logit) 

preserve

set seed 339487731

drawnorm SN_b1-SN_b11, n(10000) means(e(b)) cov(e(V)) clear

postutil clear
postfile mypost prob_hat0 lo0 hi0 prob_hat1 lo1 hi1 diff_hat diff_lo diff_hi using "McKibben_BJPS_Linking_interim", replace

noisily display "start"

local a= 0
while `a' <= 36.80302429199219 { 

{

*CALCULATION AT MIN OF DIFF-VALUED AGENDA
scalar h_diff_valued_agenda =  -.158339
scalar h_lnco2 =  9.331776
scalar h_high_level = 0
scalar h_contact_group = 0
scalar h_lngdp =  23.85371
scalar h_ln_gdp_capita =  8.342014
scalar h_bali = 0
scalar h_constant = 1 
scalar h_no_agree_costs = 23.94389       


    generate x_betahat0  = SN_b1*h_diff_valued_agenda + SN_b2*`a' + SN_b3*h_diff_valued_agenda*`a' + SN_b4*h_lnco2 + SN_b5*h_high_level + SN_b6*h_contact_group + SN_b7*h_lngdp + SN_b8*h_constant
                            
    generate x_betahat1  = SN_b1*(h_diff_valued_agenda+1) + SN_b2*`a' + SN_b3*(h_diff_valued_agenda+1)*`a' + SN_b4*h_lnco2 + SN_b5*h_high_level + SN_b6*h_contact_group + SN_b7*h_lngdp + SN_b8*h_constant
    
    gen prob0 =normal(x_betahat0)
    gen prob1=normal(x_betahat1)
    gen diff=prob1-prob0
    
    egen probhat0 =mean(prob0)
    egen probhat1=mean(prob1)
    egen diffhat=mean(diff)
    
    tempname prob_hat0 lo0 hi0 prob_hat1 lo1 hi1 diff_hat diff_lo diff_hi 

    _pctile prob0, p(2.5,97.5) 
    scalar `lo0' = r(r1)
    scalar `hi0' = r(r2)  
    
    _pctile prob1, p(2.5,97.5)
    scalar `lo1'= r(r1)
    scalar `hi1'= r(r2) 
    
    _pctile diff, p(2.5,97.5)
    scalar `diff_lo'= r(r1)
    scalar `diff_hi'= r(r2) 
   
    scalar `prob_hat0'=probhat0
    scalar `prob_hat1'=probhat1
    scalar `diff_hat'=diffhat
    
    post mypost (`prob_hat0') (`lo0') (`hi0') (`prob_hat1') (`lo1') (`hi1') (`diff_hat') (`diff_lo') (`diff_hi')
   
    }      
    
    drop x_betahat0 x_betahat1 prob0 prob1 diff probhat0 probhat1 diffhat
    
    local a=(`a')+ 1
    
    display "." _c
    
} 

display ""
postclose mypost

restore
merge using McKibben_BJPS_Linking_interim

gen yline=0
gen MV = (0 + (_n-1))
replace MV=. if (0 + (_n-1)) > 36.80302


line diff_hat MV, clwidth(thick) ///
         || line diff_lo MV, lpattern(dash) lcolor(gray) clwidth(thin) ///
         || line diff_hi MV, lpattern(dash) lcolor(gray) clwidth(thin)  ///
	     || line yline MV, lpattern(solid) lcolor(black) clwidth(thin) ///
         || ,   ///
           xlabel(0 36.80302, nogrid labsize(2)) ///
           ylabel(, axis(1) nogrid labsize(2)) ///
		   legend(off) ///
           xtitle("Costliness of No Agreement", size(2.5)) ///
           ytitle("Marginal Effect of an Increase in Differently-Valued Nature of the Agenda", size(2.5)) ///
		   title("Bargaining Over More Differently Valued Issues") ///
		   subtitle("The Marginal Effect at Different BATNA Values") ///
		   note("*Note: 95% Confidence Intervals Indicated", size(2)) ///
           xsca(titlegap(4)) ysca(titlegap(4)) ///
           scheme(s2mono) graphregion(fcolor(white) ilcolor(white) lcolor(white))
		   
		   
*Model 1C:
*RUN ISSUE ADDITION MODEL, CALCULATE MARGINAL EFFECTS OF INTERACTION AT MAX OF "DIFFERENTLY-VALUED ISSUES" VARIABLE
use McKibben_BJPS_Linking-paper.dta, replace

**Create an index using "EVI" and gdp per capita 
generate evi_inv=evi*-1+4.8125
generate interim=evi_inv*ln_gdp_capita
generate no_agree_costs=interim*-1+36.80302429199219

**Create interaction terms
generate interaction1=no_agree_costs*diff_valued_agenda
generate interaction2=lnco2*diff_valued_agenda

gllamm il_add diff_valued_agenda no_agree_costs interaction1 lnco2 high_level contact_group lngdp, i(i_bargain bali) link(logit) 

preserve

set seed 339487731

drawnorm SN_b1-SN_b11, n(10000) means(e(b)) cov(e(V)) clear

postutil clear
postfile mypost prob_hat0 lo0 hi0 prob_hat1 lo1 hi1 diff_hat diff_lo diff_hi using "McKibben_BJPS_Linking_interim", replace

noisily display "start"

local a= 0
while `a' <= 36.80302429199219 { 

{

*CALCULATION AT MAX OF OF DIFF-VALUED AGENDA
scalar h_diff_valued_agenda =  .029419
scalar h_lnco2 =  9.331776
scalar h_high_level = 0
scalar h_contact_group = 0
scalar h_lngdp =  23.85371
scalar h_ln_gdp_capita =  8.342014
scalar h_bali = 0
scalar h_constant = 1 
scalar h_no_agree_costs = 23.94389       


    generate x_betahat0  = SN_b1*h_diff_valued_agenda + SN_b2*`a' + SN_b3*h_diff_valued_agenda*`a' + SN_b4*h_lnco2 + SN_b5*h_high_level + SN_b6*h_contact_group + SN_b7*h_lngdp + SN_b8*h_constant
                            
    generate x_betahat1  = SN_b1*(h_diff_valued_agenda+1) + SN_b2*`a' + SN_b3*(h_diff_valued_agenda+1)*`a' + SN_b4*h_lnco2 + SN_b5*h_high_level + SN_b6*h_contact_group + SN_b7*h_lngdp + SN_b8*h_constant
    
    gen prob0 =normal(x_betahat0)
    gen prob1=normal(x_betahat1)
    gen diff=prob1-prob0
    
    egen probhat0 =mean(prob0)
    egen probhat1=mean(prob1)
    egen diffhat=mean(diff)
    
    tempname prob_hat0 lo0 hi0 prob_hat1 lo1 hi1 diff_hat diff_lo diff_hi 

    _pctile prob0, p(2.5,97.5) 
    scalar `lo0' = r(r1)
    scalar `hi0' = r(r2)  
    
    _pctile prob1, p(2.5,97.5)
    scalar `lo1'= r(r1)
    scalar `hi1'= r(r2) 
    
    _pctile diff, p(2.5,97.5)
    scalar `diff_lo'= r(r1)
    scalar `diff_hi'= r(r2) 
   
    scalar `prob_hat0'=probhat0
    scalar `prob_hat1'=probhat1
    scalar `diff_hat'=diffhat
    
    post mypost (`prob_hat0') (`lo0') (`hi0') (`prob_hat1') (`lo1') (`hi1') (`diff_hat') (`diff_lo') (`diff_hi')
   
    }      
    
    drop x_betahat0 x_betahat1 prob0 prob1 diff probhat0 probhat1 diffhat
    
    local a=(`a')+ 1
    
    display "." _c
    
} 

display ""
postclose mypost

restore
merge using McKibben_BJPS_Linking_interim

gen yline=0
gen MV = (0 + (_n-1))
replace MV=. if (0 + (_n-1)) > 36.80302



line diff_hat MV, clwidth(thick) ///
         || line diff_lo MV, lpattern(dash) lcolor(gray) clwidth(thin) ///
         || line diff_hi MV, lpattern(dash) lcolor(gray) clwidth(thin)  ///
	   ||  line yline MV, lpattern(solid) lcolor(black) clwidth(thin) ///
       ||  ,   ///
           xlabel(0 36.80302, nogrid labsize(2)) ///
           ylabel(, axis(1) nogrid labsize(2)) ///
		   legend(off) ///
           xtitle("Costliness of No Agreement", size(2.5)) ///
           ytitle("Marginal Effect of an Increase in Differently-Valued Nature of the Agenda", size(2.5)) ///
		   title("Bargaining Over More Differently Valued Issues") ///
		   subtitle("The Marginal Effect at Different BATNA Values") ///
		   note("*Note: 95% Confidence Intervals Indicated", size(2)) ///
           xsca(titlegap(4)) ysca(titlegap(4))  ///
           scheme(s2mono) graphregion(fcolor(white) ilcolor(white) lcolor(white))
		   
		   

		   
****************************************************
********************Model 2*************************
**************Issue Subtraction*********************
****************************************************


use McKibben_BJPS_Linking-paper.dta, replace

*Model 2A:
*RUN ISSUE SUBTRACTION MODEL, CALCULATE MARGINAL EFFECTS OF INTERACTION AT MEAN OF "DIFFERENTLY-VALUED ISSUES" VARIABLE

**Create an index using "EVI" and gdp per capita 
generate evi_inv=evi*-1+4.8125
generate interim=evi_inv*ln_gdp_capita
generate no_agree_costs=interim*-1+36.80302429199219

**Create interaction terms
generate interaction1=no_agree_costs*diff_valued_agenda
generate interaction2=lnco2*diff_valued_agenda


gllamm il_subtract diff_valued_agenda lnco2 interaction2 no_agree_costs high_level contact_group lngdp, i(i_bargain bali) link(logit) 

preserve

set seed 339487731

drawnorm SN_b1-SN_b11, n(10000) means(e(b)) cov(e(V)) clear

postutil clear
postfile mypost prob_hat0 lo0 hi0 prob_hat1 lo1 hi1 diff_hat diff_lo diff_hi using "McKibben_BJPS_Linking_interim", replace

noisily display "start"

local a= 0
while `a' <= 15.8551 { 

{

*CALCULATE USING MEAN OF DIFF-VALUED AGENDA
scalar h_diff_valued_agenda =  -.0320825
scalar h_lnco2 =  9.331776
scalar h_high_level = 1
scalar h_contact_group = 0
scalar h_lngdp =  23.85371
scalar h_ln_gdp_capita =  8.342014
scalar h_bali = 1
scalar h_constant = 1 
scalar h_no_agree_costs = 23.94389       


    generate x_betahat0  = SN_b1*h_diff_valued_agenda + SN_b2*`a' + SN_b3*h_diff_valued_agenda*`a' + SN_b4*h_no_agree_costs + SN_b5*h_high_level + SN_b6*h_contact_group + SN_b7*h_lngdp + SN_b8*h_constant
                            
    generate x_betahat1  = SN_b1*(h_diff_valued_agenda+1) + SN_b2*`a' + SN_b3*(h_diff_valued_agenda+1)*`a' + SN_b4*h_no_agree_costs + SN_b5*h_high_level + SN_b6*h_contact_group + SN_b7*h_lngdp + SN_b8*h_constant
    
    gen prob0 =normal(x_betahat0)
    gen prob1=normal(x_betahat1)
    gen diff=prob1-prob0
    
    egen probhat0 =mean(prob0)
    egen probhat1=mean(prob1)
    egen diffhat=mean(diff)
    
    tempname prob_hat0 lo0 hi0 prob_hat1 lo1 hi1 diff_hat diff_lo diff_hi 

    _pctile prob0, p(2.5,97.5) 
    scalar `lo0' = r(r1)
    scalar `hi0' = r(r2)  
    
    _pctile prob1, p(2.5,97.5)
    scalar `lo1'= r(r1)
    scalar `hi1'= r(r2) 
    
    _pctile diff, p(2.5,97.5)
    scalar `diff_lo'= r(r1)
    scalar `diff_hi'= r(r2) 
   
    scalar `prob_hat0'=probhat0
    scalar `prob_hat1'=probhat1
    scalar `diff_hat'=diffhat
    
    post mypost (`prob_hat0') (`lo0') (`hi0') (`prob_hat1') (`lo1') (`hi1') (`diff_hat') (`diff_lo') (`diff_hi')
   
    }      
    
    drop x_betahat0 x_betahat1 prob0 prob1 diff probhat0 probhat1 diffhat
    
    local a=(`a')+ 1
    
    display "." _c
    
} 

display ""
postclose mypost

restore
merge using McKibben_BJPS_Linking_interim

gen yline=0
gen MV = (0 + (_n-1))
replace MV=. if (0 + (_n-1)) > 15.8551

line diff_hat MV, clwidth(thick) ///
         || line diff_lo MV, lpattern(dash) lcolor(gray) clwidth(thin) ///
         || line diff_hi MV, lpattern(dash) lcolor(gray) clwidth(thin)  ///
	   ||  line yline MV, lpattern(solid) lcolor(black) clwidth(thin) ///
       ||  , ///
           xlabel(3.93843 15.8551, nogrid labsize(2)) ///
           ylabel(, axis(1) nogrid labsize(2)) ///
		   legend(off) ///
           xtitle("CO2 Emissions", size(2.5)) ///
           ytitle("Marginal Effect of an Increase in Differently-Valued Nature of the Agenda", size(2.5)) ///
		   title("Bargaining Over More Differently Valued Issues") ///
		   subtitle("The Marginal Effect at Different Levels of CO2 Emissions") ///
		   note("*Note: 95% Confidence Intervals Indicated", size(2)) ///
           xsca(titlegap(4)) ysca(titlegap(4))  ///
           scheme(s2mono) graphregion(fcolor(white) ilcolor(white) lcolor(white)) ///



		   
*Model 2B:
*RUN ISSUE SUBTRACTION MODEL, CALCULATE MARGINAL EFFECTS OF INTERACTION AT MIN OF "DIFFERENTLY-VALUED ISSUES" VARIABLE
use McKibben_BJPS_Linking-paper.dta, replace


**Create an index using "EVI" and gdp per capita 
generate evi_inv=evi*-1+4.8125
generate interim=evi_inv*ln_gdp_capita
generate no_agree_costs=interim*-1+36.80302429199219

**Create interaction terms
generate interaction1=no_agree_costs*diff_valued_agenda
generate interaction2=lnco2*diff_valued_agenda
		   
gllamm il_subtract diff_valued_agenda lnco2 interaction2 no_agree_costs high_level contact_group lngdp, i(i_bargain bali) link(logit) 


preserve

set seed 339487731

drawnorm SN_b1-SN_b11, n(10000) means(e(b)) cov(e(V)) clear

postutil clear
postfile mypost prob_hat0 lo0 hi0 prob_hat1 lo1 hi1 diff_hat diff_lo diff_hi using "McKibben_BJPS_Linking_interim", replace

noisily display "start"

local a= 0
while `a' <= 15.8551 { 

{

*CALCULATE USING MIN OF DIFF-VALUED AGENDA
scalar h_diff_valued_agenda =  -.158339
scalar h_lnco2 =  9.331776
scalar h_high_level = 1
scalar h_contact_group = 0
scalar h_lngdp =  23.85371
scalar h_ln_gdp_capita =  8.342014
scalar h_bali = 0
scalar h_constant = 1 
scalar h_no_agree_costs = 23.94389       


    generate x_betahat0  = SN_b1*h_diff_valued_agenda + SN_b2*`a' + SN_b3*h_diff_valued_agenda*`a' + SN_b4*h_no_agree_costs + SN_b5*h_high_level + SN_b6*h_contact_group + SN_b7*h_lngdp + SN_b8*h_constant
                            
    generate x_betahat1  = SN_b1*(h_diff_valued_agenda+1) + SN_b2*`a' + SN_b3*(h_diff_valued_agenda+1)*`a' + SN_b4*h_no_agree_costs + SN_b5*h_high_level + SN_b6*h_contact_group + SN_b7*h_lngdp + SN_b8*h_constant
    
    gen prob0 =normal(x_betahat0)
    gen prob1=normal(x_betahat1)
    gen diff=prob1-prob0
    
    egen probhat0 =mean(prob0)
    egen probhat1=mean(prob1)
    egen diffhat=mean(diff)
    
    tempname prob_hat0 lo0 hi0 prob_hat1 lo1 hi1 diff_hat diff_lo diff_hi 

    _pctile prob0, p(2.5,97.5) 
    scalar `lo0' = r(r1)
    scalar `hi0' = r(r2)  
    
    _pctile prob1, p(2.5,97.5)
    scalar `lo1'= r(r1)
    scalar `hi1'= r(r2) 
    
    _pctile diff, p(2.5,97.5)
    scalar `diff_lo'= r(r1)
    scalar `diff_hi'= r(r2) 
   
    scalar `prob_hat0'=probhat0
    scalar `prob_hat1'=probhat1
    scalar `diff_hat'=diffhat
    
    post mypost (`prob_hat0') (`lo0') (`hi0') (`prob_hat1') (`lo1') (`hi1') (`diff_hat') (`diff_lo') (`diff_hi')
   
    }      
    
    drop x_betahat0 x_betahat1 prob0 prob1 diff probhat0 probhat1 diffhat
    
    local a=(`a')+ 1
    
    display "." _c
    
} 

display ""
postclose mypost

restore
merge using McKibben_BJPS_Linking_interim

gen yline=0
gen MV = (0 + (_n-1))
replace MV=. if (0 + (_n-1)) > 15.8551

line diff_hat MV, clwidth(thick) ///
         || line diff_lo MV, lpattern(dash) lcolor(gray) clwidth(thin) ///
         || line diff_hi MV, lpattern(dash) lcolor(gray) clwidth(thin)  ///
	   ||  line yline MV, lpattern(solid) lcolor(black) clwidth(thin) ///
       ||  , ///
           xlabel(3.93843 15.8551, nogrid labsize(2)) ///
           ylabel(, axis(1) nogrid labsize(2)) ///
		   legend(off) ///
           xtitle("CO2 Emissions", size(2.5)) ///
           ytitle("Marginal Effect of an Increase in Differently-Valued Nature of the Agenda", size(2.5)) ///
		   title("Bargaining Over More Differently Valued Issues") ///
		   subtitle("The Marginal Effect at Different Levels of CO2 Emissions") ///
		   note("*Note: 95% Confidence Intervals Indicated", size(2)) ///
           xsca(titlegap(4)) ysca(titlegap(4))  ///
           scheme(s2mono) graphregion(fcolor(white) ilcolor(white) lcolor(white)) ///


		   
*Model 2C:
*RUN ISSUE SUBTRACTION MODEL, CALCULATE MARGINAL EFFECTS OF INTERACTION AT MAX OF "DIFFERENTLY-VALUED ISSUES" VARIABLE
use McKibben_BJPS_Linking-paper.dta, replace


**Create an index using "EVI" and gdp per capita 
generate evi_inv=evi*-1+4.8125
generate interim=evi_inv*ln_gdp_capita
generate no_agree_costs=interim*-1+36.80302429199219

**Create interaction terms
generate interaction1=no_agree_costs*diff_valued_agenda
generate interaction2=lnco2*diff_valued_agenda
		   
*RUN MODEL FOR INTERACTION2 MARGINAL EFFECTS INTERPRETATION AT MAX OF IL STRUCTURE VARIABLE
gllamm il_subtract diff_valued_agenda lnco2 interaction2 no_agree_costs high_level contact_group lngdp, i(i_bargain bali) link(logit) 


preserve

set seed 339487731

drawnorm SN_b1-SN_b11, n(10000) means(e(b)) cov(e(V)) clear

postutil clear
postfile mypost prob_hat0 lo0 hi0 prob_hat1 lo1 hi1 diff_hat diff_lo diff_hi using "McKibben_BJPS_Linking_interim", replace

noisily display "start"

local a= 0
while `a' <= 15.8551 { 

{

*CALCULATE USING MAX OF DIFF-VALUED AGENDA
scalar h_diff_valued_agenda =  .029419
scalar h_lnco2 =  9.331776
scalar h_high_level = 1
scalar h_contact_group = 0
scalar h_lngdp =  23.85371
scalar h_ln_gdp_capita =  8.342014
scalar h_bali = 0
scalar h_constant = 1 
scalar h_no_agree_costs = 23.94389       


    generate x_betahat0  = SN_b1*h_diff_valued_agenda + SN_b2*`a' + SN_b3*h_diff_valued_agenda*`a' + SN_b4*h_no_agree_costs + SN_b5*h_high_level + SN_b6*h_contact_group + SN_b7*h_lngdp + SN_b8*h_constant
                            
    generate x_betahat1  = SN_b1*(h_diff_valued_agenda+1) + SN_b2*`a' + SN_b3*(h_diff_valued_agenda+1)*`a' + SN_b4*h_no_agree_costs + SN_b5*h_high_level + SN_b6*h_contact_group + SN_b7*h_lngdp + SN_b8*h_constant
    
    gen prob0 =normal(x_betahat0)
    gen prob1=normal(x_betahat1)
    gen diff=prob1-prob0
    
    egen probhat0 =mean(prob0)
    egen probhat1=mean(prob1)
    egen diffhat=mean(diff)
    
    tempname prob_hat0 lo0 hi0 prob_hat1 lo1 hi1 diff_hat diff_lo diff_hi 

    _pctile prob0, p(2.5,97.5) 
    scalar `lo0' = r(r1)
    scalar `hi0' = r(r2)  
    
    _pctile prob1, p(2.5,97.5)
    scalar `lo1'= r(r1)
    scalar `hi1'= r(r2) 
    
    _pctile diff, p(2.5,97.5)
    scalar `diff_lo'= r(r1)
    scalar `diff_hi'= r(r2) 
   
    scalar `prob_hat0'=probhat0
    scalar `prob_hat1'=probhat1
    scalar `diff_hat'=diffhat
    
    post mypost (`prob_hat0') (`lo0') (`hi0') (`prob_hat1') (`lo1') (`hi1') (`diff_hat') (`diff_lo') (`diff_hi')
   
    }      
    
    drop x_betahat0 x_betahat1 prob0 prob1 diff probhat0 probhat1 diffhat
    
    local a=(`a')+ 1
    
    display "." _c
    
} 

display ""
postclose mypost

restore
merge using McKibben_BJPS_Linking_interim

gen yline=0
gen MV = (0 + (_n-1))
replace MV=. if (0 + (_n-1)) > 15.8551

line diff_hat MV, clwidth(thick) ///
         || line diff_lo MV, lpattern(dash) lcolor(gray) clwidth(thin) ///
         || line diff_hi MV, lpattern(dash) lcolor(gray) clwidth(thin)  ///
	     ||  line yline MV, lpattern(solid) lcolor(black) clwidth(thin) ///
         ||  , ///
           xlabel(3.93843 15.8551, nogrid labsize(2)) ///
           ylabel(, axis(1) nogrid labsize(2)) ///
		   legend(off) ///
           xtitle("CO2 Emissions", size(2.5)) ///
           ytitle("Marginal Effect of a 1-unit increase in the IL Structure", size(2.5)) ///
		   title("Bargaining Over More Differently Valued Issues") ///
		   subtitle("The Marginal Effect at Different Levels of CO2 Emissions") ///
		   note("*Note: 95% Confidence Intervals Indicated", size(2)) ///
           xsca(titlegap(4)) ysca(titlegap(4))  ///
           scheme(s2mono) graphregion(fcolor(white) ilcolor(white) lcolor(white)) ///