开始使用免费开始使用

假设检验:b 值是否不同?

请完成上一练习中勾勒的假设检验。变量 mags_premags_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)
编辑并运行代码