9

Допустим, Вы исследовали некоторую выборку из 8 аллелей из популяции, и первое событие коалесценции произошло 25 поколений назад. Какова вероятность этого события ($P(t_8 = 25)$), если популяция диплоидна и $N_e=1000$?

Округлите ответ до тысячных.

import pandas as pd
import numpy as np
from matplotlib import pyplot as plt

from math import factorial
def C(n, m):
    return (factorial (n)) / ((factorial (m)) * (factorial (n-m)))

N = 1000
x = 25
k = 8
p = (1-(C(k,2))/(2*N))**(x-1)*((C(k,2))/(2*N))
print(p)

0.009980998590710544