假設檢定: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)