MulaiMulai sekarang secara gratis

Nilai b untuk Parkfield

ECDF efektif untuk memperlihatkan penurunan (roll-off), seperti yang dapat Anda lihat di bawah magnitudo 1. Karena terdapat banyak gempa di atas magnitudo 3, Anda dapat menggunakan mt = 3 sebagai ambang kelengkapan (completeness threshold). Dengan ambang kelengkapan ini, hitung nilai b untuk wilayah Parkfield dari 1950 hingga 2016, beserta interval kepercayaan 95%. Cetak hasilnya ke layar. Variabel mags yang berisi semua magnitudo sudah tersedia dalam namespace Anda.

Tumpangtindihkan CDF Eksponensial teoretis untuk memverifikasi bahwa wilayah Parkfield mengikuti Hukum Gutenberg–Richter.

Latihan ini adalah bagian dari kursus

Studi Kasus dalam Pemikiran Statistik

Lihat Kursus

Petunjuk latihan

  • Hitung nilai b dan interval kepercayaan 95% menggunakan fungsi b_value() Anda. Gunakan 10.000 replikasi bootstrap.
  • Gunakan np.random.exponential() untuk mengambil 100.000 sampel dari distribusi teoretis. Rata-rata untuk distribusi tersebut adalah b/np.log(10), dan Anda perlu menambahkan mt ke sampel Anda untuk menangani parameter lokasi secara tepat. Simpan hasilnya dalam m_theor.
  • Plot ECDF dari m_theor sebagai garis.
  • Plot ECDF dari semua magnitudo di atas mt sebagai titik.
  • Tekan 'Submit Answer' untuk menampilkan plot dan mencetak nilai b serta interval kepercayaan ke layar.

Latihan interaktif praktis

Cobalah latihan ini dengan menyelesaikan kode contoh berikut.

# Compute b-value and 95% confidence interval
b, conf_int = ____(____, ____, ____=[____, ____], n_reps=____)

# Generate samples to for theoretical ECDF
m_theor = ____(____, size=____) + ____

# Plot the theoretical CDF
_ = ____(*____)

# Plot the ECDF (slicing mags >= mt)
_ = plt.plot(*____(____[____ >= ____]), marker='.', linestyle='none')

# Pretty up and show the plot
_ = plt.xlabel('magnitude')
_ = plt.ylabel('ECDF')
_ = plt.xlim(2.8, 6.2)
plt.show()

# Report the results
print("""
b-value: {0:.2f}
95% conf int: [{1:.2f}, {2:.2f}]""".format(b, *conf_int))
Edit dan Jalankan Kode