Schattingen van de gemiddelde tijd tussen aardbevingen
De grafische EDA in de vorige oefening laat een duidelijke verandering in de frequentie van aardbevingen rond 2010 zien. Ter vergelijking: bereken de gemiddelde tijd tussen aardbevingen van magnitude 3 en groter van 1980 t/m 2009 en ook van 2010 t/m medio 2017. Neem ook 95%-betrouwbaarheidsintervallen van het gemiddelde op. De variabelen dt_pre en dt_post bevatten respectievelijk de tijdsintervallen tussen alle aardbevingen met een magnitude van minstens 3 van vóór 2010 en na 2010, in dagen.
Deze oefening maakt deel uit van de cursus
Casestudies in statistisch denken
Oefeninstructies
- Bereken de gemiddelde tijd tussen aardbevingen voor vóór (
dt_pre) en na 2010 (dt_post). - Trek 10.000 bootstrap-replicaties van het gemiddelde voor de datasets van vóór en na 2010.
- Gebruik
np.percentile()om het 95%-betrouwbaarheidsinterval van het gemiddelde voor beide datasets te berekenen. - Klik op 'Antwoord verzenden' om de resultaten op het scherm te printen.
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# Compute mean interearthquake time
mean_dt_pre = ____
mean_dt_post = ____
# Draw 10,000 bootstrap replicates of the mean
bs_reps_pre = ____
bs_reps_post = ____
# Compute the confidence interval
conf_int_pre = ____
conf_int_post = ____
# Print the results
print("""1980 through 2009
mean time gap: {0:.2f} days
95% conf int: [{1:.2f}, {2:.2f}] days""".format(mean_dt_pre, *conf_int_pre))
print("""
2010 through mid-2017
mean time gap: {0:.2f} days
95% conf int: [{1:.2f}, {2:.2f}] days""".format(mean_dt_post, *conf_int_post))