MulaiMulai sekarang secara gratis

Menghitung nilai b

Nilai b adalah metrik umum untuk mengukur seismisitas suatu wilayah. Anda bisa membayangkan bahwa Anda akan sering menghitungnya saat bekerja dengan data gempa bumi. Untuk tugas semacam ini yang sering dilakukan, sebaiknya tulis sebuah fungsi! Jadi, tulislah fungsi dengan tanda tangan b_value(mags, mt, perc=[2.5, 97.5], n_reps=None) yang mengembalikan nilai b dan (opsional, jika n_reps bukan None) interval kepercayaannya untuk sekumpulan magnitudo, mags. Ambang kelengkapan diberikan oleh mt. Argumen kata kunci perc memberikan persentil untuk batas bawah dan atas interval kepercayaan, dan n_reps adalah jumlah replikasi bootstrap yang digunakan dalam menghitung interval kepercayaan.

Latihan ini adalah bagian dari kursus

Studi Kasus dalam Pemikiran Statistik

Lihat Kursus

Petunjuk latihan

  • Definisikan fungsi dengan tanda tangan b_value(mags, mt, perc=[2.5, 97.5], n_reps=None) yang melakukan hal-hal berikut:
    • Potong magnitudo dari mags pada dan di atas ambang kelengkapan mt menggunakan pengindeksan Boolean. Simpan hasilnya dalam variabel m.
    • Hitung estimasi terbaik dari nilai b. Ingat, estimasi terbaik untuk nilai b adalah b = (m - mt)·ln(10). Simpan hasilnya dalam variabel b.
    • jika n_reps bukan None, lakukan hal berikut.
      • Ambil n_reps replikasi bootstrap dari mean m. Simpan hasilnya dalam variabel m_bs_reps.
      • Konversi replikasi bootstrap dari mean m menjadi replikasi nilai b. Simpan hasilnya dalam b_bs_reps.
      • Hitung interval kepercayaan dari replikasi bootstrap nilai b. Simpan hasilnya dalam conf_int.
    • Kembalikan b dan conf_int, atau hanya b jika n_reps adalah None.

Latihan interaktif praktis

Cobalah latihan ini dengan menyelesaikan kode contoh berikut.

____ ____:
    """Compute the b-value and optionally its confidence interval."""
    # Extract magnitudes above completeness threshold: m
    m = ____[____ >= ____]

    # Compute b-value: b
    b = ____

    # Draw bootstrap replicates
    if n_reps is None:
        return b
    else:
        m_bs_reps = ____

        # Compute b-value from replicates: b_bs_reps
        b_bs_reps = (____ - ____) * ____

        # Compute confidence interval: conf_int
        conf_int = ____
    
        return b, conf_int
Edit dan Jalankan Kode