Аббат, не смог прикрепить два выходных файла по Нибелунгам (из ТНТ и моей программы) к сообщению в личку. Попытаюсь - здесь, поскольку Ваш мэйл не знаю.
Там с расчётом возрастов такой напряг - моя прога считает их полагая все гаплотипы равной длины - 67, в то время как у Вас есть и 37-маркерные. Так в ветвях, где они присутствуют, нужно скорректировать возраст по формуле Внов=Встар*(1+N37/(N67+N37)), где N37 и N67 - количество соответственно 37 и 67-маркерных гаплотипов в ветке, к которой указанный возраст относится. Для всего дерева поправка получится небольшая, порядка 10% в плюс, поскольку 37-маркерных у Вас мало. Но в некоторых терминальных веточках возраст может увеличиваться больше, в пределе, до двух раз (там где одни 37-маркерные).
Причина неудобства - ТНТ не выдаёт мне действующее к-во маркеров в каждом гаплотипе - только число мутаций. Позже попытаюсь обойти это ограничение, но пока - так.
И еще - смотрите файлы в Блокноте, предварительно установив там шрифт Терминал и убрав перенос, иначе увидите кашу. Выложить его в текстовом виде на форум не получится, так как в нашем движке такого шрифта нет. Придётся Вам действовать через Принтскрин и Пэйнт, чтобы получить JPEG файл, который и выкладывается.
Успеха!
Тема: I1 разное (Прочитано 6648 раз)