Beräkning av b-värden
Utifrån det diagram du skapade i föregående övning kan du tryggt använda en fullständighetsgräns på mt = 3. Beräkna b-värden för perioden 1980–2009 samt för 2010 till mitten av 2017 med detta tröskelvärde. Funktionen b_value() som du skrev i förra kapitlet – och som beräknar b-värdet och konfidensintervallet från en uppsättning magnituder och en fullständighetsgräns – finns tillgänglig i ditt namnrum, liksom numpy-arrayerna mags_pre och mags_post från föregående övning samt mt.
Den här övningen är en del av kursen
Fallstudier i statistiskt tänkande
Övningsinstruktioner
- Beräkna b-värdet och 95-procentigt konfidensintervall för jordbävningar från 1980 till och med 2009 med hjälp av 10 000 bootstrap-replikat.
- Beräkna b-värdet och 95-procentigt konfidensintervall för jordbävningar från 2010 till mitten av 2017 med hjälp av 10 000 bootstrap-replikat.
- Klicka på 'Skicka in svar' för att skriva ut resultaten på skärmen.
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
# 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))