前言

在繪制熱圖之后，仍會(huì)保留每個(gè)熱圖組件的繪圖區(qū)域。因此，可以通過獲取組件的繪圖區(qū)域來添加或修改圖形

可以使用 list_components() 獲取熱圖的所有組件名稱或 viewport 名稱，例如

set.seed(123)
mat <- matrix(rnorm(80, 2), 8, 10)
mat <- rbind(mat, matrix(rnorm(40, -2), 4, 10))
rownames(mat) <- paste0("R", 1:12)
colnames(mat) <- paste0("C", 1:10)

ha_column1 <- HeatmapAnnotation(
  points = anno_points(rnorm(10)), 
  annotation_name_side = "left"
)
ht1 <- Heatmap(
  mat, name = "ht1", km = 2, 
  column_title = "Heatmap 1", 
  top_annotation = ha_column1, 
  row_names_side = "left"
)
ha_column2 <- HeatmapAnnotation(
  type = c(rep("a", 5), rep("b", 5)),
  col = list(type = c("a" = "red", "b" = "blue"))
)
ht2 <- Heatmap(
  mat, name = "ht2", row_title = "Heatmap 2", 
  column_title = "Heatmap 2",
  bottom_annotation = ha_column2, column_km = 2
)

ht_list <- ht1 + ht2 + 
  rowAnnotation(
    bar = anno_barplot(
      rowMeans(mat), width = unit(2, "cm"))
)
draw(
  ht_list, row_title = "Heatmap list", 
  column_title = "Heatmap list",
)

列出所有組件

> list_components()
 [1] "ROOT"                           "global"                         "global_layout"                 
 [4] "global-heatmaplist"             "main_heatmap_list"              "heatmap_ht1"                   
 [7] "ht1_heatmap_body_wrap"          "ht1_heatmap_body_1_1"           "ht1_heatmap_body_2_1"          
[10] "ht1_column_title_1"             "ht1_row_title_1"                "ht1_row_title_2"               
[13] "ht1_dend_row_1"                 "ht1_dend_row_2"                 "ht1_dend_column_1"             
[16] "ht1_row_names_1"                "ht1_row_names_2"                "ht1_column_names_1"            
[19] "annotation_points_1"            "heatmap_ht2"                    "ht2_heatmap_body_wrap"         
[22] "ht2_heatmap_body_1_1"           "ht2_heatmap_body_1_2"           "ht2_heatmap_body_2_1"          
[25] "ht2_heatmap_body_2_2"           "ht2_column_title_1"             "ht2_dend_column_1"             
[28] "ht2_dend_column_2"              "ht2_column_names_1"             "ht2_column_names_2"            
[31] "annotation_type_1"              "annotation_type_2"              "heatmap_heatmap_annotation_2"  
[34] "annotation_bar_1"               "annotation_bar_2"               "global-column_title_top"       
[37] "global_column_title"            "global-row_title_left"          "global_row_title"              
[40] "global-heatmap_legend_right"    "heatmap_legend"                 "global-annotation_legend_right"
[43] "annotation_legend"             

1. 裝飾函數(shù)

獲取了所有組件的 viewport 名稱之后，我們可以使用 seekViewport() 來獲取對(duì)應(yīng)的 viewport。

為了避免使用繁瑣的 viewport 名稱，ComplexHeatmap 包提供了 decorate_*() 形式的簡(jiǎn)便函數(shù)

• decorate_heatmap_body()：裝飾熱圖主體
• decorate_annotation()：注釋
• decorate_dend()：樹狀圖
• decorate_title()：標(biāo)題
• decorate_dimnames()：行列名稱
• decorate_row_names()：行名，相當(dāng)于 decorate_dimnames(..., which = "row")
• decorate_column_names()：列名，相當(dāng)于 decorate_dimnames(..., which = "column")
• decorate_row_dend()：行樹狀圖，相當(dāng)于 decorate_dend(..., which = "row")
• decorate_column_dend()：列樹狀圖，相當(dāng)于 decorate_dend(..., which = "column")
• decorate_row_title()：行標(biāo)題，相當(dāng)于 decorate_title(..., which = "row")
• decorate_column_title()：列標(biāo)題，相當(dāng)于 to decorate_title(..., which = "column")

這些函數(shù)需要傳入熱圖或注釋的名稱、繪圖塊，如果熱圖分塊了，需要設(shè)置分塊索引

例如，對(duì)于下面的熱圖

ht_list <- draw(ht_list, row_title = "Heatmap list", column_title = "Heatmap list", 
               heatmap_legend_side = "right", annotation_legend_side = "left")

在第一個(gè)熱圖的第二個(gè)分塊中，文本和線條注釋

decorate_heatmap_body("ht1", {
  grid.text("outlier", 1.5/10, 2.5/4, default.units = "npc")
  grid.lines(
    c(0.5, 0.5), c(0, 1), 
    gp = gpar(lty = 2, lwd = 2, col = "green"))
}, slice = 2)

為第一個(gè)熱圖的樹形圖的不同類別添加不同的背景顏色

decorate_column_dend("ht1", {
  # 獲取樹狀圖
  tree = column_dend(ht_list)$ht1[[1]]
  # 將樹狀圖分為兩類
  ind = cutree(as.hclust(tree), k = 2)[order.dendrogram(tree)]
  # 獲取每類的起始索引
  first_index = function(l) which(l)[1]
  # 獲取沒類的終止索引
  last_index = function(l) { x = which(l); x[length(x)] }
  x1 = c(first_index(ind == 1), first_index(ind == 2)) - 1
  x2 = c(last_index(ind == 1), last_index(ind == 2))
  # 為沒類設(shè)置矩形背景色
  grid.rect(x = x1/length(ind), width = (x2 - x1)/length(ind), just = "left",
            default.units = "npc", gp = gpar(fill = c("#FF000040", "#00FF0040"), col = NA))
})

為第一個(gè)熱圖的第二個(gè)分塊的行名，設(shè)置背景填充色

decorate_row_names("ht1", {
  grid.rect(gp = gpar(fill = "#FF000040"))
}, slice = 2)

為第一個(gè)熱圖的第一個(gè)分塊的行標(biāo)題設(shè)置背景填充色

decorate_row_title("ht1", {
  grid.rect(gp = gpar(fill = "#00FF0040"))
}, slice = 1)

在點(diǎn)圖注釋中添加一條水平線

decorate_annotation("points", {
  grid.lines(
    c(0, 1), unit(c(0, 0), "native"), 
    gp = gpar(col = "red"))
})