特别声明：本部分（系列）内容均来自顾祖光博士对ComplexHeatmap的介绍，仅为学习交流，尊重原创。

闲话少说，先说顾老师：
顾祖光老师，现在德国癌症中心从事与NGS相关的生物信息学研究，其本科和研究生均在南京大学生命科学学院就读：

• 2001-09 to 2005-07 | Undergraduate (School of Life Science, Nanjing University)
• 2005-09 to 2008-07 | Master (School of Life Science, Nanjing University)
• 2008-09 to 2012-07 | PhD (School of Life Science, Nanjing University)

为什么要介绍顾祖光博士，因为很有可能你所需要的工具（尤其是可视化），很多你都能从他那儿找到答案，看看顾祖光老师有多高产：

R

安装

• 通过 bioconductor 安装：

if (!requireNamespace("BiocManager", quietly = TRUE))
    install.packages("BiocManager")
BiocManager::install("ComplexHeatmap")

• 通过 github 安装最新版

library(devtools)
install_github("jokergoo/ComplexHeatmap")

基础——单个热图绘制

library(ComplexHeatmap)
library(dplyr)
library(tidyverse)
library(circlize)

#生成矩阵
mat = matrix(rnorm(100, mean = 0, sd = 1), nrow = 10)
rownames(mat) = paste('row', seq_len(10), sep = '_')
colnames(mat) = paste('col', seq_len(10), sep = '_')

#绘制基础热图
Heatmap(mat)

• 指定热图颜色（包括连续型变量、离散型变量、NA值配色以及热图边界配色四部分内容）

热图中的数据一般都是连续型变量，ComplexHeatmap推荐大家使用circlize包中的colorRamp2进行颜色映射设定，函数参数包括断点和断点对应的颜色，这个包也是顾博士的。

col = colorRamp2(breaks = c(-4, 0, 4), colors = c('green', 'white', 'red'))
Heatmap(mat, name = 'mat', col = col)

这个颜色还有另外的妙用，我们常常想比较不同热图中的值的大小，但是又担心不同热图的颜色与数值大小的映射标准不一，ComplexHeatmap为我们提供了解决方案：只要几个热图共用一个相同的配色方案，那么它们之间就是可比的。

当然你也可以通过直接输入含有一系列颜色的向量来对热图进行颜色指定，不过需要注意：你所指定的颜色ComplexHeatmap是默认将其映射到整个数据矩阵的数据范围内，这意味着异常值会影响到整个热图的颜色配置。举例来说，这两种配色是一样的：

col = rainbow(3)
col = colorRamp2(breaks = seq(min(mat), max(mat), length = 3), colors = rainbow(3))

如果热图中的数据是离散型变量，比如是不同的字符，我们应该如何绘制呢？

#离散型变量
#生成矩阵
mat = matrix(sample(letters[1:4], 100, replace = T), ncol = 10)
#指定颜色（structure()函数）
col = structure(rainbow(4), names = letters[1:4])
#绘制热图
Heatmap(mat, name = 'mat', col = col)

col = rainbow(4)
names(col) <- letters[1:4]

那么我们如何处理NA值的配色？（在热图的绘制当中是允许存在NA值的）

#NA值配色
#生成矩阵
index = sample(c(TRUE, FALSE), 100, replace = T, prob = c(1, 9))
mat[index] = NA
#绘制热图
Heatmap(mat, name = 'mat', col = c('blue', 'white', 'red'))

Heatmap(mat, name = 'mat', col = c('blue', 'white', 'red'), na_col = 'black')

mat = matrix(sample(letters[1:4], 100, replace = T), ncol = 10)
col = rainbow(4)
names(col) <- letters[1:4]
Heatmap(matrix = mat, col = col, rect_gp = gpar(col = 'white', lwd = 2))

Heatmap(matrix = mat, col = col, border_gp = gpar(col = 'black', lty = 2))

• 自定义标题（这部分主要包括如何给热图的行和列添加标题）

这部分内容概括为以下代码框架：

Heatmap(matrix, #绘图矩阵 [matrix]
        col, #矩阵配色 [vector]
        name, #矩阵图例标题 [character]
        column_title, #矩阵列标题 [character]
        column_title_gp, #指定矩阵列标题相关属性 
        column_title_side, #指定矩阵列标题位置（顶部和底部）[character; 'top' or 'bottom']
        column_title_rot #指定列标题旋转角度，仅支持水平和竖直两种呈现形式 [numeric]
        )

其余参数没什么好讲的，主要介绍一下 column_title_gp 参数。这个参数主要是通过gpar函数指定标题的性质。

column_title_gp = gpar(fontsize = 20,  #指定字体大小
                       fontface = "bold",  #指定字体风格
                       fill = "red",  #指定标题背景填充色
                       col = "white",  #指定字体颜色
                       border = "blue" #指定标题背景框边框色
                       )

例如：

mat = matrix(rnorm(100, mean = 0, sd = 1), nrow = 10)
Heatmap(matrix = mat,
        col = col,
        name = 'mat',
        column_title = 'The heatmap for mat',
        column_title_side = 'top', 
        column_title_rot = 0,
        column_title_gp = gpar(fontsize = 20,
                               fontface = 'bold', 
                               fill = 'red',
                               col = 'white',
                               border = 'blue'))

• 聚类

默认情况下，绘制热图时会对行和列进行聚类，我们可以通过cluster_rows和cluster_cols关闭聚类功能：

mat = matrix(rnorm(100, mean = 0, sd = 1), nrow = 10)
rownames(mat) = paste('row', seq_len(10), sep = '_')
colnames(mat) = paste('col', seq_len(10), sep = '_')
Heatmap(matrix = mat, cluster_rows = F, name = 'mat')

Heatmap(matrix = mat, show_row_dend = F, name = 'mat')

当然有的时候，树是很重要的，你可能不能进行隐藏，但是你可以把树的size变小：

Heatmap(mat, name = "mat", 
        column_dend_height = unit(5, "mm"), 
        row_dend_width = unit(5, "mm"))

如你所想，聚类的树状结构你可以放在任何位置，这是通过row_dend_side和col_dend_side实现的：

Heatmap(matrix = mat, row_dend_side = 'right', name = 'mat')

本期到此为止，下期内容：行（列）顺序、行（列）名。