教你畫一個掰彎的熱圖(Heatmap)省撑,展示更多的基因表達量


難道赌蔑,這不像孔雀開屏嗎?

寫在前面

組學數據已經泛濫竟秫,但是信息的挖掘仍任重道遠娃惯。順手的工具,可以節(jié)省使用者盡可能多的時間肥败,或者將看起來很復雜渐逃,很難以完成的事情等浊,變得非常簡單犁罩。
圖形的掰彎届良,這個是我很久以前已經實現的功能哄啄,當然這個功能的靈感(而不是實現邏輯)來至《The Grammer of Graphics》(這本書也是R包ggplot2的實現基礎)蟆融。圖形的掰彎姐仅,正如本公眾號早前的推文介紹刊头,事實是笛卡爾坐標系像極坐標系轉換的過程仗颈。坐標的轉換是一個事情佛舱,而圖形的變換是另一個事情,比如矩形要辦成中空的扇形等挨决。
熱圖请祖,我已經不止一次談及,如果早前有關注生信札記公眾號的朋友脖祈,那么應該知道肆捕,TBtools有兩個熱圖功能,Amazing HeatmapFancy Heatmap Browser盖高,各有所長慎陵。一般用過的人眼虱,大部分時候都會覺得真的是Amazing。后者目前似乎用的人不多席纽,不過我們課題組有一些相關成果會應用上捏悬。或許有時候润梯,過于非常規(guī)的出圖邏輯(往往也是更優(yōu)秀的發(fā)表級圖片)并不能為多數人所接受过牙。所以沒有人用,也是意料之內纺铭。
優(yōu)秀的人寇钉,總是會突破已有局限的。
我寫Amazing Heatmap在先彤蔽,而后續(xù)才實現了掰彎邏輯摧莽。掰彎邏輯的實現,本身是為了Circos圖工具的實現顿痪。但是我并沒有想到镊辕,會有TBtools用戶在我之前挖掘出了掰彎的熱圖的使用邏輯。
掰彎的熱圖蚁袭,目前已可在數個已發(fā)表工作中見到征懈。往往,使用這一功能的目的揩悄,可能是盡可能展示更多的基因卖哎。這個需求其實很容易產生。有時候删性,我們展示十幾個基因亏娜,那么常規(guī)的矩形熱圖效果就很不錯;但是蹬挺,當基因數目上升到數十上百维贺,那么可能圖版就不好放下。此時巴帮,圓環(huán)狀的優(yōu)勢就會體現出來溯泣。

工具的使用

首先,應該升級TBtools到最新版本榕茧,然后垃沦,還是打開Amazing Heatmap這個工具


然后輸入數據,我直接用一個好基友的數據(屏蔽了樣品和基因名字)用押,

于是得到下圖肢簿,

恩,基因數太多,一般可以縮小文字译仗,當然抬虽,更多的時候,是調整圖片的大小纵菌,同時也調一下配色阐污,得到下圖

很明顯,此時圖片太長咱圆。

掰彎這個熱圖和聚類樹

第一步是笛辟,大概調整一些文字與圖形的距離,主要原因是文字總歸是文字序苏,使用的空間與我們圖形元件使用的空間不同手幢,所以距離需要自己多試試..., 我同時也調整了樹的寬度


首先是一個提示操作,TBtools的所有圖片是交互的忱详,原本圖是這樣的


在圖形外的空白處围来,鼠標右擊兩次,可顯示隱藏元件


可以看出這個圖片是之后包含幾個Panel匈睁,如聚類樹的Panel监透,和熱力圖矩形的Panel,大體了解了Panel的位置之后航唆,鼠標右擊兩次胀蛮,恢復隱藏模式
開始掰彎熱圖
以下關鍵操作,
鍵盤摁住Shift鍵不放手糯钙,鼠標左鍵分別單擊熱圖的TilePanel粪狼,和TreePanel,大概是下述兩個位置任岸。這樣就同時選擇了兩個Panel再榄,

摁住Shift鍵不放,鼠標右鍵TilePanel享潜,此時會出現菜單不跟,選擇



如此,即可合并兩個Panel(用于后續(xù)統(tǒng)一轉換坐標系)
接下來的操作就比較簡單米碰,對合并后的Panel,鼠標右鍵购城,調出菜單吕座,



隨后在彈出菜單中,設置自己需要的極坐標的參數瘪板,

這些參數吴趴,可以依據自己需要,不斷細調侮攀。具體就看個人需要了
比如你看都不想看到一個gap



只是锣枝,這個有什么用嗎厢拭?
當然,可能不要忘了撇叁,Amazing HeatMaps供鸠,是可以自定義樹結構的。

一些無關緊要的事情陨闹。

其實楞捂,如果你確實花時間去了解一個工具,那么并不會有那么多使用問題趋厉。如果有寨闹,那么一般是你想要使用別人甚至是開發(fā)者的時間來實現你的夢想。既然如此君账,其實你需要做的是考慮支付一定額度的酬勞繁堡。
真正做事情的,往往遇到問題乡数,會自己找到解決辦法椭蹄。比如這位朋友,自己搞定一個圖瞳脓,最后整到微信頭像上


比如塑娇,我并沒有寫這個教程,但是看到朋友圈劫侧,別人自己自然就可以使用TBtools繪制出來

很多時候埋酬,完成一件事情,使用一個工具烧栋,達成一個目的写妥,需要的或許是解決問題的思維,而不是問到的就是學會的审姓。

寫在后面

近日珍特,事情挺多的∧拢或許生物學問題值得做的事情已經很多了扎筒,寫工具已經放在第二甚至第三位。正如我們所看到的酬姆,過得好的嗜桌,更多的還是做生物學問題的,而不是寫工具的辞色。
當然骨宠,寫工具,本身也是為了解決生物學問題,只是為誰解決生物學問題层亿,這是一個需要問自己的問題桦卒。

?著作權歸作者所有,轉載或內容合作請聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市匿又,隨后出現的幾起案子方灾,更是在濱河造成了極大的恐慌,老刑警劉巖琳省,帶你破解...
    沈念sama閱讀 206,013評論 6 481
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件迎吵,死亡現場離奇詭異,居然都是意外死亡针贬,警方通過查閱死者的電腦和手機击费,發(fā)現死者居然都...
    沈念sama閱讀 88,205評論 2 382
  • 文/潘曉璐 我一進店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來桦他,“玉大人蔫巩,你說我怎么就攤上這事】煅梗” “怎么了圆仔?”我有些...
    開封第一講書人閱讀 152,370評論 0 342
  • 文/不壞的土叔 我叫張陵,是天一觀的道長蔫劣。 經常有香客問我坪郭,道長,這世上最難降的妖魔是什么脉幢? 我笑而不...
    開封第一講書人閱讀 55,168評論 1 278
  • 正文 為了忘掉前任歪沃,我火速辦了婚禮,結果婚禮上嫌松,老公的妹妹穿的比我還像新娘沪曙。我一直安慰自己,他們只是感情好萎羔,可當我...
    茶點故事閱讀 64,153評論 5 371
  • 文/花漫 我一把揭開白布液走。 她就那樣靜靜地躺著,像睡著了一般贾陷。 火紅的嫁衣襯著肌膚如雪缘眶。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 48,954評論 1 283
  • 那天髓废,我揣著相機與錄音巷懈,去河邊找鬼。 笑死瓦哎,一個胖子當著我的面吹牛,可吹牛的內容都是我干的。 我是一名探鬼主播蒋譬,決...
    沈念sama閱讀 38,271評論 3 399
  • 文/蒼蘭香墨 我猛地睜開眼割岛,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了犯助?” 一聲冷哼從身側響起癣漆,我...
    開封第一講書人閱讀 36,916評論 0 259
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎剂买,沒想到半個月后惠爽,有當地人在樹林里發(fā)現了一具尸體,經...
    沈念sama閱讀 43,382評論 1 300
  • 正文 獨居荒郊野嶺守林人離奇死亡瞬哼,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內容為張勛視角 年9月15日...
    茶點故事閱讀 35,877評論 2 323
  • 正文 我和宋清朗相戀三年婚肆,在試婚紗的時候發(fā)現自己被綠了。 大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片坐慰。...
    茶點故事閱讀 37,989評論 1 333
  • 序言:一個原本活蹦亂跳的男人離奇死亡较性,死狀恐怖,靈堂內的尸體忽然破棺而出结胀,到底是詐尸還是另有隱情赞咙,我是刑警寧澤,帶...
    沈念sama閱讀 33,624評論 4 322
  • 正文 年R本政府宣布糟港,位于F島的核電站攀操,受9級特大地震影響,放射性物質發(fā)生泄漏秸抚。R本人自食惡果不足惜速和,卻給世界環(huán)境...
    茶點故事閱讀 39,209評論 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望耸别。 院中可真熱鬧健芭,春花似錦、人聲如沸秀姐。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,199評論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽省有。三九已至痒留,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間蠢沿,已是汗流浹背伸头。 一陣腳步聲響...
    開封第一講書人閱讀 31,418評論 1 260
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留舷蟀,地道東北人恤磷。 一個月前我還...
    沈念sama閱讀 45,401評論 2 352
  • 正文 我出身青樓面哼,卻偏偏與公主長得像,于是被迫代替她去往敵國和親扫步。 傳聞我的和親對象是個殘疾皇子魔策,可洞房花燭夜當晚...
    茶點故事閱讀 42,700評論 2 345

推薦閱讀更多精彩內容