Игорь,
Можно ли модифицировать Ваш скрипт таким образом, чтобы он на "базальных" предковых ветвях (тех ветвях, которые соединяют узлы MRCA гаплогрупп, субкладов и кластеров) вместо возраста ветви по ро-статистике, производил бы оценку возраста по кросс-дисперсии по формула Нордтведта? Формула Нордтведта:
G(AB) = Var(AB) / 2M,
где Var(AB) -кросс-дисперсия пар гаплотипов из популяции A (например, I1) и популяции B (например, I2b1), М- скорость мутации.
Думаю, что самая трудная задача здесь, как запрограммировать разделение ветвей
на те, оценка возраста которых дается по ро-статистке, и те, рассчет возраста которых идет по интеркладовой дисперсии.