基因家族鉴定-ggtree-大神专用进化树美化R package
2022-09-14 本文已影响0人
小陈生信日记
我是个屁的大神,我不配用这个东西。一句两句三句报错连成线哈哈。
打开R编辑器,也就这步比较顺。
library('ggtree')
library('treeio')
library('ggplot2')
setwd('C:/Users/25799/Desktop')
inF='88and91.aln.contree'
tree=read.iqtree(inF)
data=read.table('groups.txt', header=T)
tree=as_tibble(tree)
get.fields(tree) ### "SH_aLRT", "UFboot"
#inF='88and91.nwk'
#tree=read.mega(inF)
#tree=read.newick(inF)
treeDat=full_join(tree, data, by='label')
treeDat=as.treedata(treeDat)
### show all node number in the tree graph
p<-ggtree(treeDat, layout='circular', branch.length='none', linesize=0.5) ## linetype="dotted")
p<-p+geom_text(aes(label=node), hjust=.5, size=2)
print(p)
到这里画出树来。