在 R 中规范化系统发育树

Question

在 R 中处理系统发育树数据时（特别是在处理“phylo”或“phylo4”对象时），标准化分支长度会很有用，这样某些分类群（进化更快的分类群）不会贡献不成比例的数量树的分支长度。这似乎在计算 UniFrac 值时很常见，可以在此处的讨论中找到：http://bmf2.colorado.edu/unifrac/help.psp。 （不过，我需要的不仅仅是 UniFrac 值）。

但是，我找不到执行此标准化步骤的函数。我看过ape、picante、adephylo和phylobase。有人可以指导我找到包含此功能的包，或者可以直接编写此类功能的包吗？

Answer 1

您是否正在寻找一个函数来缩放树的分支长度？如果是这样，猿中的compute.brlen() 会这样做。 Grafen 的 rho 和 all = 1 有内置选项。您也可以提供自己的函数。

我不知道 UniFrac 是否会进行其他类型的分支长度缩放。但如果是这样，你可以编写你的函数并传递它。