1. 学ぶ
    2. /
  2. コース
    3. /
  3. Rで学ぶデータサイエンスのための線形代数
Connected

演習

対立遺伝子頻度のマルコフモデル

講義では、マルコフ行列 \(M\) の最大固有値に対応する固有ベクトルが、対立遺伝子が等しく表現される状況（それぞれ確率 0.25）を与えることを確認しました。R の出力は次のとおりです。

      [,1]  [,2]  [,3]  [,4]
[1,] 0.980 0.005 0.005 0.010
[2,] 0.005 0.980 0.010 0.005
[3,] 0.005 0.010 0.980 0.005
[4,] 0.010 0.005 0.005 0.980

この固有ベクトルは、対立遺伝子が等確率で表れる状況をモデル化します（各 0.25）。

この演習では、初期の対立遺伝子分布：

[1] 1 0 0 0

から突然変異の過程を for ループで反復し、実際にそのようになること、すなわち for ループを使わずとも固有ベクトルが正しい情報を与えることを示します。

マルコフ過程の詳細は、このリンクをご覧ください。

指示

100 XP
  • 1000 回の突然変異後の対立遺伝子分布 x を print() で表示してください。
  • M（すでに読み込まれています）の第1固有ベクトルを取り出し、和が 1 になるようにスケーリングして v1 に代入してください。
  • スケーリングした第1固有ベクトル v1 を表示し、x と比較してください。