假设检验:地震频率是否发生了变化?
显然,随着废水回注开始,地震频率大幅上升。不过,出于练习目的,您仍将进行一次假设检验。您不会检验 2010 年前后地震间隔时间是否具有相同分布,因为废水回注可能会影响分布。相反,您将假设两者具有相同的均值。因此,请计算以下假设对应的 p 值:2010 年前后的地震间隔时间均值相同。以"2010 年前的平均时间间隔减去 2010 年后的平均时间间隔"作为检验统计量。
本练习是课程的一部分
统计思维案例研习
练习说明
- 计算观测到的检验统计量。前面练习中的变量
mean_dt_pre和mean_dt_post已在您的命名空间中。 - 平移 2010 年后的数据,使其均值与 2010 年前的数据相同。将结果存为
dt_post_shift。 - 分别对
dt_pre和dt_post_shift的均值抽取各 10,000 个自举重复样本。分别将结果存入bs_reps_pre和bs_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)