两棵进化树一致性的校对
首先, 要保证获得的进化树是有根树, 设定好外类群,并且,分支出现的次序经过排序reorder,这可以在Figtree软件中, 对newick格式的进化树进行操作, export as currently displayed 即可。
其次, 用R软件ape程序包的 cophyloplot函数, 可以实现, 将两棵进化树面对面绘制, 并将名称相同的分类单元用线段联系起来, 此时, 如果两个进化树结构有差别, 就一目了然了。
举例:
library(ape)
labs1 <- paste("t", 1:10, sep = "")
labs2 <- paste("t", 1:10, sep = "")
cophyloplot(x = rtree(20), y = rtree(10), space = 100, length.line = 0.1, gap = 10,assoc = as.matrix(data.frame(labs1, labs2)))
https://blog.sciencenet.cn/blog-255662-639595.html
上一篇:
R绘图一则 (2017年1月9日修订)下一篇:
从植物采集记录生成标本记录标签