ПочатиПочніть безкоштовно

Обчислення значення b

Значення b — поширений показник сейсмічності регіону. Його часто доводиться обчислювати під час роботи з даними про землетруси. Для таких рутинних завдань найкраще написати функцію. Тож створіть функцію з сигнатурою b_value(mags, mt, perc=[2.5, 97.5], n_reps=None), яка повертає значення b і (за потреби, якщо n_reps не None) його довірчий інтервал для набору магнітуд mags. Поріг повноти задано як mt. Параметр perc задає перцентилі для нижньої та верхньої меж довірчого інтервалу, а n_reps — кількість бутстрепових реплік, що використовуються для його обчислення.

Ця вправа є частиною курсу

Кейс-стаді зі статистичного мислення

Переглянути курс

Інструкції до вправи

  • Створіть функцію з сигнатурою b_value(mags, mt, perc=[2.5, 97.5], n_reps=None), яка виконує таке:
    • Виділіть магнітуди з mags, що дорівнюють або перевищують поріг повноти mt, використовуючи булеву індексацію. Збережіть результат у змінній m.
    • Обчисліть найкращу оцінку значення b. Нагадаємо, найкраща оцінка для b — це b = (m - mt)·ln(10). Збережіть результат у змінній b.
    • якщо n_reps не None, виконайте таке.
      • Згенеруйте n_reps бутстрепових реплік середнього m. Збережіть результат у змінній m_bs_reps.
      • Перетворіть бутстрепові репліки середнього m на репліки значення b. Збережіть результат у b_bs_reps.
      • Обчисліть довірчий інтервал з бутстрепових реплік значення b. Збережіть результат у conf_int.
    • Поверніть b і conf_int, або лише b, якщо n_reps дорівнює None.

Інтерактивна практична вправа

Спробуйте виконати цю вправу, доповнивши цей зразок коду.

____ ____:
    """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
Редагувати та запускати код