АвторТема: yMurMur - Web-обертка над Murka  (Прочитано 2485 раз)

0 Пользователей и 1 Гость просматривают эту тему.

Оффлайн Daemon2017Автор темы

  • Сообщений: 2146
  • Страна: ru
  • Рейтинг +1041/-18
  • Y-ДНК: R1a-Y35174
  • мтДНК: V7-a2a2a2b*
Re: yMurMur - Web-обертка над Murka
« Ответ #75 : 06 Ноябрь 2023, 19:22:51 »
нет, почему, они осмыслены
просто стрелочка указывает на узел, который датирует подпись ко всему дереву в текстовых таблицах
эту величину можно усреднять по разным деревьям (ессно, она распределена нормально вокруг "ро", хотя ро может быть лажей))))

Т.е., например, узел A датируется с помощью 2 нижестоящих узлов (которые были датированы шагом ранее) и 5 нижестоящих нод. Но т.к. расстояния между А и нодами 3-7 оказались сильно меньше, чем между А и узлом 1, то они уменьшили общую датировку узла А?

Оффлайн Valery

  • Сообщений: 10098
  • Страна: 00
  • Рейтинг +1372/-7
  • Ultimate Matriarchy
Re: yMurMur - Web-обертка над Murka
« Ответ #76 : 06 Ноябрь 2023, 20:01:01 »
безусловно
Ро это же просто { сумма_по_i (частота узла i) * (мутаций от корня до i)  } / (общая частота)

Оффлайн Daemon2017Автор темы

  • Сообщений: 2146
  • Страна: ru
  • Рейтинг +1041/-18
  • Y-ДНК: R1a-Y35174
  • мтДНК: V7-a2a2a2b*
Re: yMurMur - Web-обертка над Murka
« Ответ #77 : 06 Ноябрь 2023, 20:23:12 »
безусловно
Ро это же просто { сумма_по_i (частота узла i) * (мутаций от корня до i)  } / (общая частота)

Спасибо! Теперь дошло :D

Оффлайн Semargl

  • "El sueño de la razón produce monstruos" ©
  • Сообщений: 5981
  • Страна: ru
  • Рейтинг +4154/-5
  • www.semargl.me
    • www.semargl.me
  • Y-ДНК: R1a [CTS3402+]
  • мтДНК: U4a2g
Re: yMurMur - Web-обертка над Murka
« Ответ #78 : 08 Ноябрь 2023, 10:26:33 »
нет, почему, они осмыслены
просто стрелочка указывает на узел, который датирует подпись ко всему дереву в текстовых таблицах
эту величину можно усреднять по разным деревьям (ессно, она распределена нормально вокруг "ро", хотя ро может быть лажей))))

Т.е., например, узел A датируется с помощью 2 нижестоящих узлов (которые были датированы шагом ранее) и 5 нижестоящих нод. Но т.к. расстояния между А и нодами 3-7 оказались сильно меньше, чем между А и узлом 1, то они уменьшили общую датировку узла А?
Теоретически, если режет глаз, то можно пройтись по дереву скриптом и на основании весов изменить возраста в ту или иную сторону. В указанном примере возраст узла "A" замещается возрастом узла потомка с наибольшим значением, т.е. возрастом "1", или возраст узла "1", замещается возрастом узла "A". Это если в случаях со снипами и точной топологией дерева. С STR мы не совсем уверены в топологии и плюс чем выше возраст, тем он более занижен. Но наверное все равно можно применить указанный метод коррекции.

Оффлайн Daemon2017Автор темы

  • Сообщений: 2146
  • Страна: ru
  • Рейтинг +1041/-18
  • Y-ДНК: R1a-Y35174
  • мтДНК: V7-a2a2a2b*
Re: yMurMur - Web-обертка над Murka
« Ответ #79 : 08 Ноябрь 2023, 11:14:13 »
Теоретически, если режет глаз, то можно пройтись по дереву скриптом и на основании весов изменить возраста в ту или иную сторону. В указанном примере возраст узла "A" замещается возрастом узла потомка с наибольшим значением, т.е. возрастом "1", или возраст узла "1", замещается возрастом узла "A". Это если в случаях со снипами и точной топологией дерева. С STR мы не совсем уверены в топологии и плюс чем выше возраст, тем он более занижен. Но наверное все равно можно применить указанный метод коррекции.

Спасибо! Боюсь вмешиваться в результаты работы Мурки, но идея опираться на нижестоящие возраста (мб с привлечением доверительных интервалов) мне нравится - мб сделаю какой-нибудь такой включаемый/отключаемый функционал)

Оффлайн Valery

  • Сообщений: 10098
  • Страна: 00
  • Рейтинг +1372/-7
  • Ultimate Matriarchy
Re: yMurMur - Web-обертка над Murka
« Ответ #80 : 08 Ноябрь 2023, 11:23:05 »
В новую версию мурки я бы включил исправленные доверительные интервалы - если yFull опубликует хотя бы упрощенную версию своего алгоритма. Собственно, если не будет новой версии сигмы, то не буду включать и старую версию ро-сигмы тк оригинальный метод Сейлард основанный на дисперсии ро - никудышний.


Оффлайн Valery

  • Сообщений: 10098
  • Страна: 00
  • Рейтинг +1372/-7
  • Ultimate Matriarchy
Re: yMurMur - Web-обертка над Murka
« Ответ #81 : 08 Ноябрь 2023, 11:25:17 »
Боюсь вмешиваться в результаты работы Мурки

в возраста вполне можно вмешиваться, тк остальное от них не зависит.

Оффлайн Semargl

  • "El sueño de la razón produce monstruos" ©
  • Сообщений: 5981
  • Страна: ru
  • Рейтинг +4154/-5
  • www.semargl.me
    • www.semargl.me
  • Y-ДНК: R1a [CTS3402+]
  • мтДНК: U4a2g
Re: yMurMur - Web-обертка над Murka
« Ответ #82 : 08 Ноябрь 2023, 11:28:42 »
Вероятно не совсем алгоритм YFull, но думаю так:

Для примера рассчитаем возраст узла "A":
1) Для начала рассчитаем веса дочерних нод и листьев:
нода 1: 2*(2+1+1+1+1+1) = 14
нода 2: 2*(2+1+1+1+1+1) = 14
лист 3: 1*(2+2+1+1+1+1) = 8
лист 4: 1*(2+2+1+1+1+1) = 8
лист 5: 1*(2+2+1+1+1+1) = 8
лист 6: 1*(2+2+1+1+1+1) = 8
лист 7: 1*(2+2+1+1+1+1) = 8
2) рассчитаем возраст узла с учетом возраста дочерних нод/листьев и полученных весов "A":
(309.68*14 + 61.94*14 + возраст3*8 + возраст4*8 + возраст5*8 + возраст6*8 + возраст7*8 )/(14 + 14 + 8 + 8 + 8 + 8 + 8 )
3) Если полученный таким методом возраст для "A" выше чем полученный просто расчетом ро, то корректируем его вверх с учетом возраста "1" или даже с учетом родителя для "A" (если расширить алгоритм).

Аналогично считаем взвешенные возраста для всех узлов. Расчет лучше делать от листьев к корню, а последующую корректировку от корня к листьям.

PS. В данных расчетах понятно что под возрастом листьев и узлов подразумеваются возраста ребер.
PPS. Расшифровка формулы расчета веса ноды
нода 1: 2(2 семпла в ноде 1)*(2(2 семпла в ноде 2)+1(1 семпл в листе 3)+1(1 семпл в листе 4)+1(1 семпл в листе 5)+1(1 семпл в листе 6)+1(1 семпл в листе 7)) = 14

Оффлайн Daemon2017Автор темы

  • Сообщений: 2146
  • Страна: ru
  • Рейтинг +1041/-18
  • Y-ДНК: R1a-Y35174
  • мтДНК: V7-a2a2a2b*
Re: yMurMur - Web-обертка над Murka
« Ответ #83 : 11 Ноябрь 2023, 16:37:04 »
Добавил поле "Средний возраст протестированных:", чтобы можно было принимать в расчет возраст протестированных. Средний возраст населения в РФ =39 годам, так что по умолчанию значение =40. Если хочется оставить вычисления Мурки неизменными, то можно сделать это поле =0.

 

© 2007 Молекулярная Генеалогия (МолГен)

Внимание! Все сообщения отражают только мнения их авторов.
Все права на материалы принадлежат их авторам (владельцам) и сетевым изданиям, с которых они взяты.