假设检验:b 值是否不同?
请完成上一练习中勾勒的假设检验。变量 mags_pre 和 mags_post 已加载到您的命名空间,且有 mt = 3。
本练习是课程的一部分
统计思维案例研习
练习说明
- 筛选 2010 年之前地震中震级大于(或等于)完备性阈值的记录,并用结果覆盖
mags_pre。对mags_post做同样处理。 - 计算观测到的均值差,用 2010 年之后地震的震级均值减去 2010 年之前地震的震级均值。
- 使用
dcst.draw_perm_reps()生成 10,000 个置换重复。func参数使用dcst.diff_of_means。 - 计算并打印 p 值,其中"至少同样极端"指检验统计量比观测值更小。
交互式实操练习
通过完成这段示例代码来试试这个练习。
# Only magnitudes above completeness threshold
mags_pre = ____[____ >= ____]
mags_post = ____[____ >= ____]
# Observed difference in mean magnitudes: diff_obs
diff_obs = ____ - ____
# Generate permutation replicates: perm_reps
perm_reps = ____
# Compute and print p-value
p_val = ____(____ < ____) / 10000
print('p =', p_val)