我想用3种不同的线宽在系统树上标记不同的分支。如何指定分支2具有width= 2,行10具有width= 4,其余行具有width=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))
谢谢!
发布于 2015-04-22 16:44:01
首先,您需要弄清楚边的顺序,当plot.phylo()
的参数方向设置为direction='leftwards'
时,它显然遵循以下规则
<>G211
一旦获得了边缘的映射,就必须为每个分支构建一个具有所需宽度的向量。向量位置与边缘位置匹配。使用此问题中的示例,此代码将使分支1到5具有width=4,其余分支将采用默认宽度(1):
plot(bird.orders, edge.width = c(4,4,4,4,4))
https://stackoverflow.com/questions/29731077
复制相似问题