Как рассчитать изменение частоты аллелей в гаплоидной популяции при отборе?

Из этой книги

Для простоты рассмотрим гаплоидный организм и предположим, что частоты аллелей А 1 а также А 2 даны Икс а также у знак равно 1 Икс , соответственно. Мы также предполагаем, что пригодность А 1 а также А 2 находятся ж 1 знак равно 1 а также ж 2 знак равно 1 с , соответственно. В этом случае средняя пригодность ж ¯ дан кем-то Икс + ( 1 Икс ) ( 1 с ) знак равно 1 с у , а изменение частоты аллеля за поколение становится

Δ Икс знак равно д Икс д т знак равно с Икс у 1 с у

Если бы мне пришлось найти то, что д Икс д т равно, я бы использовал уравнение Райта-Фишера и обнаружил, что:

д Икс д т знак равно ж 1 Икс ж ¯ знак равно Икс 1 с у

, что, очевидно, не совпадает с тем, что нашел автор...

Что мне не хватает? Как автор узнал этот результат д Икс д т знак равно с Икс у 1 с у ?

Перепроверьте все знаки равенства во втором уравнении (выглядит странно).
Ну, это просто переформулировать одно и то же уравнение по-другому. Я избавился от всего, что было в скобках, в надежде, что все покажется более ясным. Я отредактировал некоторые другие вещи, надеясь, что все станет яснее. Спасибо за помощь @fileunderwater
Привет, Реми, учитывая твою информацию, мне удалось получить ответ. Проверьте мой ответ ниже и спросите меня, нужно ли вам какое-либо дополнительное объяснение, это важно, поскольку я дал вам ответ, но я все еще не уверен, как вы ошиблись.
Вот простая проверка согласованности, чтобы показать, что Δ Икс знак равно Икс 1 с у не может быть правильным ответом: в чистом мономорфном А 1 популяции, очевидно, не может быть никаких изменений в частотах аллелей посредством отбора, поэтому Икс знак равно 1 должно подразумевать Δ Икс знак равно 0 . Ваша формула, однако, дает Δ Икс знак равно 1 1 знак равно 1 в этом случае.
Пс. Как математик я хотел бы заявить о своем возражении против использования обозначения д Икс д т для «изменения частоты аллеля за поколение»; он правильно обозначает скорость изменения частоты аллелей с течением времени (предположительно в постоянно размножающейся популяции, чтобы такая скорость была четко определена). Даже если мы измеряем время в (средних) поколениях, нетрудно привести примеры, в которых они не равны.

Ответы (1)

Вот мой полный вывод из приведенного вами примера книги, надеюсь, это поможет вам понять, что пошло не так:

введите описание изображения здесь

Вы должны помнить, что после отбора, действующего на популяцию, у вас больше нет 1 после отбора. Думайте об отборе как об «убийстве» особей, что означает, что общее количество теперь равно 1 минус то, что было «отобрано». s y — это то, что выбирается, поэтому новая сумма равна 1-s y. Это означает, что теперь вы делите каждую из ваших частот на 1-s*y (см. рисунок).

Пожалуйста, прокомментируйте, если вам нужно больше объяснений.

Хорошо, я полностью понял свою ошибку. Это было очень просто. Большое спасибо за это хорошее объяснение.
Самое замечательное, что вы поняли «самую сложную часть» в моделировании популяционной генетики. Процедура разбивки шагов по таблице и нормализации аналогична процедуре для диплоидов, но вместо A1 у вас есть A1B1, вместо x у вас есть 2*p q, где p — это f(A1), а q — это f(B1) . )), после выбора A1B1 это просто 2*p q - s(2*p*q) и т. д. Даже если смоделировать несколько отдельных локусов и где каждый локус может иметь> 2 аллелей, просто: (1) занесите в таблицу все возможные генотипы в виде столбцов, (2) сведите в таблицу частоты на каждом этапе жизненного цикла, (3) и обязательно нормализуйте, если есть выбор. :)