开始使用免费开始使用

假设检验:地震频率是否发生了变化?

显然,随着废水回注开始,地震频率大幅上升。不过,出于练习目的,您仍将进行一次假设检验。您不会检验 2010 年前后地震间隔时间是否具有相同分布,因为废水回注可能会影响分布。相反,您将假设两者具有相同的均值。因此,请计算以下假设对应的 p 值:2010 年前后的地震间隔时间均值相同。以"2010 年前的平均时间间隔减去 2010 年后的平均时间间隔"作为检验统计量。

本练习是课程的一部分

统计思维案例研习

查看课程

练习说明

  • 计算观测到的检验统计量。前面练习中的变量 mean_dt_premean_dt_post 已在您的命名空间中。
  • 平移 2010 年后的数据,使其均值与 2010 年前的数据相同。将结果存为 dt_post_shift
  • 分别对 dt_predt_post_shift 的均值抽取各 10,000 个自举重复样本。分别将结果存入 bs_reps_prebs_reps_post
  • 通过用 bs_reps_pre 减去 bs_reps_post 来计算"均值之差"的重复样本。
  • 计算并打印 p 值。将"至少同样极端"理解为:检验统计量大于或等于观测到的值。

交互式实操练习

通过完成这段示例代码来试试这个练习。

# Compute the observed test statistic
mean_dt_diff = ____ - ____

# Shift the post-2010 data to have the same mean as the pre-2010 data
dt_post_shift = ____ - ____ + ____

# Compute 10,000 bootstrap replicates from arrays
bs_reps_pre = ____
bs_reps_post = ____

# Get replicates of difference of means
bs_reps = ____ - ____

# Compute and print the p-value
p_val = ____(____ >= ____) / 10000
print('p =', p_val)
编辑并运行代码