Математические модели в естествознании
Вопросы - Математика и статистика
Другие вопросы по предмету Математика и статистика
>
. (29)
где -номера субпопуляций. В (28), (29) средние частоты и заданы формулами (26), а нормировочные коэффициенты -формулами (27). Отметим, что отображение (28), (29) не меняет вид после замены и . В связи с этим положим и , где . Если , то аллель имеет преимущество в борьбе за существование перед аллелем . Если же , то преимущество - за аллелем . Коэффициент легко преобразовать к виду:
.
В результате формулы (28) приобретают вид:
, (30)
Пусть во всех субпопуляциях отбор не происходит, т.е. . Используя (30) для средней частоты аллеля получаем:
Средние частоты аллелей сохраняются. Из (30) получаем, что с ростом номера поколения частоты . Результат естественен: миграция выравнивает частоты.
Предположим, что во всех субпопуляциях отбор действует против аллеля , т.е. . Предположим, что в -ом поколении хотя бы для одной субпопуляции частота аллеля меньше единицы. Очевидно, что и средняя частота аллеля . Тогда для всех субпопуляций получаем: . В этом случае из (30) следует, что . Для средней частоты выполнено неравенство , которое является строгим, пока по крайней мере в одной субпопуляции частота аллеля не равна единице. Последовательность монотонно растет. Она ограничена, а, следовательно, имеет предел. В предельной точке приращения нет, а это возможно только, если для всех субпопуляций . Тем самым, происходит повсеместное вытеснение аллеля . В случае, когда отбор действует против аллелей (т.е. ), совершенно аналогичные рассуждения показывают, что из популяции вытесняется аллель .
Направление отбора в разных субпопуляциях может быть различным -числа имеют разные знаки. В этом случае проверка сходимости итерационного процесса (30) представляет собой не простую задачу. Для равновесных частот и ( неподвижных точек отображения) получаем систему связанных между собой уравнений:
. (31)
Рассмотрим два частных случая. Пусть параметры малы. Смысл допущения заключается в том, что влияние миграции более существенно, чем процесс отбора. Уравнения (31) перепишем в виде:
. (32)
Отбрасывая малые слагаемые, приближенно получим: либо , либо . Второй вариант невозможен т.к. . Из (32) с точностью до слагаемых порядка получаем:
.
Поскольку , получаем, что либо , либо , или же удовлетворяет соотношению:
. (33)
Заметим, что правая часть уравнения -монотонно растущая функция от . Действительно, для тех , где , сдагаемые монотонно убывают. Наоборот, если , то соответствующие слагаемые монотонно растут. Пусть , т.е. отбор интегрально действует против аллелей . Тогда при правая часть уравненния положительна. Она положительна при всех , и уравнение (33) не имеет корней для . Если же , т.е. отбор интегрально действует против аллелей , то зеркальные рассуждения показывают, что уравнение также не может иметь состояний равновесия. Тем самым, при интенсивной миграции один из аллелей , или вытесняется из популяции. Какой конкретно аллель вытесняется, определяется знаком величины (если она положительна, то аллели имеют преимущество, и, наоборот, если она отрицательна, то преимущество принадлежит аллелям ). Прогнозы модели полностью согласуются с биологическим смыслом.
Пусть теперь . Перепишем уравнения (31) в виде:
.
Каждое из имеет два корня:
.
Обратим внимание, что полученные формулы не задают в явном виде решение системы (31), поскольку в правых частях фигурирует средняя частота , которая сама выражается через . Упростим формулу для корней. Деля числитель и знаменатель дроби на и пренебрегая слагаемыми порядка и по отношению к слагаемым порядка , последовательно получим:
,
,
.
Здесь использовано то, что для . Для получаем: . В свою очередь, для получаем: . Найдем среднюю частоту аллеля . Пусть для и для . Тогда, усредняя частоты в субпопуляциях, получим:
.
Отсюда следует, что
.
Средняя частота аллеля найдена. Тем самым, найдены и частоты аллелей в субпопуляциях. В тех субпопуляциях, где отбор действует против аллелей , частоты аллелей близки к единице. Наоборот, если отбор действует против аллелей , их частоты близки к нулю. Предсказания модели полностью согласуются с биологическим смыслом.
Дрейф генов
Дрейф генов -это случайные отклонения частот аллелей от теоретически ожидаемых, возникающие в результате недостаточного объема выборки. Такие явления часто называют ошибками выборки. Дрейф генов постоянно происходит в популяциях, поскольку их численность всегда конечна. Дополнительно заметим, что правильное представление о численности популяции дает не общее число особей, а число особей дающих начало следующему поколению. Действительно, только они дают вклад в генофонд следующего поколения.
Будем рассуждать в терминах аллелей, не переходя к генотипам. Рассмотрим популяцию аллелей и . Пусть априорно их частоты суть и . Случайным образом сформируем выборку из аллелей, которые оставят потомство. Пусть -число аллелей в выборке. Согласно теореме Муавра -Лапласа вероятность события , где , стремится при к числу . Здесь -нормальное распределение. В частности, если , то . Для эмпирической частоты аллеля в выборке получаем оценку: , которая выполнена с вероятностью . Поскольку , то . Чем длиннее выборка, тем эмпирическая частота ближе к априорной. Например, при получаем . Наоборот, при эмпирическая частота аллеля может принимать лишь одно из трех значений , т.е. эмпиричес?/p>