基因家族鉴定-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)

到这里画出树来。

上一篇下一篇

猜你喜欢

热点阅读