Кількісна оцінка 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 бутстреп-реплік.
- Натисніть "Submit Answer", щоб вивести результати на екран.
Інтерактивна практична вправа
Спробуйте виконати цю вправу, доповнивши цей зразок коду.
# 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))