Вычисление b-значений
Судя по графику, построенному в предыдущем упражнении, можно смело использовать порог полноты mt = 3. Используя это значение, вычислите b-значения для периодов 1980–2009 годов и с 2010 года по середину 2017 года. Функция b_value(), написанная вами в прошлой главе, — она вычисляет b-значение и доверительный интервал по набору магнитуд и порогу полноты — уже доступна в вашем пространстве имён, как и массивы numpy mags_pre и mags_post из предыдущего упражнения, а также mt.
Это упражнение является частью курса
Примеры применения статистического мышления
Инструкции к упражнению
- Вычислите b-значение и 95%-й доверительный интервал для землетрясений за период с 1980 по 2009 год, используя 10 000 бутстрап-выборок.
- Вычислите b-значение и 95%-й доверительный интервал для землетрясений за период с 2010 года по середину 2017 года, используя 10 000 бутстрап-выборок.
- Нажмите «Отправить ответ», чтобы вывести результаты на экран.
Интерактивное практическое упражнение
Попробуйте выполнить это упражнение, дополнив этот пример кода.
# Compute b-value and confidence interval for pre-2010
b_pre, conf_int_pre = ____(____, ____, ____=[____, ____], ____=____)
# Compute b-value and confidence interval for post-2010
b_post, conf_int_post = ____
# Report the results
print("""
1980 through 2009
b-value: {0:.2f}
95% conf int: [{1:.2f}, {2:.2f}]
2010 through mid-2017
b-value: {3:.2f}
95% conf int: [{4:.2f}, {5:.2f}]
""".format(b_pre, *conf_int_pre, b_post, *conf_int_post))