寫在前面
很多時候,我們做完轉(zhuǎn)錄組煤辨,會想著看看樣品與樣品之間的表達(dá)相關(guān)性如何裳涛?依此,可以看看是否存在異常樣品众辨,也或許可以找到一些生物學(xué)相關(guān)問題端三。常常,我會使用 R 語言的 Corr 函數(shù)鹃彻,然后用 pheatmap 出個圖郊闯。當(dāng)然,這個兩行命令就解決了蛛株。但是呢 pheatmap 出的熱圖調(diào)整起來還是麻煩团赁。為了偷懶,我決定花點(diǎn)時間谨履,在 TBtools 中直接寫一個欢摄。因?yàn)檫@個功能比較簡單,其實(shí)就是讀取表達(dá)量矩陣笋粟,計算樣品間相關(guān)系數(shù)怀挠,用系數(shù)矩陣?yán)L制熱圖析蝴。
Expression Corr Calc
感覺就是十來分鐘的事情,真的太簡單绿淋。打界面麻煩了點(diǎn)闷畸,不過我也有些好的組件,所以吞滞,其實(shí)也簡單佑菩。結(jié)果如下,
從界面來看冯吓,使用更簡單倘待,只需要給一個基因表達(dá)量矩陣即可。
點(diǎn)擊 Start 之后组贺,可以看到輸出的 相關(guān)系數(shù)矩陣凸舵。
方便起見,同步輸出 Heatmap失尖,如下
接下來就是按照熱圖常見操作啊奄,做個微調(diào),得到圖稿如下
看起來不錯掀潮,當(dāng)然其實(shí)就是調(diào)用 TBtools 的 Heatmap菇夸,所以你可以做各種各樣的微調(diào).....
寫在最后
Emmm,又是一個簡單的小工具仪吧。事實(shí)上庄新,我很久以前也有想過寫這個,但是呢....用不上薯鼠。這次寫的原因簡單择诈,我覺得原有的方式出來的熱圖,不是我要的效果出皇,因?yàn)椴粔蚝每础?br> 但是羞芍,TBtools 出的熱圖,好看郊艘!而且荷科,細(xì)節(jié)容易調(diào)整。