Tính toán giá trị b

Giá trị b là một thước đo phổ biến cho mức độ hoạt động địa chấn của một khu vực. Bạn có thể sẽ cần tính nó thường xuyên khi làm việc với dữ liệu động đất. Với những tác vụ lặp lại như vậy, tốt nhất là viết một hàm! Vì vậy, hãy viết một hàm với chữ ký b_value(mags, mt, perc=[2.5, 97.5], n_reps=None) trả về giá trị b và (tùy chọn, nếu n_reps không phải None) khoảng tin cậy của nó cho một tập độ lớn, mags. Ngưỡng đầy đủ được cho bởi mt. Tham số từ khóa perc cho biết các phân vị cho cận dưới và cận trên của khoảng tin cậy, và n_reps là số bản sao bootstrap dùng để tính khoảng tin cậy.

Bài tập này là một phần của khóa học

Các nghiên cứu tình huống về tư duy thống kê

Xem khóa học

Hướng dẫn bài tập

Định nghĩa một hàm với chữ ký b_value(mags, mt, perc=[2.5, 97.5], n_reps=None) thực hiện các bước sau:
- Lọc các độ lớn trong mags tại và trên ngưỡng đầy đủ mt bằng chỉ mục Boolean. Lưu kết quả vào biến m.
- Tính ước lượng tốt nhất của giá trị b. Nhớ rằng ước lượng tốt nhất cho giá trị b là b = (m - m_t)·ln(10). Lưu kết quả vào biến b.
- nếu n_reps không phải None, thực hiện:
  - Rút n_reps bản sao bootstrap của trung bình của m. Lưu kết quả vào biến m_bs_reps.
  - Chuyển các bản sao bootstrap của trung bình m thành các bản sao của giá trị b. Lưu kết quả vào b_bs_reps.
  - Tính khoảng tin cậy từ các bản sao bootstrap của giá trị b. Lưu kết quả vào conf_int.
- Trả về b và conf_int, hoặc chỉ b nếu n_reps là None.

Bài tập tương tác thực hành trực tiếp

Hãy thử làm bài tập này bằng cách hoàn thành đoạn mã mẫu này.

____ ____:
    """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

Chỉnh sửa và Chạy Mã