1. Lära sig
    2. /
  2. Courses
    3. /
  3. Pythonで学ぶ統計思考（パート2）
Connected

exercise

遺伝率を測る

ピアソンの相関係数は、共分散を2つのデータセットの分散の幾何平均で割った比でした。これは親と子の相関を測る指標ですが、遺伝率の推定としては最適ではないかもしれません。立ち止まって考えると、遺伝率は「親と子の共分散」を「親のみの分散」で割った比として定義するほうが理にかなっています。この演習では、遺伝率を推定し、ペア・ブートストラップにより 95% 信頼区間を求めます。

この演習は非常に重要な点を強調します。統計的推測（および一般的なデータ分析）は、公式に当てはめて機械的に計算するだけの作業ではありません。データで答えたい問いをよく考え、それにふさわしい分析を行う必要があります。形質の遺伝のしやすさに関心があるなら、既製の統計量であるピアソンの相関係数よりも、ここで定義した遺伝率のほうが適切です。

データは bd_parent_scandens、bd_offspring_scandens、bd_parent_fortis、bd_offspring_fortis に格納されています。

Instruktioner

100 XP
  • 親と子の形質の共分散を親の形質の分散で割った比として定義される遺伝率を計算する関数 heritability(parents, offspring) を作成してください。ヒント: このコースの前編で扱った np.cov() 関数を思い出してください。
  • この関数を使って、G. scandens と G. fortis の遺伝率を計算してください。
  • G. scandens と G. fortis について、ペア・ブートストラップを用いて遺伝率のブートストラップ標本を 1000 個取得してください。
  • 得られたブートストラップ標本から、両方の 95% 信頼区間を計算してください。
  • 結果を出力してください。