ggtree是R語言里對進(jìn)化樹進(jìn)行可視化展示的一個功能非常強(qiáng)大的R包留潦，ggtree的作者還專門寫了一本書對ggtree的用法進(jìn)行了詳細(xì)的介紹，相關(guān)鏈接是 https://yulab-smu.top/treedata-book/勇边。最新版的ggtree還可以接受R語言里層次聚類分析的結(jié)果娩贷，畫聚類樹展示結(jié)果跋选，非常方便。我之前也錄制過視頻進(jìn)行介紹官脓。
今天的內(nèi)容是介紹一下
ggtree可視化進(jìn)化樹的一下常用操作
現(xiàn)在假設(shè)你已經(jīng)拿到了nwk格式的進(jìn)化樹文件协怒，如下
(Synergus:0.1976902387,(((((Periclistus:0.1403183720,Synophromorpha:0.0325185390)93:0.0313182375,(Xestophanes:0.0275715134,(Diastrophus:0.0456139475,Gonaspis:0.1146402107)97:0.0603746476)86:0.0275523221)91:0.0396704245,Ibalia:0.1295291852)93:0.0678466304,(((Liposthenes_ker:0.0568838340,Rhodus:0.4243267334)73:0.0825510697,Plagiotrochus:0.0778290252)71:0.0457931797,Phanacis_2:0.1416544135)42:0.0142517743)48:0.0209026386,(((Liposthenes_gle:0.1641119081,((((Antistrophus:0.1098867540,Hedickiana:0.2313789580)73:0.0566918206,Neaylax:0.1747090949)53:0.0027850349,(Isocolus:0.0980216531,Aulacidea:0.1315344980)40:0.0147148853)54:0.0123010924,((Andricus:0.0479556214,Neuroterus:0.0392025403)95:0.0395094917,Biorhiza:0.0640188941)87:0.0159496082)20:0.0000025961)50:0.0194234721,((((Panteliella:0.0792235900,Diplolepis:0.3184402599)84:0.0461941800,Phanacis_1:0.1153410113)66:0.0099961323,(Eschatocerus:0.2548694740,Parnips:0.0000022831)64:0.0802390069)34:0.0241704495,((Barbotinia:0.0731026287,Aylax:0.0957869567)87:0.0269932737,Iraella:0.0390833327)95:0.0797807340)18:0.0000021284)23:0.0095262346,Timaspis:0.0585073936)19:0.0170106400)57:0.0526944283,(Ceroptres:0.1057541047,(Pediaspis:0.1932340906,Paramblynotus:0.1711455809)28:0.0000021043)48:0.0416999011);
首先是讀取nwk格式的進(jìn)化樹文件
讀取nwk格式的進(jìn)化樹文件需要用到
treeio這個包中的read.newick()函數(shù)
library(treeio)
tree<-read.newick("ggtree_practice_aligned.fasta.treefile",
node.label = "support")
現(xiàn)在進(jìn)化樹的所有信息都存儲在了tree這個變量里
接下來是對進(jìn)化樹進(jìn)行可視化展示
最基本就是ggtree()函數(shù)，直接加讀進(jìn)來的樹文件
library(ggtree)
ggtree(tree)
image.png
添加文字標(biāo)簽
用到的的geom_tiplab()
ggtree(tree)+
geom_tiplab()
image.png
從上圖可以看到有的文字標(biāo)簽超出了繪圖邊界
可以首先加上theme_tree2()函數(shù)顯示出坐標(biāo)軸范圍卑笨，然后用xlim()函數(shù)更改坐標(biāo)軸范圍
ggtree(tree)+
geom_tiplab()+
theme_tree2()
image.png
ggtree(tree)+
geom_tiplab()+
theme_tree2()+
xlim(NA,0.8)
image.png
添加標(biāo)尺
ggtree(tree)+
geom_tiplab()+
theme_tree2()+
xlim(NA,0.8)+
geom_treescale(x=0.7,y=30,color="red")
image.png
這里遇到一個問題：geom_treescale()函數(shù)如果設(shè)置width參數(shù)孕暇，標(biāo)尺就顯示不出來，不知道是什么原因
更改樹的布局
這里布局的參數(shù)就不一一介紹了湾趾，可以參考 https://yulab-smu.top/treedata-book/chapter4.html
image.png
image.png
給指定的分組添加背景顏色
ggtree(tree)+
geom_tiplab()+
theme_tree2()+
xlim(NA,0.8)+
geom_treescale(x=0.7,y=30,color="red")+
#geom_text(aes(label=node))+
geom_highlight(node=35,fill="red")
image.png
ggtree(tree)+
geom_tiplab()+
theme_tree2()+
xlim(NA,0.8)+
geom_treescale(x=0.7,y=30,color="red")+
#geom_text(aes(label=node))+
geom_highlight(node=35,fill="red")+
geom_strip(6,11,label = "AAA",offset = 0.1,offset.text = 0.02,
color = "green",barsize = 3,fontsize = 5,angle = 90,
hjust = 0.5)
image.png
添加支持率的信息
ggtree(tree)+
geom_tiplab()+
theme_tree2()+
xlim(NA,0.8)+
geom_treescale(x=0.7,y=30,color="red")+
#geom_text(aes(label=node))+
geom_highlight(node=35,fill="red")+
geom_strip(6,11,label = "AAA",offset = 0.1,offset.text = 0.02,
color = "green",barsize = 3,fontsize = 5,angle = 90,
hjust = 0.5)+
geom_text(aes(label=support),hjust=0.5,vjust=0.5)
image.png
支持率可能會有部分重疊芭商，我暫時想不到如何用代碼把這些重疊分開，目前只能出圖后手動編輯
論文中通常只展示支持率大于某些值的搀缠，比如只顯示支持率大于75
tree@data$support1<-ifelse(tree@data$support<75,NA,tree@data$support)
ggtree(tree)+
geom_tiplab()+
theme_tree2()+
xlim(NA,0.8)+
geom_treescale(x=0.7,y=30,color="red")+
#geom_text(aes(label=node))+
geom_highlight(node=35,fill="red")+
geom_strip(6,11,label = "AAA",offset = 0.1,offset.text = 0.02,
color = "green",barsize = 3,fontsize = 5,angle = 90,
hjust = 0.5)+
geom_text(aes(label=support1),hjust=-0.5)
image.png
今天的內(nèi)容就到這里了
歡迎大家關(guān)注我的公眾號
小明的數(shù)據(jù)分析筆記本
