假設檢定：地震發生頻率是否改變？
很明顯，當開始進行廢水回注後，地震發生頻率大幅上升。不過，作為練習，你仍會進行一次假設檢定。你不會去檢定 2010 年前後「地震間隔時間」是否來自同一個分配，因為廢水回注可能會改變分配形狀。相反地，你會假設它們的「平均數」相同。因此，請計算這個假設所對應的 p 值：以前－後（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)