#安装
source("https://bioconductor.org/biocLite.R")
biocLite("ggtree")
接着，我们需要准备好树文件和分组文件。
树文件就是我们通过Treebest、PAML、RAxML等工具获得的结果。
分组文件的格式如下：
sample group
ID A
ID1 A
ID2 B
ID3 B
美化树
# 加载R包
library("ggtree")
# 读取树文件（ggtree针对不同工具生成的树文件有不同的函数进行读取，如没有对应的函数，可以使用通用的办法读取）
tree <- read.tree("file.tree")
# 读取分组信息
group_file <- read.table("group_file.txt",header = T,row.names = 1)
# 按类分组
groupInfo <- split(row.names(group_file), group_file$Group)
# 将分组信息添加到树中
tree <- groupOTU(tree, groupInfo)
# 绘制进化树
ggtree(tree, layout="fan", ladderize = FALSE, branch.length = "none",aes(color=group)) + geom_tiplab2(size=3) + theme(legend.position = "right")
