Hypothesentest: Sind die b-Werte unterschiedlich?
Führe den im vorherigen Aufgabe skizzierten Hypothesentest durch. Die Variablen mags_pre und mags_post sind bereits in deinem Namespace geladen, ebenso mt = 3.
Diese Übung ist Teil des Kurses
<Kurs>Fallstudien zum statistischen Denken</Kurs>Übungsanweisungen
- Schneide die Magnituden der Erdbeben vor 2010 heraus, die eine Magnitude über (oder gleich) der Vollständigkeitsschwelle haben, und überschreibe
mags_premit dem Ergebnis. Mache dasselbe fürmags_post. - Berechne die beobachtete Differenz der mittleren Magnituden, indem du die Magnituden der Erdbeben vor 2010 von denen nach 2010 subtrahierst.
- Erzeuge 10.000 Permutationsreplikate mit
dcst.draw_perm_reps(). Verwendedcst.diff_of_meansals Argument fürfunc. - Berechne und gib den p-Wert aus, wobei „mindestens so extrem wie“ hier bedeutet, dass die Teststatistik kleiner ist als der beobachtete Wert.
Interaktive praktische Übung
Versuche dich an dieser Übung, indem du diesen Beispielcode vervollständigst.
# Only magnitudes above completeness threshold
mags_pre = ____[____ >= ____]
mags_post = ____[____ >= ____]
# Observed difference in mean magnitudes: diff_obs
diff_obs = ____ - ____
# Generate permutation replicates: perm_reps
perm_reps = ____
# Compute and print p-value
p_val = ____(____ < ____) / 10000
print('p =', p_val)