b 值的量化
基于您在上一个练习中生成的图，可以放心使用完备性阈值 mt = 3。在该阈值下，分别计算 1980–2009 年与 2010 年至 2017 年年中这两个时段的 b 值。上一章您编写的函数 b_value()（用于基于一组震级和完备性阈值计算 b 值及其置信区间）已在您的命名空间中可用；上一个练习中的 numpy 数组 mags_pre 和 mags_post 以及 mt 也已可用。
本练习是课程的一部分
统计思维案例研习
练习说明
- 使用 10,000 次自助法重抽样，计算 1980–2009 年地震的 b 值及其 95% 置信区间。
- 使用 10,000 次自助法重抽样，计算 2010 年至 2017 年年中的地震 b 值及其 95% 置信区间。
- 点击 "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))