我希望在存储为multiPhylo对象的多个系统发育中运行相同的功能。
例如,假设我有1000棵树的multiPhylo,我想在每棵树中加总边/分支长度。我知道我可以使用一棵树:
sum(tree$edge.length)
但我无法解决如何在multiPhylo中为所有树做到这一点。我确信这很简单,但它超出了我的范围。有人可以帮忙吗?
由于
丹
答案 0 :(得分:3)
类multiPhylo
是一个列表(str(tree)
),因此提供R
用于处理列表的功能。要对各个树的边长进行求和,请使用lapply
函数。
lapply(tree, FUN = function(x) sum(x$edge.length))