******************************* *Millenials and Foreign Policy* ******************************* ******************************* ***********Alex Kirss********** ******************************* clear cd "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/Figures" use "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/CCES2016_Feb2018.dta" *************** *Data Cleaning* *************** gen millennial =. replace millennial = 0 if birthyr !=. replace millennial = 1 if birthyr ==1981 replace millennial = 1 if birthyr ==1982 replace millennial = 1 if birthyr ==1983 replace millennial = 1 if birthyr ==1984 replace millennial = 1 if birthyr ==1985 replace millennial = 1 if birthyr ==1986 replace millennial = 1 if birthyr ==1987 replace millennial = 1 if birthyr ==1988 replace millennial = 1 if birthyr ==1989 replace millennial = 1 if birthyr ==1990 replace millennial = 1 if birthyr ==1991 replace millennial = 1 if birthyr ==1992 replace millennial = 1 if birthyr ==1993 replace millennial = 1 if birthyr ==1994 replace millennial = 1 if birthyr ==1995 replace millennial = 1 if birthyr ==1996 gen millennial2 =. replace millennial2 = 0 if birthyr !=. replace millennial2 = 1 if millennial == 1 replace millennial2 = 1 if birthyr == 1997 gen female =. replace female = 1 if gender == 2 replace female = 0 if gender ==1 gen race2 = race replace race2 =. if race == 6 replace race2 =. if race == 7 gen party =. replace party =0 if pid3 ==3 replace party =1 if pid3 ==1 replace party =2 if pid3 ==2 gen bornagain =. replace bornagain=0 if pew_bornagain ==2 replace bornagain=1 if pew_bornagain ==1 gen faminc2 = faminc replace faminc2=. if faminc==97 replace faminc2=. if faminc==98 gen veteran = 0 replace veteran = 1 if milstat_1==1 *Syria Questions: "As you may know, there are on-going conflicts in Syria and Iraq led by the organization *Islamic State of Iraq and Syria (or ISIS). What do you think the United States should do * in response to ISIS? (Check all that apply)" gen syria1 =. replace syria1 = 0 if CC16_312_1==2 replace syria1 = 1 if CC16_312_1==1 gen syria2 =. replace syria2 = 0 if CC16_312_2==2 replace syria2 = 1 if CC16_312_2==1 gen syria3 =. replace syria3 = 0 if CC16_312_3==2 replace syria3 = 1 if CC16_312_3==1 gen syria4 =. replace syria4 = 0 if CC16_312_4==2 replace syria4 = 1 if CC16_312_4==1 gen syria5 =. replace syria5 = 0 if CC16_312_5==2 replace syria5 = 1 if CC16_312_5==1 gen syria6 =. replace syria6 = 0 if CC16_312_6==2 replace syria6 = 1 if CC16_312_6==1 gen syria7 =. replace syria7 = 0 if CC16_312_7==2 replace syria7 = 1 if CC16_312_7==1 *Use of milforce Questions gen milforce1 =. replace milforce1 = 0 if CC16_414_1==2 replace milforce1 = 1 if CC16_414_1==1 gen milforce2 =. replace milforce2 = 0 if CC16_414_2==2 replace milforce2 = 1 if CC16_414_2==1 gen milforce3 =. replace milforce3 = 0 if CC16_414_3==2 replace milforce3 = 1 if CC16_414_3==1 gen milforce4 =. replace milforce4 = 0 if CC16_414_4==2 replace milforce4 = 1 if CC16_414_4==1 gen milforce5 =. replace milforce5 = 0 if CC16_414_5==2 replace milforce5 = 1 if CC16_414_5==1 gen milforce6 =. replace milforce6 = 0 if CC16_414_6==2 replace milforce6 = 1 if CC16_414_6==1 gen milforce7 =. replace milforce7 = 0 if CC16_414_7==2 replace milforce7 = 1 if CC16_414_7==1 save "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/CCES2016_2.dta", replace ********************* *Simple Syria Models* ********************* clear cd "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/Figures" use "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/CCES2016_2.dta" * U.S. should "not get involved" in Syria eststo: qui logit syria1 millennial female i.educ i.race2 i.party bornagain i.faminc2 margins, dydx(millennial) post parmest, label list(parm estimate min* max* p) saving(syriamarg1, replace) * U.S. should "Send food, medicine and other aid to countries affected" in Syria eststo: qui logit syria2 millennial female i.educ i.race2 i.party bornagain i.faminc2 margins, dydx(millennial) post parmest, label list(parm estimate min* max* p) saving(syriamarg2, replace) * U.S. should "Provide arms to those opposing ISIS" in Syria eststo: qui logit syria3 millennial female i.educ i.race2 i.party bornagain i.faminc2 margins, dydx(millennial) post parmest, label list(parm estimate min* max* p) saving(syriamarg3, replace) * U.S. should "Enforce a no-fly zone" in Syria eststo: qui logit syria4 millennial female i.educ i.race2 i.party bornagain i.faminc2 margins, dydx(millennial) post parmest, label list(parm estimate min* max* p) saving(syriamarg4, replace) * U.S. should "Use drones and aircraft to bomb ISIS troops" in Syria eststo: qui logit syria5 millennial female i.educ i.race2 i.party bornagain i.faminc2 margins, dydx(millennial) post parmest, label list(parm estimate min* max* p) saving(syriamarg5, replace) * U.S. should "Send military support staff(non-combat)" in Syria eststo: qui logit syria6 millennial female i.educ i.race2 i.party bornagain i.faminc2 margins, dydx(millennial) post parmest, label list(parm estimate min* max* p) saving(syriamarg6, replace) * U.S. should "Send significant force to fight ISIS" in Syria eststo: qui logit syria7 millennial female i.educ i.race2 i.party bornagain i.faminc2 margins, dydx(millennial) post parmest, label list(parm estimate min* max* p) saving(syriamarg7, replace) esttab using syria.csv, replace eststo clear clear use "syriamarg1.dta" append using syriamarg2.dta append using syriamarg3.dta append using syriamarg4.dta append using syriamarg5.dta append using syriamarg6.dta append using syriamarg7.dta replace label = "1" in 1 replace label = "2" in 2 replace label = "3" in 3 replace label = "4" in 4 replace label = "5" in 5 replace label = "6" in 6 replace label = "7" in 7 destring label, generate(label2) eclplot estimate min95 max95 label2, /// rplottype(rspike) estopts(m(O) mc(gs0)) ciopts(lc(gs0)) /// hori xline(0, lpattern(shortdash) lc(gs8)) /// ylabel(1 "Not Involved" 2 "Send Food/Medicine" 3 "Provide Arms to Opposition" 4 "No-Fly Zone" 5 "Bomb ISIS" 6 "Send Non-Combat Staff" 7 "Send Significant Forces") /// xlab(, nogrid) ylab(, nogrid) /// ytitle("Syria Policy") xtitle(Marginal Effect (AME)) /// title("Figure 1: Millennial Support for Syria Policy") caption("Bars represent 95% confidence intervals", size(vsmall)) graph export "Figure1.png", replace graph export "Figure1.pdf", replace ************************************* *Simple Use of Military Force Models* ************************************* clear cd "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/Figures" use "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/CCES2016_2.dta" *Approve of "Ensuring Supply of Oil" eststo: qui logit milforce1 millennial female i.educ i.race2 i.party bornagain i.faminc2 margins, dydx(millennial) post parmest, label list(parm estimate min* max* p) saving(milforcemarg1, replace) *Approve of "Destroying a Terrorist Camp" eststo: qui logit milforce2 millennial female i.educ i.race2 i.party bornagain i.faminc2 margins, dydx(millennial) post parmest, label list(parm estimate min* max* p) saving(milforcemarg2, replace) *Approve of "Intervening in a Region Where There is Genocide of a Civil War" eststo: qui logit milforce3 millennial female i.educ i.race2 i.party bornagain i.faminc2 margins, dydx(millennial) post parmest, label list(parm estimate min* max* p) saving(milforcemarg3, replace) *Approve of "Assisting the Spread of Democracy" eststo: qui logit milforce4 millennial female i.educ i.race2 i.party bornagain i.faminc2 margins, dydx(millennial) post parmest, label list(parm estimate min* max* p) saving(milforcemarg4, replace) *Approve of "Protecting American Allies under attack by foreign nations" eststo: qui logit milforce5 millennial female i.educ i.race2 i.party bornagain i.faminc2 margins, dydx(millennial) post parmest, label list(parm estimate min* max* p) saving(milforcemarg5, replace) *Approve of "Helping the United Nations Upholding International Law" eststo: qui logit milforce6 millennial female i.educ i.race2 i.party bornagain i.faminc2 margins, dydx(millennial) post parmest, label list(parm estimate min* max* p) saving(milforcemarg6, replace) *Approve of "None of the Above" eststo: qui logit milforce7 millennial female i.educ i.race2 i.party bornagain i.faminc2 margins, dydx(millennial) post parmest, label list(parm estimate min* max* p) saving(milforcemarg7, replace) esttab using milforce.csv, replace eststo clear clear use "milforcemarg1.dta" append using milforcemarg2.dta append using milforcemarg3.dta append using milforcemarg4.dta append using milforcemarg5.dta append using milforcemarg6.dta append using milforcemarg7.dta replace label = "1" in 1 replace label = "2" in 2 replace label = "3" in 3 replace label = "4" in 4 replace label = "5" in 5 replace label = "6" in 6 replace label = "7" in 7 destring label, generate(label2) eclplot estimate min95 max95 label2, /// rplottype(rspike) estopts(m(O) mc(gs0)) ciopts(lc(gs0)) /// hori xline(0, lpattern(shortdash) lc(gs8)) /// ylabel(1 "Ensure Oil Supply" 2 "Destroy Terrorist Camp" 3 "Intervene Genocide/Civil War" 4 "Assist Spread of Democracy" 5 "Protect Allies" 6 "Help UN" 7 "None of The Above") /// xlab(, nogrid) ylab(, nogrid) /// ytitle("Use of Force") xtitle(Marginal Effect (AME)) /// title("Figure 2: Millennial Support for Use of Military Force") caption("Bars represent 95% confidence intervals", size(vsmall)) graph export "Figure2.png", replace graph export "Figure2.pdf", replace ************************************** *Interactive Male/Female Syria Models* ************************************** clear cd "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/Figures" use "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/CCES2016_2.dta" set matsize 800 * Q1: U.S. should "not get involved" in Syria ********************************************** eststo: qui logit syria1 i.millennial#i.female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post mplotoffset, offset(0.1) horiz recast(scatter) recastci(rspike) plotop(m(O)) /// xlab(, nogrid) ylab(, nogrid) ytitle("Millennial") xtitle("Marginal Mean") title("Not Get Involved") /// legend(order(1 "Male" 2 "Female")) name(femsyria1, replace) * Compare male non-millennial to female millennial contrast {millennial#female 1 0 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(syria1con1, replace) * Can use the following code to check the comparison groups: matrix list r(L) * Compare male non-millennial to male millennial...need to run whole model again, ugh qui logit syria1 i.millennial#i.female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 0 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(syria1con2, replace) * Compare male non-millennial to female non-millennial qui logit syria1 i.millennial#i.female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 -1 0 0}, effects post parmest, label list(parm estimate min* max* p) saving(syria1con3, replace) * Compare female non-millennial to female millennial qui logit syria1 i.millennial#i.female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(syria1con4, replace) * Compare female non-millennial to male millennial qui logit syria1 i.millennial#i.female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(syria1con5, replace) * Compare male millennial to female millennial qui logit syria1 i.millennial#i.female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 0 1 -1}, effects post parmest, label list(parm estimate min* max* p) saving(syria1con6, replace) * Q2: U.S. should "Send food, medicine and other aid to countries affected" in Syria *************************************************************************************** eststo: qui logit syria2 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post mplotoffset, offset(0.1) horiz recast(scatter) recastci(rspike) plotop(m(O)) /// xlab(, nogrid) ylab(, nogrid) ytitle("Millennial") xtitle("Marginal Mean") title("Send Food and Aid") /// legend(order(1 "Male" 2 "Female")) name(femsyria2, replace) * Compare male non-millennial to female millennial contrast {millennial#female 1 0 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(syria2con1, replace) * Compare male non-millennial to male millennial qui logit syria2 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 0 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(syria2con2, replace) * Compare male non-millennial to female non-millennial qui logit syria2 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 -1 0 0}, effects post parmest, label list(parm estimate min* max* p) saving(syria2con3, replace) * Compare female non-millennial to female millennial qui logit syria2 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(syria2con4, replace) * Compare female non-millennial to male millennial qui logit syria2 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(syria2con5, replace) * Compare male millennial to female millennial qui logit syria2 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 0 1 -1}, effects post parmest, label list(parm estimate min* max* p) saving(syria2con6, replace) * Q3: U.S. should "Provide arms to those opposing ISIS" in Syria ****************************************************************** eststo: qui logit syria3 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post mplotoffset, offset(0.1) horiz recast(scatter) recastci(rspike) plotop(m(O)) /// xlab(, nogrid) ylab(, nogrid) ytitle("Millennial") xtitle("Marginal Mean") title("Arm Opposition") /// legend(order(1 "Male" 2 "Female")) name(femsyria3, replace) * Compare male non-millennial to female millennial contrast {millennial#female 1 0 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(syria3con1, replace) * Compare male non-millennial to male millennial qui logit syria3 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 0 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(syria3con2, replace) * Compare male non-millennial to female non-millennial qui logit syria3 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 -1 0 0}, effects post parmest, label list(parm estimate min* max* p) saving(syria3con3, replace) * Compare female non-millennial to female millennial qui logit syria3 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(syria3con4, replace) * Compare female non-millennial to male millennial qui logit syria3 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(syria3con5, replace) * Compare male millennial to female millennial qui logit syria3 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 0 1 -1}, effects post parmest, label list(parm estimate min* max* p) saving(syria3con6, replace) * Q4: U.S. should "Enforce a no-fly zone" in Syria ******************************************************* eststo: qui logit syria4 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post mplotoffset, offset(0.1) horiz recast(scatter) recastci(rspike) plotop(m(O)) /// xlab(, nogrid) ylab(, nogrid) ytitle("Millennial") xtitle("Marginal Mean") title("No Fly Zone") /// legend(order(1 "Male" 2 "Female")) name(femsyria4, replace) * Compare male non-millennial to female millennial contrast {millennial#female 1 0 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(syria4con1, replace) * Compare male non-millennial to male millennial qui logit syria4 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 0 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(syria4con2, replace) * Compare male non-millennial to female non-millennial qui logit syria4 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 -1 0 0}, effects post parmest, label list(parm estimate min* max* p) saving(syria4con3, replace) * Compare female non-millennial to female millennial qui logit syria4 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(syria4con4, replace) * Compare female non-millennial to male millennial qui logit syria4 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(syria4con5, replace) * Compare male millennial to female millennial qui logit syria4 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 0 1 -1}, effects post parmest, label list(parm estimate min* max* p) saving(syria4con6, replace) * Q5: U.S. should "Use drones and aircraft to bomb ISIS troops" in Syria ************************************************************************** eststo: qui logit syria5 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post mplotoffset, offset(0.1) horiz recast(scatter) recastci(rspike) plotop(m(O)) /// xlab(, nogrid) ylab(, nogrid) ytitle("Millennial") xtitle("Marginal Mean") title("Bomb ISIS") /// legend(order(1 "Male" 2 "Female")) name(femsyria5, replace) * Compare male non-millennial to female millennial contrast {millennial#female 1 0 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(syria5con1, replace) * Compare male non-millennial to male millennial qui logit syria5 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 0 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(syria5con2, replace) * Compare male non-millennial to female non-millennial qui logit syria5 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 -1 0 0}, effects post parmest, label list(parm estimate min* max* p) saving(syria5con3, replace) * Compare female non-millennial to female millennial qui logit syria5 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(syria5con4, replace) * Compare female non-millennial to male millennial qui logit syria5 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(syria5con5, replace) * Compare male millennial to female millennial qui logit syria5 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 0 1 -1}, effects post parmest, label list(parm estimate min* max* p) saving(syria5con6, replace) * Q6: U.S. should "Send military support staff(non-combat)" in Syria ************************************************************************* eststo: qui logit syria6 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post mplotoffset, offset(0.1) horiz recast(scatter) recastci(rspike) plotop(m(O)) /// xlab(, nogrid) ylab(, nogrid) ytitle("Millennial") xtitle("Marginal Mean") title("Send Non-Combat Staff") /// legend(order(1 "Male" 2 "Female")) name(femsyria6, replace) * Compare male non-millennial to female millennial contrast {millennial#female 1 0 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(syria6con1, replace) * Compare male non-millennial to male millennial qui logit syria6 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 0 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(syria6con2, replace) * Compare male non-millennial to female non-millennial qui logit syria6 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 -1 0 0}, effects post parmest, label list(parm estimate min* max* p) saving(syria6con3, replace) * Compare female non-millennial to female millennial qui logit syria6 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(syria6con4, replace) * Compare female non-millennial to male millennial qui logit syria6 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(syria6con5, replace) * Compare male millennial to female millennial qui logit syria6 millennial#female i.race2 i.educ i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 0 1 -1}, effects post parmest, label list(parm estimate min* max* p) saving(syria6con6, replace) * Q7: U.S. should "Send significant force to fight ISIS" in Syria ******************************************************************** eststo: qui logit syria7 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post mplotoffset, offset(0.1) horiz recast(scatter) recastci(rspike) plotop(m(O)) /// xlab(, nogrid) ylab(, nogrid) ytitle("Millennial") xtitle("Marginal Mean") title("Send Combat Force") /// legend(order(1 "Male" 2 "Female")) name(femsyria7, replace) * Compare male non-millennial to female millennial contrast {millennial#female 1 0 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(syria7con1, replace) * Compare male non-millennial to male millennial qui logit syria7 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 0 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(syria7con2, replace) * Compare male non-millennial to female non-millennial qui logit syria7 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 -1 0 0}, effects post parmest, label list(parm estimate min* max* p) saving(syria7con3, replace) * Compare female non-millennial to female millennial qui logit syria7 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(syria7con4, replace) * Compare female non-millennial to male millennial qui logit syria7 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(syria7con5, replace) * Compare male millennial to female millennial qui logit syria7 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 0 1 -1}, effects post parmest, label list(parm estimate min* max* p) saving(syria7con6, replace) gr combine femsyria1 femsyria2 femsyria3 femsyria4 femsyria5 femsyria6 femsyria7, rows(3) title("Figure 4: Male vs. Female Millennials on Syria Policy ") /// caption("Bars represent 95% confidence intervals", size(vsmall)) graph export "Figure4.png", replace graph export "Figure4.pdf", replace * esttab using femsyria.csv, replace eststo clear * Build Syria Contrast Plots ****************************** clear cd "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/Figures" use "syria1con1.dta" append using syria1con2.dta append using syria1con3.dta append using syria1con4.dta append using syria1con5.dta append using syria1con6.dta replace label = "1" in 1 replace label = "2" in 2 replace label = "3" in 3 replace label = "4" in 4 replace label = "5" in 5 replace label = "6" in 6 destring label, generate(label2) eclplot estimate min95 max95 label2, /// rplottype(rspike) estopts(m(O) mc(gs0)) ciopts(lc(gs0)) /// hori xline(0, lpattern(shortdash) lc(gs8)) /// xlab(, nogrid) ylab(, nogrid) /// ylabel(1 "MNM vs. FM" 2 "MNM vs. MM" 3 "MNM vs. FNM" 4 "FNM vs. FM" 5 "FNM vs. MM" 6 "MM vs. FM") /// ytitle("Comparison") xtitle(Contrast) /// title("Not Get Involved") name(syria1con, replace) clear cd "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/Figures" use "syria2con1.dta" append using syria2con2.dta append using syria2con3.dta append using syria2con4.dta append using syria2con5.dta append using syria2con6.dta replace label = "1" in 1 replace label = "2" in 2 replace label = "3" in 3 replace label = "4" in 4 replace label = "5" in 5 replace label = "6" in 6 destring label, generate(label2) eclplot estimate min95 max95 label2, /// rplottype(rspike) estopts(m(O) mc(gs0)) ciopts(lc(gs0)) /// hori xline(0, lpattern(shortdash) lc(gs8)) /// xlab(, nogrid) ylab(, nogrid) /// ylabel(1 "MNM vs. FM" 2 "MNM vs. MM" 3 "MNM vs. FNM" 4 "FNM vs. FM" 5 "FNM vs. MM" 6 "MM vs. FM") /// ytitle("Comparison") xtitle(Contrast) /// title("Send Food and Aid") name(syria2con, replace) clear cd "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/Figures" use "syria3con1.dta" append using syria3con2.dta append using syria3con3.dta append using syria3con4.dta append using syria3con5.dta append using syria3con6.dta replace label = "1" in 1 replace label = "2" in 2 replace label = "3" in 3 replace label = "4" in 4 replace label = "5" in 5 replace label = "6" in 6 destring label, generate(label2) eclplot estimate min95 max95 label2, /// rplottype(rspike) estopts(m(O) mc(gs0)) ciopts(lc(gs0)) /// hori xline(0, lpattern(shortdash) lc(gs8)) /// xlab(, nogrid) ylab(, nogrid) /// ylabel(1 "MNM vs. FM" 2 "MNM vs. MM" 3 "MNM vs. FNM" 4 "FNM vs. FM" 5 "FNM vs. MM" 6 "FM vs. MM") /// ytitle("Comparison") xtitle(Contrast) /// title("Arm Opposition") name(syria3con, replace) clear cd "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/Figures" use "syria4con1.dta" append using syria4con2.dta append using syria4con3.dta append using syria4con4.dta append using syria4con5.dta append using syria4con6.dta replace label = "1" in 1 replace label = "2" in 2 replace label = "3" in 3 replace label = "4" in 4 replace label = "5" in 5 replace label = "6" in 6 destring label, generate(label2) eclplot estimate min95 max95 label2, /// rplottype(rspike) estopts(m(O) mc(gs0)) ciopts(lc(gs0)) /// hori xline(0, lpattern(shortdash) lc(gs8)) /// xlab(, nogrid) ylab(, nogrid) /// ylabel(1 "MNM vs. FM" 2 "MNM vs. MM" 3 "MNM vs. FNM" 4 "FNM vs. FM" 5 "FNM vs. MM" 6 "MM vs. FM") /// ytitle("Comparison") xtitle(Contrast) /// title("No Fly Zone") name(syria4con, replace) clear cd "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/Figures" use "syria5con1.dta" append using syria5con2.dta append using syria5con3.dta append using syria5con4.dta append using syria5con5.dta append using syria5con6.dta replace label = "1" in 1 replace label = "2" in 2 replace label = "3" in 3 replace label = "4" in 4 replace label = "5" in 5 replace label = "6" in 6 destring label, generate(label2) eclplot estimate min95 max95 label2, /// rplottype(rspike) estopts(m(O) mc(gs0)) ciopts(lc(gs0)) /// hori xline(0, lpattern(shortdash) lc(gs8)) /// xlab(, nogrid) ylab(, nogrid) /// ylabel(1 "MNM vs. FM" 2 "MNM vs. MM" 3 "MNM vs. FNM" 4 "FNM vs. FM" 5 "FNM vs. MM" 6 "MM vs. FM") /// ytitle("Comparison") xtitle(Contrast) /// title("Bomb ISIS") name(syria5con, replace) clear cd "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/Figures" use "syria6con1.dta" append using syria6con2.dta append using syria6con3.dta append using syria6con4.dta append using syria6con5.dta append using syria6con6.dta replace label = "1" in 1 replace label = "2" in 2 replace label = "3" in 3 replace label = "4" in 4 replace label = "5" in 5 replace label = "6" in 6 destring label, generate(label2) eclplot estimate min95 max95 label2, /// rplottype(rspike) estopts(m(O) mc(gs0)) ciopts(lc(gs0)) /// hori xline(0, lpattern(shortdash) lc(gs8)) /// xlab(, nogrid) ylab(, nogrid) /// ylabel(1 "MNM vs. FM" 2 "MNM vs. MM" 3 "MNM vs. FNM" 4 "FNM vs. FM" 5 "FNM vs. MM" 6 "MM vs. FM") /// ytitle("Comparison") xtitle(Contrast) /// title("Send Non-Combat Staff") name(syria6con, replace) clear cd "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/Figures" use "syria7con1.dta" append using syria7con2.dta append using syria7con3.dta append using syria7con4.dta append using syria7con5.dta append using syria7con6.dta replace label = "1" in 1 replace label = "2" in 2 replace label = "3" in 3 replace label = "4" in 4 replace label = "5" in 5 replace label = "6" in 6 destring label, generate(label2) eclplot estimate min95 max95 label2, /// rplottype(rspike) estopts(m(O) mc(gs0)) ciopts(lc(gs0)) /// hori xline(0, lpattern(shortdash) lc(gs8)) /// xlab(, nogrid) ylab(, nogrid) /// ylabel(1 "MNM vs. FM" 2 "MNM vs. MM" 3 "MNM vs. FNM" 4 "FNM vs. FM" 5 "FNM vs. MM" 6 "MM vs. FM") /// ytitle("Comparison") xtitle(Contrast) /// title("Send Combat Forces") name(syria7con, replace) gr combine syria1con syria2con syria3con syria4con syria5con syria6con syria7con, rows(3) /// title("Figure 5: Differences in Marginal Means Between Gender Generation Groups, Syria Policy") /// caption("MNM: Male Non-Millennial, FNM: Female Non-Millennial, MM: Male Millennial, FM: Female Millennial; Bars represent 95% confidence intervals", size(vsmall)) graph export "Figure5.png", replace graph export "Figure5.pdf", replace ****************************************************** *Interactive Male/Female Use of Military Force Models* ****************************************************** clear cd "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/Figures" use "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/CCES2016_2.dta" set matsize 800 * Q1: Approve of "Ensuring Supply of Oil" ******************************************* eststo: qui logit milforce1 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post mplotoffset, offset(0.1) horiz recast(scatter) recastci(rspike) plotop(m(O)) /// xlab(, nogrid) ylab(, nogrid) ytitle("Millennial") xtitle("Marginal Mean") title("Ensure Oil Supply") /// legend(order(1 "Male" 2 "Female")) name(femmil1, replace) * Compare male non-millennial to female millennial contrast {millennial#female 1 0 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(milforce1con1, replace) * Compare male non-millennial to male millennial qui logit milforce1 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 0 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(milforce1con2, replace) * Compare male non-millennial to female non-millennial qui logit milforce1 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 -1 0 0}, effects post parmest, label list(parm estimate min* max* p) saving(milforce1con3, replace) * Compare female non-millennial to female millennial qui logit milforce1 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(milforce1con4, replace) * Compare female non-millennial to male millennial qui logit milforce1 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(milforce1con5, replace) * Compare male millennial to female millennial qui logit milforce1 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 0 1 -1}, effects post parmest, label list(parm estimate min* max* p) saving(milforce1con6, replace) * Q2: Approve of "Destroying a Terrorist Camp" ************************************************ eststo: qui logit milforce2 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post mplotoffset, offset(0.1) horiz recast(scatter) recastci(rspike) plotop(m(O)) /// xlab(, nogrid) ylab(, nogrid) ytitle("Millennial") xtitle("Marginal Mean") title("Destroy Terrorist Camp") /// legend(order(1 "Male" 2 "Female")) name(femmil2, replace) * Compare male non-millennial to female millennial contrast {millennial#female 1 0 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(milforce2con1, replace) * Compare male non-millennial to male millennial qui logit milforce2 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 0 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(milforce2con2, replace) * Compare male non-millennial to female non-millennial qui logit milforce2 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 -1 0 0}, effects post parmest, label list(parm estimate min* max* p) saving(milforce2con3, replace) * Compare female non-millennial to female millennial qui logit milforce2 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(milforce2con4, replace) * Compare female non-millennial to male millennial qui logit milforce2 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(milforce2con5, replace) * Compare male millennial to female millennial qui logit milforce2 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 0 1 -1}, effects post parmest, label list(parm estimate min* max* p) saving(milforce2con6, replace) * Q3: Approve of "Intervening in a Region Where There is Genocide of a Civil War" *********************************************************************************** eststo: qui logit milforce3 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post mplotoffset, offset(0.1) horiz recast(scatter) recastci(rspike) plotop(m(O)) /// xlab(, nogrid) ylab(, nogrid) ytitle("Millennial") xtitle("Marginal Mean") title("Intervene in Genocide") /// legend(order(1 "Male" 2 "Female")) name(femmil3, replace) * Compare male non-millennial to female millennial contrast {millennial#female 1 0 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(milforce3con1, replace) * Compare male non-millennial to male millennial qui logit milforce3 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 0 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(milforce3con2, replace) * Compare male non-millennial to female non-millennial qui logit milforce3 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 -1 0 0}, effects post parmest, label list(parm estimate min* max* p) saving(milforce3con3, replace) * Compare female non-millennial to female millennial qui logit milforce3 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(milforce3con4, replace) * Compare female non-millennial to male millennial qui logit milforce3 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(milforce3con5, replace) * Compare male millennial to female millennial qui logit milforce3 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 0 1 -1}, effects post parmest, label list(parm estimate min* max* p) saving(milforce3con6, replace) * Q4: Approve of "Assisting the Spread of Democracy" ******************************************************* eststo: qui logit milforce4 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post mplotoffset, offset(0.1) horiz recast(scatter) recastci(rspike) plotop(m(O)) /// xlab(, nogrid) ylab(, nogrid) ytitle("Millennial") xtitle("Marginal Mean") title("Assist Spread of Democracy") /// legend(order(1 "Male" 2 "Female")) name(femmil4, replace) * Compare male non-millennial to female millennial contrast {millennial#female 1 0 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(milforce4con1, replace) * Compare male non-millennial to male millennial qui logit milforce4 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 0 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(milforce4con2, replace) * Compare male non-millennial to female non-millennial qui logit milforce4 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 -1 0 0}, effects post parmest, label list(parm estimate min* max* p) saving(milforce4con3, replace) * Compare female non-millennial to female millennial qui logit milforce4 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(milforce4con4, replace) * Compare female non-millennial to male millennial qui logit milforce4 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(milforce4con5, replace) * Compare male millennial to female millennial qui logit milforce4 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 0 1 -1}, effects post parmest, label list(parm estimate min* max* p) saving(milforce4con6, replace) * Q5: Approve of "Protecting American Allies under attack by foreign nations" ******************************************************************************* eststo: qui logit milforce5 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post mplotoffset, offset(0.1) horiz recast(scatter) recastci(rspike) plotop(m(O)) /// xlab(, nogrid) ylab(, nogrid) ytitle("Millennial") xtitle("Marginal Mean") title("Protect Allies") /// legend(order(1 "Male" 2 "Female")) name(femmil5, replace) * Compare male non-millennial to female millennial contrast {millennial#female 1 0 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(milforce5con1, replace) * Compare male non-millennial to male millennial qui logit milforce5 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 0 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(milforce5con2, replace) * Compare male non-millennial to female non-millennial qui logit milforce5 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 -1 0 0}, effects post parmest, label list(parm estimate min* max* p) saving(milforce5con3, replace) * Compare female non-millennial to female millennial qui logit milforce5 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(milforce5con4, replace) * Compare female non-millennial to male millennial qui logit milforce5 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(milforce5con5, replace) * Compare male millennial to female millennial qui logit milforce5 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 0 1 -1}, effects post parmest, label list(parm estimate min* max* p) saving(milforce5con6, replace) * Q6: Approve of "Helping the United Nations Upholding International Law" ************************************************************************** eststo: qui logit milforce6 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post mplotoffset, offset(0.1) horiz recast(scatter) recastci(rspike) plotop(m(O)) /// xlab(, nogrid) ylab(, nogrid) ytitle("Millennial") xtitle("Marginal Mean") title("Help UN") /// legend(order(1 "Male" 2 "Female")) name(femmil6, replace) * Compare male non-millennial to female millennial contrast {millennial#female 1 0 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(milforce6con1, replace) * Compare male non-millennial to male millennial qui logit milforce6 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 0 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(milforce6con2, replace) * Compare male non-millennial to female non-millennial qui logit milforce6 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 -1 0 0}, effects post parmest, label list(parm estimate min* max* p) saving(milforce6con3, replace) * Compare female non-millennial to female millennial qui logit milforce6 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(milforce6con4, replace) * Compare female non-millennial to male millennial qui logit milforce6 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(milforce6con5, replace) * Compare male millennial to female millennial qui logit milforce6 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 0 1 -1}, effects post parmest, label list(parm estimate min* max* p) saving(milforce6con6, replace) * Q7: Approve of "None of the Above" ************************************** eststo: qui logit milforce7 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post mplotoffset, offset(0.1) horiz recast(scatter) recastci(rspike) plotop(m(O)) /// xlab(, nogrid) ylab(, nogrid) ytitle("Millennial") xtitle("Marginal Mean") title("None of the Above") /// legend(order(1 "Male" 2 "Female")) name(femmil7, replace) * Compare male non-millennial to female millennial contrast {millennial#female 1 0 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(milforce7con1, replace) * Compare male non-millennial to male millennial qui logit milforce7 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 0 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(milforce7con2, replace) * Compare male non-millennial to female non-millennial qui logit milforce7 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 1 -1 0 0}, effects post parmest, label list(parm estimate min* max* p) saving(milforce7con3, replace) * Compare female non-millennial to female millennial qui logit milforce7 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 0 -1}, effects post parmest, label list(parm estimate min* max* p) saving(milforce7con4, replace) * Compare female non-millennial to male millennial qui logit milforce7 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 1 -1 0}, effects post parmest, label list(parm estimate min* max* p) saving(milforce7con5, replace) * Compare male millennial to female millennial qui logit milforce7 millennial#female i.educ i.race2 i.party bornagain i.faminc2 margins millennial#female, post contrast {millennial#female 0 0 1 -1}, effects post parmest, label list(parm estimate min* max* p) saving(milforce7con6, replace) gr combine femmil1 femmil2 femmil3 femmil4 femmil5 femmil6 femmil7, rows(3) title("Figure 6: Male vs. Female Millennials on Use of Military Force") /// caption("Bars represent 95% confidence intervals", size(vsmall)) graph export "Figure6.png", replace graph export "Figure6.pdf", replace esttab using femmilforce.csv, replace eststo clear * Build Military Force Contrast Plots ****************************** clear cd "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/Figures" use "milforce1con1.dta" append using milforce1con2.dta append using milforce1con3.dta append using milforce1con4.dta append using milforce1con5.dta append using milforce1con6.dta replace label = "1" in 1 replace label = "2" in 2 replace label = "3" in 3 replace label = "4" in 4 replace label = "5" in 5 replace label = "6" in 6 destring label, generate(label2) eclplot estimate min95 max95 label2, /// rplottype(rspike) estopts(m(O) mc(gs0)) ciopts(lc(gs0)) /// hori xline(0, lpattern(shortdash) lc(gs8)) /// xlab(, nogrid) ylab(, nogrid) /// ylabel(1 "MNM vs. FM" 2 "MNM vs. MM" 3 "MNM vs. FNM" 4 "FNM vs. FM" 5 "FNM vs. MM" 6 "MM vs. FM") /// ytitle("Comparison") xtitle(Contrast) /// title("Ensure Oil Supply") name(milforcecon1, replace) clear cd "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/Figures" use "milforce2con1.dta" append using milforce2con2.dta append using milforce2con3.dta append using milforce2con4.dta append using milforce2con5.dta append using milforce2con6.dta replace label = "1" in 1 replace label = "2" in 2 replace label = "3" in 3 replace label = "4" in 4 replace label = "5" in 5 replace label = "6" in 6 destring label, generate(label2) eclplot estimate min95 max95 label2, /// rplottype(rspike) estopts(m(O) mc(gs0)) ciopts(lc(gs0)) /// hori xline(0, lpattern(shortdash) lc(gs8)) /// xlab(, nogrid) ylab(, nogrid) /// ylabel(1 "MNM vs. FM" 2 "MNM vs. MM" 3 "MNM vs. FNM" 4 "FNM vs. FM" 5 "FNM vs. MM" 6 "MM vs. FM") /// ytitle("Comparison") xtitle(Contrast) /// title("Destroy Terrorist Camp") name(milforcecon2, replace) clear cd "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/Figures" use "milforce3con1.dta" append using milforce3con2.dta append using milforce3con3.dta append using milforce3con4.dta append using milforce3con5.dta append using milforce3con6.dta replace label = "1" in 1 replace label = "2" in 2 replace label = "3" in 3 replace label = "4" in 4 replace label = "5" in 5 replace label = "6" in 6 destring label, generate(label2) eclplot estimate min95 max95 label2, /// rplottype(rspike) estopts(m(O) mc(gs0)) ciopts(lc(gs0)) /// hori xline(0, lpattern(shortdash) lc(gs8)) /// xlab(, nogrid) ylab(, nogrid) /// ylabel(1 "MNM vs. FM" 2 "MNM vs. MM" 3 "MNM vs. FNM" 4 "FNM vs. FM" 5 "FNM vs. MM" 6 "MM vs. FM") /// ytitle("Comparison") xtitle(Contrast) /// title("Intervene Genocide/Civil War") name(milforcecon3, replace) clear cd "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/Figures" use "milforce4con1.dta" append using milforce4con2.dta append using milforce4con3.dta append using milforce4con4.dta append using milforce4con5.dta append using milforce4con6.dta replace label = "1" in 1 replace label = "2" in 2 replace label = "3" in 3 replace label = "4" in 4 replace label = "5" in 5 replace label = "6" in 6 destring label, generate(label2) eclplot estimate min95 max95 label2, /// rplottype(rspike) estopts(m(O) mc(gs0)) ciopts(lc(gs0)) /// hori xline(0, lpattern(shortdash) lc(gs8)) /// xlab(, nogrid) ylab(, nogrid) /// ylabel(1 "MNM vs. FM" 2 "MNM vs. MM" 3 "MNM vs. FNM" 4 "FNM vs. FM" 5 "FNM vs. MM" 6 "MM vs. FM") /// ytitle("Comparison") xtitle(Contrast) /// title("Assist Spread of Democracy") name(milforcecon4, replace) clear cd "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/Figures" use "milforce5con1.dta" append using milforce5con2.dta append using milforce5con3.dta append using milforce5con4.dta append using milforce5con5.dta append using milforce5con6.dta replace label = "1" in 1 replace label = "2" in 2 replace label = "3" in 3 replace label = "4" in 4 replace label = "5" in 5 replace label = "6" in 6 destring label, generate(label2) eclplot estimate min95 max95 label2, /// rplottype(rspike) estopts(m(O) mc(gs0)) ciopts(lc(gs0)) /// hori xline(0, lpattern(shortdash) lc(gs8)) /// xlab(, nogrid) ylab(, nogrid) /// ylabel(1 "MNM vs. FM" 2 "MNM vs. MM" 3 "MNM vs. FNM" 4 "FNM vs. FM" 5 "FNM vs. MM" 6 "MM vs. FM") /// ytitle("Comparison") xtitle(Contrast) /// title("Protect Allies") name(milforcecon5, replace) clear cd "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/Figures" use "milforce6con1.dta" append using milforce6con2.dta append using milforce6con3.dta append using milforce6con4.dta append using milforce6con5.dta append using milforce6con6.dta replace label = "1" in 1 replace label = "2" in 2 replace label = "3" in 3 replace label = "4" in 4 replace label = "5" in 5 replace label = "6" in 6 destring label, generate(label2) eclplot estimate min95 max95 label2, /// rplottype(rspike) estopts(m(O) mc(gs0)) ciopts(lc(gs0)) /// hori xline(0, lpattern(shortdash) lc(gs8)) /// xlab(, nogrid) ylab(, nogrid) /// ylabel(1 "MNM vs. FM" 2 "MNM vs. MM" 3 "MNM vs. FNM" 4 "FNM vs. FM" 5 "FNM vs. MM" 6 "MM vs. FM") /// ytitle("Comparison") xtitle(Contrast) /// title("Help UN") name(milforcecon6, replace) clear cd "/Users/alexanderkirss/Desktop/Ongoing Projects/MilennialsFP/CCES2016/Figures" use "milforce7con1.dta" append using milforce7con2.dta append using milforce7con3.dta append using milforce7con4.dta append using milforce7con5.dta append using milforce7con6.dta replace label = "1" in 1 replace label = "2" in 2 replace label = "3" in 3 replace label = "4" in 4 replace label = "5" in 5 replace label = "6" in 6 destring label, generate(label2) eclplot estimate min95 max95 label2, /// rplottype(rspike) estopts(m(O) mc(gs0)) ciopts(lc(gs0)) /// hori xline(0, lpattern(shortdash) lc(gs8)) /// xlab(, nogrid) ylab(, nogrid) /// ylabel(1 "MNM vs. FM" 2 "MNM vs. MM" 3 "MNM vs. FNM" 4 "FNM vs. FM" 5 "FNM vs. MM" 6 "FM vs. MM") /// ytitle("Comparison") xtitle(Contrast) /// title("None of the Above") name(milforcecon7, replace) gr combine milforcecon1 milforcecon2 milforcecon3 milforcecon4 milforcecon5 milforcecon6 milforcecon7, rows(3) /// title("Figure 7: Differences in Marginal Means Between Gender Generation Groups, Military Force") /// caption("MNM: Male Non-Millennial, FNM: Female Non-Millennial, MM: Male Millennial, FM: Female Millennial; Bars represent 95% confidence intervals", size(vsmall)) graph export "Figure7.png", replace graph export "Figure7.pdf", replace *Appendix *************** *INCOMPLETE *Other Preferences *Defense Spending *NEED TO RECODE gen cutdefense= . replace cutdefense= 0 if CC16_337_1 == *Iran Sanctions: Vote For (1) or Against (0) Iran Sanctions Act gen iransanctions =. replace iransanctions = 0 if CC16_351G==2 replace iransanctions = 1 if CC16_351G==1 logit iransanctions millennial female i.educ i.race2 i.party bornagain i.faminc2