На лекции вы видели, что ведущее собственное значение марковской матрицы \(M\), вывод которой в 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).

В этом упражнении мы используем цикл for, чтобы смоделировать процесс мутаций, начиная с исходного распределения аллелей:

[1] 1 0 0 0

и убедиться, что это действительно так — собственный вектор даёт ту же информацию, что и цикл.

Подробнее о марковских процессах можно узнать по этой ссылке.

Вычисление собственных значений и собственных векторов в R