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
