Tabellen und Grafiken in R (2024)

2.1 Einfache Häufigkeitstabelle

Wir berechnen eine einfache Häufigkeitstabelle für die Variable geschlecht. Dies geschieht mit den Funktionen group_by und summarise. Die Funktion n übernimmt das eigentliche Auszählen der Häufigkeit einer Variablenausprägung.

geschlechterverteilung <- sample_klein %>%  group_by(geschlecht) %>%  summarise(anzahl = n())geschlechterverteilung

# A tibble: 2 × 2 geschlecht anzahl <fct> <int>1 MANN 152 FRAU 5

Im nächsten Schritt fügen wir den absoluten Zahlen relative Anteile hinzu, um das Geschlechterverhältnis besser zu verstehen.

geschlechterverteilung <- sample_klein %>%  group_by(geschlecht) %>%  summarise(anzahl = n()) %>%  mutate(anteil = anzahl/sum(anzahl))geschlechterverteilung

# A tibble: 2 × 3 geschlecht anzahl anteil <fct> <int> <dbl>1 MANN 15 0.752 FRAU 5 0.25

Nun fügen wir den Anteilen von 1 noch Prozentanteile hinzu, indem with mit 100 multiplizieren und runden (was in erster Linie kosmetischer Natur ist).

geschlechterverteilung <- sample_klein %>%  group_by(geschlecht) %>%  summarise(anzahl = n()) %>%  mutate(anteil = anzahl/sum(anzahl), prozent = round(anteil * 100))geschlechterverteilung

# A tibble: 2 × 4 geschlecht anzahl anteil prozent <fct> <int> <dbl> <dbl>1 MANN 15 0.75 752 FRAU 5 0.25 25

2.2 Tabellen speichern

Jetzt schreiben wir die Daten in eine CSV-Datei, die wir bspw. später mit Excel oder einer anderen Tabellenkalkulation öffnen können.

write_excel_csv2(geschlechterverteilung, file = "geschlechterverteilung.csv")

2.3 Kreuztabellen

Was, wenn wir zwei Variablen mit Blick auf ihre Werte in Beziehung setzen wollen? Das bezeichnet man als Kreuz- oder Kontigenztabelle. Man erzeugt eine Kreuztabelle mit dplyr, indem man nach mehreren Variablen gruppiert.

bildung_und_geschlecht <- sample_klein %>%  group_by(bildung, geschlecht) %>%  summarise(anzahl = n()) %>%  mutate(anteil = anzahl/sum(anzahl), prozent = round(anteil * 100))bildung_und_geschlecht

# A tibble: 9 × 5# Groups: bildung [6] bildung geschlecht anzahl anteil prozent <fct> <fct> <int> <dbl> <dbl>1 VOLKS-,HAUPTSCHULE MANN 1 0.5 502 VOLKS-,HAUPTSCHULE FRAU 1 0.5 503 MITTLERE REIFE MANN 4 0.571 574 MITTLERE REIFE FRAU 3 0.429 435 FACHHOCHSCHULREIFE MANN 3 1 1006 HOCHSCHULREIFE MANN 5 0.833 837 HOCHSCHULREIFE FRAU 1 0.167 178 NOCH SCHUELER MANN 1 1 1009 <NA> MANN 1 1 100

Möglicherweise möchte ich den Anteil anders berechnen und nicht die relative Geschlechterverteilung innerhalb eines Bildungsabschlusses in den Blick nehmen, sondern etwa die Verteilung der Bildungsabschlüsse jeweils für männliche und weibliche Studienteilnehmer anschauen. Dies lässt sich mit einer Neugruppierung der Daten durch group_by() erreichen.

geschlecht_und_bildung <- sample_klein %>%  group_by(geschlecht, bildung) %>%  summarise(anzahl = n()) %>%  mutate(anteil = anzahl/sum(anzahl), prozent = round(anteil * 100))geschlecht_und_bildung

# A tibble: 9 × 5# Groups: geschlecht [2] geschlecht bildung anzahl anteil prozent <fct> <fct> <int> <dbl> <dbl>1 MANN VOLKS-,HAUPTSCHULE 1 0.0667 72 MANN MITTLERE REIFE 4 0.267 273 MANN FACHHOCHSCHULREIFE 3 0.2 204 MANN HOCHSCHULREIFE 5 0.333 335 MANN NOCH SCHUELER 1 0.0667 76 MANN <NA> 1 0.0667 77 FRAU VOLKS-,HAUPTSCHULE 1 0.2 208 FRAU MITTLERE REIFE 3 0.6 609 FRAU HOCHSCHULREIFE 1 0.2 20

Lag also im ersten Anlauf das Geschlecht der Gruppierung zugrunde, ist es in der zweiten Version der Bildungsabschluss.

3.1 Balkendiagramme

Das vermutlich einfachste Diagramm, das man mit dem ggplot2-Paket erstellen kann, ist ein Balken- oder Säulendiagramm (barplot). Es zeigt die Häufigkeitsverteilung einer diskreten Variablen, indem es Säulen auf der x-Achse darstellt. Balkendiagramme sind besonders geeignet, um wenige Ausprägungen (bis ca. 15) zu veranschaulichen. Wenn es mehr Kategorien gibt, wird die Anschaulichkeit beeinträchtigt, und es empfiehlt sich, auf Liniendiagramme auszuweichen.

Das nachstehende Beispiel zeigt die Häufigkeitsverteilung der Variable geschlecht im Datensatz (bzw. dem kleinen Sample).

ggplot(sample_klein, aes(geschlecht)) +  geom_bar()

In einem nächsten Schritt fügen wir eine Überschrift hinzu und formatieren die Balken und Achsen so, dass das Plot lesbarer ist.

ggplot(sample_klein, aes(bildung)) +  geom_bar() +  ggtitle("Verteilung der Bildungsabschlüsse der Befragten") +  theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1)) +  xlab("") + ylab("")

Farben lassen sich in ggplot vielseitig einsetzen, um Kategorienunterschiede anzuzeigen. Dies geschieht mit den Argumenten ‘fill’ bzw. ‘color’.

ggplot(sample_klein, aes(bildung, fill = bildung)) +  geom_bar() + ggtitle("Verteilung der Bildungsabschlüsse der Befragten") +  theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1)) +  xlab("") + ylab("")

Es existieren zahlreiche Farbpaletten für ggplot, um unterschiedliche Arten von Beziehungen darzustellen. Das nachstehende Palette unterscheidet verschiedenen Kategorien. Für gradierte Variablen (‘viel’ - ‘wenig’) sind andere Paletten z.T. besser geeignet. Eine gute Auswahl an Palette enthält u.a. das Paket RColorBrewer.

ggplot(sample_klein, aes(bildung, fill = bildung)) +  geom_bar() + coord_flip() +  scale_fill_brewer(palette = "Set1") +  ggtitle("Verteilung der Bildungsabschlüsse der Befragten") +  xlab("") + ylab("")

Schließlich lassen sich auch die Beschriftung und weitere Aspekte eines Plots anpassen (Achsenorientierung, Legende etc).

ggplot(sample_klein, aes(as_factor(fernsehkonsum))) +  geom_bar() + theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1)) +  ggtitle("Fernsehkonsum pro Woche in Tagen") +  xlab("Fernsehkonsum pro Woche in Tagen") + ylab("Anzahl der Respondenten")

3.2 Histogramme

Für ein besseres Verständnis einer Verteilung sind oftmals die Häufigkeitsausprägungen einer kontinuierlichen Variable in gleich großen Gruppen interessant (sog. “bins”). Dabei hilft der Visualisierungstyp Histogramm. Ein Histogramm ist eine grafische Darstellung der Häufigkeitsverteilung kardinal skalierter Merkmale. Dabei werden Daten in Klassen unterteilt, die unterschiedliche Breiten haben können. Rechtecke mit der Breite der Klassen werden nebeneinander gezeichnet, wobei ihre Flächen die (relativen oder absoluten) Klassenhäufigkeiten repräsentieren. Die Höhe jedes Rechtecks zeigt die (relative oder absolute) Häufigkeitsdichte an, also die Häufigkeit geteilt durch die Breite der jeweiligen Klasse.

ggplot(sample_mittel, aes(alter)) +  geom_histogram()

Auch hier lassen sich relevante Aspekte anpassen, etwa die Anzahl und Breite der Flächen.

ggplot(sample_mittel, aes(alter)) +  geom_histogram(bins = 40) +  ggtitle("Altersverteilung der Respondenten") +  xlab("Alter") + ylab("Anzahl der RespondentInnen")

Hier können wir erstmalig die Farbe eines Elements gezielt einsetzen, um eine zusätzliche (also nach Vetrauen auf der x-Achse und der Anzahl der Respondenten auf der y-Achse eine dritte Variable) darzustellen, nämlich das Geschlecht der RespondentInnen.

ggplot(sample_gross, aes(vertrauen_polizei, fill = geschlecht)) +  geom_histogram(binwidth = 1, position = "dodge") +  ggtitle("Vertrauen in die Polizei nach Geschlecht") +  xlab("Vetrauen (1-7)") + ylab("Anzahl der RespondentInnen") + labs(fill = "Geschlecht") 

Neben Blaken und Flächen beherrscht ggplot auch zahlreiche weitere Darstellungsformen (sog. geoms). Eine interessante Alternative zum klassischen Histogramm ist etwa das Dichte-Plot (density plot). Nachstehend verwenden wir zwei Flächen und einen Trasparenz-Effekt für die Darstellung.

3.3 Dichte-Plots

polizei <- sample_gross %>% select(vertrauen_polizei, geschlecht) %>% filter(!is.na(geschlecht))ggplot(polizei, aes(vertrauen_polizei, fill = geschlecht)) +  geom_density(alpha = 0.5) +  ggtitle("Vertrauen in die Polizei nach Geschlecht") +  xlab("Vetrauen (1-7)") + ylab("Anteil der Respondenten") + labs(fill = "Geschlecht") 

3.4 Liniendiagramme

Zu den klassischen Plot-Typen gehören neben Barplots und Histogrammen auch Linien-, Punkt- und Streudiagramme, sowie Boxplots.

Liniendiagramme zeigen den Zusammenhang von zwei Variablen, in diesem Beispiel einer nominalen und einer ordinalen (oder, wenn wir großzügig sind [pseudo]metrischen) Variable, nämlich Wahlabsicht nach Partei und Vertrauen in die Presse.

vertrauen_nach_partei <- sample_gross %>%  rename(Partei = wahlabsicht_partei) %>%  group_by(Partei) %>%  summarise(Vertrauenswürdigkeit = mean(vertrauen_zeitungswesen, na.rm = T))ggplot(vertrauen_nach_partei, aes(Partei, Vertrauenswürdigkeit, group = 1)) +  geom_line() + geom_point(size = 3) +  theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1)) + ggtitle("Vertrauen in die Presse nach Partei")

Es lassen sich auch problemlos mehrere Geoms kominieren (hier: Linie und Punkte). Im folgenden Beispiel lässt sich durch eine absteigende Sortierung der Ergebnisse ein klareres Resultat erzielen.

vertrauen_nach_partei_sortiert <- vertrauen_nach_partei %>%  arrange(desc(Vertrauenswürdigkeit)) %>%  mutate(Rang = row_number())ggplot(vertrauen_nach_partei_sortiert, aes(reorder(Partei, Rang), Vertrauenswürdigkeit, group = 1)) +  geom_line() + geom_point(size = 3) +  theme(axis.text.x = element_text(angle = 90, vjust = 1, hjust = 1)) + ggtitle("Verrauen in die Presse nach Partei") + xlab("")

3.5 Streudiagramme

Ein weiterer Diagrammtyp, der häufig zum Einsatz kommt, ist das sog. Streudiagramm (scatter plot). Ein Streudiagramm, auch als Punktwolkebekannt, ist eine visuelle Darstellung von beobachteten Wertepaaren zweier statistischer Merkmale. Diese Wertepaare werden in ein kartesisches Koordinatensystem eingetragen, was eine Ansammlung von Punkten ergibt. Die Darstellung der Punkte kann mit verschiedenen Symbolen erfolgen.

Mit diesem Plottypen lassen sich die Beziehung mehrerer Variablen (i.d.R. zwei, jeweils auf der x/y-Achse) anschaulich darstellen und sowohl (lineare oder nicht-lineare) Zusammenhänge indetifizieren als auch Cluster bilder.

Wir erzeugen zunächst einen Data Frame, der Informationen zu Geschlecht, Alter, Einkommen und Bildung der StudienteilnehmerInnen enthält.

einkommen <- daten %>%  select(age, sex, educ, di01a) %>%  rename(alter = age, geschlecht = sex, bildung = educ, einkommen = di01a) %>%  replace_with_na_all(condition = ~.x < 0) %>%  mutate(geschlecht = as_factor(geschlecht), bildung = as_factor(bildung)) %>%  drop_na() %>%  slice_sample(n = 80)einkommen

# A tibble: 80 × 4 alter geschlecht bildung einkommen <dbl+lbl> <fct> <fct> <dbl+lbl> 1 47 MANN HOCHSCHULREIFE 9000 2 30 MANN HOCHSCHULREIFE 3300 3 52 MANN OHNE ABSCHLUSS 3000 4 34 FRAU HOCHSCHULREIFE 3100 5 49 MANN MITTLERE REIFE 2100 6 34 FRAU HOCHSCHULREIFE 1800 7 52 MANN VOLKS-,HAUPTSCHULE 4500 8 69 MANN MITTLERE REIFE 5000 9 56 MANN HOCHSCHULREIFE 4700 10 58 MANN HOCHSCHULREIFE 5000 # ℹ 70 more rows

Welche Beziehung lässt sich zwischen Alter (x-Achse) und dem Einkommen (y-Achse) der RespondentInnen feststellen?

ggplot(einkommen, aes(alter, einkommen)) + geom_point()

Wir entwickeln dieses Beispiel nun noch etwas weiter.

ggplot(einkommen, aes(alter, einkommen)) + geom_jitter(width = 1) + geom_smooth(method = 'lm', formula = 'y ~ x') +  ggtitle("Zusammenhang zwischen Alter und Nettoeinkommen") +  xlab("Alter") + ylab("Nettoeinkommen in Euro")

Bei der Linie, die wir mit dem Befehl geom_smooth gezeichnet haben, handelt es sich um eine Regressionsgerade. Der graue Bereich um die Gerade zeigt den lokalen Standardfehler an. Auf Regressionmodelle gehen wir zum Abschluss des Moduls noch intensiv ein. Schon jetzt kann man aber erahnen, was uns die Regressionsgerade illustriert.

Was, wenn wir mehr als drei Variablen (bzw. unterschiedliche Ausprägungen einer kategorialen Variable) darstellen wollen? Neben der Positionierung auf der x- und y-Achse und der Farbe können wir hier zusätzlich auch noch mit unterschiedlichen Formen arbeiten.

Das nachstehende Beispiel ist nicht unbedingt besonders informativ, zeigt aber das Prinzip nachvollziehbar auf.

ggplot(einkommen, aes(alter, einkommen, color = bildung, shape = geschlecht)) + geom_jitter(width = .5, height = .5, size = 2) + ggtitle("Zusammenhang zwischen Alter und Nettoeinkommen ") +  xlab("Alter") + ylab("Nettoeinkommen") +  labs(color = "Bildungsabschluss") + labs(shape = "Geschlecht")

3.6 Boxplots

Ein Box-Plot, auch als Box-Whisker-Plot oder Kastengrafik bezeichnet, ist ein grafisches Darstellungsinstrument, das verwendet wird, um die Verteilung von mindestens ordinalskalierten Merkmalen übersichtlich darzustellen. Es fasst robuste Streuungs- und Lagemaße in einer einzigen Grafik zusammen. Diese Darstellung bietet auf einen Blick Informationen darüber, wo die Datenwerte liegen und wie sie über diesen Bereich verteilt sind. Dies wird erreicht, indem die sogenannte Fünf-Punkte-Zusammenfassung, bestehend aus dem Median, den beiden Quartilen und den beiden Extremwerten, in einer grafischen Form präsentiert wird.

Zunächst legen wir noch einen weiteren Data Frame an, der Infortationen zur Demokratiezufriedenheit der RespondentInnen enthält.

demokratiezefriedenheit <- sample_gross %>%  mutate(altersgruppe = cut(alter,  breaks = c(0, 24, 34, 44, 54, 64, Inf), labels = c("18-24", "25-34", "35-44", "45-54", "55-64", "65+")), zufriedenheit_demokratie_zusammengefasst = as_factor(case_when( zufriedenheit_demokratie %in% c("SEHR ZUFRIEDEN",  "ZIEMLICH ZUFRIEDEN",  "ETWAS ZUFRIEDEN") ~ "eher zufrieden", zufriedenheit_demokratie %in% c("ETWAS UNZUFRIEDEN", "ZIEML. UNZUFRIEDEN", "SEHR UNZUFRIEDEN") ~ "eher unzufrieden"))) %>%  select(alter, altersgruppe, geschlecht, zufriedenheit_demokratie_zusammengefasst) %>%  drop_na()

Dann erstellen wir ein einfaches Boxplot, welches das Alter der eher zufriedenen Personen mit dem der eher unzufriedenen vergleicht.

ggplot(demokratiezefriedenheit, aes(zufriedenheit_demokratie_zusammengefasst, alter)) +  geom_boxplot()

Der Einsatz von Boxplots für den Gruppenverleich bietet sich oftmals an, wenn man verschiedene Lageparameter vergleichen möchte.

Nun erstellen wir einen anderen Data Frame, jetzt zum Vertrauen in gesellschaftliche Institutionen, um ein weiteres Beispiel für die Verwendung von Boxplots in den Blick zu nehmen.

vertrauen_summiert <- sample_gross %>%  rowwise() %>%  mutate(vertrauen_gesamt = sum(across(starts_with("vertrauen_")))) %>%  select(alter, geschlecht, entwicklung_kriminalitaet, vertrauen_gesamt) %>%  ungroup() %>% drop_na()

ggplot(vertrauen_summiert, aes(entwicklung_kriminalitaet, vertrauen_gesamt)) + geom_boxplot() + geom_jitter(alpha = 0.3) +  theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1)) + ggtitle("Kriminalitätseinschätzung und Vertrauen in gesellschaftliche Institutionen") +  xlab("Kriminalitätseinschätzung") + ylab("Vertrauen in gesellschaftliche Institutionen")

Wir speichern unser Plot zudem als PDF.

ggsave("Kriminalität_und_Vertrauen.pdf")

3.7 Facettierte Plots

Im letzten Abschnitt behandeln wir eine weitere nützliche Funktion von ggplot, um eine größere Zahl von Variablen zu visualisieren – die sog. Facettierung.

Dazu zählen wir zunächst wie zu Beginn des Kapitels mehrere Variablen aus, hier die Altersgruppe, das Geschlecht und die zuvor ermittelte Demokratiezufriedenheit als binäre Variable (eher zufrieden/eher unzufrieden).

demokratiezefriedenheit_gruppen <- demokratiezefriedenheit %>%  count(altersgruppe, geschlecht, zufriedenheit_demokratie_zusammengefasst)

Das nun folgenden Plot hat die Altersgruppe auf der X-Achse, die Anzahl der RespondentInnen auf der Y-Achse und das Geschlecht als Füllfarbe. Zu diesen drei Variablen kommt noch eine weitere dazu (Demokratiezufriedenheit) nach der facettiert wird.

ggplot(demokratiezefriedenheit_gruppen, aes(altersgruppe, n, fill = geschlecht)) + geom_bar(stat = "identity") +  facet_grid(cols = vars(zufriedenheit_demokratie_zusammengefasst)) +  theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1)) + ggtitle("Demokratiezufriedenheit nach Alter und Geschlecht") +  xlab("Altersgruppe") + ylab("StudienteilnehmerInnen")