clear all
sysuse stimmenanteile_data

gen float date = date(datum, "DMY")
format date %td
gen datesif = date

rename land bundesland

replace ebene = "1_Europa" if ebene == "Europa"
replace ebene = "2_Bund" if ebene == "Bund"
replace ebene = "3_Land" if ebene == "Land"
replace ebene = "4_Kreis" if ebene == "Kreis"

gen kreis = 0
replace kreis = 1 if ebene == "4_Kreis"

gen land = 0
replace land = 1 if ebene == "3_Land"

gen bund = 0
replace bund = 1 if ebene == "2_Bund"

gen europa = 0
replace europa = 1 if ebene == "1_Europa"


gen neue_bl = 0
replace neue_bl = 1 if bundesland == "Brandenburg" | bundesland == "Mecklenburg-Vorpommern" | ///
bundesland == "Sachsen" | bundesland == "Sachsen-Anhalt" | bundesland == "Thüringen"

*WS-Variablen
gen huerde = 0
replace huerde = 1 if sperrklausel > 0

gen pv = 0
replace pv = 1 if wahlsystem != "rv"


*Gruppen nach Niedermayer 2014
gen anteil_große = cdu_share + spd_share

gen anteil_mittel = grne_share + fdp_share + linke_share
**replace anteil_mittel = grne_share + linke_share if datesif > 19623

gen anteil_mittel_ol =  anteil_mittel - linke_share

gen anteil_klein = 1 - anteil_große - anteil_mittel

twoway (scatter anteil_klein jahr if kreis == 1, legend(label(1 "Kreis"))) ///
(scatter anteil_klein jahr if land == 1, legend(label(2 "Land"))) ///
(scatter anteil_klein jahr if bund == 1, legend(label(3 "Bund"))) ///
(scatter anteil_klein jahr if europa == 1, legend(label(4 "Europa")))


*Deskriptiv
tabulate ebene, summarize(anteil_große)
tabulate ebene, summarize(anteil_mittel)
tabulate ebene, summarize(anteil_klein)

*Zeit
pwcorr anteil_große anteil_mittel anteil_klein jahr, sig

pwcorr anteil_große anteil_mittel anteil_klein jahr if kreis == 1, sig
pwcorr anteil_große anteil_mittel anteil_klein jahr if land == 1, sig
pwcorr anteil_große anteil_mittel anteil_klein jahr if bund == 1, sig
pwcorr anteil_große anteil_mittel anteil_klein jahr if europa == 1, sig


*Abbildungen
*Ebenen und Zeittrends
twoway (scatter anteil_klein date if europa == 1, legend(label(1 "Europa"))) (scatter anteil_klein date if bund == 1, legend(label(2 "Bund"))) (scatter anteil_klein date if land == 1, legend(label(3 "Land"))) (scatter anteil_klein date if kreis == 1, legend(label(4 "Kreis"))) (fpfit anteil_klein date if europa == 1) (fpfit anteil_klein date if bund == 1) (fpfit anteil_klein date if land == 1) (fpfit anteil_klein date if kreis == 1)
twoway (fpfit anteil_klein date if europa == 1, legend(label(1 "Europa"))) (fpfit anteil_klein date if bund == 1, legend(label(2 "Bund"))) (fpfit anteil_klein date if land == 1, legend(label(3 "Land"))) (fpfit anteil_klein date if kreis == 1, legend(label(4 "Kreis")))

twoway (scatter anteil_mittel date if europa == 1, legend(label(1 "Europa"))) (scatter anteil_mittel date if bund == 1, legend(label(2 "Bund"))) (scatter anteil_mittel date if land == 1, legend(label(3 "Land"))) (scatter anteil_mittel date if kreis == 1, legend(label(4 "Kreis")))
twoway (fpfit anteil_mittel date if europa == 1, legend(label(1 "Europa"))) (fpfit anteil_mittel date if bund == 1, legend(label(2 "Bund"))) (fpfit anteil_mittel date if land == 1, legend(label(3 "Land"))) (fpfit anteil_mittel date if kreis == 1, legend(label(4 "Kreis")))

twoway (scatter anteil_mittel_ol date if europa == 1, legend(label(1 "Europa"))) (scatter anteil_mittel_ol date if bund == 1, legend(label(2 "Bund"))) (scatter anteil_mittel_ol date if land == 1, legend(label(3 "Land"))) (scatter anteil_mittel_ol date if kreis == 1, legend(label(4 "Kreis")))
twoway (fpfit anteil_mittel_ol date if europa == 1, legend(label(1 "Europa"))) (fpfit anteil_mittel_ol date if bund == 1, legend(label(2 "Bund"))) (fpfit anteil_mittel_ol date if land == 1, legend(label(3 "Land"))) (fpfit anteil_mittel_ol date if kreis == 1, legend(label(4 "Kreis")))

twoway (scatter anteil_große date if europa == 1, legend(label(1 "Europa"))) (scatter anteil_große date if bund == 1, legend(label(2 "Bund"))) (scatter anteil_große date if land == 1, legend(label(3 "Land"))) (scatter anteil_große date if kreis == 1, legend(label(4 "Kreis")))
twoway (fpfit anteil_große date if europa == 1, legend(label(1 "Europa"))) (fpfit anteil_große date if bund == 1, legend(label(2 "Bund"))) (fpfit anteil_große date if land == 1, legend(label(3 "Land"))) (fpfit anteil_große date if kreis == 1, legend(label(4 "Kreis")))


xi: reg anteil_klein jahr kreis land europa i.bundesland
xi: reg anteil_mittel jahr kreis land europa i.bundesland
xi: reg anteil_mittel_ol jahr kreis land europa i.bundesland
xi: reg anteil_große jahr kreis land europa i.bundesland


*Hürde
twoway (scatter anteil_klein date if huerde == 1, legend(label(1 "Hürde"))) (scatter anteil_klein date if huerde == 0, legend(label(2 "ohne")))
twoway (fpfit anteil_klein date if huerde == 1, legend(label(1 "Hürde"))) (fpfit anteil_klein date if huerde == 0, legend(label(2 "ohne")))

*Nur Kreis
twoway (scatter anteil_klein date if huerde == 1 & kreis == 1, legend(label(1 "Hürde"))) (scatter anteil_klein date if huerde == 0 & kreis == 1, legend(label(2 "ohne"))) (fpfit anteil_klein date if huerde == 1 & kreis == 1) (fpfit anteil_klein date if huerde == 0 & kreis == 1)
twoway (scatter anteil_klein date if huerde == 1 & kreis == 1 & jahr <= 2005, legend(label(1 "Hürde"))) (scatter anteil_klein date if huerde == 0 & kreis == 1 & jahr <= 2005, legend(label(2 "ohne"))) (fpfit anteil_klein date if huerde == 1 & kreis == 1 & jahr <= 2005) (fpfit anteil_klein date if huerde == 0 & kreis == 1 & jahr <= 2005)
twoway (fpfitci anteil_klein date if huerde == 1 & kreis == 1) (fpfitci anteil_klein date if huerde == 0 & kreis == 1)

sum anteil_große if huerde == 1 & kreis == 1
sum anteil_große if huerde == 0 & kreis == 1
ttest anteil_große if kreis == 1, by(huerde)

sum anteil_mittel if huerde == 1 & kreis == 1
sum anteil_mittel if huerde == 0 & kreis == 1
ttest anteil_mittel if kreis == 1, by(huerde)

sum anteil_klein if huerde == 1 & kreis == 1
sum anteil_klein if huerde == 0 & kreis == 1
ttest anteil_klein if kreis == 1, by(huerde)

sum anteil_klein if huerde == 1 & kreis == 1 & jahr < 2005
sum anteil_klein if huerde == 0 & kreis == 1 & jahr < 2005
sum anteil_klein if huerde == 1 & kreis == 1 & jahr < 2000
sum anteil_klein if huerde == 0 & kreis == 1 & jahr < 2000
sum anteil_klein if huerde == 1 & kreis == 1 & jahr < 1995
sum anteil_klein if huerde == 0 & kreis == 1 & jahr < 1995

sum anteil_klein if huerde == 1 & kreis == 1 & neue_bl == 0
sum anteil_klein if huerde == 0 & kreis == 1 & neue_bl == 0
sum anteil_klein if huerde == 1 & kreis == 1 & neue_bl == 1
sum anteil_klein if huerde == 0 & kreis == 1 & neue_bl == 1

twoway (scatter anteil_mittel date if huerde == 1 & kreis == 1, legend(label(1 "Hürde"))) (scatter anteil_mittel date if huerde == 0 & kreis == 1, legend(label(2 "ohne")))
twoway (fpfit anteil_mittel date if huerde == 1 & kreis == 1, legend(label(1 "Hürde"))) (fpfit anteil_mittel date if huerde == 0 & kreis == 1, legend(label(2 "ohne")))

sum anteil_mittel if huerde == 1 & kreis == 1
sum anteil_mittel if huerde == 0 & kreis == 1
sum anteil_mittel if huerde == 1 & kreis == 1 & jahr < 2005
sum anteil_mittel if huerde == 0 & kreis == 1 & jahr < 2005
sum anteil_mittel if huerde == 1 & kreis == 1 & jahr < 2000
sum anteil_mittel if huerde == 0 & kreis == 1 & jahr < 2000
sum anteil_mittel if huerde == 1 & kreis == 1 & jahr < 1995
sum anteil_mittel if huerde == 0 & kreis == 1 & jahr < 1995

sum anteil_mittel if huerde == 1 & kreis == 1 & neue_bl == 0
sum anteil_mittel if huerde == 0 & kreis == 1 & neue_bl == 0
sum anteil_mittel if huerde == 1 & kreis == 1 & neue_bl == 1
sum anteil_mittel if huerde == 0 & kreis == 1 & neue_bl == 1

sum anteil_große if huerde == 1 & kreis == 1 & neue_bl == 0
sum anteil_große if huerde == 0 & kreis == 1 & neue_bl == 0
sum anteil_große if huerde == 1 & kreis == 1 & neue_bl == 1
sum anteil_große if huerde == 0 & kreis == 1 & neue_bl == 1


reg anteil_große huerde jahr if kreis == 1 
reg anteil_mittel huerde jahr if kreis == 1 
reg anteil_klein huerde jahr if kreis == 1 

reg anteil_große huerde jahr neue_bl if kreis == 1 
reg anteil_mittel huerde jahr neue_bl if kreis == 1 
reg anteil_klein huerde jahr neue_bl if kreis == 1


*Vorher-Nachher Länder
gen reform_land = 0
replace reform_land = 1 if bundesland == "Hessen" | bundesland == "Mecklenburg-Vorpommern" | ///
bundesland == "Nordrhein-Westfalen" | bundesland == "Rheinland-Pfalz" | bundesland == "Schleswig-Holstein" | ///
bundesland == "Saarland" | bundesland == "Thüringen"

ttest anteil_klein if reform_land == 1 & kreis == 1, by(huerde)
ttest anteil_mittel if reform_land == 1 & kreis == 1 & neue_bl == 0, by(huerde)
ttest anteil_große if reform_land == 1 & kreis == 1 & neue_bl == 0, by(huerde)

xi: reg anteil_große huerde i.bundesland if kreis == 1 & reform_land == 1
xi: reg anteil_mittel huerde i.bundesland if kreis == 1  & reform_land == 1
xi: reg anteil_klein huerde i.bundesland if kreis == 1 & reform_land == 1

reg anteil_große huerde jahr neue_bl if kreis == 1 & reform_land == 1
reg anteil_mittel huerde jahr neue_bl if kreis == 1  & reform_land == 1
reg anteil_klein huerde jahr neue_bl if kreis == 1 & reform_land == 1


***
gen huerdenzeit = 0
replace huerdenzeit = 1 if bundesland == "Schleswig-Holstein" & datesif <= 16721.5
replace huerdenzeit = 1 if bundesland == "Hessen" & datesif <= 14313.5
replace huerdenzeit = 1 if bundesland == "Mecklenburg-Vorpommern" & datesif <= 15321.5
replace huerdenzeit = 1 if bundesland == "Nordrhein-Westfalen" & datesif <= 13603
replace huerdenzeit = 1 if bundesland == "Rheinland-Pfalz" & datesif <= 17145
replace huerdenzeit = 1 if bundesland == "Saarland" & datesif <= 17145
replace huerdenzeit = 1 if bundesland == "Thüringen" & datesif <= 17152

sum anteil_klein if kreis == 1 & reform_land == 1 & huerdenzeit == 0
sum anteil_klein if kreis == 1 & reform_land == 1 & huerdenzeit == 1
sum anteil_klein if land == 1 & reform_land == 1 & huerdenzeit == 0
sum anteil_klein if land == 1 & reform_land == 1 & huerdenzeit == 1
sum anteil_klein if bund == 1 & reform_land == 1 & huerdenzeit == 0
sum anteil_klein if bund == 1 & reform_land == 1 & huerdenzeit == 1

sum anteil_große if kreis == 1 & reform_land == 1 & huerdenzeit == 0
sum anteil_große if kreis == 1 & reform_land == 1 & huerdenzeit == 1
sum anteil_große if land == 1 & reform_land == 1 & huerdenzeit == 0
sum anteil_große if land == 1 & reform_land == 1 & huerdenzeit == 1
sum anteil_große if bund == 1 & reform_land == 1 & huerdenzeit == 0
sum anteil_große if bund == 1 & reform_land == 1 & huerdenzeit == 1


*pV vs. rV
twoway (scatter anteil_klein date if pv == 1, legend(label(1 "pV"))) (scatter anteil_klein date if pv == 0, legend(label(2 "rV")))
twoway (fpfit anteil_klein date if pv == 1, legend(label(1 "pV"))) (fpfit anteil_klein date if pv == 0, legend(label(2 "rV")))

twoway (scatter anteil_klein date if pv == 1 & huerde == 1, legend(label(1 "pV + Hürde"))) (scatter anteil_klein date if pv == 0 & huerde == 1, legend(label(2 "rV + Hürde"))) (scatter anteil_klein date if pv == 0 & huerde == 0, legend(label(3 "rV ohne"))) (scatter anteil_klein date if pv == 1 & huerde == 0, legend(label(4 "pV ohne")))
twoway (fpfit anteil_klein date if pv == 1 & huerde == 1, legend(label(1 "pV + Hürde"))) (fpfit anteil_klein date if pv == 0 & huerde == 1, legend(label(2 "rV + Hürde"))) (fpfit anteil_klein date if pv == 0 & huerde == 0, legend(label(3 "rV ohne"))) (fpfit anteil_klein date if pv == 1 & huerde == 0, legend(label(4 "pV ohne")))

*Nur Land
twoway (scatter anteil_klein date if pv == 1 & land == 1, legend(label(1 "pV"))) (scatter anteil_klein date if pv == 0 & land == 1, legend(label(2 "rV")))
twoway (fpfit anteil_klein date if pv == 1 & land == 1, legend(label(1 "pV"))) (fpfit anteil_klein date if pv == 0 & land == 1, legend(label(2 "rV")))

*Nur Kreis
twoway (scatter anteil_klein date if pv == 1 & huerde == 1 & kreis == 1, legend(label(1 "pV + Hürde"))) (scatter anteil_klein date if pv == 0 & huerde == 1 & kreis == 1, legend(label(2 "rV + Hürde"))) (scatter anteil_klein date if pv == 0 & huerde == 0 & kreis == 1, legend(label(3 "rV ohne"))) (scatter anteil_klein date if pv == 1 & huerde == 0 & kreis == 1, legend(label(4 "pV ohne")))
twoway (fpfit anteil_klein date if pv == 1 & huerde == 1 & kreis == 1, legend(label(1 "pV + Hürde"))) (fpfit anteil_klein date if pv == 0 & huerde == 1 & kreis == 1, legend(label(2 "rV + Hürde"))) (fpfit anteil_klein date if pv == 0 & huerde == 0 & kreis == 1, legend(label(3 "rV ohne"))) (fpfit anteil_klein date if pv == 1 & huerde == 0 & kreis == 1, legend(label(4 "pV ohne")))


sum anteil_klein if pv == 0 & kreis == 1
sum anteil_klein if pv == 1 & kreis == 1

sum anteil_klein if pv == 0 & kreis == 1 & huerde == 0
sum anteil_klein if pv == 1 & kreis == 1 & huerde == 0

sum anteil_klein if pv == 0 & kreis == 1 & huerde == 1
sum anteil_klein if pv == 1 & kreis == 1 & huerde == 1

sum anteil_klein if pv == 0 & land == 1
sum anteil_klein if pv == 1 & land == 1

sum anteil_mittel if pv == 0 & kreis == 1 & huerde == 0
sum anteil_mittel if pv == 1 & kreis == 1 & huerde == 0

sum anteil_mittel if pv == 0 & kreis == 1 & huerde == 1
sum anteil_mittel if pv == 1 & kreis == 1 & huerde == 1

sum anteil_mittel if pv == 0 & land == 1
sum anteil_mittel if pv == 1 & land == 1

sum anteil_große if pv == 0 & kreis == 1 & huerde == 0
sum anteil_große if pv == 1 & kreis == 1 & huerde == 0

sum anteil_große if pv == 0 & kreis == 1 & huerde == 1
sum anteil_große if pv == 1 & kreis == 1 & huerde == 1

sum anteil_große if pv == 0 & land == 1
sum anteil_große if pv == 1 & land == 1


*Regressionsmodelle
*Große Parteien
xi: reg anteil_große kreis land europa

xi: reg anteil_große huerde pv

xi: reg anteil_große huerde pv kreis land europa

xi: reg anteil_große huerde pv kreis land europa jahr i.bundesland


*Mittlere Parteien
xi: reg anteil_mittel kreis land europa

xi: reg anteil_mittel huerde pv

xi: reg anteil_mittel huerde pv kreis land europa

xi: reg anteil_mittel huerde pv kreis land europa jahr i.bundesland


*Mittlere Parteien ohne Linke
xi: reg anteil_mittel_ol kreis land europa

xi: reg anteil_mittel_ol huerde pv

xi: reg anteil_mittel_ol huerde pv kreis land europa

xi: reg anteil_mittel_ol huerde pv kreis land europa jahr i.bundesland


*Kleine und Kleinstparteien
xi: reg anteil_klein kreis land europa

xi: reg anteil_klein huerde pv

xi: reg anteil_klein huerde pv kreis land europa

xi: reg anteil_klein huerde pv kreis land europa jahr i.bundesland