в R (эволюция белка в пределах одного вида)? Напишите в личку, пожалуйста
Сам по себе Р не строит деревьев, вам нужна сторонняя программа. Если дерево у вас уже есть, воспользуйтесь пакетом ape для визуализации
как это не строит? https://fuzzyatelin.github.io/bioanth-stats/module-24/module-24.html
Первая ссылка в гугле https://fuzzyatelin.github.io/bioanth-stats/module-24/module-24.html
Сразу дисклеймер: я не владею R, есть филогенетическое дерево из UGene, но мне нужно его построить не по родству, а по временной динамике накопления мутаций в белке
Вот еще вариант https://brouwern.github.io/lbrb/worked-example-building-a-phylogeny-in-r.html
Это просто реализация, причем не очень понятно чего. Я не нашла ссылки на описание модели, по крайней мере для maximum likelihood
Это же вообще другая задача. И зачем Вам для этого Р?
Мне посоветовали для этого R, но если Вы знаете другой способ построения таких деревьев, пожалуйста, поделитесь. В UGene на основе выравнивания у меня получается только дерево по родству
Я не знаю. Но я бы обратилась к человеку, который это посоветовал))
это реализация построения филогенетического дерева
Дело в том, что это было довольно давно, контакт потерян. И после неудачной попытки разобраться с этим самостоятельно я вновь прошу о помощи у более широкой аудитории
построить филогенетическое дерево так, чтобы оно не было филогенетическим? :)
А алгоритм-то какой? Я знаю как минимум три реализации maximum likelihood, и у них к всех получаются разные результаты)) Не говоря уже о том, что в каждой есть дюжина эволюционных моделей
Неточно выразилась, хотела чтобы было понятнее. Наверное это называется timetree или эволюционное дерево
Что у вас есть на входе? Выровненные последовательности?
Да, и дерево, построенное на их основе в UGene
То есть дерево уже есть и его нужно визуализировать?
Это не то дерево)) Это простое филогенетическое (хз каким алгоритмом юджин строит, но что-то стандартное должно быть)
There are three different phylogenetic trees building methods based on different algorithms: Neighbor-joining from the PHYLIP toolset, Bayesian inference from the MrBayes software and maximum likelihood from PhyML.
Просто для timetree нужны даты, а если их нет то как его строить то
Возможно даты в названии образцов в выравнивании стоят? Посмотрите
Я выбирала метод максимального подобия
*правдо*подобия?
Да, даты есть в названии образцов, но мне нужно сделать отдельно временную шкалу, и с помощью UGene это не удаётся
Timetree попробуйте
Нет, их тех поддержка ответила, что построить эволюционное дерево с помощью этого инструмента нельзя
Очень странно, тогда Mr Bayes, BEAST
Попробую, спасибо
Конечно, не достаточно! Люди десятки статей написали на эту тему, а тут такие - давайте максимизируем P(D|M)😂😂😂
А про какие три разные реализации вы говорили?
Ну как минимум RAxML, iqtree, fasttree. Есть и другие
Тут просто штука в том, что, по-сути, все инструменты врут, но кто-то врёт больше, кто-то меньше) Вот в статьях Вы не увидите NJ методов, там минимум нужен ML
Все, что я перечислила - это ML. А в статьях можно увидеть что угодно, зависит от задачи)
Обсуждают сегодня