在代謝組學(xué)的數(shù)據(jù)分析中，通常根據(jù)代謝物的含量進(jìn)行層次聚類(lèi)分析，以反映各個(gè)樣本之間代謝物含量的差異。

 

常用的數(shù)據(jù)分析工具中，R語(yǔ)言中的程序包可以靈活繪制出非常美觀的層次聚類(lèi)熱圖形。本文將利用R語(yǔ)言中的pheatmap包對(duì)繪制層次聚類(lèi)熱圖進(jìn)行詳細(xì)的介紹（詳細(xì)參數(shù)見(jiàn)文末彩蛋�。�

一、pheatmap( )函數(shù)常用參數(shù)簡(jiǎn)介

 

pheatmap(mat,  # 熱圖的輸入數(shù)據(jù)，數(shù)據(jù)類(lèi)型為數(shù)值型data.frame或matrix

color = colorRampPalette(rev(brewer.pal(n = 7, name = "RdYlBu")))(100),  # 熱圖顏色設(shè)置

cellwidth = NA,  # 設(shè)置熱圖單元格的寬度

cellheight = NA,  # 設(shè)置熱圖單元格的高度

treeheight_row = ifelse((class(cluster_rows) == "hclust") || cluster_rows, 50, 0),  # 設(shè)置行聚類(lèi)樹(shù)的高度

treeheight_col = ifelse((class(cluster_cols) == "hclust") || cluster_cols, 50, 0),  # 設(shè)置列聚類(lèi)樹(shù)的高度

cluster_rows = TRUE,  # 是否對(duì)行進(jìn)行聚類(lèi)

cluster_cols = TRUE,  # 是否對(duì)列進(jìn)行聚類(lèi)

cutree_rows = NA,  # 設(shè)置將行聚類(lèi)的結(jié)果分割成多少個(gè)集群，分割依據(jù)時(shí)基于層次聚類(lèi)結(jié)果，如果未對(duì)行進(jìn)行聚類(lèi)，則該參數(shù)被直接忽略

cutree_cols = NA,  # 設(shè)置將列聚類(lèi)的結(jié)果分割成多少個(gè)集群，分割依據(jù)時(shí)基于層次聚類(lèi)結(jié)果

annotation_row = NA,  # 在熱圖左側(cè)添加注釋?zhuān)�可用于顯示代謝物的分類(lèi)信息，數(shù)據(jù)為data.frame格式

annotation_col = NA,  # 在熱圖頂部添加注釋?zhuān)�可用于顯示分組信息，數(shù)據(jù)為data.frame類(lèi)型

border_color = "grey60",  # 熱圖內(nèi)部單元格的邊界顏色，可設(shè)置為NA無(wú)邊界

display_numbers = F,  # 是否在熱圖單元格中顯示相應(yīng)的數(shù)值

fontsize = 10,  # 基礎(chǔ)字體大小設(shè)置

fontsize_row = fontsize,  # 行名字體大小設(shè)置

fontsize_col = fontsize,  # 列名字體大小設(shè)置

filename = NA,  # 保存熱圖的全文件路徑，包含路徑和輸出的熱圖文件名（當(dāng)設(shè)置了filename時(shí)，圖形顯示器將不再顯示熱圖

…）

 

 

二、pheatmap( )函數(shù)繪制熱圖實(shí)戰(zhàn)

 第一步：安裝所需的R包及原始數(shù)據(jù)導(dǎo)入

install.packages("openxlsx")   # 導(dǎo)入excel數(shù)據(jù)所需的R包

install.packages("pheatmap")   # 繪制熱圖所需的R包

library(openxlsx)

library(pheatmap)

data <- read.xlsx(‘E: /R_TEST.xlsx’)  # 導(dǎo)入數(shù)據(jù)（案例中的數(shù)據(jù)非真實(shí)數(shù)據(jù)）

rownames(data) <- data[, 1]  # 設(shè)置行名


第二步：基礎(chǔ)熱圖-默認(rèn)參數(shù)

pheatmap(data[, c(3:8)])

第三步：設(shè)置顏色

pheatmap(data[, c(3:8)], color = colorRampPalette(c('red', 'white', 'blue'))(100))  

# 其中('red', 'white', 'blue')和(100)可根據(jù)需求自行修改

第四步：設(shè)置單元格大小

pheatmap(data[, c(3:8)], cellwidth = 18, cellheight = 18)

第五步：設(shè)置聚類(lèi)樹(shù)的高度

pheatmap(data[, c(3:8)], treeheight_row = 100)  

# 設(shè)置行（代謝物）聚類(lèi)樹(shù)的高度

pheatmap(data[, c(3:8)], treeheight_col = 100)  

# 設(shè)置列（樣本）聚類(lèi)樹(shù)的高度

第六步：取消聚類(lèi)

pheatmap(data[, c(3:8)], cluster_rows = FALSE)  

# 取消行聚類(lèi)

pheatmap(data[, c(3:8)], cluster_cols = FALSE)  

# 取消列聚類(lèi)

第七步：依據(jù)聚類(lèi)結(jié)果對(duì)熱圖進(jìn)行分割

pheatmap(data[, c(3:8)], cluster_rows = TRUE, cluster_cols = FALSE, cutree_rows = 3) 

# cutree_rows = 3 是根據(jù)聚類(lèi)結(jié)果進(jìn)行的具體值設(shè)置

第八步：添加樣本分組信息

annotation_col = data.frame(Group = factor(rep(c('A','B'), c(3,3))), row.names = colnames(data[, c(3:8)]))   

# 對(duì)每個(gè)樣本設(shè)置相對(duì)于的分組，(A、B為組名，3為每個(gè)組的樣本數(shù))

注意：原始數(shù)據(jù)表中每個(gè)分組的樣本放在一起，否則factor( )需一一對(duì)應(yīng)樣本名設(shè)置相對(duì)應(yīng)的分組名。
pheatmap(data[, c(3:8)], cluster_rows = TRUE, cluster_cols = FALSE, cutree_rows = 3, annotation_col = annotation_col)


第九步：添加代謝物分類(lèi)信息

annotation_row = data.frame(Class = data$Class, row.names = rownames(data))

pheatmap(data[, c(3:8)], cluster_rows = TRUE, cluster_cols = FALSE, cutree_rows = 3, annotation_col = annotation_col, annotation_row = annotation_row)

第十步：設(shè)置單元格邊框顏色，默認(rèn)顏色為灰色（grey30）

pheatmap(data[, c(3:8)], cluster_rows = TRUE, cluster_cols = FALSE, cutree_rows = 3, annotation_col = annotation_col, annotation_row = annotation_row, border_color = NA)

# border_color = NA 表示不設(shè)置單元格邊框顏色

第十一步：?jiǎn)卧�格�?nèi)顯示數(shù)值

pheatmap(data[, c(3:8)], cluster_rows = TRUE, cluster_cols = FALSE, cutree_rows = 3, annotation_col = annotation_col, annotation_row = annotation_row, border_color = NA, display_numbers = T)

第十二步：設(shè)置字體大小

pheatmap(data[, c(3:8)], cluster_rows = TRUE, cluster_cols = FALSE, cutree_rows = 3, annotation_col = annotation_col, annotation_row = annotation_row, border_color = NA, display_numbers = T, fontsize = 11)
也可使用fontsize_row和fontsize_col分別對(duì)行標(biāo)簽和列標(biāo)簽進(jìn)行字體大小設(shè)置。
pheatmap(data[, c(3:8)], cluster_rows = TRUE, cluster_cols = FALSE, cutree_rows = 3, annotation_col = annotation_col, annotation_row = annotation_row, border_color = NA, display_numbers = T, fontsize = 11, fontsize_row = 13, fontsize_col = 15)

第十三步：保存熱圖

pheatmap(data[, c(3:8)], cluster_rows = TRUE, cluster_cols = FALSE, cutree_rows = 3, annotation_col = annotation_col, annotation_row = annotation_row, border_color = NA, display_numbers = T, filename = 'E:/R_TEST.png')   
# filename為保存路徑和文件名
注意：設(shè)置filename后，RStudio圖形界面將不再顯示出熱圖。

 

彩蛋：pheatmap( )函數(shù)的參數(shù)詳解

網(wǎng)盤(pán)下載鏈接：

https://pan.baidu.com/s/1ao690Lpk2upH57_pu5AElw 

提取碼：98e4