heatmap.2() 中的颜色零值

Question

我正在使用gplots 和Rcolorbrewer 使用heatmap2() 函数制作热图。
我的数据为零，表示此样本中没有记录任何数据。为了创建热图，我希望将这些零着色为白色，其余部分以三色图例显示。

这是目前为止的代码：

my_palette <- colorRampPalette(c("blue", "green", "red"))(n = 299)
col_breaks = c(seq(0.001,0.2,length=100), seq(0.3,0.4,length=100), seq(0.5,0.6,length=100))

library(gplots)
heatmap.2(deco.hmBB, main = "BCC at day 0 and day 10", margins = c(6,9), trace = "none", 
density.info = "none", dendrogram = "none", col = my_palette,  breaks =  col_breaks, 
Colv = "NA")

Answer 1

如评论所述，将零替换为 NA，如下所示：

deco.hmBB[deco.hmBB == 0] <- NA

现在重新运行热图，添加na.color。

heatmap.2(deco.hmBB, main = "BCC at day 0 and day 10", margins = c(6,9),
trace = "none", density.info = "none", dendrogram = "none", 
col = my_palette, breaks = col_breaks, Colv = "NA", na.color = "White")