Kapan gempa besar Parkfield berikutnya akan terjadi?
Gempa besar terakhir di wilayah Parkfield terjadi pada malam 27 September 2004 waktu setempat. Tugas Anda adalah memperkirakan kapan gempa Parkfield berikutnya akan terjadi, dengan mengasumsikan model Eksponensial dan juga model Gaussian. Pada kedua kasus, estimasi terbaik diberikan oleh selang waktu rata-rata, yang Anda hitung pada latihan sebelumnya sebesar 24,62 tahun, yang berarti gempa berikutnya diperkirakan pada 2029. Hitung interval kepercayaan 95% untuk waktu terjadinya gempa berikutnya dengan mengasumsikan sebaran Eksponensial yang diparameterkan oleh mean_time_gap yang Anda hitung pada latihan sebelumnya. Lakukan hal yang sama dengan mengasumsikan sebaran Normal yang diparameterkan oleh mean_time_gap dan std_time_gap.
Latihan ini adalah bagian dari kursus
Studi Kasus dalam Pemikiran Statistik
Petunjuk latihan
- Ambil 100.000 sampel dari sebaran Eksponensial dengan mean yang diberikan oleh
mean_time_gap. Simpan hasilnya dalamexp_samples. - Ambil 100.000 sampel dari sebaran Normal dengan mean yang diberikan oleh
mean_time_gapdan simpangan baku yang diberikan olehstd_time_gap. Simpan hasilnya dalamnorm_samples. - Karena hingga hari ini belum ada gempa Parkfield, lakukan slicing untuk mengambil sampel yang lebih besar dari
today - last_quake, di mana tahun desimal hari ini disimpan sebagaitoday, danlast_quake = 2004.74, tahun desimal dari gempa Parkfield terakhir. Timpa variabelexp_samplesdannorm_samplesmasing-masing dengan array hasil slicing tersebut. - Gunakan
np.percentile()untuk menghitung interval kepercayaan 95% untuk waktu terjadinya gempa Parkfield berikutnya. Dalam pemanggilan fungsi yang sama, Anda juga dapat menghitung median dengan menyertakan persentil ke-50.
Latihan interaktif praktis
Cobalah latihan ini dengan menyelesaikan kode contoh berikut.
# Draw samples from the Exponential distribution: exp_samples
exp_samples = ____
# Draw samples from the Normal distribution: norm_samples
norm_samples = ____
# No earthquake as of today, so only keep samples that are long enough
exp_samples = ____[____ > ____ - ____]
norm_samples = ____[____ > ____ - ____]
# Compute the confidence intervals with medians
conf_int_exp = ____(____, [____, ____, ____]) + last_quake
conf_int_norm = ____(____, [____, ____, ____]) + last_quake
# Print the results
print('Exponential:', conf_int_exp)
print(' Normal:', conf_int_norm)