Стоп. Способ довольно простой, давайте в три действия проверим мои расчёты на конкретном примере.
Допустим, у меня два 67-маркёрных гаплотипа и между ними 7 мутаций. То есть М=7, N=2, L=67.
Считаем по
Вашей рекомендации общую лямбду:
lambda_obs =lambda_obs = M/ (L N)7/(2*67)=0,052
Отсюда просто лямда:
lambda = lambda_obs (1 + exp(lambda_obs))/2=0,052(1+exp(0,052))/2=0,054
ВБОП:
TMRCA = lambda/ mu=0,054/0,00216=25 поколений
где мю=0,00216 для 67-маркёрного гаплотипа
где мю=0,00243 для 37-маркёрного гаплотипа.
где мю=0,00184 для 25-маркёрного гаплотипа.
где мю=0,00183 для 12-маркёрного гаплотипа.