開始使用免費開始

假設檢定:b 值是否不同?

請依照前一個練習所示意的步驟執行假設檢定。變數 mags_premags_post 已載入到你的命名空間,且 mt = 3 也已設定。

本練習屬於課程

統計思維個案研究

檢視課程

練習說明

  • 切出 2010 年前、震級大於或等於完備性門檻的地震,並以結果覆寫 mags_premags_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)
編輯並執行程式碼