首先利用mega进行多序列(clustalw)比对然后用最大似然法建树 树文件另存为*.nwk clustalw结果另存为fas结尾的fasta格式文件 利用ggtree进行可视化 library(ggtree) tree <- read.tree("tree.nwk") # 读入树文件 p1 <- ggtree(tree) + geom_tiplab(size = 4) + # 加物种名称 geom_treescale(fontsize=4, ...