clear all
capture log close
set more off
set logtype text


//install the following packages if not already done so
//fitstat.ado
//spost.ado
//pciplot.ado

log using ...

*****************
*** Kapitel 3 ***
*****************

//Use the EPRG MEP Survey 2010 dataset available online under http://www.lse.ac.uk/government/research/resgroups/EPRG/MEPsurveyData.aspx 

//modify variable names and labels

rename q1_1 MS
label var MS "Member State"
rename q1_2 PG
label var PG "EP Party Group"
rename q1_3 NP
label var NP "National Party"

rename q8_1_1 VoteInstructNP
rename q8_1_2 VoteInstructPG
rename q8_1_3 VoteInstructNDel
rename q8_1_4 VoteInstructComm
rename q8_1_5 VoteInstructCOM
rename q8_1_6 VoteInstructGovt
rename q8_1_7 VoteInstructEUIGs
rename q8_1_9 VoteInstructCitiz

rename q8_6_1 persInterest
rename q8_6_2 ConstInterest
rename q8_6_3 ProfExp
rename q8_6_4 CommCompet
rename q8_6_5 EPGDecision
rename q8_6_6 NPDecision
rename q8_6_7 PrevCommMemb


*** Abb. 3.1 Beweggründe für die Mitgliedschaft in einem Ausschuss ***

//select only PGs represented in my analysis of EP6

tab PG persInterest if PG<5 | PG==6 | PG==9 | PG==10
tab PG ProfExp if PG<5 | PG==6 | PG==9 | PG==10
tab PG PrevCommMemb if PG<5 | PG==6 | PG==9 | PG==10
tab PG CommCompet if PG<5 | PG==6 | PG==9 | PG==10
tab PG ConstInterest if PG<5 | PG==6 | PG==9 | PG==10
tab PG NPDecision if PG<5 | PG==6 | PG==9 | PG==10
tab PG EPGDecision if PG<5 | PG==6 | PG==9 | PG==10

//use trick to plot bars for frequencies of categories
//use count function for ID (b_1)

graph bar (count) b_1 if PG<5 | PG==6 | PG==9 | PG==10, over(persInterest, label(angle(45) labsize(small))) asyvars percent showyvars bargap(20) bar(1, color(gs6)) bar(2, color(gs6)) bar(3, color(gs6)) bar(4, color(gs6)) bar(5, color(gs6)) b1title("persönliches Interesse") ytitle("Prozent") legend(off) scheme(s2manual)
graph bar (count) b_1 if PG<5 | PG==6 | PG==9 | PG==10, over(ProfExp, label(angle(45) labsize(small))) asyvars percent showyvars bargap(20) bar(1, color(gs6)) bar(2, color(gs6)) bar(3, color(gs6)) bar(4, color(gs6)) bar(5, color(gs6)) b1title("fachliche Expertise") ytitle("Prozent") legend(off) scheme(s2manual)
graph bar (count) b_1 if PG<5 | PG==6 | PG==9 | PG==10, over(PrevCommMemb, label(angle(45) labsize(small))) asyvars percent showyvars bargap(20) bar(1, color(gs6)) bar(2, color(gs6)) bar(3, color(gs6)) bar(4, color(gs6)) bar(5, color(gs6)) b1title("Arbeitserfahrung in Ausschuss") ytitle("Prozent") legend(off) scheme(s2manual)
graph bar (count) b_1 if PG<5 | PG==6 | PG==9 | PG==10, over(CommCompet, label(angle(45) labsize(small))) asyvars percent showyvars bargap(20) bar(1, color(gs6)) bar(2, color(gs6)) bar(3, color(gs6)) bar(4, color(gs6)) bar(5, color(gs6)) b1title("Ausschusskompetenz") ytitle("Prozent") legend(off) scheme(s2manual)
graph bar (count) b_1 if PG<5 | PG==6 | PG==9 | PG==10, over(ConstInterest, label(angle(45) labsize(small))) asyvars percent showyvars bargap(20) bar(1, color(gs6)) bar(2, color(gs6)) bar(3, color(gs6)) bar(4, color(gs6)) bar(5, color(gs6)) b1title("Wahlkreisinteressen") ytitle("Prozent") legend(off) scheme(s2manual)
graph bar (count) b_1 if PG<5 | PG==6 | PG==9 | PG==10, over(NPDecision, label(angle(45) labsize(small))) asyvars percent showyvars bargap(20) bar(1, color(gs6)) bar(2, color(gs6)) bar(3, color(gs6)) bar(4, color(gs6)) bar(5, color(gs6)) b1title("Entscheidung nationale Partei") ytitle("Prozent") legend(off) scheme(s2manual)
graph bar (count) b_1 if PG<5 | PG==6 | PG==9 | PG==10, over(EPGDecision, label(angle(45) labsize(small))) asyvars percent showyvars bargap(20) bar(1, color(gs6)) bar(2, color(gs6)) bar(3, color(gs6)) bar(4, color(gs6)) bar(5, color(gs6)) b1title("Entscheidung EP-Fraktion") ytitle("Prozent") legend(off) scheme(s2manual)


*** Abb. 3.2 Abstimmungsempfehlungen nach Häufigkeit und Quelle ***


//select only PGs represented in my analysis of EP6

//Party Group
tab PG VoteInstructPG if PG<5 | PG==6 | PG==9 | PG==10

//National Delegation
tab PG VoteInstructNDel if PG<5 | PG==6 | PG==9 | PG==10

//Committee
tab PG VoteInstructComm if PG<5 | PG==6 | PG==9 | PG==10

//National Party
tab PG VoteInstructNP if PG<5 | PG==6 | PG==9 | PG==10

//Commission
tab PG VoteInstructCOM if PG<5 | PG==6 | PG==9 | PG==10

//National government
tab PG VoteInstructGovt if PG<5 | PG==6 | PG==9 | PG==10

//EU interest groups
tab PG VoteInstructEUIGs if PG<5 | PG==6 | PG==9 | PG==10

//citizens
tab PG VoteInstructCitiz if PG<5 | PG==6 | PG==9 | PG==10


//use same trick as above to plot bars for frequencies of categories
//use count function for ID (b_1)

graph bar (count) b_1 if PG<5 | PG==6 | PG==9 | PG==10, over(VoteInstructPG) asyvars percent showyvars bargap(20) bar(1, color(gs6)) bar(2, color(gs6)) bar(3, color(gs6)) bar(4, color(gs6)) bar(5, color(gs6)) ytitle("Prozent") b1title("EP-Fraktion", placement(south)) legend(off)  scheme(s2manual)
graph bar (count) b_1 if PG<5 | PG==6 | PG==9 | PG==10, over(VoteInstructNDel) asyvars percent showyvars bargap(20) bar(1, color(gs6)) bar(2, color(gs6)) bar(3, color(gs6)) bar(4, color(gs6)) bar(5, color(gs6)) ytitle("Prozent") b1title("nationale Delegation", placement(south)) legend(off)  scheme(s2manual)
graph bar (count) b_1 if PG<5 | PG==6 | PG==9 | PG==10, over(VoteInstructComm) asyvars percent showyvars bargap(20) bar(1, color(gs6)) bar(2, color(gs6)) bar(3, color(gs6)) bar(4, color(gs6)) bar(5, color(gs6)) ytitle("Prozent") b1title("Ausschuss", placement(south)) legend(off)  scheme(s2manual)
graph bar (count) b_1 if PG<5 | PG==6 | PG==9 | PG==10, over(VoteInstructNP) asyvars percent showyvars bargap(20) bar(1, color(gs6)) bar(2, color(gs6)) bar(3, color(gs6)) bar(4, color(gs6)) bar(5, color(gs6)) ytitle("Prozent") b1title("nationale Partei", placement(south)) legend(off)  scheme(s2manual)
graph bar (count) b_1 if PG<5 | PG==6 | PG==9 | PG==10, over(VoteInstructCOM) asyvars percent showyvars bargap(20) bar(1, color(gs6)) bar(2, color(gs6)) bar(3, color(gs6)) bar(4, color(gs6)) bar(5, color(gs6)) ytitle("Prozent") b1title("Kommission", placement(south)) legend(off)  scheme(s2manual)
graph bar (count) b_1 if PG<5 | PG==6 | PG==9 | PG==10, over(VoteInstructGovt) asyvars percent showyvars bargap(20) bar(1, color(gs6)) bar(2, color(gs6)) bar(3, color(gs6)) bar(4, color(gs6)) bar(5, color(gs6)) ytitle("Prozent") b1title("nationale Regierung", placement(south)) legend(off)  scheme(s2manual)
graph bar (count) b_1 if PG<5 | PG==6 | PG==9 | PG==10, over(VoteInstructEUIGs) asyvars percent showyvars bargap(20) bar(1, color(gs6)) bar(2, color(gs6)) bar(3, color(gs6)) bar(4, color(gs6)) bar(5, color(gs6)) ytitle("Prozent") b1title("EU-Interessengruppen", placement(south)) legend(off)  scheme(s2manual)
graph bar (count) b_1 if PG<5 | PG==6 | PG==9 | PG==10, over(VoteInstructCitiz) asyvars percent showyvars bargap(20) bar(1, color(gs6)) bar(2, color(gs6)) bar(3, color(gs6)) bar(4, color(gs6)) bar(5, color(gs6)) ytitle("Prozent") b1title("Bürger", placement(south)) legend(off)  scheme(s2manual)

*****************
*** Kapitel 5 ***
*****************

use ...DataDissWIDE.dta


***Abb. 5.1***
graph bar (sum) rcvReqPPE (sum) rcvReqPSE (sum) rcvReqALDE (sum) rcvReqVerts (sum) rcvReqGUE (sum) rcvReqUEN (sum) rcvReqIND (sum) rcvReqMEP, blabel(total)


***Tab. 5.3***

//für erste Lesung
tab rappOwnPPE result_dummy if pgAmPPE==1 & reading==0 & result_dummy!=2
tab rappOwnPSE result_dummy if pgAmPSE==1 & reading==0 & result_dummy!=2
tab rappOwnALDE result_dummy if pgAmALDE==1 & reading==0 & result_dummy!=2
tab rappOwnVerts result_dummy if pgAmVerts==1 & reading==0 & result_dummy!=2
tab rappOwnGUE result_dummy if pgAmGUE==1 & reading==0 & result_dummy!=2
tab rappOwnUEN result_dummy if pgAmUEN==1 & reading==0 & result_dummy!=2
tab rappOwnIND result_dummy if pgAmIND==1 & reading==0 & result_dummy!=2

//für zweite Lesung
tab rappOwnPPE result_dummy if pgAmPPE==1 & reading==1 & result_dummy!=2
tab rappOwnPSE result_dummy if pgAmPSE==1 & reading==1 & result_dummy!=2
tab rappOwnALDE result_dummy if pgAmALDE==1 & reading==1 & result_dummy!=2
tab rappOwnVerts result_dummy if pgAmVerts==1 & reading==1 & result_dummy!=2
tab rappOwnGUE result_dummy if pgAmGUE==1 & reading==1 & result_dummy!=2
tab rappOwnUEN result_dummy if pgAmUEN==1 & reading==1 & result_dummy!=2
tab rappOwnIND result_dummy if pgAmIND==1 & reading==1 & result_dummy!=2



*** Tab. 5.4. Determinanten von Plenaränderungsanträgen *** 
//
use ...COUNTS.dta


eststo, title("EVP-ED"): quietly nbreg copgAmPPE ipp reading amComStage absDisMedianD1 absDisMedianD2 rappOwnPPE margWinPPE policy_posPPE policy_salPPE if singleVote!=1, robust nolog
eststo, title("SPE"): quietly nbreg copgAmPSE ipp reading amComStage absDisMedianD1 absDisMedianD2 rappOwnPSE margWinPSE policy_posPSE policy_salPSE if singleVote!=1, robust nolog
eststo, title("ALDE"): quietly nbreg copgAmALDE ipp reading amComStage absDisMedianD1 absDisMedianD2 rappOwnALDE margWinALDE policy_posALDE policy_salALDE if singleVote!=1, robust nolog
eststo, title("Grüne/EFA"): quietly nbreg copgAmVerts ipp reading amComStage absDisMedianD1 absDisMedianD2 rappOwnVerts margWinVerts policy_posVerts policy_salVerts if singleVote!=1, robust nolog
eststo, title("GUE/NGL"): quietly nbreg copgAmGUE ipp reading amComStage absDisMedianD1 absDisMedianD2 rappOwnGUE margWinGUE policy_posGUE policy_salGUE if singleVote!=1, robust nolog
eststo, title("UEN"): quietly nbreg copgAmUEN ipp reading amComStage absDisMedianD1 absDisMedianD2 rappOwnUEN margWinUEN policy_posUEN policy_salUEN if singleVote!=1, robust nolog
eststo, title("IND/DEM"): quietly nbreg copgAmIND ipp reading amComStage absDisMedianD1 absDisMedianD2 rappOwnIND margWinIND policy_posIND policy_salIND if singleVote!=1, robust nolog 
esttab using tableESTCountsAbsDisMedian.booktabs, cells(b(star fmt(3)) se(par fmt(3))) stats(r2 bic N) starlevels(* 0.1 ** 0.05 *** 0.01) legend label booktabs title(Ergebnisse der negativ binomialen Count-Modelle\label{tab1}) alignment(D{.}{.}{-1}) compress replace

*** Abb. 5.2 und 5.3 *** 

//EVP-ED
quietly nbreg copgAmPPE ipp reading amComStage absDisMedianD1 absDisMedianD2 rappOwnPPE margWinPPE policy_posPPE policy_salPPE if singleVote!=1, robust nolog
prgen ipp, from(0) to (50) generate(pPPE) x(reading=0 rappOwnPPE=0) rest(mean) ci
prgen ipp, from(0) to (50) generate(pPPE1) x(reading=0 rappOwnPPE=1) rest(mean) ci
label var pPPEx "ipp"
label var pPPEp0 "nicht Berichterstatter"
label var pPPE1p0 "Berichterstatter"
label var pPPEp0lb "95%-Konfidenzintervall untere Grenze"
label var pPPEp0ub "95%-Konfidenzintervall obere Grenze"
label var pPPE1p0lb "95%-Konfidenzintervall untere Grenze"
label var pPPE1p0ub "95%-Konfidenzintervall obere Grenze"

//SPE
quietly nbreg copgAmPSE ipp reading amComStage absDisMedianD1 absDisMedianD2 rappOwnPSE margWinPSE policy_posPSE policy_salPSE if singleVote!=1, robust nolog
prgen ipp, from(0) to (50) generate(pPSE) x(reading=0 rappOwnPSE=0) rest(mean) ci
prgen ipp, from(0) to (50) generate(pPSE1) x(reading=0 rappOwnPSE=1) rest(mean) ci
label var pPSEx "ipp"
label var pPSEp0 "nicht Berichterstatter"
label var pPSE1p0 "Berichterstatter"
label var pPSEp0lb "95%-Konfidenzintervall untere Grenze"
label var pPSEp0ub "95%-Konfidenzintervall obere Grenze"
label var pPSE1p0lb "95%-Konfidenzintervall untere Grenze"
label var pPSE1p0ub "95%-Konfidenzintervall obere Grenze"

//ALDE
quietly nbreg copgAmALDE ipp reading amComStage absDisMedianD1 absDisMedianD2 rappOwnALDE margWinALDE policy_posALDE policy_salALDE if singleVote!=1, robust nolog
prgen ipp, from(0) to (50) generate(pALDE) x(reading=0 rappOwnALDE=0) rest(mean) ci
prgen ipp, from(0) to (50) generate(pALDE1) x(reading=0 rappOwnALDE=1) rest(mean) ci
label var pALDEx "ipp"
label var pALDEp0 "nicht Berichterstatter"
label var pALDE1p0 "Berichterstatter"
label var pALDEp0lb "95%-Konfidenzintervall untere Grenze"
label var pALDEp0ub "95%-Konfidenzintervall obere Grenze"
label var pALDE1p0lb "95%-Konfidenzintervall untere Grenze"
label var pALDE1p0ub "95%-Konfidenzintervall obere Grenze"

//Grüne/EFA
quietly nbreg copgAmVerts ipp reading amComStage absDisMedianD1 absDisMedianD2 rappOwnVerts margWinVerts policy_posVerts policy_salVerts if singleVote!=1, robust nolog
prgen ipp, from(0) to (50) generate(pVerts) x(reading=0 rappOwnVerts=0) rest(mean) ci
prgen ipp, from(0) to (50) generate(pVerts1) x(reading=0 rappOwnVerts=1) rest(mean) ci
label var pVertsx "ipp"
label var pVertsp0 "nicht Berichterstatter"
label var pVerts1p0 "Berichterstatter"
label var pVertsp0lb "95%-Konfidenzintervall untere Grenze"
label var pVertsp0ub "95%-Konfidenzintervall obere Grenze"
label var pVerts1p0lb "95%-Konfidenzintervall untere Grenze"
label var pVerts1p0ub "95%-Konfidenzintervall obere Grenze"

//GUE/NGL
quietly nbreg copgAmGUE ipp reading amComStage absDisMedianD1 absDisMedianD2 rappOwnGUE margWinGUE policy_posGUE policy_salGUE if singleVote!=1, robust nolog
prgen ipp, from(0) to (50) generate(pGUE) x(reading=0 rappOwnGUE=0) rest(mean) ci
prgen ipp, from(0) to (50) generate(pGUE1) x(reading=0 rappOwnGUE=1) rest(mean) ci
label var pGUEx "ipp"
label var pGUEp0 "nicht Berichterstatter"
label var pGUEp0 "Berichterstatter"
label var pGUEp0lb "95%-Konfidenzintervall untere Grenze"
label var pGUEp0ub "95%-Konfidenzintervall obere Grenze"
label var pGUE1p0lb "95%-Konfidenzintervall untere Grenze"
label var pGUE1p0ub "95%-Konfidenzintervall obere Grenze"


//UEN
quietly nbreg copgAmUEN ipp reading amComStage absDisMedianD1 absDisMedianD2 rappOwnUEN margWinUEN policy_posUEN policy_salUEN if singleVote!=1, robust nolog
prgen ipp, from(0) to (50) generate(pUEN) x(reading=0 rappOwnUEN=0) rest(mean) ci
prgen ipp, from(0) to (50) generate(pUEN1) x(reading=0 rappOwnUEN=1) rest(mean) ci
label var pUENx "ipp"
label var pUENp0 "nicht Berichterstatter"
label var pUENp0 "Berichterstatter"
label var pUENp0lb "95%-Konfidenzintervall untere Grenze"
label var pUENp0ub "95%-Konfidenzintervall obere Grenze"
label var pUEN1p0lb "95%-Konfidenzintervall untere Grenze"
label var pUEN1p0ub "95%-Konfidenzintervall obere Grenze"


//IND/DEM
quietly nbreg copgAmIND ipp reading amComStage absDisMedianD1 absDisMedianD2 rappOwnIND margWinIND policy_posIND policy_salIND if singleVote!=1, robust nolog
prgen ipp, from(0) to (50) generate(pIND) x(reading=0 rappOwnIND=0) rest(mean) ci
prgen ipp, from(0) to (50) generate(pIND1) x(reading=0 rappOwnIND=1) rest(mean) ci
label var pINDx "ipp"
label var pINDp0 "nicht Berichterstatter"
label var pINDp0 "Berichterstatter"
label var pINDp0lb "95%-Konfidenzintervall untere Grenze"
label var pINDp0ub "95%-Konfidenzintervall obere Grenze"
label var pIND1p0lb "95%-Konfidenzintervall untere Grenze"
label var pIND1p0ub "95%-Konfidenzintervall obere Grenze"



*** Tab. 5.5 Wahrscheinlichkeit eines Antrags auf namentliche Abstimmung für unterschiedliche Antragstypen und in Abhängigkeit von der Position als Berichterstatter 95%-Konfidenzintervalle *** 

use "\...DataDissWIDE.dta"


//generate interaction effects
gen pgAmrappOwnPPE=pgAmPPE*rappOwnPPE
gen finalVoterappOwnPPE=finalVote*rappOwnPPE
gen pgAmrappOwnPSE=pgAmPSE*rappOwnPSE
gen finalVoterappOwnPSE=finalVote*rappOwnPSE
gen pgAmrappOwnALDE=pgAmPSE*rappOwnALDE
gen finalVoterappOwnALDE=finalVote*rappOwnALDE
gen pgAmrappOwnVerts=pgAmVerts*rappOwnVerts
gen finalVoterappOwnVerts=finalVote*rappOwnVerts
gen pgAmrappOwnGUE=pgAmGUE*rappOwnGUE
gen finalVoterappOwnGUE=finalVote*rappOwnGUE
gen pgAmrappOwnUEN=pgAmUEN*rappOwnUEN
gen finalVoterappOwnUEN=finalVote*rappOwnUEN
gen pgAmrappOwnIND=pgAmIND*rappOwnIND
gen finalVoterappOwnIND=finalVote*rappOwnIND


eststo: logit rcvReqPPE pgAmPPE jAmPPE commAm finalVote pgAmPSE pgAmALDE pgAmVerts pgAmGUE pgAmUEN pgAmIND ipp margWinPPE rappOwnPPE pgAmrappOwnPPE finalVoterappOwnPPE policy_posPPE policy_salPPE reading MediaCov amComStage, cluster(Order) nolog
eststo: logit rcvReqPSE pgAmPSE jAmPSE commAm finalVote pgAmPPE pgAmALDE pgAmVerts pgAmGUE pgAmUEN pgAmIND ipp margWinPSE rappOwnPSE pgAmrappOwnPSE finalVoterappOwnPSE policy_posPSE policy_salPSE reading MediaCov amComStage, cluster(Order) nolog
eststo: logit rcvReqALDE pgAmALDE jAmALDE commAm finalVote pgAmPPE pgAmPSE pgAmVert pgAmGUE pgAmUEN pgAmIND ipp margWinALDE rappOwnALDE pgAmrappOwnALDE finalVoterappOwnALDE policy_posALDE policy_salALDE reading MediaCov amComStage, cluster(Order) nolog
eststo: logit rcvReqVerts pgAmVerts jAmVerts commAm finalVote pgAmPPE pgAmPSE pgAmALDE pgAmGUE pgAmUEN pgAmIND ipp margWinVerts rappOwnVerts pgAmrappOwnVerts finalVoterappOwnVerts policy_posVerts policy_salVerts reading MediaCov amComStage, cluster(Order) nolog
eststo: logit rcvReqGUE pgAmGUE jAmGUE commAm finalVote pgAmPPE pgAmPSE pgAmALDE pgAmVerts pgAmUEN pgAmIND ipp margWinGUE rappOwnGUE pgAmrappOwnGUE finalVoterappOwnGUE policy_posGUE policy_salGUE reading MediaCov amComStage, cluster(Order) nolog 
eststo: logit rcvReqUEN pgAmUEN jAmUEN commAm finalVote pgAmPPE pgAmPSE pgAmALDE pgAmVerts pgAmGUE pgAmIND ipp margWinUEN rappOwnUEN pgAmrappOwnUEN finalVoterappOwnUEN policy_posUEN policy_salUEN reading MediaCov amComStage, cluster(Order) nolog 
eststo: logit rcvReqIND pgAmIND jAmIND commAm finalVote pgAmPPE pgAmPSE pgAmALDE pgAmVerts pgAmGUE pgAmUEN ipp margWinIND rappOwnIND pgAmrappOwnIND finalVoterappOwnIND policy_posIND policy_salIND reading MediaCov amComStage, cluster(Order) nolog  
estout using Model3, replace cells(b(star fmt(3)) se(par fmt(3))) stats(r2 N) starlevels(* 0.1 ** 0.05 *** 0.01) style(tex) legend varlabels(_cons Constant commAm "Committee Am." finalVote "Final Vote" ipp "IPP" amComStage "# Amendments" MediaCov "Media Coverage")



***
//First, graph predicted probabilities
///Caution: only include the covariates that Stata kept!! (i.e. omit covariates that are collinear)
///only some examples shown for illustrative purposes
//ALDE
quietly logit rcvReqALDE pgAmALDE pgAmPPE pgAmPSE pgAmVerts jAmALDE commAm ipp reading MediaCov rappOwnALDE margWinALDE policy_posALDE policy_salALDE finalVote, cluster(Order) nolog
prgen ipp, from(0) to (50) generate(pALDE) x(reading=0 rappOwnALDE=0 pgAmALDE=0 jAmALDE=0 finalVote=1 commAm=0) rest(mean) ci
prgen ipp, from(0) to (50) generate(pALDE1) x(reading=0 rappOwnALDE=0 pgAmALDE=1 jAmALDE=0 finalVote=0 commAm=0) rest(mean) ci
prgen ipp, from(0) to (50) generate(pALDE2) x(reading=0 rappOwnALDE=0 pgAmALDE=0 jAmALDE=1 finalVote=0 commAm=0) rest(mean) ci
prgen ipp, from(0) to (50) generate(pALDE3) x(reading=0 rappOwnALDE=0 pgAmALDE=0 jAmALDE=0 finalVote=0 commAm=1) rest(mean) ci
label var pALDEx "ipp"
label var pALDEp1 "Schlussabstimmung"
label var pALDEp1ub "95% upper bound"
label var pALDEp1lb "95% lower bound"
label var pALDE1p1 "Einzelantrag"
label var pALDE1p1ub "95% upper bound"
label var pALDE1p1lb "95% lower bound"
label var pALDE2p1 "Gruppenantrag"
label var pALDE3p1 "Ausschussantrag"


//predicted values for ALDE holding rapporteurship (rappOwnALDE==1)
prgen ipp, from(0) to (50) generate(pALDErapp) x(reading=0 rappOwnALDE=1 pgAmALDE=0 jAmALDE=0 finalVote=1 commAm=0) rest(mean) ci
prgen ipp, from(0) to (50) generate(pALDErapp1) x(reading=0 rappOwnALDE=1 pgAmALDE=1 jAmALDE=0 finalVote=0 commAm=0) rest(mean) ci
prgen ipp, from(0) to (50) generate(pALDErapp2) x(reading=0 rappOwnALDE=1 pgAmALDE=0 jAmALDE=1 finalVote=0 commAm=0) rest(mean) ci
prgen ipp, from(0) to (50) generate(pALDErapp3) x(reading=0 rappOwnALDE=1 pgAmALDE=0 jAmALDE=0 finalVote=0 commAm=1) rest(mean) ci

//ciplot of predicted probability of RCV request for different categories of amendments, varying ALDE rapporteurship 
ciplot pALDEp1 pALDErappp1 pALDE1p1 pALDErapp1p1 pALDE2p1 pALDErapp2p1 pALDE3p1 pALDErapp3p1, msymbol(square) mcolor(black) title("ALDE")  ylabel(0 (.05) .25) ytitle(Wahrscheinlichkeit eines Antrags auf NA) scheme(s2manual) title("ALDE") ylabel(0 (.05) .25) ytitle(Wahrscheinlichkeit eines Antrags auf NA) scheme(s2manual) saving("xxx.gph")
 
 
//to combine graphs (side-by-side), use "graph combine" command with the "ycommon" option to get a single legend (y-axis)

*****************
*** Kapitel 6 ***
*****************


*** Abb. 6.1 Beobachtete und vorhergesagte Anzahl individueller Änderungsanträge *** 

use "\...COUNTS.dta"

//transform dataset in long format 
reshape long margWin CatmargWin strMargWin rappOwn pgAm jAm rcvReq policy_pos policy_sal corcvReq copgAm, i(VoteID) j(party) string
//create a factor variable denominating the PGs under var party
gen PG=1 if party=="PPE"
replace PG=2 if party=="PSE"
replace PG=3 if party=="ALDE"
replace PG=4 if party=="Verts"
replace PG=5 if party=="GUE"
replace PG=6 if party=="UEN"
replace PG=7 if party=="IND"
replace PG=8 if party=="MEP"
label define PG 1"PPE" 2"PSE" 3"ALDE" 4"Verts" 5"GUE" 6"UEN" 7"IND" 8"MEP"
label var PG "Party group"
label values PG PG
order PG, after(votetype)
order party, first
order CatmargWin, after(PG)

quietly nbreg copgAm ipp reading amComStage rappOwn margWin absDisMedianD1 absDisMedianD2 policy_pos policy_sal if singleVote!=1 & PG!=8, robust //exclude single votes (no amendmnets) and MEPs

tw (scatter c copgAm, mlab(Order)) (scatter c copgAm if Order==39, mlab(party) mlabposition(1) mlabgap(2) mlcolor(black)), ylabel(0(50)100) xlabel(0(50)100) legend(off) scheme(s2manual) ytitle("vorhergesagte Anzahl an Änderungsanträgen") xtitle("beobachtete Anzahl an Änderungsanträgen") || line copgAm copgAm

*** QCA *** 

use "...\REACH Amendments Selection.dta" 

//recode compCommAm as a dichotomous variable capturing only whether there were amendments from other committess
drop compCommAm
gen compCommAm=1 if strpos(CommSponsorCompAm, "IMCO") | strpos(CommSponsorCompAm, "ITRE") | strpos(CommSponsorCompAm, "INTA") | strpos(CommSponsorCompAm, "EMPL") ///
| strpos(CommSponsorCompAm, "FEMM") | strpos(CommSponsorCompAm, "JURI") | strpos(CommSponsorCompAm, "ECON") | strpos(CommSponsorCompAm, "BUDG")
replace compCommAm=0 if compCommAm!=1
order compCommAm, after(CommSponsorCompAm)


replace friendlyCommittee=1 if ID==5 | ID==9 | ID==13 | ID==34 | ID==41 | ID==56 | ID==90 | ID==102 | ID==122 | ID==116 | ID==125

//correct observations
replace Reintro=1 if ID==22
replace Reintro=0 if ID==15
replace compCommAm=1 if ID==15
replace CommSponsor=" " if ID==15
replace CommSponsorID=7 if ID==15
replace Intro4DraftReport=0 if ID==15
replace Amendment=" " if ID==15
replace CompAmDraftReport=0 if ID==15
replace AdoptENVI=1 if ID==15
replace compCommAm=0 if ID==15
replace AdoptENVI=1 if CommSponsorID==7

//csQCA
//Analysis with 128 Amendments

fuzzy Reintro CompAmDraftReport AdoptENVI friendlyCommittee compCommAm
fuzzy Reintro CompAmDraftReport AdoptENVI friendlyCommittee compCommAm, necessity
fuzzy Reintro CompAmDraftReport AdoptENVI friendlyCommittee compCommAm, matx(coincid suffnec) standardized
fuzzy Reintro CompAmDraftReport AdoptENVI friendlyCommittee compCommAm, settest(yvo yvv) greater(col1) conval(.80) reduce

//exclude amendments tabled by committee, MEPs and several authors
//for some reason, excluding leads to empty matrices (bestfit variables are not created)

//I am keeping only individual groups and MEPs -> approximate minority rights

// Sample with 106 Amendmens
drop if PlenSponsorID==7
drop if PlenSponsorID==9

fuzzy Reintro CompAmDraftReport AdoptENVI friendlyCommittee compCommAm
fuzzy Reintro CompAmDraftReport AdoptENVI friendlyCommittee compCommAm, necessity //this command is bull, because it tests the necessity of logical AND configurations - which is superfluous (Schneider/Wagemann SMSS: 72f.)
fuzzy Reintro CompAmDraftReport AdoptENVI friendlyCommittee compCommAm,  matx(coincid suffnec) standardized
fuzzy Reintro CompAmDraftReport AdoptENVI friendlyCommittee compCommAm, settest(yvo yvv) greater(col1) conval(.80) reduce //this command yields the same output as that of fsQCA
//in other words, fsQCA does not apply significance thresholds and the like in computing the reduced set
//AbD = AbcD + AbCD 
l VoteResult ID bestfit Reintro PlenSponsorID CommSponsorID Subject Am_no RCV if bestfit==3 
l VoteResult ID bestfit Reintro PlenSponsorID CommSponsorID Subject Am_no RCV if bestfit==5 
//AbC = AbCd + AbCD
l VoteResult ID bestfit Reintro PlenSponsorID CommSponsorID Subject Am_no RCV if bestfit==4