爾云間生信代碼|R可視化:復(fù)雜熱圖繪制技巧之熱圖中添加柱狀圖
科研有捷徑,輸入代碼,一鍵獲取科研成果!就是這么省事,來具體看下有多方便!
搜索http://985.so/a9kb查看全部代碼(目前共計(jì)50+持續(xù)新增中),也可以點(diǎn)擊右側(cè)【目錄】,可以看到更多有趣的代碼;真香提示:文末可以知道如何獲取代碼~
熱圖(Heatmap)一直是生信分析數(shù)據(jù)結(jié)果展示的好幫手,能很好的將結(jié)果清晰呈現(xiàn),在組學(xué)研究的相關(guān)文章中,熱圖作為一種對實(shí)驗(yàn)數(shù)據(jù)及其分析結(jié)果的直觀的表達(dá)方式,在很多文章中都有著不可或缺的地位。它是一種將規(guī)則化矩陣數(shù)據(jù)信息轉(zhuǎn)換成顏色色調(diào)的常用的可視化方法,其中每個(gè)單元對應(yīng)數(shù)據(jù)的某些屬性,屬性的值通過顏色映射轉(zhuǎn)換為不同色調(diào)并填充規(guī)則單元。

熱圖是非常常見的圖形,在R語言中有非常多的R包可以畫熱圖,比如pheatmap等,但complexheatmap包無疑是其中的老大哥,其功能之全面遠(yuǎn)超其他同類R包。ComplexHeatmap 不僅可以畫單個(gè)的熱圖,也可以制作多個(gè)熱圖及注釋信息相互組合起來的更為復(fù)雜的圖形。復(fù)雜熱圖可用于展示同一個(gè)數(shù)據(jù)集或不同數(shù)據(jù)集之間的關(guān)系或揭示內(nèi)部規(guī)律。ComplexHeatmap包可提供靈活的熱圖展示及高度自定義的注釋圖形。
高分期刊文獻(xiàn)往往不僅創(chuàng)新性好,工作量大,其中的配圖同樣美觀且能很好地反映數(shù)據(jù)內(nèi)容。下面我們來看下在熱圖的基礎(chǔ)上怎么增加柱狀圖來將圖形變得更高級。代碼相關(guān)文件見如下3個(gè)文件夾,文末可知如何領(lǐng)取。

1.????? 加載相應(yīng)R依賴包:
library("ComplexHeatmap")
2.????? 數(shù)據(jù)導(dǎo)入
輸入數(shù)據(jù)格式如下所示,我們用一個(gè)簡單的矩陣作為例子。

3.????? 基本熱圖繪制
?????? 先來看下這個(gè)數(shù)據(jù)的熱圖結(jié)果。


?????? 然后我們把在右邊的標(biāo)簽換到左邊,使用參數(shù)row_names_side = "left"。


4.????? 添加柱狀圖
接下來就是添加柱狀圖,也是本文最關(guān)鍵的一步,對熱圖每個(gè)基因的數(shù)據(jù)之和進(jìn)行統(tǒng)計(jì)并在熱圖正上方繪制柱狀圖,對每行數(shù)據(jù)之和進(jìn)行統(tǒng)計(jì)并在熱圖右側(cè)繪制柱狀圖。

到這里圖形就完成了,效果如下圖所示,大家可以嘗試用我們的示例數(shù)據(jù)跑起來,也可以用自己的數(shù)據(jù)進(jìn)行嘗試。

如需代碼及示例數(shù)據(jù)等文件,請掃碼聊天框回復(fù) “B31”領(lǐng)??!?

寫在文末:
如果您近期想做生信方面的文章而苦于沒有思路,或者不知道如何來入手生信分析,或者兌具體的某一個(gè)圖有作圖需求,都可以掃碼咨詢小云,我們有專業(yè)的技術(shù)團(tuán)隊(duì),生信熱點(diǎn)思路設(shè)計(jì)、生信分析、熱點(diǎn)方向生信挖掘等,如有需要,可掃碼下方二維碼了解詳情:

?