1. 学ぶ
    2. /
  2. コース
    3. /
  3. 統計的思考 ケーススタディ
Connected

演習

ブートストラップ仮説検定

順列検定は、ヘテロ接合体と野生型のブート長が同一の分布に従うという、かなり制約の強い仮説に基づいています。次に、ブートストラップ仮説検定を使って、分布に関する仮定を設けずに「平均値が等しい」という仮説を検証してみましょう。

指示

100 XP
  • np.concatenate() を使い、野生型（bout_lengths_wt）とヘテロ接合体（bout_lengths_het）のブート長をすべて結合した配列 bout_lengths_concat を作成します。
  • 結合後の配列（bout_lengths_concat）から全ブート長の平均を計算し、結果を変数 mean_bout_length に格納します。
  • 両データセットの平均がいずれも mean_bout_length になるよう、それぞれをシフトします。シフト後の配列をそれぞれ wt_shifted と het_shifted に格納します。
  • dcst.draw_bs_reps() を使い、シフト後の各データセットの平均についてブートストラップ複製を10,000個生成します。それぞれの複製を bs_reps_wt と bs_reps_het に格納します。
  • bs_reps_het から bs_reps_wt を引いて平均値の差のブートストラップ複製を求め、結果を変数 bs_reps に格納します。
  • p値を計算します。「少なくとも同程度に極端」の定義は、帰無仮説のもとでの平均値の差が実験で観測された値以上であることとします。前の演習の変数 diff_means_exp はすでに名前空間に存在しています。