在R语言中如何可视化大矩阵中的10×10对角分组结构？

哦，这个问题我太熟了！原生的heatmap()函数之所以搞不定，核心原因是它默认会对行和列做层次聚类，直接把你精心安排的分组顺序打乱了。要画出你想要的10×10对角块热图，关键就是保留原行列顺序+可视化分组边界，给你几个实用的解决方案：

方法一：用pheatmap包（最推荐，简洁直观） #

pheatmap是专门用于热图绘制的工具，支持自定义间隙和分组注释，完美适配你的需求。

步骤： #

1. 先构造一个模拟矩阵（方便你测试，替换成自己的矩阵即可）：

set.seed(123) # 保证结果可重复
mat <- matrix(0, nrow = 1000, ncol = 1000)
# 给每个对角块填充较高的值，模拟真实分组的信号
for (i in 1:100) {
  start <- (i - 1) * 10 + 1
  end <- i * 10
  mat[start:end, start:end] <- rnorm(100, mean = 5, sd = 1)
}

2. 安装并加载pheatmap，然后绘制带间隙的热图：

install.packages("pheatmap") # 首次使用需要安装
library(pheatmap)

# 定义分组间隙：在每10行/列的末尾插入间隙
gaps <- seq(10, 990, by = 10)

pheatmap(mat,
         cluster_rows = FALSE, # 关闭行聚类，保留原顺序
         cluster_cols = FALSE, # 关闭列聚类，保留原顺序
         gaps_row = gaps,      # 行之间的间隙位置
         gaps_col = gaps,      # 列之间的间隙位置
         border_color = "white", # 块内单元格边框设为白色，突出间隙
         show_rownames = FALSE,  # 1000行太多，隐藏行名
         show_colnames = FALSE,  # 隐藏列名
         main = "10×10对角块分组热图" # 添加标题
)

3. 进阶：添加分组注释（侧边/顶部的颜色条）
   如果想更清晰地看到每个分组的对应关系，可以给行和列添加分组标签：

# 创建分组注释数据框
group_labels <- rep(1:100, each = 10)
row_annotation <- data.frame(Group = factor(group_labels))
col_annotation <- data.frame(Group = factor(group_labels))

pheatmap(mat,
         cluster_rows = FALSE,
         cluster_cols = FALSE,
         gaps_row = gaps,
         gaps_col = gaps,
         annotation_row = row_annotation, # 行分组注释
         annotation_col = col_annotation, # 列分组注释
         show_rownames = FALSE,
         show_colnames = FALSE
)

方法二：用ggplot2绘制（灵活定制） #

如果你习惯用ggplot2的语法，也可以通过转换矩阵格式+添加分隔线来实现：

install.packages(c("ggplot2", "reshape2"))
library(ggplot2)
library(reshape2)

# 将宽格式矩阵转换为长格式
mat_long <- melt(mat)

# 给每个单元格添加行/列分组标签
mat_long$row_group <- rep(1:100, each = 10 * 1000)
mat_long$col_group <- rep(rep(1:100, each = 10), 1000)

# 绘制热图并添加分组分隔线
ggplot(mat_long, aes(x = Var2, y = Var1, fill = value)) +
  geom_tile() +
  scale_fill_viridis_c(option = "plasma") + # 用更美观的颜色映射
  # 添加白色分隔线区分分组
  geom_vline(xintercept = seq(10, 990, by = 10) + 0.5, color = "white", linewidth = 0.5) +
  geom_hline(yintercept = seq(10, 990, by = 10) + 0.5, color = "white", linewidth = 0.5) +
  theme_minimal() +
  theme(
    axis.text = element_blank(), # 隐藏轴标签
    axis.title = element_blank(),
    panel.grid = element_blank(), # 隐藏网格线
    plot.title = element_text(hjust = 0.5)
  ) +
  labs(title = "分组矩阵对角块热图", fill = "数值")

关键注意事项 #

• 无论用哪种方法，必须关闭聚类或保留原行列顺序，否则分组结构会被聚类算法打乱
• 间隙/分隔线是可视化分组的核心，能让10×10的对角块一目了然
• 如果你的矩阵数值范围差异大，可以先做标准化（比如scale(mat)），让热图颜色更有区分度

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