Kuantifikasi nilai b
Berdasarkan plot yang Anda buat pada latihan sebelumnya, Anda dapat dengan aman menggunakan ambang kelengkapan mt = 3. Dengan ambang ini, hitung nilai b untuk periode antara 1980 hingga 2009 dan untuk 2010 hingga pertengahan 2017. Fungsi b_value() yang Anda tulis pada bab sebelumnya—yang menghitung nilai b dan interval kepercayaan dari sekumpulan magnitudo dan ambang kelengkapan—tersedia dalam namespace Anda, begitu pula array numpy mags_pre dan mags_post dari latihan terakhir, serta mt.
Latihan ini adalah bagian dari kursus
Studi Kasus dalam Pemikiran Statistik
Petunjuk latihan
- Hitung nilai b dan interval kepercayaan 95% untuk gempa bumi dari 1980 hingga 2009 menggunakan 10.000 replikasi bootstrap.
- Hitung nilai b dan interval kepercayaan 95% untuk gempa bumi dari 2010 hingga pertengahan 2017 menggunakan 10.000 replikasi bootstrap.
- Tekan 'Kirim Jawaban' untuk mencetak hasil ke layar.
Latihan interaktif praktis
Cobalah latihan ini dengan menyelesaikan kode contoh berikut.
# 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))