LoslegenKostenlos loslegen

Hypothesentest: Hat sich die Erdbebenhäufigkeit verändert?

Offensichtlich gab es einen massiven Anstieg der Erdbebenhäufigkeit, sobald die Einleitung von Abwasser begann. Zur Übung machst du trotzdem einen Hypothesentest. Du testest nicht die Hypothese, dass die Zeiten zwischen Erdbeben vor und nach 2010 dieselbe Verteilung haben, da die Abwassereinleitung die Verteilung beeinflussen kann. Stattdessen nimmst du an, dass sie denselben Mittelwert haben. Berechne also den p-Wert zur Hypothese, dass die Zeiten zwischen Erdbeben vor und nach 2010 denselben Mittelwert haben, und verwende dabei als Teststatistik den Mittelwert der Zeitabstände vor 2010 minus den Mittelwert der Zeitabstände nach 2010.

Diese Übung ist Teil des Kurses

Fallstudien zum statistischen Denken

Kurs anzeigen

Anleitung zur Übung

  • Berechne die beobachtete Teststatistik. Die Variablen mean_dt_pre und mean_dt_post aus den vorherigen Übungen sind in deinem Namespace vorhanden.
  • Verschiebe die Daten nach 2010 so, dass sie denselben Mittelwert wie die Daten vor 2010 haben. Speichere das Ergebnis als dt_post_shift.
  • Ziehe jeweils 10.000 Bootstrap-Replikate des Mittelwerts von dt_pre und dt_post_shift. Speichere die jeweiligen Ergebnisse in bs_reps_pre und bs_reps_post.
  • Berechne Replikate der Mittelwertsdifferenz, indem du bs_reps_post von bs_reps_pre subtrahierst.
  • Berechne und gib den p-Wert aus. Interpretiere „mindestens genauso extrem wie“ so, dass die Teststatistik größer oder gleich dem beobachteten Wert ist.

Interaktive Übung

Vervollständige den Beispielcode, um diese Übung erfolgreich abzuschließen.

# Compute the observed test statistic
mean_dt_diff = ____ - ____

# Shift the post-2010 data to have the same mean as the pre-2010 data
dt_post_shift = ____ - ____ + ____

# Compute 10,000 bootstrap replicates from arrays
bs_reps_pre = ____
bs_reps_post = ____

# Get replicates of difference of means
bs_reps = ____ - ____

# Compute and print the p-value
p_val = ____(____ >= ____) / 10000
print('p =', p_val)
Code bearbeiten und ausführen