# ***************************** # Forschungsseminar Mastertrack # Politischer Datenjournalismus # # # Garret Binding # 08-711-343 # März 2014 # **************************** # Workspace räumen und working directory setzen rm(list=ls(all=TRUE)) setwd("/home/ddj/gbinding/Datenjournalismus FS 14/Aufgabe 1") # Pakete laden require(foreign) require(MASS) require(sandwich) require(car) require(robust) require(texreg) require(sjPlot) # Laden der Datensätze. data beinhaltet die Werte, data.info die Informationen zu den Variablen. data <- read.csv("~/Datenjournalismus FS 14/Aufgabe 1/data.csv") data$X <- NULL data.info <- read.csv ("~/Datenjournalismus FS 14/Aufgabe 1/data_info.csv", sep="\t") # Schlussendlich interessieren zwei Fragestellungen: # Hat der Ausländeranteil auch noch Kontrolle verschiedener Faktoren einen Effekt auf die Zustimmung zur MEI? (Tagi Datenblog) # Hat die Verschiebung zuungunsten der Freizügigkeit mit dem Raumtyp zu tun? (Moser NZZ) # Ausländeranteil inkl. aller Raumtypen summary(ausl <- lm(Ant.Ja ~ Ant.Ausl + Ant.Ausl^2 + log(Einwo) + SekII + Tert + Eink.Med + Eink.Med^2 + Arbeitslos + Soz.Hilf + Pendl.Einw + Ant.SVP + Bilanz.CH + Bilanz.not.CH + Bilanz.int + Raum.2 + Raum.4 + Ant.Ausl*Raum.2 + Ant.Ausl*Raum.4 , data=data)) # Mit einer rlm() sind die Resultate im Vergleich mit lm() robust. # summary(rlm(Ant.Ja ~ Ant.Ausl + Ant.Ausl^2 + log(Einwo) + SekII + Tert + # Eink.Med + Eink.Med^2 + Arbeitslos + Soz.Hilf + Pendl.Einw + Ant.SVP + # Bilanz.CH + Bilanz.not.CH + Bilanz.int + Raum.2 + Raum.4 + Ant.Ausl*Raum.2 + Ant.Ausl*Raum.4 # , data=data)) lab <- c("Ausländeranteil [%, 2013]", "Einwohnerzahl zivilr. [Pers., log., 2013]", "Schüler Sekundarstufe II [Pers. pro 1000 Einw, 2012]", "Schüler Tertiärstufe [Pers. pro 1000 Einw, 2012]", "Steuerb. Einkommen natürl. Pers. Median [Fr., 2010]", "Arbeitslosenanteil an Bevölkerung 15-64 Jahre [%, 2012]", "Sozialhilfequote [%, 2012]", "Pendleranteil pro Einwohner [2000 / 2013]", "Anteil SVP [%, 2011]", "Wand.Bilanz Schweizer pro 1000 [2012]", "Wand.Bilanz Ausländer pro 1000 [2012]","Wand.Bilanz international pro 1000 [2012]", "Agglomeration [2014]", "Land [2014]", "Agglomeration * Ausländeranteil", "Land * Ausländeranteil") pl.Ja <- sjp.lm(ausl, title="Einfluss auf Ja-Anteil zur MEI (alle Raumtypen)", axisLabels.y=lab, axisLabelSize=1, breakLabelsAt=30, axisTitle.x="Koeffizienten (95% Konfidenzintervall)", theme="minimal") jpeg("pl.Ja.jpeg", width=30, height=30, units='cm', res=300) pl.Ja dev.off() # Welche Faktoren beeinflussten die Zustimmung zur MEI in der Agglomeration? summary(ausl.Raum.2 <- lm(Ant.Ja ~ Ant.Ausl + Ant.Ausl^2 + log(Einwo) + SekII + Tert + Eink.Med + Eink.Med^2 + Arbeitslos + Soz.Hilf + Pendl.Einw + Ant.SVP + Bilanz.CH + Bilanz.not.CH + Bilanz.int , data=subset(data, Raum.2==1))) summary(ausl.Raum.4 <- lm(Ant.Ja ~ Ant.Ausl + Ant.Ausl^2 + log(Einwo) + SekII + Tert + Eink.Med + Eink.Med^2 + Arbeitslos + Soz.Hilf + Pendl.Einw + Ant.SVP + Bilanz.CH + Bilanz.not.CH + Bilanz.int , data=subset(data, Raum.4==1))) lab <- c("Ausländeranteil [%, 2013]", "Einwohnerzahl zivilr. [Pers., log., 2013]", "Schüler Sekundarstufe II [Pers. pro 1000 Einw, 2012]", "Schüler Tertiärstufe [Pers. pro 1000 Einw, 2012]", "Steuerb. Einkommen natürl. Pers. Median [Fr., 2010]", "Arbeitslosenanteil an Bevölkerung 15-64 Jahre [%, 2012]", "Sozialhilfequote [%, 2012]", "Pendleranteil pro Einwohner [2000 / 2013]", "Anteil SVP [%, 2011]", "Wand.Bilanz Schweizer pro 1000 [2012]", "Wand.Bilanz Ausländer pro 1000 [2012]","Wand.Bilanz international pro 1000 [2012]") pl.Ja.agglo <- sjp.lm(ausl.Raum.2, title="Einfluss auf Ja-Anteil zur MEI (Agglomeration)", axisLabels.y=lab, axisLabelSize=1, breakLabelsAt=30, axisTitle.x="Koeffizienten (95% Konfidenzintervall)", theme="minimal") jpeg("pl.Ja.agglo.jpeg", width=30, height=30, units='cm', res=300) pl.Ja.agglo dev.off() pl.Ja.land <- sjp.lm(ausl.Raum.4, title="Einfluss auf Ja-Anteil zur MEI (Land)", axisLabels.y=lab, axisLabelSize=1, breakLabelsAt=30, axisTitle.x="Koeffizienten (95% Konfidenzintervall)", theme="minimal") jpeg("pl.Ja.Land.jpeg", width=30, height=30, units='cm', res=300) pl.Ja.land dev.off() l <- list(ausl, ausl.Raum.2, ausl.Raum.4) c <- c("Konstante", "Ausländeranteil [%, 2013]", "Einwohnerzahl zivilr. [Pers., log., 2013]", "Schüler Sekundarstufe II [Pers. pro 1000 Einw, 2012]", "Schüler Tertiärstufe [Pers. pro 1000 Einw, 2012]", "Steuerb. Einkommen natürl. Pers. Median [Fr., 2010]", "Arbeitslosenanteil an Bevölkerung 15-64 Jahre [%, 2012]", "Sozialhilfequote [%, 2012]", "Pendleranteil pro Einwohner [%, 2000 / 2013]", "Anteil SVP [%, 2011]", "Wand.Bilanz Schweizer pro 1000 [2012]", "Wand.Bilanz Ausländer pro 1000 [2012]","Wand.Bilanz international pro 1000 [2012]", "Agglomeration", "Land", "Agglomeration * Ausländeranteil", "Land * Ausländeranteil") htmlreg(l, file = "MEI.html", custom.model.names = c("Alle Raumtypen", "Agglomeration", "Land"), custom.coef.names=c, caption.above=TRUE, caption = "Tabelle 1: Einfluss auf Ja-Anteil zur MEI", bold = 0.05) # Wenn wir die Raumtypen "Andere Agglomerationsgemeinde" und "Ländliche Gemeinde" separat berechnen, verschwindet dieser # Effekt. Für diese Raumtypen ist der Effekt des Ausländeranteils in den Gemeinden des Kt. Zürich nicht signifikant. # Kommen wir zur Veränderung der Zustimmung zur Personenfreizügigkeit. Dazu brauchen wir erst mal das delta(2000-2014). data$Bil.Ant.Nein <- 100-(100 * data$Bil.Ant.Ja) data$delta <- (data$Bil.Ant.Nein-data$Ant.Ja) summary(delta <- lm(delta ~ Ant.Ausl + Ant.Ausl^2 + log(Einwo) + SekII + Tert + Eink.Med + Eink.Med^2 + Arbeitslos + Soz.Hilf + Pendl.Einw + Ant.SVP + Bilanz.CH + Bilanz.not.CH + Bilanz.int + Raum.2 + Raum.4 + Ant.Ausl*Raum.2 + Ant.Ausl*Raum.4 , data=data)) lab <- c("Ausländeranteil [%, 2013]", "Einwohnerzahl zivilr. [Pers., log., 2013]", "Schüler Sekundarstufe II [Pers. pro 1000 Einw, 2012]", "Schüler Tertiärstufe [Pers. pro 1000 Einw, 2012]", "Steuerb. Einkommen natürl. Pers. Median [Fr., 2010]", "Arbeitslosenanteil an Bevölkerung 15-64 Jahre [%, 2012]", "Sozialhilfequote [%, 2012]", "Pendleranteil pro Einwohner [2000 / 2013]", "Anteil SVP [%, 2011]", "Wand.Bilanz Schweizer pro 1000 [2012]", "Wand.Bilanz Ausländer pro 1000 [2012]","Wand.Bilanz international pro 1000 [2012]", "Agglomeration [2014]", "Land [2014]", "Agglomeration * Ausländeranteil", "Land * Ausländeranteil") pl.delta <- sjp.lm(delta, title="Einfluss auf Zustimmungsänderung (alle Raumtypen)", axisLabels.y=lab, axisLabelSize=1, breakLabelsAt=30, axisTitle.x="Koeffizienten (95% Konfidenzintervall)", theme="minimal") jpeg("pl.delta.jpeg", width=30, height=30, units='cm', res=300) pl.delta dev.off() # Der Raumtyp der Gemeinden im Kt. Zürich hat keinen Einfluss auf die Veränderung der Zustimmung zur Personenfreizügigkeit. # Was geschieht, wenn wir nur den Vergleich Stadt - Agglo und Stadt - Land machen? # Stadt vs. Agglo summary(delta.Raum.2 <- lm(delta ~ Ant.Ausl + Ant.Ausl^2 + log(Einwo) + SekII + Tert + Eink.Med + Eink.Med^2 + Arbeitslos + Soz.Hilf + Pendl.Einw + Ant.SVP + Bilanz.CH + Bilanz.not.CH + Bilanz.int + Raum.2 + Ant.Ausl*Raum.2 , data=subset(data, Raum.4==0))) # Stadt vs. Land summary(delta.Raum.4 <- lm(delta ~ Ant.Ausl + Ant.Ausl^2 + log(Einwo) + SekII + Tert + Eink.Med + Eink.Med^2 + Arbeitslos + Soz.Hilf + Pendl.Einw + Ant.SVP + Bilanz.CH + Bilanz.not.CH + Bilanz.int + Raum.4 + Ant.Ausl*Raum.4 , data=subset(data, Raum.2==0))) lab <- c("Ausländeranteil [%, 2013]", "Einwohnerzahl zivilr. [Pers., log., 2013]", "Schüler Sekundarstufe II [Pers. pro 1000 Einw, 2012]", "Schüler Tertiärstufe [Pers. pro 1000 Einw, 2012]", "Steuerb. Einkommen natürl. Pers. Median [Fr., 2010]", "Arbeitslosenanteil an Bevölkerung 15-64 Jahre [%, 2012]", "Sozialhilfequote [%, 2012]", "Pendleranteil pro Einwohner [2000 / 2013]", "Anteil SVP [%, 2011]", "Wand.Bilanz Schweizer pro 1000 [2012]", "Wand.Bilanz Ausländer pro 1000 [2012]","Wand.Bilanz international pro 1000 [2012]", "Agglomeration [2014]", "Agglomeration * Ausländeranteil") pl.delta.agglo <- sjp.lm(delta.Raum.2, title="Einfluss auf Zustimmungsänderung (Stadt vs. Agglomeration)", axisLabels.y=lab, axisLabelSize=1, breakLabelsAt=30, axisTitle.x="Koeffizienten (95% Konfidenzintervall)", theme="minimal") jpeg("pl.delta.agglo.jpeg", width=30, height=30, units='cm', res=300) pl.delta.agglo dev.off() lab <- c("Ausländeranteil [%, 2013]", "Einwohnerzahl zivilr. [Pers., log., 2013]", "Schüler Sekundarstufe II [Pers. pro 1000 Einw, 2012]", "Schüler Tertiärstufe [Pers. pro 1000 Einw, 2012]", "Steuerb. Einkommen natürl. Pers. Median [Fr., 2010]", "Arbeitslosenanteil an Bevölkerung 15-64 Jahre [%, 2012]", "Sozialhilfequote [%, 2012]", "Pendleranteil pro Einwohner [2000 / 2013]", "Anteil SVP [%, 2011]", "Wand.Bilanz Schweizer pro 1000 [2012]", "Wand.Bilanz Ausländer pro 1000 [2012]","Wand.Bilanz international pro 1000 [2012]", "Land [2014]", "Land * Ausländeranteil") pl.delta.land <- sjp.lm(delta.Raum.4, title="Einfluss auf Zustimmungsänderung (Stadt vs. Land)", axisLabels.y=lab, axisLabelSize=1, breakLabelsAt=30, axisTitle.x="Koeffizienten (95% Konfidenzintervall)", theme="minimal") jpeg("pl.delta.land.jpeg", width=30, height=30, units='cm', res=300) pl.delta.land dev.off() l <- list(delta, delta.Raum.2, delta.Raum.4) c <- c("Konstante", "Ausländeranteil [%, 2013]", "Einwohnerzahl zivilr. [Pers., log., 2013]", "Schüler Sekundarstufe II [Pers. pro 1000 Einw, 2012]", "Schüler Tertiärstufe [Pers. pro 1000 Einw, 2012]", "Steuerb. Einkommen natürl. Pers. Median [Fr., 2010]", "Arbeitslosenanteil an Bevölkerung 15-64 Jahre [%, 2012]", "Sozialhilfequote [%, 2012]", "Pendleranteil pro Einwohner [%, 2000 / 2013]", "Anteil SVP [%, 2011]", "Wand.Bilanz Schweizer pro 1000 [2012]", "Wand.Bilanz Ausländer pro 1000 [2012]","Wand.Bilanz international pro 1000 [2012]", "Agglomeration", "Land", "Agglomeration * Ausländeranteil", "Land * Ausländeranteil") htmlreg(l, file = "delta.html", custom.model.names = c("Alle Raumtypen", "Stadt vs. Agglomeration", "Stadt vs. Land"), custom.coef.names=c, caption.above=TRUE, caption = "Tabelle 2: Einfluss auf Zustimmungsveränderung", bold = 0.05)