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
Fallstudien zum statistischen Denken
Anleitung zur Übung
- 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 Übung
Vervollständige den Beispielcode, um diese Übung erfolgreich abzuschließen.
# 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)