5  Daten verbildlichen

5.1 Lernsteuerung

5.1.1 Standort im Lernpfad

Abb. Abbildung 1.2 zeigt den Standort dieses Kapitels im Lernpfad und gibt damit einen Überblick über das Thema dieses Kapitels im Kontext aller Kapitel.

5.1.2 Lernziele

• Sie können erläutern, wann und wozu das Visualisieren statistischer Inhalte sinnvoll ist.
• Sie kennen typische Arte von Datendiagrammen.
• Sie können typische Datendiagramme mit R visualisieren.
• Sie können zentrale Ergebnisse aus Datendiagrammen herauslesen.

5.1.3 Benötigte R-Pakete

library(tidyverse)
library(easystats)
library(DataExplorer)  # nicht vergessen zu installieren

5.1.4 Benötigte Daten

mariokart <- read.csv("https://vincentarelbundock.github.io/Rdatasets/csv/openintro/mariokart.csv")

5.1.5 Quiz zum Einstieg

Vielleicht fordert Sie die Lehrkraft zu einem Einstiegsquiz auf, etwas mittels der Plattform antworte.jetzt.¹ Alternativ überlegen Sie sich selber 10 Quiz-Aufgaben zum Stoff des letzten Kapitels.

5.1.6 Wozu das alles?

Quelle: GIPHY

🥷 Wir müssen die Galaxis retten, Kermit.

🐸 Schlock

5.2 Ein Dino sagt mehr als 1000 Worte

Es heißt, ein Bild sage mehr als 1000 Worte. Schon richtig, aber ein Dinosaurier sagt auch mehr als 1000 Worte, s. Abbildung 5.1.²

Abbildung 5.1: Dinosaurier und Kreis: Gleiche statistischen Kennwerte

In Abbildung 5.1 sieht man verschiedene “Bilder”, also Datensätze: etwa einen Dino und einmal einen Kreis. Obwohl die Bilder grundverschiedene sind, sind die zentralen statistischen Kennwerte (praktisch) identisch.

In die gleiche Bresche schlägt “Anscombes Quartett” (Anscombe, 1973), s. Abbildung 5.2: Es zeigt vier Datensätze, in denen die zentralen Statistiken fast identisch sind,
also Mittelwerte, Streuungen, Korrelationen. Aber die Streudiagramme sind grundverschieden.³

Anscombes Beispiel zeigt (zugespitzt): Eine Visualisierung enthüllt, was der Statistik (als Kennzahl) verhüllt bleibt.

Wichtig

Statistische Diagramme können Einblicke geben, die sich nicht (leicht) in grundlegenden Statistiken (Kennwerten) abbilden. \(\square\)

Abbildung 5.2: Anscombes Quartet: Gleiche statistischen Kennwerte in vier Datensätzen

Unter visueller Cortex ist sehr leistungsfähig. Wir können ohne Mühe eine große Anzahl an Informationen aufnehmen und parallel verarbeiten. Aus diesem Grund sind Datendiagramme eine effektive und einfache Art, aus Daten Erkenntnisse zu ziehen.

Tipp

Nutzen Sie Datendiagramme umfassend; sie sind einfach zu verstehen und doch sehr mächtig.

5.2.1 Datendiagramm

Ein Datendiagramm (kurz: Diagramm) ist ein Diagramm, das Daten und Statistiken zeigt, mit dem Zweck, Erkenntnisse daraus zu ziehen.

Beispiel 5.1 (Aus der Forschung: Ein aufwändiges (und ansprechendes) Datendiagramm) Hier finden Sie ein Beispiel für ein Datendiagramm, das mit R erzeugt wurde (Scherer et al., 2019).

5.2.2 Ein Bild hat nicht so viele Dimensionen

Abbildung 5.3 zeigt ein Bild mit mehreren (5) Variablen, die jeweils einer “Dimension” entsprechen. Wie man (nicht) sieht, wird es langsam unübersichtlich. Offenbar kann man in einem Bild nicht beliebig viele Variablen sinnvoll reinquetschen. Die “Dimensionalität” eines Diagramms hat ihre Grenzen, vielleicht bei 4-6 Variablen.

Abbildung 5.3: Ein Diagramm kann nur eine begrenzte Anzahl von Variablen zeigen. Wenn Sie dieses Bild nicht checken: Prima. Genau das soll das Bild zeigen.

Möchten wir den Zusammenhang von vielen Variablen, z.B. mehr als 5, verstehen, kommen wir mit Bildern nicht weiter. Dann brauchen wir andere Werkzeuge: statistics to the rescue.

Hinweis

Bei klaren Zusammenhängen und wenig Variablen braucht man keine (aufwändige) Statistik. Ein Bild (Datendiagramm) ist dann (oft) ausreichend. Man könnte sagen, dass es Statistik nur deshalb gibt, weil unser Auge mit mehr als ca. 4-6 Variablen nicht gleichzeitig umgehen kann.

Übungsaufgabe 5.1 Wie viele Variablen sind in Abbildung 5.3 dargestellt?⁴

Eine weitere Möglichkeit, mehr Variablen in einem Diagramm unterzubringen, ist, die “Flatlands” zu verlassen, also von 2D auf 3D zu wechseln, s. hier, s. Abbildung 5.4.⁵

Abbildung 5.4: Eine 3D-Karte der Erde

5.3 Nomenklatur von Datendiagrammen

Tabelle 5.1 zeigt eine - sehr kurze Nomenklatur - an Datendigrammen.⁶

Tabelle 5.1: Ein (sehr kurze) Nomenklatur von Datendiagrammen

Erkenntnisziel	qualitativ	quantitativ
Verteilung	Balkendiagramm	Histogramm und Dichtediagramm
Zusammenhang	gefülltes Balkendiagramm	Streudiagramm
Unterschied	gefülltes Balkendiagramm	Boxplot

Hinweis

Wir arbeiten hier mit dem Datensatz mariokart. Hilfe bzw. ein Data-Dictionary (Codebook) finden Sie hier.

5.4 Verteilungen verbildlichen

5.4.1 Verteilung: nominale Variable

Definition 5.1 (Verteilung) Eine (Häufigkeits-)Verteilung einer Variablen \(X\) schlüsselt auf, wie häufig jede Ausprägung von \(X\) ist.\(\square\)

Beispiel 5.2 Tabelle 5.2 zeigt die Häufigkeitsverteilung von cond aus dem Datensatz mariokart. Die Variable hat 5 Ausprägungen; z.B. kommt die Ausprägung new 59 mal vor.\(\square\)

Tabelle 5.2: Häufigkeitsverteilung von cond aus dem Datensatz mariokart

cond	n
new	59
used	84

Zugegeben, das Datendiagramm von cond ist nicht so aufregend, s. Abbildung 5.5. Wie man sieht, besteht so ein Diagramm als Balken, daher heißt es Balkendiagramm⁷. Man kann so ein Diagramm um 90° drehen; keine Ausrichtung ist unbedingt besser als die andere.

Definition 5.2 (Balkendiagramm) Ein Balkendiagramm eignet sich, um Häufigkeiten darzustellen

Abbildung 5.5: Häufigkeitsverteilung der Variable cond

Abbildung 5.6: Balkendiagramm mit dem R-Paket DataExplorer

Es gibt viele Methoden, sich mit R ein Balkendiagramm ausgeben zu lassen. Eine einfache, komfortable ist die mit dem Paket DataExplorer, s. Abbildung 5.6.

Zuerst importieren wir die Daten, s. Listing 5.1.

Listing 5.1: Mariokart-Daten importieren von einer Webseite

mariokart <- read.csv("https://vincentarelbundock.github.io/Rdatasets/csv/openintro/mariokart.csv")

Außerdem nicht vergessen, das Paket DataExplorer zu starten, s. Listing 5.2.⁸ In diesem Paket “wohnen” die Befehle, die wir zum Erstellen der Datendiagramme nutzen werden. Listing 5.3 zeigt die Syntax, um ein Balkendiagramm zu erstellen.⁹

Listing 5.2: Wir starten das R-Paket DataExplorer

library(DataExplorer)

Listing 5.3: Syntax zur Erstellung eines Balkendiagramms

mariokart %>% 
  select(cond) %>% 
  plot_bar()

Die Syntax ist in Listing 5.3 abgedruckt¹⁰ . Übersetzen wir die Syntax ins Deutsche:

Nimm den Datensatz mariokart *und dann*
  wähle die Spalte cond *und dann*
  zeichne ein Balkendiagramm.

Übungsaufgabe 5.2 (Spalten wählen für das Balkendiagramm) Hätten wir andere Spalten ausgewählt, so würde das Balkendiagramm die Verteilung jener Variablen zeigen. Ja, Sie können auch mehrere Variablen auf einmal auswählen. Probieren Sie das doch mal aus!

Übungsaufgabe 5.3  

Aufgabe

Visualisieren Sie die Verteilung von stock_photo!

Lösung

mariokart |> 
  select(stock_photo) |> 
  plot_bar()

5.4.2 Verteilung: quantitative Variable

5.4.2.1 Histogramm

Bei einer quantitativen Variablen mit vielen Ausprägungen wäre ein Balkendiagramm nicht so aussagekräftig, s. Abbildung 5.7. Es gibt einfach zu viele Ausprägungen.

Abbildung 5.7: Balkendiagramm für total_pr

Die Lösung: Wir reduzieren die Anzahl der Ausprägungen, in dem wir auf ganze Dollar runden. Oder, um noch weniger Ausprägungen zu bekommen, können wir einfach Gruppen definieren, z.B.

• Gruppe 1: 0-5 Dollar
• Gruppe 2: 6-10 Dollar
• Gruppe 2: 11-15 Dollar …

In Abbildung 5.8 sind z.B. die Ausprägungen des Verkaufspreis (total_pr) in in Gruppen der Breite von 5 Dollar aufgeteilt worden. Zusätzlich sind noch die einzelnen Werte als schwarze Punkte gezeigt.

Abbildung 5.8: Balkendiagramm für total_pr

Definition 5.3 (Histogramm) Ein Histogramm ist ein Diagramm zur Darstellung der Häufigkeitsverteilung einer quantitativen Variablen. Die Daten werden in Gruppen (Klassen) eingeteilt, die dann durch einen Balken (pro Klasse) dargestellt sind. Die Höhe der Balken zeigt die Häufigkeit der Daten in dieser Gruppe/in diesem Balken¹¹.

Es gibt keine klare Regel, in wie viele Balken ein Histogramm gegliedert sein sollte. Nur: Es sollten nicht sehr viele und nicht sehr wenig sein, s. Abbildung 5.9 links bzw. Abbildung 5.9, rechts.

(a) Zu viele Gruppen (Balken)

(b) Zu wenige Gruppen (Balken)

Abbildung 5.9: Nicht zu wenig und nicht zu viele Balken im Balkendiagramm

Zur Erstellung eines Histogramms können Sie die Syntax Listing 5.4 nützen, vgl. Abbildung 5.10, links.

Listing 5.4: Syntax zur Erstellung eines Histogramms

mariokart %>% 
  select(total_pr) %>% 
  filter(total_pr < 100) %>%  # ohne Extremwerte
  plot_histogram()

(a) Histogramm

(b) Dichtediagramm

Abbildung 5.10: Eine stetige Verteilung verbildlichen

Übungsaufgabe 5.4  

Aufgabe

Visualisieren Sie die Verteilung von ship_pr anhand eines Histogramms!

Lösung

mariokart |> 
  select(ship_pr) |> 
  plot_histogram()

5.4.2.2 Dichtediagramm

Abbildung 5.11 fügt zu Abbildung 5.8 ein Dichtediagramm hinzu (rote Linie). Ein Dichtediagramm ähnelt einem “glattgeschmirgeltem” Histogramm.

Definition 5.4 (Dichtediagramm) Ein Dichtediagramm visualisiert die Verteilung einer stetigen Variablen. Im Gegensatz zum Histogramm wird der Verlauf der Kurve geglättet, so kann Rauschen (Zufallsschwankung) besser ausgeblendet werden.¹²

Abbildung 5.11: Histogramm (graue Balken) und Dichtediagramm (orange Linie) für total_pr

Übungsaufgabe 5.5 Erstellen Sie das Diagramm Abbildung 5.10, rechtes Teildiagramm!¹³\(\square\)

5.4.2.3 Eigenschaften von Verteilungen

Verteilungen unterscheiden sich z.B. einerseits in ihrem “typischen” oder “mittleren” Wert¹⁴ und anderseits in ihrer Streuung¹⁵

(Diagramme von) Verteilungen können symmetrisch oder schief (nicht symmetrisch) sein, s. Abbildung 5.12.

(a) Symmetrisch (Normal)

(b) Schief

Abbildung 5.12: Symmetrische vs. schiefe Verteilung, verbildlicht

Abbildung 5.13 zeigt verschiedene Formen von Verteilungen. “Bimodal” meint “zweigipflig” und “multimodal” entsprechend “mehrgipflig”.¹⁶

Abbildung 5.13: Verschiedene Verteilungsformen

Übungsaufgabe 5.6  

Aufgabe

Bennenen Sie die am besten passende Verteilungsform für die Variable total_pr.

Lösung

Die Verteilung ist rechtschief.

mariokart |> 
  select(total_pr) |> 
  plot_density()

5.4.3 Normalverteilung

Eine Normalverteilung ist eine bestimmte Art von Verteilung einer quantitativen Variablen. Aber sie ist besonders wichtig, und ist daher hier herausgestellt.

Eine Normalverteilung sehen Sie in Abbildung 5.12, links. Sie hat u.a. folgende Eigenschaften:

• symmetrisch
• glockenförmig
• stetig
• eingipflig (unimodal)
• Mittelwert, Median und Modus sind identisch

Beispiel 5.3 Beispiele für normalverteilte Variablen sind Körpergröße von Männern oder Frauen, IQ-Werte, Prüfungsergebnisse, Messfehler, Lebensdauer von Glühbirnen, Gewichte von Brotlaiben, Milchproduktion von Kühen, Brustumfang schottischer Soldaten (Lyon, 2014).\(\square\)

Die Normalverteilung ist von hoher Bedeutung, da sich diese Verteilung unter (recht häufigen) Bedingungen zwangsläufig ergeben muss.

Definition 5.5 (Entstehung einer Normalverteilung) Wenn sich eine Variable \(X\) als Summe mehrerer, unabhängiger, etwa gleich starker Summanden, dann kann man erwarten, dass sich diese Variable \(X\) tendenziell normalverteilt. \(\square\)

Dieses Phänomen kann man gut anhand des Galton-Bretts veranschaulichen.

Parameter der Normalverteilung

Eine Normalverteilung lässt sich exakt beschreiben anhand zweier Parameter: ihres zentralen Werts (Mittelwerts), \(\mu\), und ihrer Streuung (Standardabweichung), \(\sigma\). \(\square\)

Abbildung 5.14 zeigt interaktive Beispiele für Normalverteilung. Wählen Sie einfach Mittelwert (\(\mu\)) und Streuung (\(\sigma\)) anhand der Schieberegler.¹⁷

sliders = {
  let div = d3.create("div");

  let m0 = d3.mean(pts);
  let s0 = d3.deviation(pts);
  let mu = Inputs.range([1, 8], {
    value: m0,
    step: 0.001,
    format: d3.format("0.3f"),
    label: tex`\large\pmb\mu:`
  });
  let sigma = Inputs.range([0.2, 4], {
    value: s0,
    step: 0.001,
    format: d3.format("0.3f"),
    label: tex`\large\pmb\sigma:`
  });

  d3.select(mu).on("input", redraw).select("label").style("width", "30px");
  d3.select(sigma).on("input", redraw).select("label").style("width", "30px");

  div.append(() => mu);
  div.append(() => sigma);

  return div.node();

  function redraw() {
    let m = mu.value;
    let s = sigma.value;
    d3.select(normal_model).select("svg").remove();
    let standardized = pts.map((x) => (x - m0) / s0);
    let new_pts = standardized.map((z) => z * s + m);
    let new_plot = create_plot(new_pts);
    d3.select(normal_model).append(() => new_plot);
  }
}

viewof steely_dan_says = Inputs.button("Neuer Zufallsversuch")

normal_model = {
  let div = d3.create("div");
  let plot = create_plot(pts);

  d3.select(plot).selectAll("circle").attr("opacity", 0);

  let initials = d3
    .select(plot)
    .selectAll("rect")
    .nodes()
    .map((r) => ({ height: r.getAttribute("height"), y: r.getAttribute("y") }));
  let y_scale = plot.scale("y");
  d3.select(plot)
    .selectAll("rect")
    .attr("height", 0)
    .attr("y", y_scale.apply(0));
  d3.select(plot).select("path").attr("opacity", 0);
  Promises.delay(500).then(function () {
    d3.select(plot)
      .selectAll("circle")
      .attr("opacity", 0)
      .transition()
      .duration(1000)
      .attr("opacity", 0.0);
  });
  Promises.delay(1500).then(function () {
    d3.select(plot)
      .selectAll("rect")
      .attr("height", 0)
      .attr("y", y_scale.apply(0))
      .transition()
      .duration(850)
      .attr("height", (d, i) => initials[i].height)
      .attr("y", (d, i) => initials[i].y);
  });
  if (show_curve) {
    Promises.delay(1500).then(function () {
      d3.select(plot)
        .selectAll("path")
        .attr("opacity", 0)
        .transition()
        .duration(1000)
        .attr("opacity", 0.8);
    });
  }

  div.append(() => plot);

  return div.node();
}

pts = {
  steely_dan_says;
  let n = 1000;
  let m0 = d3.randomUniform(1, 8)();
  let s0 = d3.randomUniform(1 / 2, 2)();
  let pts = d3.range(n).map(d3.randomNormal(m0, s0));

  return pts;
}

create_plot = function (pts) {
  let m = d3.mean(pts);
  let s = d3.deviation(pts);

  let w = 800;
  let h = 0.4 * w;

  let f = (x) =>
    Math.exp((-(x - m) * (x - m)) / (2 * s * s)) / (Math.sqrt(2 * Math.PI) * s);

  let marks = [
    Plot.rectY(
      pts,
      Plot.binX(
        {
          y: (a, bin) => {
            return a.length / pts.length / (bin.x2 - bin.x1);
          },
          title: "proportion"
        },
        { x: (pt) => pt, fill: "#b00" }
      )
    ),
    Plot.dot(pts, {
      x: (x) => x,
      y: (_) => 0,
      stroke: "black",
      fill: "black",
      opacity: 0.2
    }),
    Plot.ruleX([0]),
    Plot.ruleY([0])
  ];
  if (show_curve) {
    marks.push(
      Plot.line(build_samples(f, -1, 12, { N: 100 }), {
        strokeWidth: 5,
        stroke: "#111",
        opacity: 0
      })
    );
  }

  let plot = Plot.plot({
    x: { domain: [0, 11] },
    y: { domain: [0, 1] },
    width: w,
    height: h,
    marks: marks
  });

  d3.select(plot)
    .selectAll("rect")
    .on("pointerenter", function () {
      d3.select(this).attr("opacity", 0.5);
    })
    .on("pointerleave", function () {
      d3.select(this).attr("stroke", null).attr("opacity", null);
    })
    .nodes()
    .forEach((bar) =>
      tippy(bar, { content: d3.select(bar).select("title").text() })
    );
  d3.select(plot).selectAll("rect").select("title").remove();
  return plot;
}

show_curve = true

import { build_samples } from '@mcmcclur/adaptive-plotter'

tippy = require("tippy.js@6")

Abbildung 5.14: Interaktives Beispiel für Normalverteilungen.

Kennt man diese beiden Parameter, so kann man einfach angeben, welcher Anteil der Fläche sich in einem bestimmten Bereich befindet, s. Abbildung 5.15.¹⁸

Davon leitet sich die “68-95-99-Prozentregel” ab:

• \(68\,\%\) der Werte im Bereich \(\mu\pm 1 \cdot \sigma\)
• \(95\,\%\) der Werte im Bereich \(\mu\pm 2 \cdot \sigma\)
• \(99{,}7\,\%\) der Werte im Bereich \(\mu\pm 3 \cdot \sigma\)

Abbildung 5.15: Die Flächeninhalte (Wahrscheinlichkeitsmasse) einer Normalverteilung in Abhängigkeit der SD-Einheiten

5.5 Zusammenhänge verbildlichen

5.5.1 Zusammenhang: nominale Variablen

Beispiel 5.4 (Beispiele für Zusammenhänge bei nominalen Variablen)  

• Hängt Berufserfolg (Führungskraft ja/nein) mit dem Geschlecht zusammen?
• Hängt der Beruf des Vaters mit dem Schulabschluss des Kindes (Abitur, Realschule, Mittelschule) zusammen?
• Gibt es einen Zusammenhang zwischen Automarke und politische Präferenz einer Partei? \(\square\)

Sagen wir, Sie arbeiten immer noch beim Online-Auktionshaus und Sie fragen sich, ob ein Produktfoto wohl primär bei neuwertigen Produkten beiliegt, aber nicht bei gebrauchten? Dazu betrachten Sie wieder die mariokart-Daten, s. Abbildung 5.16.

(a) Es findet sich ein Zusammenhang von Foto und Zustand in den Daten

(b) Es findet sich (fast) kein Zusammenhang von wheel und Foto in den Daten

Abbildung 5.16: Zusammenhang zwischen nominalskalierten Variablen verbildlichen

Tatsächlich: Es findet sich ein Zusammenhang zwischen der Tatsache, ob dem versteigerten Produkt ein Foto bei lag und ob es neuwertig oder gebraucht war (Abbildung 5.16, links). Bei neuen Spielen war fast immer (ca. 90%) ein Foto dabei; bei gebrauchten Spielen immerhin bei gut der Hälfte der Fälle.

Anders sieht es aus für die Frage, ob ein (oder mehrere) Lenkräder dem Spiel beilagen (oder nicht) in Zusammenhang mit der Fotofrage Hier gab es fast keinen Unterschied zwischen neuen und alten Spielen, was die Frage nach “Foto des Produkts dabei” betraf (Abbildung 5.16, rechts), der Anteil betrug jeweils ca. 70%. Das zeigt, dass es keinen Zusammenhang zwischen Foto und Neuwertigkeit des Spiels gibt (laut unseren Daten).

Anders gesagt: Unterscheiden sich die “Füllhöhe” in den Diagrammen, so gibt es einen Unterschied hinsichtlich “Foto ist dabei” zwischen den beiden Gruppen (linker vs. rechter Balken). Unterscheiden sich die Anteile in den Gruppen (neuwertige vs. gebrauchte Spiele), so spielt z.B. die Variable “Foto dabei” offenbar eine Rolle. Dann hängen Neuwertigkeit und “Foto dabei” also zusammen!

So können Sie sich in R ein gefülltes Balkendiagramm ausgeben lassen, s. Abbildung 5.17. Diese Darstellung eignet sich, um Zusammenhänge zwischen zwei zweistufigen nominal skalierten Variablen zu verbildlichen. Die verschiedenen Werte der Füllfarbe werden den Stufen der Variablen cond zugewiesen, s. Listing 5.5.

Listing 5.5: R-Syntax für ein gefülltes Balkendiagramm

mariokart %>% 
  select(cond, stock_photo) %>% 
  plot_bar(by = "cond")  # aus dem Paket DataExplorer

Abbildung 5.17: Ein gefülltes Balkendiagramm zur Untersuchung eines Zusammenhangs zwischen nominalskalierter Variablen

Hinweis

Gefüllte Balkendiagramme eignen sich zur Analyse eines Zusammenhangs zwischen nominalskalierten Variablen. Allerdings sollte eine der beiden Variablen nur zwei Ausprägungen aufweisen, sonst sind die Zusammenhänge nicht mehr so gut zu erkennen.\(\square\)

Übungsaufgabe 5.7  

Aufgabe

Visualisieren Sie den Zusammenhang der beiden nominalen Variablen cond und wheels!

Lösung

wheels ist als metrische Variable (int: Integer, d.h. Ganzzahl) formatiert im Datensatz mariokart. Wir müssen Sie zunächst als Faktorvariable umformatieren, damit R sie als nominal skalierte Variable erkennt.

mariokart |> 
  # Mache aus einer metrischen eine nominale Variable: 
  mutate(wheels = factor(wheels)) |> 
  select(cond, wheels) |> 
  plot_bar(by = "cond")

5.5.2 Zusammenhang: metrisch

Den (etwaigen) Zusammenhang zweier metrischer Variablen kann man mit einem Streudiagramm visualisieren (engl. scatterplot). Abbildung 5.18 links untersucht den Zusammenhang des Einstiegpreises (X-Achse) und Abschlusspreises (Y-Achse) von Geboten bei Versteigerungen des Computerspiels Mariokart. In dem Diagramm ist eine “Trendgerade” (Regressionsgerade), um die Art des Zusammenhangs besser zu verdeutlichen. Die Trendgerade steigt an (von links nach recht). Daraus kann man schließen: Es handelt sich um einen gleichsinnigen (positiven) Zusammenhang: Je höher der Startpreis, desto höher der Abschlusspreis, zumindest tendenziell. Diese Gerade liegt “mittig” in den Daten (wir definieren dies später genauer). Diese Trendgerade gibt Aufschluss über “typische” Werte: Welcher Y-Wert ist “typisch” für einen bestimmten X-Wert?

Abbildung 5.18 rechts untersucht den Zusammenhang zwischen Anzahl der Gebote (X-Achse) und Abschlusspreises (Y-Achse). Es handelt sich um einen negativen Zusammenhang: Je mehr Gebote, desto geringer der Abschlusspreis. Das erkennt man an der sinkenden Trendgeraden.

Die Ellipse zeigt an, wie eng die Daten um die Trendgerade streuen. Daraus kann man ableiten, wie stark der Absolutwert des Zusammenhangs ist, vgl. Abbildung 5.20.

(a) positiver, mittelstarker Zusammenhang

(b) negativer, eher schwacher Zusammenhang

Abbildung 5.18: Streudiagramm zur Darstellung eines Zusammenhangs zweier metrischer Variablen

Definition 5.6 (Linearer Zusammenhang) Lässt sich die Beziehung zwischen zwei Variablen mit einer Gerade visualisieren, so spricht man von einem linearen Zusammenhang. Ändert man eine der beiden Variablen um einen bestimmten Wert (z.B. 1), so ändert sich die andere um einen proportionalen Wert (z.B. 0.5). \(\square\)

Natürlich könnte man auch nicht-lineare Zusammenhänge untersuchen, aber der Einfachheit halber konzentrieren wir uns hier mit linearen; Beispiele für nicht-lineare Zusammenhänge sind in Abbildung 5.19 zu sehen.

Abbildung 5.19: Beispiele nichtlinearer Zusammenhänge

Definition 5.7 (Richtig und Stärke eines Zusammenhang) Gleichsinnige (positive) Zusammenhänge erkennt man an aufsteigenden Trendgeraden; gegensinnigen (negative) Zusammenhänge an absteigenden Trendgeraden:

• ➕ : ⬆️
• ➖ : ⬇️

Starke Zusammenhänge erkennt man an schmalen Ellipsen (“Baguette”); schwache Zusammenhänge an breiten Ellipsen (“Torte”):

• schwach: 🥮
• stark: 🥖

\(\square\)

Abbildung 5.20 bietet einen Überblick über verschiedene Beispiele von Richtung und Stärke von Zusammenhängen.¹⁹

Abbildung 5.20: Lineare Zusammenhänge verschiedener Stärke und Richtung

In Abbildung 5.20 ist für jedes Teildiagramm eine Zahl angegeben: der Korrelationskoeffizient. Diese Statistik quantifiziert Richtung und Stärke des Zusammenhangs (mehr dazu in Kap. Kapitel 8). Ein positives Vorzeichen steht für einen positiven Zusammenhang, ein negatives Vorzeichen für einen negativen Zusammenhang. Der (Absolut-)Wert gibt die Stärke des linearen Zusammenhangs an (Cohen, 1992):

• ±0: Kein Zusammenhang
• ±0.1: schwacher Zusammenhang
• ±0.3: mittlerer Zusammenhang
• ±0.5: starker Zusammenhang
• ±1: perfekter Zusammenhang

Abbildung 5.21 hat die gleiche Aussage, ist aber plakativer, indem Stärke (schwach, stark) und Richtung (positiv, negativ) gegenübergestellt sind.

Abbildung 5.21: Überblick über starke vs. schwache bzw. positive vs. negative Zusammenhänge

Man sieht in Abbildung 5.20 und Abbildung 5.21, dass ein negativer Korrelationskoeffizient mit einer absinkenden Trendgerade ²⁰ (blaue Linie) einhergeht. Umgekehrt geht ein positiver Trend mit einer ansteigenden Trendgerade einher. Zweitens erkennt man, dass starke Zusammenhänge mit einer schmaler Ellipse einhergehen und schwache Zusammenhänge mit einer breiten Ellipse einhergehen.

Abbildung 5.22 zeigt interaktive Beispiele für (lineare) Zusammenhänge.²¹

viewof cor_type = select({
  title: "Correlation type",
  options: [
    "A perfect linear relationship",
    "A close to linear relationship",
    "A close to linear, but negative, relationship",
    "A weaker relationship",
    "A nonlinear relationship",
    "No relationship"
  ],
  value: "A perfect linear relationship"
})

viewof redo = Inputs.button("Redo")

pic = (redo, graph_from_type(cor_type))

function graph_from_type(s) {
  if (s == "A perfect linear relationship") {
    let a = jstat.uniform.sample(1 / 3, 2 / 3);
    let b = jstat.uniform.sample(-6, 6);
    return make_perturbed_graph(
      (x) => a * x + b,
      (x) => 0,
      -10,
      10,
      100
    );
  } else if (s == "A close to linear relationship") {
    let a = jstat.uniform.sample(1 / 3, 2 / 3);
    let b = jstat.uniform.sample(-6, 6);
    return make_perturbed_graph(
      (x) => a * x + b,
      (x) => jstat.randn(),
      -10,
      10,
      100
    );
  } else if (s == "A close to linear, but negative, relationship") {
    let a = jstat.uniform.sample(1 / 3, 2 / 3);
    let b = jstat.uniform.sample(-6, 6);
    return make_perturbed_graph(
      (x) => -a * x + b,
      (x) => jstat.randn(),
      -10,
      10,
      100
    );
  } else if (s == "A weaker relationship") {
    let a = jstat.uniform.sample(-2 / 3, 2 / 3);
    let b = jstat.uniform.sample(-6, 6);
    return make_perturbed_graph(
      (x) => a * x + b,
      (x) => 4 * jstat.randn(),
      -10,
      10,
      100
    );
  } else if (s == "A nonlinear relationship") {
    let a = jstat.uniform.sample(-3, 1);
    let b = jstat.uniform.sample(1, 3);
    return make_perturbed_graph(
      (x) => (x - a) * (x - b),
      (x) => 0.6 * jstat.randn(),
      -10,
      10,
      100
    );
  } else if (s == "No relationship") {
    let a = jstat.uniform.sample(-3, 1);
    let b = jstat.uniform.sample(1, 3);
    return make_perturbed_graph(
      (x) => 0,
      (x) => jstat.uniform.sample(-10, 10),
      -10,
      10,
      100
    );
  }
}

function make_perturbed_graph(f, r, a, b, n) {
  let xs = jstat.arange(n).map(() => jstat.uniform.sample(a, b));
  let ys = xs.map((x) => f(x) + r());
  let plot = plotter({ width: 500, height: 400, grid: false });
  jstat.arange(n).forEach((_, i) => plot.point(xs[i], ys[i]));

  let R = jstat.corrcoeff(xs, ys);

  return html`<div style="text-align:center; width:500px">R = ${d3.format(
    "0.4f"
  )(R)}</div>${plot.node}`;
}

import {select} from "@jashkenas/inputs"
jstat = require('jstat')
import { plotter } from '50dadfdec01c15a8'
import { rk4 } from '@mcmcclur/runge-kutta-for-systems-of-odes'
import { slider } from "@jashkenas/inputs"

Abbildung 5.22: Interaktives Beispiel für Zusammenhangsdiagramme.

Beispiel 5.5 Sie arbeiten nach wie vor bei einem Online-Auktionshaus, und machmal gehört Datenanalyse zu Ihren Aufgaben. Daher interessiert Sie, ob welche Variablen mit dem Abschlusspreis (total_pr) im Datensatz mariokart zusammenhängen. Sie verbildlichen die Daten mit R, und zwar nutzen Sie das Paket DataExplorer. Starten Sie dieses Paket, s. Listing 5.2. Außerdem müssen wir noch die Daten importieren, falls noch nicht getan, s. Listing 5.1.

So, jetzt kann die eigentliche Arbeit losgehen. Da Sie sich nur auf metrische Variablen konzentrieren wollen, wählen Sie (mit select) nur diese Variablen aus. Dann weisen Sie R an, einen Scatterplot zu malen (plot_scatterplot) und zwar jeweils den Zusammenhang einer der gewählten Variablen mit dem Abschlusspreis (total_pr), da das die Variable ist, die Sie primär interessiert. Das Ergebnis sieht man in Abbildung 5.23 bzw. Listing 5.6.

Listing 5.6: Streudiagramm erstellen mit dem R-Paket ‘DataExplorer’

mariokart %>% 
  select(duration, n_bids, start_pr, ship_pr, total_pr, seller_rate, wheels) %>% 
  plot_scatterplot(by = "total_pr")

Abbildung 5.23: Der Zusammenhang metrischer Variablen mit Abschlusspreis

Aha… Was sagt uns das Bild? Hm. Es scheint einige Extremwerte zu geben, die dafür sorgen, dass der Rest der Daten recht zusammengequetscht auf dem Bild erscheint. Vielleicht sollten Sie solche Extremwerte lieber entfernen? Sie entscheiden sich, nur Verkäufe mit einem Abschlusspreis von weniger als 100 Dollar anzuschauen (total_pr < 100). Das Ergebnis ist in Abbildung 5.24 zu sehen.

mariokart_no_extreme <-
  mariokart %>% 
  filter(total_pr < 100)

mariokart_no_extreme %>% 
  select(duration, n_bids, start_pr, ship_pr, total_pr, seller_rate, wheels) %>% 
  plot_scatterplot(by = "total_pr")

Abbildung 5.24: Der Zusammenhang metrischer Variablen mit Abschlusspreis

Ohne Extremwerte schält sich ein deutlicheres Bild (Abbildung 5.24) hervor: Startpreis (start_pr) und Anzahl der Räder (wheels) scheinen am stärksten mit dem Abschlusspreis zusammenzuhängen.

Das Argument by = "total_pr" bei plot_scatterplot weist R an, als Y-Variable stets total_pr zu verwenden. Alle übrigen Variablen kommen jeweils einmal als X-Variable vor.\(\square\)

Übungsaufgabe 5.8  

Aufgabe

Visualisieren Sie den Zusammenhang der beiden metrischen Variablen start_pr und total_pr. Verwenden Sie den Datensatz ohne Extremwerte wie oben definiert.

Lösung

mariokart_no_extreme |> 
  select(start_pr, total_pr) |> 
  plot_scatterplot(by = "total_pr")

5.6 Unterschiede verbildlichen

5.6.1 Unterschied: nominale Variablen

Gute Nachrichten: Für nominale Variablen bieten sich Balkendiagramme sowohl zur Darstellung von Zusammenhängen als auch von Unterschieden an. Genau genommen zeigt ja Abbildung 5.16 (links) den Unterschied zwischen neuen und gebrauchten Spielen hinsichtlich der Frage, ob Photos beiliegen. Und wie man in Abbildung 5.16 sieht, ist der Anteil der Spiele mit Foto bei den neuen Spielen höher als bei gebrauchten Spielen.²²

5.6.2 Unterschied: quantitative Variablen

Eine typische Analysefrage ist, ob sich zwei Gruppen hinsichtlich einer metrischen Zielvariablen deutlich unterscheiden. Genauer gesagt untersucht man z.B. oft, ob sich die Mittelwerte der beiden Gruppen zwischen der Zielvariablen deutlich unterscheiden. Das hört sich abstrakt an? Am besten wir schauen uns einige Beispiele an, s. Abbildung 5.25.

(a) Histogramm pro Gruppe

(b) Boxplot pro Gruppe

Abbildung 5.25: Unterschiede zwischen zwei Gruppen: Metrische Y-Variable, nominale X-Variable

Das linke Teildiagramm von Abbildung 5.25 zeigt das Histogramm von total_pr, getrennt für neue und gebrauchte Spiele, vgl. Abbildung 5.10. Das rechte Teildiagramm zeigt die gleichen Verteilungen, aber mit einer vereinfachten, groberen Darstellungsfrom, den Boxplot.²³

(a) Y: Abschlusspreis, X: Zustand

(b) Y: Abschlusspreis, X: Photo dabei?

Abbildung 5.26: Unterschiede zwischen zwei Gruppen: Metrische Y-Variable, nominale X-Variable

Das linke Teildiagramm von Abbildung 5.26 zeigt den Unterschied in den Verteilungen von total_pr, einmal für die neuen Computerspiele (cond == new) und einmal für gebrauchte Spiele (cond == used).

Was ein “deutlicher”²⁴ Zusammenhang ist, ist keine statistische, sondern inhaltliche Frage, die man mit Sachverstand zum Forschungsgegenstand beantworten muss.

Definition 5.8 (Boxplot) Der Boxplot ist eine Vereinfachung bzw. eine Zusammenfassung eines Histograms.²⁵ Damit stellt der Boxplot auch eine Verteilung (einer metrischen Variablen) dar.\(\square\)

In Abbildung 5.27 sieht man die “Übersetzung” von Histogramm (oben) zu einem Boxplot (unten).

Abbildung 5.27: Übersetzung eines Histogramms zu einem Boxplot

Schauen wir uns die “Anatomie” des Boxplots näher an:

1. Der dicke Strich in der Box zeigt den Median der Verteilung, vgl. Kapitel 6.3.
2. Die Enden der Box zeigen das 1. Quartil (41) bzw. das 3. Quartil (54). Damit zeigt die Breite der Box die Streuung der Verteilung an, genauer gesagt die Streuung der inneren 50% der Beobachtungen. Je breiter die Box, desto größer die Streuung. Die Breite der Box nennt man auch den Interquartilsabstand (IQR).
3. Die “Antennen” des Boxplots zeigen die Streuung in den kleinsten 25% der Werte (linke Antenne) bzw. die Streuung der größten 25% der Werte (rechte Antennen). Je länger die Antenne, desto größer die Streuung.
4. Falls es aber Extremwerte gibt, so sollten die lieber einzeln, separat, außerhalb der Antennen gezeigt werden. Daher ist die Antennenlänge auf die 1,5-fache Länge der Box beschränkt. Werte die außerhalb dieses Bereichs liegen (also mehr als das 1,5-fache der Boxlänge von Q3 entfernt sind) werden mittels eines Punktes dargestellt.
5. Liegt der Median-Strich in der Mitte der Box, so ist die Verteilung symmetrisch (bezogen auf die inneren 50% der Werte), liegt der Median-Strich nicht in der Mitte der Box, so ist die Verteilung nicht symmetrisch (d.h. sie ist schief). Gleiches gilt für die Antennenlängen: Sind die Antennen gleich lang, so ist der äußere Teil der Verteilung symmetrisch, andernfalls schief.

Beispiel 5.6 In einer vorherigen Analyse haben Sie den Zusammenhang von Abschlusspreis und der Anzahl der Lenkräder untersucht. Jetzt möchten Sie eine sehr ähnliche Fragestellung betrachten: Wie unterscheiden sich die Verkaufspreise je nach Anzahl der beigelegten Lenkräder? Flink erstellen Sie dazu folgendes Diagramm, Abbildung 5.28, links. Es zeigt die Verteilung des Abschlusspreises, aufgebrochen nach Anzahl Lenkräder (by = "wheels). \(\square\)

Aber ganz glücklich sind Sie mit dem Diagramm nicht: R hat die Variable wheels komisch aufgeteilt. Es wäre eigentlich ganz einfach, wenn R die Gruppen 0, 1, 2, 3 und 4 aufteilen würde. Aber schaut man sich die Y-Achse (im linken Teildiagramm von Abbildung 5.28) an, so erkennt man, dass R wheels als stetige Zahl betrachtet und nicht in ganze Zahlen gruppiert.²⁶ Aber wir möchten jeden einzelnen Wert von wheels (0, 1, 2, 3, 4) als Gruppe verstehen. Mit anderen Worten, wir möchten wheels als nominale Variable definieren. Das kann man mit dem Befehle factor(wheels) erreichen (verpackt in mutate), s. Abbildung 5.28 rechts.

mariokart_no_extreme %>% 
  select(total_pr, wheels) %>% 
  plot_boxplot(by = "wheels")

mariokart_no_extreme %>% 
  select(total_pr, wheels) %>% 
  mutate(wheels = factor(wheels)) %>% 
  plot_boxplot(by = "wheels")

(a) wheels als metrische Variable

(b) wheels als nominale Variable

Abbildung 5.28: Abschlusspreis nach Anzahl von beigelegten Lenkrädern

Sie schließen aus dem Bild, dass Lenkräder und Preis (positiv) zusammenhängen. Allerdings scheint es wenig Daten für wheels == 4 zu geben. Das prüfen Sie nach:

mariokart_no_extreme %>% 
  count(wheels)

Tatsächlich gibt es (in mariokart_no_extreme) auch für 3 Lenkräder schon wenig Daten, so dass wir die Belastbarkeit dieses Ergebnisses skeptisch betrachten sollten.

Übrigens bezeichnet Sie Ihre Chefin nur noch als “Datengott”.

Beispiel 5.7  

Aufgabe

Visualisieren Sie den Unterschied im Verkaufspreis zwischen gebrauchten und neuen Spielen.

Lösung

mariokart_no_extreme |> 
  select(cond, total_pr) |> 
  plot_boxplot(by = "cond")

Beispiel 5.8  

Aufgabe

Visualisieren Sie den Unterschied im Verkaufspreis abhängig von ship_pr; betrachten Sie ship_pr als ein Gruppierungsvariable. Interpretieren Sie das Ergebnis.

Lösung

mariokart_no_extreme |> 
  select(ship_pr, total_pr) |> 
  plot_boxplot(by = "ship_pr")

plot_boxplot gruppiert metrische Variablen, wie ship_pr automatisch in fünf Gruppen (mit gleichen Ranges). Wir müssen also nichts tun, um die metrische Variable ship_pr in eine Gruppierungsvariable (Faktorvariable) umzuwandeln.

Es sieht so aus, als würde der Median zwischen den Gruppen leicht steigen, mit Ausnahme der mittleren Gruppe.

5.7 So lügt man mit Statistik

Diagramme werden häufig eingesetzt, um die Wahrheit “aufzuhübschen”.

5.7.1 Achsen manipulieren

Achsen zu stauchen ist ein einfacher Trick, s. Abbildung 5.29.

(a) Oh nein, dramatischer Einbruch des Umsatzes!

(b) Kaum der Rede wert, ist nur ein bisschen Schwankung!

Abbildung 5.29: Stauchen der Y-Achse, um mit Statistik zu lügen

Natürlich kann man auch durch “Abschneiden” der Y-Achse einen eindrucksvollen Effekt erzielen, s. Abbildung 5.30.

(a) Oh nein, dramatischer Einbruch des Umsatzes!

(b) Kaum der Rede wert, ist nur ein bisschen Schwankung!

Abbildung 5.30: Abschneiden der Y-Achse, um mit Statistik zu lügen

5.7.2 Scheinkorrelation

Messerli (2012) berichtet von einem Zusammenhang von Schokoladenkonsum und Anzahl von Nobelpreisen (Beobachtungseinheit: Länder), s. Abbildung 5.31. Das ist doch ganz klar: Schoki futtern macht schlau und Nobelpreise! (?)

Abbildung 5.31: Schokolodenkonsum und Nobelpreise

Leider ist hier von einer Scheinkorrelation auszugehen: Auch wenn die beiden Variablen Schokoladenkonsum und Nobelpreise zusammenhängen, heißt das nicht, dass die Variable die Ursache und die andere die Wirkung sein muss. So könnte auch eine Drittvariable im Hintergrund die gleichzeitige Ursache von Schokoladenkonsum und Nobelpreise sein, etwa der allgemeine Entwicklungsstand des Landes: In höher entwickelten Ländern wird mehr Schokolade konsumiert und es werden mehr Nobelpreise gewonnen im Vergleich zu Ländern mit geringerem Entwicklungsstand.

5.8 Praxisbezug

Ein, wie ich finde schlagendes Beispiel zur Stärke von Datendiagrammen ist Abbildung 5.32. Das Diagramm zeigt die Häufigkeit von Masern, vor und nach der Einführung der Impfung. Die Daten und die Idee zur Visualisierung gehen auf van Panhuis et al. (2013) zurück. Das Diagramm und weitere finden sich in ähnlicher Form im Wall Street Journal.

Abbildung 5.32: Häufigkeit von Masern und Impfung in den USA, Lizenz: MIT

Quellcode²⁷

In der “freien Wildbahn” findet man häufig sog. “Tortendiagramme”. Zwar sind sie beliebt, doch ist von ihrer Verwendung zumeist abzuraten; vgl. auch hier.²⁸

5.9 Vertiefung

Mehr Informationen zu {DataExplorer} finden Sie hier.

5.9.1 Animation

Eine weitere nützliche Art von Visualisierung sind Karten, 3D-Bilder und Animationen. So zeigt z.B. Abbildung 5.33 die Veränderung der Lebenserwartung (in Jahren) über die letzten Dekaden.²⁹

Abbildung 5.33: Animation zur Veränderung der Lebenserwartung

In einigen Situation können Animationen zweckdienlich sein. Außerdem sind sie mitunter nett anzuschauen, s. Abbildung 5.34.

Abbildung 5.34: Veränderung des Zusammenhangs von Lebenswertung und Bruttosozialprodukt pro Land, gegliedert nach Kontinenten

Weniger spektakulär als Abbildung 5.4, aber näher an der Datenanalyse ist Abbildung 5.35.

Abbildung 5.35: 3D-Punktediagramm zum Datensatz mariokart

Leider ist Abbildung 5.35 nicht sehr aufschlussreich.

Natürlich sind der Fantasie keine Grenzen beim Visualisieren gesetzt, so ist etwa diese Animationen ziemlich atemberaubend.³⁰

5.9.2 Schicke Diagramme

Ein Teil der Diagramm dieses Kapitels wurden mit dem R-Paket ggpubr erstellt. Mit diesem Paket lassen sich einfach ansprechende Datendiagramme erstellen, so lautet die etwa die Syntax von Abbildung 5.26 wie folgt.

library(ggpubr)  # einmalig instalieren nicht vergessen
mariokart %>% 
  filter(total_pr < 100) %>% 
  ggboxplot(x = "cond", y = "total_pr")

Möchte man Mittelwerte vergleichen, so sind Boxplots nicht ideal, da diese ja nicht den Mittelwert, sondern den Median heraustellen. Eine Abhilfe (also eine Darstellung des Mittelwerts) schafft man (z.B.) mit ggpubr, s. Abbildung 5.36.

ggviolin(mariokart_no_extreme, x = "cond", y = "total_pr",
         add = "mean_sd") 

Abbildung 5.36: Vergleich der Verteilungen zweier Gruppen mit Mittelwert und Standardabweichung pro Gruppe hervorgehoben

Ein “Violinenplot” hat die gleiche Aussage wie ein Dichtediagramm: Je breiter die “Violine”, desto mehr Beobachtungen gibt es an dieser Stelle. Weitere Varianten zum Violinenplot mit ggpubr finden sich hier.

Übrigens sind Modelle - und Diagramme sind Modelle - immer eine Vereinfachung, lassen also Informationen weg. Manchmal auch wichtige Informationen. Dieses Beispiel zeigt, wie etwa Histogramme wichtige Informationen unter den Tisch fallen lassen.³¹

Ein weiteres R-Paket zur Erstellung ansprechender Datenvisualisierung heißt ggstatsplot.³²

Abbildung 5.37 zeigt ein Histogramm, das mit ggstatsplot erstellt wurde.³³

library(ggstatsplot)

gghistostats(
  data       = mariokart_no_extreme,
  x          = total_pr,
  xlab       = "Verkaufspreis" 
  # results.subtitle = FALSE   # unterdrückt statistische Kennzahlen
)

Abbildung 5.37: Ein Histogramm mit ggstatsplot

Die Menge der statistischen Kennzahlen bei ggstatsplot schindet ordentlich Eindruck. Möchte man auf die Kennzahlen verzichten, so nutzt man den Schalter results.subtitle = FALSE.^[Weitere Hinweise finden sich [auf der Hilfeseite der Funktion: https://indrajeetpatil.github.io/ggstatsplot/articles/web_only/gghistostats.html].

🧑‍🎓 Ich würde gerne mal Beispiele von schlechten Datendiagrammen sehen.

👨‍🏫 Auf der Seite von Flowingdata findet sich eine nette Liste mit schlechten Datendiagrammen.³⁴

5.9.3 Farbwahl

Einige Überlegungen zur Farbwahl findet sich bei Wilke (2019), s. Kap. 4.³⁵ Die Farbpalette von Okabe und Ito ist (vgl. Ichihara et al., 2008) empfehlenswert, da sie auch bei Schwarz-Weiß-Druck und bei Sehschwächen die Farben noch recht gut unterscheiden lässt, s. Abbildung 5.38.

mariokart %>% 
  filter(total_pr < 100) %>% 
  ggboxplot(x = "cond", y = "total_pr", fill = "cond") +
  scale_fill_okabeito()

Abbildung 5.38: Die Farbskala von Okabe und Ito: Geeignet bei Farbseh-Schwächen und für Schwarz-Weiß-Druck. Außerdem nett anzuschauen.

5.9.4 Literaturhinweise

Sowohl ggpubr als auch DataExplorer (und viele andere R-Pakete) bauen auf dem R-Paket ggplot2 auf. ggplot2 ist eines der am weitesten ausgearbeiteten Softwarepakete zur Erstellung von Datendiagrammen. Das Buch zur Software (vom Autor von ggplot2) ist empfehlenswert (Wickham, 2009). Eine neue, gute Einführung in Datenvisualisierung findet sich bei Wilke (2019). Beide Bücher sind kostenfrei online lesbar.

5.10 Aufgaben

Die Webseite datenwerk.netlify.app stellt eine Reihe von einschlägigen Übungsaufgaben bereit. Sie können die Suchfunktion der Webseite nutzen, um die Aufgaben mit den folgenden Namen zu suchen:

1. boxhist
2. max-corr1
3. max-corr2
4. Histogramm-in-Boxplot
5. Diamonds-Histogramm-Vergleich2
6. Boxplot-Aussagen
7. boxplots-de1a
8. movies-vis1
9. movies-vis2
10. vis-gapminder
11. boxplots-de1a
12. diamonds-histogramm-vergleich
13. wozu-balkendiagramm
14. diamonds-histogram
15. n-vars-diagram

Noch mehr Aufgaben zum Thema Datenvisualisierung finden Sie im Datenwerk unter dem Tag vis.

5.11 Literaturhinweise

Wilke (2019) gibt einen hervorragenden Überblick über praktische Aspekte der Datenvisualisierung; gut geeignet, wenn man mit R arbeitet. In ähnlicher Richtung geht Fisher & Meyer (2018). Hier ist eine Liste von Büchern zum Thema; dort können Sie bei Interesse tiefer suchen.

5.12 Literatur

Anscombe, F. J. (1973). Graphs in Statistical Analysis. The American Statistician, 27(1), 17–21.

Cohen, J. (1992). A Power Primer. Psychological Bulletin, 112(1), 155–159.

Fisher, D., & Meyer, M. (2018). Making Data Visual: A Practical Guide to Using Visualization for Insight (First edition). O’Reilly.

Ichihara, Y. G., Okabe, M., Iga, K., Tanaka, Y., Musha, K., & Ito, K. (2008). Color Universal Design: The Selection of Four Easily Distinguishable Colors for All Color Vision Types. Color Imaging XIII: Processing, Hardcopy, and Applications, 6807, 206–213. https://doi.org/10.1117/12.765420

Lyon, A. (2014). Why Are Normal Distributions Normal? The British Journal for the Philosophy of Science, 65(3), 621–649. https://doi.org/10.1093/bjps/axs046

Messerli, F. H. (2012). Chocolate Consumption, Cognitive Function, and Nobel Laureates. New England Journal of Medicine, 367(16), 1562–1564. https://doi.org/10.1056/NEJMon1211064

Scherer, C., Radchuk, V., Staubach, C., Müller, S., Blaum, N., Thulke, H., & Kramer‐Schadt, S. (2019). Seasonal Host Life‐history Processes Fuel Disease Dynamics at Different Spatial Scales. Journal of Animal Ecology, 88(11), 1812–1824. https://doi.org/10.1111/1365-2656.13070

van Panhuis, W. G., Grefenstette, J., Jung, S. Y., Chok, N. S., Cross, A., Eng, H., Lee, B. Y., Zadorozhny, V., Brown, S., Cummings, D., & Burke, D. S. (2013). Contagious Diseases in the United States from 1888 to the Present. New England Journal of Medicine, 369(22), 2152–2158. https://doi.org/10.1056/NEJMms1215400

Wickham, H. (2009). Ggplot2: Elegant Graphics for Data Analysis. Springer. https://doi.org/10.1007/978-0-387-98141-3

Wilke, C. (2019). Fundamentals of Data Visualization: A Primer on Making Informative and Compelling Figures (First edition). O’Reilly Media. https://clauswilke.com/dataviz/

---

1. https://antworte.jetzt/

2. Quelle: https://towardsdatascience.com/how-to-turn-a-dinosaur-dataset-into-a-circle-dataset-with-the-same-statistics-64136c2e2ca0

3. Quelle: https://en.wikipedia.org/wiki/Anscombe%27s_quartet

4. 5

5. Quelle: Plotly, https://plotly.com/r/trisurf/

6. Weitere Nomenklaturen sind möglich, aber wir halten hier die Sache einfach. Wer an Vertiefung interessiert ist, findet bei data-to-vis einen Überblick über verschiedene Typen an Diagrammen, sogar in Form einer systematischen Nomenklatur: https://www.data-to-viz.com/.

7. synonym: Säulendiagramm

8. Natürlich müssen Sie das Paket einmalig installiert haben, bevor Sie es starten können.

9. Auf der Hilfeseite der Funktion finden Sie weitere Details zur Funktion.

10. Zur Erinnerung: %>% nennt man die “Pfeife und lässt sich als”und dann” übersetzen, vgl. Kapitel 4.4.

11. bei konstanter Balkenbreite

12. Mit Dichte ist die Anzahl der Beobachtungen pro Einheit der Variablen auf der X-Achse gemeint.

13. Grob gesagt: mariokart %>% plot_density().

14. vgl. Kapitel 6.5

15. vgl. Kapitel 7.4.

16. Quelle: ifes/FOM Hochschule, https://github.com/FOM-ifes/VL-Vorlesungsfolien

17. Quelle: https://observablehq.com/@mcmcclur/the-normal-model

18. Quelle: Ainali - Own work, CC BY-SA 3.0, https://commons.wikimedia.org/w/index.php?curid=3141713

19. Quelle: Aufbauend auf FOM/ifes, Autor: Norman Markgraf

20. synonym: Regressionsgerade

21. Quelle: https://observablehq.com/d/bb7ad3ecfb1ac2a6

22. Aber Freunde lassen Freunde keine Tortendiagramme verwenden: https://github.com/cxli233/FriendsDontLetFriends#10-friends-dont-let-friends-make-pie-chart.

23. Übrigens: Freunde lassen Freunde nicht Balkendiagramme verwenden, um Mittelwerte darzustellen: https://github.com/cxli233/FriendsDontLetFriends#1-friends-dont-let-friends-make-bar-plots-for-means-separation.

24. “substanzieller”, “bedeutsamer”, “relevanter” oder “(inhaltlich) signifikanter”

25. Ob der Boxplot horizontal oder vertikal steht, ist Ihrem Geschmack überlassen.

26. Vielleicht so, dass in jeder Gruppe gleich viele Wert sind?

27. Datenquelle: https://www.tycho.pitt.edu

28. https://www.data-to-viz.com/caveat/pie.html; https://github.com/cxli233/FriendsDontLetFriends#10-friends-dont-let-friends-make-pie-chart

29. Der Quellcode der Animation ist hier zu finden: https://gist.github.com/rafapereirabr/0d68f7ccfc3af1680c4c8353cf9ab345.

30. https://www.tylermw.com/wp-content/uploads/2019/06/featuredmeasles.mp4

31. https://www.autodesk.com/research/publications/same-stats-different-graphs

32. https://github.com/IndrajeetPatil/ggstatsplot/blob/main/README.md

33. https://github.com/IndrajeetPatil/ggstatsplot/blob/main/README.md#gghistostats

34. https://flowingdata.com/category/visualization/ugly-visualization/

35. Siehe auch: https://data-se.netlify.app/2023/06/30/farbpaletten/