1. 学ぶ
    2. /
  2. コース
    3. /
  3. Pythonで学ぶ統計思考（パート2）
Connected

演習

仮説検定：2012年のほうがクチバシは深いですか？

ECDF のプロットと信頼区間の結果から、Daphne Major にいる G. scandens のクチバシが深くなっていることはかなり明確に見えてきました。ですが、これは単なる偶然によるものという可能性はないでしょうか。つまり、平均が同じだとした場合に、観測された平均の差が得られる確率はどのくらいでしょうか？

注意してください。ここで検定している仮説は、クチバシの深さが同じ分布から来ているというものではありません。それには置換検定を使えます。ここでの仮説は、平均が等しいというものです。この仮説検定を行うには、2つのデータセットの平均が等しくなるようにシフトし、その後ブートストラップサンプリングを使って平均の差を計算します。

指示

100 XP
  • 1975年と2012年のクチバシの深さを連結した配列を作成し、その平均を計算して保存します。
  • 直前に計算した結合データの平均に等しくなるよう、bd_1975 と bd_2012 をシフトします。
  • 1975年と2012年のクチバシの深さについて、それぞれ平均のブートストラップ推定量を10,000個取得します。
  • 差のブートストラップ推定量を得るために、2012年の推定量から1975年の推定量を引きます。
  • p値を計算して表示します。前の演習で計算した観測された平均の差は、mean_diff として名前空間に残っています。