Technischer Hinweis: Links in den Übungen öffnen Sie, indem sie darauf rechtsklicken und “Link in neuem Tab/Fenster öffnen” wählen.
Zum Aufwärmen ein paar einfache Übungen zu Wahrscheinlichkeiten
Sie werden in jedem Fall eine Zahl zwischen 1 und 6 werfen, etwas
anderes gibt es nicht. Die Gesamtwahrscheinlichkeit addiert sich
zu
\[P = 1/6 (für 1) + 1/6 (für 2) + 1/6 (für 3)
+ 1/6 (für 4) + 1/6 (für 5) + 1/6 (für 6) = 1\]
Die Wahrscheinlichkeit, keine 6 zu würfeln beträgt also
\[P(\neq 6) = 1/6 (für 1) + 1/6 (für 2) + 1/6 (für 3) + 1/6 (für 4) + 1/6 (für 5) = 5/6 \approx 83\%\]
In den folgenden Übungen arbeiten wir mit dem bereits bekannten
Datensatz physio.csv, der Daten von 228 Studierenden der
Physiotherapie enthält. Wenn Sie den Datensatz bei den Übungen zur
deskriptiven Statistik als physio.omv gespeichert haben,
sind die Variablen bereits kategorisiert und Sie arbeiten mit diesem
Datensatz, andernfalls laden Sie physio.csvin
jamovi und kategorisieren als erstes die Variablen wie
in den Übungen zur deskriptiven Statistik beschrieben.
Als erstes verschaffen wir uns einen Überblick zu den Kennzahlen der Körpergrösse der Studentinnen und stellen diese in einem Histogramm dar.
| n | m | median | s | min | max |
|---|---|---|---|---|---|
| 183 | 166.92 | 167 | 5.66 | 148 | 183 |
Anleitung für jamovi:
Wir arbeiten mit einem Teildatensatz für Studentinnen. Dazu müssen
wir unsere Daten Filtern. Das geht folgendermassen:
jamovi > Register Data > Menüband
Filters > im Fenster geben sie = Geschlecht ==
“w” ein und schliessen mit der Enter/Return-Taste ab. Im
Tabellenfenster links ist eine neue Spalte mit der Bezeichnung
Filter 1 erschienen. Bei den Frauen
(Geschlecht == “w”) sollte in dieser Spalte ein grünes
Häckchen und bei den Männern (Geschlecht== “m”) ein rotes X
eingetragen sein.
Hilfe zum Filtern von Daten bieten diese Video-Tutorials
Wenn der Filter gesetzt ist, können Sie die Kennzahlen und das Histogramm (ohne Normalverteilungskurve) in jamovi erstellen, wie wir das bei den Übungen zur deskriptiven Statistik gelernt haben.
Die Standardnormalverteilung ist gekennzeichnet durch den Mittelwert 0 und die Standardabweichung 1. Die Umwandlung einer beliebigen Verteilung erfolgt nach der Formel:
\[z_n=\frac{x_n-\bar{x}}{s}\]
Beispiel: (Mittelwert und Standardabweichung siehe Tabelle oben)
\[z_n=\frac{x_n-166.92}{5.66}\]
Wir sehen, dass das standardisierte Histogramm und die standardisierte Normalverteilungskurve genau gleich aussehen, wie mit den nicht transformierten Originaldaten. Was bringt uns das also?
Wenn wir beliebige, normalverteilte Merkmale messen, werden wir immer wieder neue Normalverteilungen mit ihrem eigenen Mittelwert und ihrer eigenen Standardabweichung erhalten. D.h. die Normalverteilungskurve wird manchmal breiter, manchmal schmaler, manchmal höher, manchmal tiefer sein. Jedes beliebige Merkmal weist seine charakteristische Verteilung auf. In der Statistik interessieren wir uns oft für die Wahrscheinlichkeit, dass ein bestimmtes Ereignis eintritt (z.B. das Ereignis, dass eine Studentin 168 cm gross ist), und effektiv ist die Normalverteilung eine Wahrscheinlichkeitsverteilung. Die Berechnung von Wahrscheinlichkeiten erfolgt über die Berechnung von Flächeninhalten unter der Normalverteilungskurve wobei die gesamte Fläche unter der Kurve stets gleich 1 ist (Tönt kompliziert, ist es aber nicht. Wir werden das unten üben). Die Berechnung von Flächen unter Kurven ist tatsächlich eine nicht ganz triviale Angelegenheit. Um diese Berechnungen zu erleichtern, wurde die Standardnormalverteilung “erfunden”. Die Standardisierung erfolgt über die Berechnung von sog. z-Werten. Die z-Werte-Tabelle ist ein Hilfsmittel, das es uns erlaubt, für beliebige Normalverteilungen beliebige Flächeninhalte und somit Wahrscheinlichkeiten zu berechnen.
Uns interessiert jetzt die Frage Wie gross ist die Wahrscheinlichkeit, dass eine Studentin eine bestimmte Körpergrösse aufweist?. Aus hier nicht näher erläuterten Gründen ist es nicht möglich, diese Wahrscheinlichkeit exakt zu berechnen (sie wäre gleich Null). Was jedoch möglich ist, ist die Beantwortung der Frage Wie gross ist die Wahrscheinlichkeit dass eine Studentin gleich gross oder grösser, bzw. kleiner ist?.
Beispiel: Wie gross ist die Wahrscheinlichkeit, dass eine Studientin 163 cm gross oder kleiner ist?
\[z_{163} = \frac{163 - 166.92}{5.66} = -0.693\]
Die Wahrscheinlichkeit, dass eine Studentin 163 cm gross oder kleiner ist, entspricht der Grösse der Fläche unter der Kurve links von der violetten Linie. Wir können schon von blossem Auge sagen, dass die Wahrscheinlichkeit kleiner als 50% sein muss, da die violette Linie bei 163 cm die Fläche in eine kleinere linke und eine grössere rechte Fläche teilt.
Aber wie genau berechnet man jetzt die Grösse der blauen Fläche? Bevor Computerprogramme die Berechnung übernommen haben, hat man sich mit sogenannten z-Werte-Tabellen (siehe https://rpubs.com/BFH-PHY-WA/1007683) beholfen. Weil es unmöglich ist für alle möglichen Verteilungen Tabellen zu berechnen, wurde die z-Transformation entwickelt. Damit lassen sich beliebige Normalverteilungen in eine standardisierte Verteilung umwandeln. Heute übernimmt diese Arbeit der PC, wir machen es hier aber von Hand, um zu verstehen, was da vorgeht.
Übung: Welche Wahrscheinlichkeit gibt die Z-Wertetabelle für den z-Wert -0.6926956 an?
Der Computer berechnet die Wahrscheinlichkeit, dass eine Studentin in PHY13-PHY17 163 cm oder kleiner ist mit p = 0.2442503.
Übung: Wie gross ist die Wahrscheinlichkeit, dass eine Studentin 163 cm oder grösser ist?
Dies ist die Umkehrung der Frage. Da die Wahrscheinlichkeit dafür, dass die Studentin 163 cm oder kleiner ist jetzt bekannt ist, lässt sich die Frage ganz einfach beantworten, da bekanntlicherweise die Fläche unter der Kurve gleich 1 ist.
Lösung: Die Wahrscheinlichkeit, dass eine Studentin 163 cm oder grösser ist beträgt
\[1-(Wahrscheinlichkeit für \leq 163 cm) = 1-0.2403793\] also p = 0.7557497.
Ziel: \(z\)-Werte für Körpergrössen von Studentinnen berechnen
Tipp: Berechnen Sie die \(z\)-Werte von Hand mit einem Taschenrechner.
\(z = \frac{176 - 166.92}{5.66} = 1.660424 \approx 1.66\)
\(z = \frac{166.92 - 166.92}{5.66} = 0\)
\(z = \frac{195 - 166.92}{5.66} = 4.961131 \approx 4.96\)
Ziel: Wahrscheinlichkeiten für Körpergrössen von Studentinnen berechnen.
Erstellen Sie jeweils eine Skizze für jede Frage und überlegen Sie, ob der z-Wert jeweils links oder rechts vom Mittelwert liegt und ob p kleiner oder grösser 50% ist.
Vorgehen: Sie können die Wahrscheinlichkeiten für ihre berechneten z-Werte entweder mit der z-Werte-Tabelle bestimmen oder mit jamovi. Wir beschreiben hier das Vorgehen für jamovi.
pnorm().pnorm() die Wahrscheinlichkeit für einen
bestimmten z-Wert berechnen kann, benötigt die Funktion die Kennzahlen
der Verteilung Mittelwert und Standardabweichung,
daraus ergibt sich die Syntax von pnorm()\[pnorm(z-Wert, mean = Mittelwert,~ sd = Standardabweichung, ~ lower.tail = FALSE/TRUE)\]
lower.tail gibt an, ob die Fläche unter
der Kurve links vom z-Wert berechnet werden soll (p für \(\leq\) z-Wert, lower.tail = TRUE), oder die
Fläche unter der Kurve rechts vom Wert (p für \(\geq\) Wert, lower.tail = FALSE).lower.tail = TRUE eingestellt. Statt
lower.tail = FALSE kann die Fläche rechts vom z-Wert auch
als 1 - pnorm() berechnet werden.Codebeispiel:
pnorm(z-Wert, mean = 0, sd = 1, lower.tail = FALSE) # Fläche links vom z-Wert
1 - pnorm(z-Wert, mean = 0, sd = 1, lower.tail = FALSE) # Fläche rechts vom z-WertZur Erinnerung: Mittelwert und Standardabweichung für die Grösse der Studentinnen sind
| Mittelwert | Standardabweichung |
|---|---|
| 166.92 | 5.66 |
Skizze:
pnorm(176, mean = 166.92, sd = 5.66, lower.tail = TRUE)## [1] 0.9456695
Die Wahrscheinlichkeit, dass eine Studentin 176 cm oder kleiner ist beträgt p = 0.946 (94.6%). Mit anderen Worten, 94.6% der Studentinnen sind 176 cm gross oder kleiner.
Skizze
pnorm(176, mean = 166.92, sd = 5.66, lower.tail = FALSE)## [1] 0.05433055
1 - pnorm(176, mean = 166.92, sd = 5.66) # ergibt das gleiche Resultat## [1] 0.05433055
Die Wahrscheinlichkeit, dass eine Studentin 176 cm oder grösser ist beträgt p = 0.054 (5.4%). Mit anderen Worten, 5.4% der Studentinnen sind 176 cm gross oder grösser.
pnorm(166.92, mean = 166.92, sd = 5.66, lower.tail = TRUE)## [1] 0.5
Die Wahrscheinlichkeit, dass eine Studentin 166.92 cm oder kleiner ist beträgt p = 0.5 (50%). Mit anderen Worten, 50% der Studentinnen sind 176 cm gross oder kleiner.
pnorm(166.92, mean = 166.92, sd = 5.66, lower.tail = FALSE) ## [1] 0.5
1 - pnorm(166.92, mean = 166.92, sd = 5.66) # ergibt das gleiche Resultat## [1] 0.5
Die Wahrscheinlichkeit, dass eine Studentin 166.924 cm oder grösser ist beträgt p = 0.5 (50%). Mit anderen Worten, 50% der Studentinnen sind 166.924 cm gross oder grösser.
pnorm(195, mean = 166.92, sd = 5.66, lower.tail = FALSE)## [1] 3.5042e-07
Die Wahrscheinlichkeit, dass eine Studentin 195 cm oder grösser ist beträgt p = 0.00000035 (0.000035%). Mit anderen Worten, die Wahrscheinlichkeit, dass eine Studentin grösser als 195 cm ist, ist extrem klein.
Lösungsweg: Berechnen Sie die Wahrscheinlichkeit, dass eine Studentin kleiner als 176 cm ist (grüne Fläche) und ziehen Sie von diesem Wert die Wahrscheinlichkeit ab, dass eine Studentin kleiner als 163 cm gross ist (rote Fläche).
pnorm(176, mean = 166.92, sd = 5.66, lower.tail = TRUE) - pnorm(163, mean = 166.92, sd = 5.66, lower.tail = TRUE)## [1] 0.7013827
Die Wahrscheinlichkeit, dass eine Studentin grösser als 163 cm und kleiner als 176 cm ist ist p = 0.701 (70.1%), d.h. ca. 70% der Studentinnen sind zwischen 163 cm und 176 cm gross.
Wie gross ist eine “normalgrosse” Studentin in PHY13 - PHY17 (Tipp: Gesucht ist der 95%-Normbereich)?
Wir erinnern uns, dass unter “normal” die zentralen 95% der Merkmalsausprägungen definiert sind. Die zentralen 95% bedeutet, die Fläche unter der Normalverteilungskurve zwischen Mittelwert minus 2 (eigentlich 1.96) Standardabweichungen und Mittelwert plus 2 (eigentlich 1.96) Standardabweichungen.D.h.
\(Normbereich_{untere Grenze} = \bar{x}-2s\)
\(Normbereich_{obere Grenze} = \bar{x}+2s\)
Daraus folgt: “Normalgrosse” Studentinnen sind zwischen
\(Normbereich_{untere Grenze} = 166.92 - 2 \times 5.66 = 166.92 - 11.34 = 155.62\)
\(Normbereich_{obere Grenze} = 166.92 + 2 \times 5.66 = 166.92 + 11.34 = 178.26\)
cm gross.
Wie gross ist die Wahrscheinlichkeit, dass die Körpergrösse einer Studentin ausserhalb (also oberhalb oder unterhalb) des 95%-Normbereichs liegt?
Die beiden 2.5%-Flächen in der Abbildung zu Aufgabe 3 liegen ausserhalb des 95%-Normbereichs: 100% - 95% = 5%.
\[p = 0.025 + 0.025 = 0.05 = 5\%\]
Wie gross muss eine Studentin mindestens sein, damit sie zu den 10% grössten Studentinnen gehört?
qnorm() in Rj-Editor bestimmen (s. unten)z-Werte-Tabelle
\[ x_i =\bar{x} + z \times s = 166.92 + 1.28 \times 5.66 = 174.16 \]
## z-Wert für die 90. Perzentile in Rj-Editor bestimmen
qnorm(.9)
## Groesse auf der 90. Perzentile aus einer Normalverteilung bestimmen
qnorm(.9, mean = mean(data$Groesse), sd = sd(data$Groesse))## [1] 1.281552
## [1] 174.1823
Eine Studentin muss mindestens 174.18 cm (jamovi gibt 174.0 cm an) gross sein, damit sie zu den 10% der grössten Studentinnen gehört (es entsteht eine kleine Rundungsdifferenz gegenüber dem von Hand berechneten Wert).
qnorm() in Rj-Editor bestimmen (s.
unten)## z-Wert für die 25. Perzentile bestimmen in Rj-Editor
qnorm(.25)
## Groesse auf der 90. Perzentile aus einer Normalverteilung bestimmen
qnorm(.25, mean = mean(data$Groesse), sd = sd(data$Groesse))## [1] -0.6744898
## [1] 163.1031
Eine Studentin, die auf der 25. Perzentile liegt, ist 163.10 cm gross (es entsteht eine kleine Rundungsdifferenz gegenüber dem von Hand berechneten Wert, jamovi gibt 163.0 cm aus).
Wir betrachten die Standardnormalverteilung mit dem Mittelwert \(\mu=0\) und der Standardabweichung \(\sigma= 1\).
Wie gross ist die Wahrscheinlichkeit für für ein Ereignis mit \(z > 2.30\)?
Wie gross ist die Wahrscheinlichkeit, dass \(z < 3.10\) ?
Wie gross ist die Wahrscheinlichkeit dass \(z\) zwischen -1.60 und 3.10 liegt??
Welcher Wert von \(z\) markiert die Grenze zu den oberen 15% der Verteilung?
Welcher Wert von \(z\) markiert die Grenze zu den unteren 20% der Verteilung?
Tipp: Arbeiten Sie mit den R-Funktionen
pnorm() und qnorm im Rj-Editor von
jamovi.
Die Funktion pnorm() gibt die Wahrscheinlichkeit
links von \(z\) an. Die
Wahrscheinlichkeit rechts von \(z\)
berechnen wir als 1 - pnorm().
1 - pnorm(2.3)## [1] 0.01072411
# auf 3 Nachkommastellen gerundet
round(1 - pnorm(2.3), 3)## [1] 0.011
pnorm(3.1)## [1] 0.9990324
# auf 3 Nachkommastellen gerundet
round(pnorm(3.1), 3)## [1] 0.999
pnorm(3.1) - pnorm(-1.6)## [1] 0.9442331
# auf 3 Nachkommastellen gerundet
round(pnorm(3.1) - pnorm(-1.6), 3)## [1] 0.944
Die Grenze zu den oberen 15% der Verteilung ist die 85.
Perzentile. Mit der Funktion qnorm() berechnen wir die
Werte, die auf einer bestimmten Perzentile liegen.
qnorm(.85)## [1] 1.036433
#auf 3 Nachkommastellen gerundet
round(qnorm(.85), 3)## [1] 1.036
qnorm(0.2)## [1] -0.8416212
# auf 3. Nachkommastelle runden
round(qnorm(.2), 3)## [1] -0.842
Die WHO (Weltgesundheitsorganisation) definiert Osteoporose bei jungen Erwachsenen als gemessene Knochendichte, die 2.5 Standardabweichungen unter dem Mittelwert für junge Erwachsene liegt. Wie hoch ist der prozentuale Anteil an jungen Erwachsenen mit Osteoporose, unter der Voraussetzung, dass die Knochendichte bei jungen Erwachsenen normalverteilt ist?
# Berechne die Wahrscheinlichkeit für z = -2.5
pnorm(-2.5)## [1] 0.006209665
# in Prozent
100 * pnorm(-2.5)## [1] 0.6209665
# auf 2. Nachkommastelle gerundet
round(100 * pnorm(-2.5), 2)## [1] 0.62
Gemäss den WHO-Kriterien leiden 0.62% der jungen Erwachsenen an Osteoporose
Ein hoher Cholesterinspiegel im Blut ist ein Risikofaktor für eine koronare Herzkrankheit. Junge Frauen leiden seltener an einem hohen Cholesterinspiegel als andere Gruppen. Der Cholesterinspiegel von Frauen im Alter von 20 bis 34 Jahre ist annähernd normalverteilt mit einem Mittelwert von 185 mg/dl und einer Standardabweichung von 39 mg/dl Blut.
Cholesterinspiegel über 240 mg/dl müssen mediznisch überwacht werden. Wie hoch ist der prozentuale Anteil an jungen Frauen, die einen Cholesterinspiegel von 240 mg/dl oder höher haben?
Cholesterinspiegel über 200 mg/dl gelten als grenzwertig erhöht. Wie gross ist der prozentuale Anteil an jungen Frauen mit einem Cholesterinspiegel zwischen 200 mg/dl und 240 mg/dl?
Variante 1: Berechnung des \(z\)-Werts
Skizze erstellen
# z-Wert für 240 berechnen, m = 185, s = 39
z <- (240 - 185) / 39
# Wahrscheinlichkeit rechts vom z-Wert berechnen (gesucht ist >= 240)
1 - pnorm(z)## [1] 0.07923199
pnorm() geht in ihrer
Grundeinstellung von einer Standardnormalverteilung aus. Wenn eine
andere Normalverteilung als Referenz dient, können ihre Parameter in die
Funktion eingetragen werden
pnorm(Wert, mean = Mittelwert, sd = Standardabweichung)1 - pnorm(240, mean = 185, sd = 39)## [1] 0.07923199
Der Anteil junger Frauen mit einem Cholesterinspiegel über 240 mg/dl beträgt durchschnittlich 7.92%.
# Häufigkeit für <= 240 berechnen (aus Teilaufgabe a)
p_240 <- pnorm(240, mean = 185, sd = 39)
# Häufigkeit für <= 200 berechnen
p_200 <- pnorm(200, mean = 185, sd = 39)
# Differenz berechnen
p_240 - p_200## [1] 0.2710292
Der prozentuale Anteil junger Frauen mit einem Cholesterinspiegel zwischen 200 mg/dl und 240 mg/dl beträgt 27.10%.