如何在R中显示ComplexHeatmap中的批注

Question

我有31个样本，得到了31种类型的样本，比如"a"，"b"，"c"，"d"，但当我绘制热图时，它说的是下标。我检查了矩阵的列是31，这与类型的数量相同。

ha1 = HeatmapAnnotation( bar = c("b","d","d","b","b","c","b","a","b"
                     ,"b","a","d","b","b","c","c","a","b"
                     ,"d","c","c","d","b","b","b","d","a"
                     ,"c","c","b","c"),
col = list(bar = c("a" = "red", "b" = "green", "c" = "blue","d"="black"))
 )
ha = rowAnnotation(foo = anno_mark(at = seq(0,7000,1000), labels = seq(0,7000,1000)))
hist_major <- Heatmap(major,
                  column_title = "Statistics for a list of SNP.Frequency.major",
                  column_title_gp = gpar(fontsize = 15, fontface = "bold"),
                  column_names_gp = gpar(fontsize = 10),
                  row_title = "Genome Postion(bp)",
                  row_title_gp = gpar(fontsize = 10),
                  name = "Variant.Frequency", 
                  cluster_rows = FALSE,
                  column_order=order(as.numeric(gsub("SAMN041259", "", colnames(major)))),
                  show_column_dend = FALSE,
                  right_annotation = ha,
                  bottom_annotation = ha1,
                  show_row_names = FALSE,
                  use_raster= TRUE,
                  raster_resize_mat = max,
                  col = cols,
                  border_gp = gpar(col = "black", lty = 2)
)

Answer 1

这是一个适用于您的代码的简单示例。也许它能帮到你。

library(ComplexHeatmap)

# Data matrix
major <- matrix(rnorm(3100),ncol=31)
colnames(major) <- paste0("Sample",1:31)

# Your code
ha1 = HeatmapAnnotation( bar = c("b","d","d","b","b","c","b","a","b"
                                 ,"b","a","d","b","b","c","c","a","b"
                                 ,"d","c","c","d","b","b","b","d","a"
                                 ,"c","c","b","c"),
                         col = list(bar = c("a"="red", "b"="green", "c"="blue", "d"="black"))
)
ticks <- seq(0,100,10)
ticks[1] <- 1
ha = rowAnnotation(foo = anno_mark(at = ticks, labels = ticks))
cols <- heat.colors(20)

hist_major <- Heatmap(major,
                      column_title = "Statistics for a list of SNP.Frequency.major",
                      column_title_gp = gpar(fontsize = 15, fontface = "bold"),
                      column_names_gp = gpar(fontsize = 10),
                      row_title = "Genome Postion(bp)",
                      row_title_gp = gpar(fontsize = 10),
                      name = "Variant\nFrequency", 
                      cluster_rows = FALSE,
                      column_order=order(as.numeric(gsub("SAMN041259", "", colnames(major)))),
                      show_column_dend = FALSE,
                      right_annotation = ha,
                      bottom_annotation = ha1,
                      show_row_names = FALSE,
                      use_raster= TRUE,
                      raster_resize_mat = max,
                      col = cols,
                      border_gp = gpar(col = "black", lty = 2)
)
print(hist_major)

​

​