難道赌蔑,這不像孔雀開屏嗎?
寫在前面
組學數據已經泛濫竟秫,但是信息的挖掘仍任重道遠娃惯。順手的工具,可以節(jié)省使用者盡可能多的時間肥败,或者將看起來很復雜渐逃,很難以完成的事情等浊,變得非常簡單犁罩。
圖形的掰彎届良,這個是我很久以前已經實現的功能哄啄,當然這個功能的靈感(而不是實現邏輯)來至《The Grammer of Graphics》(這本書也是R包ggplot2的實現基礎)蟆融。圖形的掰彎姐仅,正如本公眾號早前的推文介紹刊头,事實是笛卡爾坐標系像極坐標系轉換的過程仗颈。坐標的轉換是一個事情佛舱,而圖形的變換是另一個事情,比如矩形要辦成中空的扇形等挨决。
熱圖请祖,我已經不止一次談及,如果早前有關注生信札記公眾號的朋友脖祈,那么應該知道肆捕,TBtools有兩個熱圖功能,Amazing Heatmap和Fancy 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繪制出來
很多時候埋酬,完成一件事情,使用一個工具烧栋,達成一個目的写妥,需要的或許是解決問題的思維,而不是問到的就是學會的审姓。
寫在后面
近日珍特,事情挺多的∧拢或許生物學問題值得做的事情已經很多了扎筒,寫工具已經放在第二甚至第三位。正如我們所看到的酬姆,過得好的嗜桌,更多的還是做生物學問題的,而不是寫工具的辞色。
當然骨宠,寫工具,本身也是為了解決生物學問題,只是為誰解決生物學問題层亿,這是一個需要問自己的問題桦卒。