Kwantisering van de b-waarden
Op basis van de plot die je in de vorige oefening hebt gemaakt, kun je veilig een volledigheidsdrempel van mt = 3 gebruiken. Gebruik deze drempel om b-waarden te berekenen voor de periode 1980–2009 en voor 2010 tot en met medio 2017. De functie b_value() die je in het vorige hoofdstuk hebt geschreven, waarmee de b-waarde en het betrouwbaarheidsinterval worden berekend op basis van een set magnitudes en een volledigheidsdrempel, is beschikbaar in je namespace. Dat geldt ook voor de numpy-arrays mags_pre en mags_post uit de vorige oefening, en mt.
Deze oefening maakt deel uit van de cursus
Casestudies in statistisch denken
Oefeninstructies
- Bereken de b-waarde en het 95%-betrouwbaarheidsinterval voor aardbevingen van 1980 tot en met 2009 met 10.000 bootstrap-replicaties.
- Bereken de b-waarde en het 95%-betrouwbaarheidsinterval voor aardbevingen van 2010 tot en met medio 2017 met 10.000 bootstrap-replicaties.
- Klik op "Antwoord verzenden" om de resultaten op het scherm af te drukken.
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# Compute b-value and confidence interval for pre-2010
b_pre, conf_int_pre = ____(____, ____, ____=[____, ____], ____=____)
# Compute b-value and confidence interval for post-2010
b_post, conf_int_post = ____
# Report the results
print("""
1980 through 2009
b-value: {0:.2f}
95% conf int: [{1:.2f}, {2:.2f}]
2010 through mid-2017
b-value: {3:.2f}
95% conf int: [{4:.2f}, {5:.2f}]
""".format(b_pre, *conf_int_pre, b_post, *conf_int_post))