0

我目前正在使用猿包。我的想法是提取与每个节点相对应的提示标签以及节点之间的距离。我怎样才能做到这一点?

非常感谢你的帮助

4

1 回答 1

0

以下是如何解释phylo对象中节点/提示之间的距离(请注意,dispRity::tree.age您需要dispRity从 GitHub https://github.com/TGuillerme/dispRity安装):

set.seed(1)
tree <- rtree(5)
tree$node.label <- paste0("n", 6:9)
plot(tree)
nodelabels(tree$node.labe)
axisPhylo()

## The distance between each tip
cophenetic(tree)

## The height of each tip and node
library(dispRity)
dispRity::tree.age(tree, order = "present")

## The distance between each nodes/tips
distances <- cbind(tree$edge, tree$edge.length)
colnames(distances) <- c("from", "to", "distance")
distances

该距离表读取为从元素(尖端/节点)n 到元素(尖端/节点)m。编号的工作方式如下:1 是树中的第一个提示 ( tree$tip.label[1]),2 是第二个等等... 6 (提示数 + 1) 是树中的第一个节点 ( tree$node.label[1]),等等。

于 2018-01-31T00:44:14.050 回答