群体遗传学

ggtree绘制进化树

2019-04-20  本文已影响18人  271828182845904
library(ggplot2)
library(treeio)
library("ggtree")
tree<-read.tree('tree.tre')
tax<-read.table('group.info',row.names = 1)
names(tax)=c('area')
groupInfo <- split(row.names(tax),  tax$area) ## 按照area分组
tree<-groupOTU(tree,groupInfo)
ggtree(tree, aes(color=group))+  theme(legend.position = "right")+geom_tiplab(size=3)

文件说明:
tree.tre : 进化树文件


图片.png

group.info : 分组信息


图片.png
ggtree(tree, layout="circular", ladderize = FALSE, branch.length = "none",aes(color=group))+geom_tiplab2(size=3)+ theme(legend.position = "right")
图片.png
ggtree(tree, layout="circular", ladderize = FALSE,aes(color=group))+geom_tiplab2(size=3)+ theme(legend.position = "right")
图片.png
ggtree(tree, aes(color=group))+  theme(legend.position = "right")+geom_tiplab(size=3)
图片.png
上一篇下一篇

猜你喜欢

热点阅读