如何在系统发育树中显示分支的长度

Question

这里我有从newick格式绘制简单系统发育树的代码：

library(ape)
t<-read.tree(text="(F:4,(  (D:2,E:2):1,(C:2,(B:1,A:1):1):1):1);")
plot(t,use.egde.length=TRUE)

我是＆＃34;显示＆＃34;正确的分支长度，但我希望所有分支都有labal。

编辑： 我希望我的情节看起来像这样： 我正在搜索文档，但我找不到在R中显示分支长度的方法。我该怎么做？

Answer 1

您可以通过提取边长并使用edgelabels()来完成。

# Load package
library(ape)

# Create data
t <- read.tree(text="(F:4,((D:2,E:2):1,(C:2,(B:1,A:1):1):1):1);")
plot(t)
edgelabels(t$edge.length, bg="black", col="white", font=2)

Answer 2

这是你可以得到你想要的情节的方法：

t$tip.label <- c("F\n4", "D\n2", "E\n2", "C\n2", "B\n1", "A\n1")
plot(t,show.node.label=TRUE, show.tip.label=TRUE)

然而，我并不知道如何在不手动完成长度的情况下提取长度。