R 指定绘图中混凝土分支的宽度。phylo

我想用3种不同的线宽在一个系统发育树上标记不同的分支。 如何区分分支2的宽度为2，行10的宽度为4，其余行的宽度为1？ 例如，bird.orders数据集中有88个分支

data(bird.orders)
bird.orders@edge

这显然没有任何作用

plot(bird.orders, edge.width = c(bird.orders$edge[2]==2, bird.orders$edge[10]==4, bird.orders$edge[-c(2,10)]==1))

---

谢谢大家!

首先，您需要计算边缘顺序，当设置

plot.phylo（）

的参数方向时，它显然遵循这些规则

direction='leftwards'

：

向左走到最远的枝尖（叶）

然后向右走

然后下去

每次分支分裂时，重复1到3次

一旦获得了边的映射，就必须为每个分支构建一个具有所需宽度的向量。Wector位置与边缘位置匹配。使用此问题中的示例，此代码将使分支1到5的宽度为4，其余将采用默认宽度（1）：

plot(bird.orders, edge.width = c(4,4,4,4,4))