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