pheatmap的安裝非常簡(jiǎn)單，只需要在R軟件中執(zhí)行一行安裝代碼即可
install.packages('pheatmap')
安裝完成后，我們來(lái)看如何使用pheatmap來(lái)繪制聚類熱圖
# 加載軟件包
library('pheatmap')
# 生成繪圖用的數(shù)據(jù)
test = matrix(rnorm(200), 20, 10)
test[1:10, seq(1, 10, 2)] = test[1:10, seq(1, 10, 2)] + 3
test[11:20, seq(2, 10, 2)] = test[11:20, seq(2, 10, 2)] + 2
test[15:20, seq(2, 10, 2)] = test[15:20, seq(2, 10, 2)] + 4
colnames(test) = paste("Test", 1:10, sep = "")
rownames(test) = paste("Gene", 1:20, sep = "")
這段代碼實(shí)際上是利用隨機(jī)數(shù)生成了一個(gè)20 X 10的矩陣。
為了模擬不同樣品和基因之間的差異，我們將第一行到第10行中奇數(shù)列的數(shù)值全部加3，將第11行到第20行的偶數(shù)列數(shù)值全部加2，將15行到20行的偶數(shù)列全部加4。最后將列名命名為T(mén)est1 ～ Test10，將行名命名為Gene1～Gene 20，最終生成的數(shù)據(jù)格式如下圖
默認(rèn)參數(shù)繪制圖形只需要執(zhí)行以下代碼。
pheatmap(test)
生成的圖片如下圖：
當(dāng)默認(rèn)參數(shù)不能滿足我們的需求時(shí)，我們可以根據(jù)自己的需要在此基礎(chǔ)上修改這個(gè)圖形。常見(jiàn)的一些參數(shù)設(shè)置如下：
pheatmap(test,?scale?=?"row")
# scale = “row”的含義是繪圖時(shí)按行進(jìn)行均一化。進(jìn)行均一化可以降低個(gè)別特殊樣品與其它樣品間的差異，這會(huì)使得其它樣品間的差異在圖形中更加顯著。一般我們?cè)诨诒磉_(dá)量進(jìn)行聚類分析時(shí)，均是常用的參數(shù)。繪圖結(jié)果如下圖所示：
如果需要對(duì)配色方案進(jìn)行修改，可以修改color參數(shù)，
# colorRampPalette函數(shù)可以設(shè)置3種顏色(只能是3種)，它可以根據(jù)給定的向量生成漸變色，這三個(gè)參數(shù)分別指定了最大值，中間值和最小值的顏色。繪圖結(jié)果如下：pheatmap(test,?color?=?colorRampPalette(c("navy",?"white",?"firebrick3"))(50))
同時(shí)，通過(guò)設(shè)置cluster_col和cluster_row參數(shù)可以控制是否取消對(duì)行或列進(jìn)行聚類分析，具體代碼及結(jié)果如下：
pheatmap(test,?cluster_col?=?FALSE)
show_rownames和show_colnames參數(shù)來(lái)控制是否顯示行名和列名，如下：
pheatmap(test,show_rownames=F,show_colnames=F)
display_numbers 和number_color 參數(shù)可以控制是否在圖中顯示數(shù)字及設(shè)置數(shù)字的顏色。
pheatmap(test,?display_numbers?=?TRUE,number_color?=?"blue")
cellwidth和cellheight兩個(gè)參數(shù)可以控制每個(gè)單元的長(zhǎng)度和寬度。參數(shù)main可以在圖片中添加標(biāo)題。
pheatmap(test,?cellwidth?=?15,?cellheight?=?12,?main?=?"Example?heatmap")
以上呢，就是在使用pheatmap繪制聚類熱圖時(shí)常用的一些參數(shù)?？梢钥吹?，使用pheatmap繪制聚類熱圖是非常簡(jiǎn)單快速的。通過(guò)組合不同的參數(shù)，我們可以控制最終生成的圖片的樣式與效果。更多的功能和參數(shù)可以通過(guò)執(zhí)行pheatmap命令查看pheatmap自帶的幫助文檔來(lái)獵??！本期聚類熱圖的繪制我們就分享完啦，敬請(qǐng)關(guān)注其他圖形繪制。
截止目前，百邁客成功舉辦了三場(chǎng)零基礎(chǔ)R語(yǔ)言繪圖培訓(xùn)班.而由于我們培訓(xùn)班主要是小班教學(xué)（20人左右），以及老師行程問(wèn)題，還有部分老師未趕上前三期培訓(xùn)，因此，現(xiàn)特征集R語(yǔ)言培訓(xùn)班東北、福建、西安、甘肅、北京場(chǎng)意向?qū)W員。如果您想學(xué)編程，想自己可以隨時(shí)繪制出編輯部要求的圖片，歡迎點(diǎn)擊下方按鈕聯(lián)系我們，我們會(huì)根據(jù)大家填寫(xiě)的意向情況，安排后期的生信培訓(xùn)班場(chǎng)次，說(shuō)不定下一場(chǎng)就在你的城市哦。
關(guān)注我們，隨時(shí)隨地學(xué)習(xí)