如何从R语言pheatmap中获取颜色刻度值（行z-score值）

嘿，这个问题我熟！当你用pheatmap()生成带行z-score标准化的热图后，想要从返回的对象里提取颜色刻度值和对应的行z-score数据，可以按以下步骤来操作：

1. 先补全你的热图代码（确保开启行标准化） #

首先得确认你在画热图时开启了行z-score标准化，也就是加上scale = "row"参数，不然没法获取行z-score值。补全后的代码大概是这样：

library(pheatmap)
library(RColorBrewer)

# 你的原始数据框
dat <- structure(c(36.28516213, -9.184523335, 5.716034994, 32.82180488, -18.5467127, 22.79930301, -24.79075287, -49.23768715), .Dim = c(2L, 4L), .Dimnames = list(c("motif_a", "motif_b"), c("cow", "dog", "snake", "cat")))

# 生成带行z-score标准化的热图
p <- pheatmap(
  dat,
  color = brewer.pal(10, "RdBu"),  # 用RColorBrewer的颜色示例
  scale = "row"  # 关键：开启行z-score标准化
)

2. 获取行z-score标准化后的数值 #

pheatmap返回的对象里，scaled_matrix字段直接存储了行标准化后的z-score矩阵，你可以直接调用：

# 提取行z-score数据
row_zscores <- p$scaled_matrix
print(row_zscores)

运行后就能得到每个单元格对应的行z-score值，和热图里的颜色一一对应。

3. 获取颜色刻度的数值（色标对应的取值范围） #

颜色刻度的分割点（也就是色标上的刻度值）存在breaks字段里，直接提取即可：

# 提取颜色刻度的分割点
color_scale_breaks <- p$breaks
print(color_scale_breaks)

这些数值就是颜色映射的边界，比如你用了10种颜色，就会得到11个breaks值，每个值对应色标上的一个刻度点。

如果想要直接获取色标上显示的文本标签，可以从热图的gtable对象里提取：

# 提取色标上的文本标签
legend_text_labels <- p$gtable$grobs[[which(sapply(p$gtable$grobs, function(x) x$name) == "legend")]]$children$labels$label
print(legend_text_labels)

内容的提问来源于stack exchange，提问作者littleworth