Нормализация филогенетического дерева в R

При работе с данными филогенетического дерева в R (особенно при работе с объектами "phylo" или "phylo4") было бы полезно нормализовать длину ветвей, чтобы определенные таксоны (те, которые развиваются быстрее) не вносили непропорционально большую длину ветвей. к дереву. Похоже, что это обычное дело при вычислении значений UniFrac, как можно найти в обсуждении здесь: http://bmf2.colorado.edu/unifrac/help.psp. (Однако мне нужно больше, чем просто значения UniFrac).

Однако я не могу найти функцию, которая выполняет этот шаг нормализации. Я искал обезьяны, пиканте, адефило и филобазу. Может ли кто-нибудь направить меня к пакету, который включает эту функцию, или к пакету, который упрощает написание такой функции?


person Krisrs1128    schedule 27.07.2011    source источник


Ответы (1)


Вы ищете функцию для простого масштабирования ветвей дерева? Если так, compute.brlen() in ape сделает это. Есть встроенные опции для графена rho и all = 1. Вы также можете указать свою собственную функцию.

Я не знаю, выполняет ли UniFrac какое-либо другое масштабирование длины ветки. Но если это так, вы можете написать свою функцию и передать ее.

person kmm    schedule 07.08.2011
comment
Это делает то, о чем я просил, спасибо! Я не думаю, что UniFrac требует какого-либо другого масштабирования длины ветки. Кроме того, я обнаружил, что еще один способ сделать это - разрешить branch.sum ‹- [email protected], а затем установить [email protected]‹ - [email protected]/branch.sum. - person Krisrs1128; 08.08.2011