仮説検定：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)