如何为pheatmap绘制的细菌频率热图添加样本分类标签？

• 第一步：对齐样本名并预处理分类数据
  确保表格2（样本类型表）的行名和热图数据data.for.heatmap2的行名完全匹配，同时将类型列转为因子类型（让分类逻辑更清晰）：

  # 假设你的样本类型表名为sample_type，列名为type
  rownames(sample_type) <- rownames(data.for.heatmap2)
  sample_type$type <- factor(sample_type$type)
  

• 第二步：修改pheatmap代码添加注释
  在原有代码基础上，加入annotation_row参数绑定分类数据，还可以通过annotation_colors自定义分类配色：

  library(pheatmap)
  library(RColorBrewer)
  
  # 自定义分类颜色（可选，不设置则用默认配色）
  anno_colors <- list(
    type = c(type1 = "#E41A1C", type2 = "#377EB8", type3 = "#4DAF4A")
  )
  
  heatmap <- pheatmap(data.for.heatmap2,
                      color = brewer.pal(9, "Blues"),
                      show_rownames = F,
                      cluster_cols = F,
                      annotation_row = sample_type,  # 绑定行（样本）的分类注释
                      annotation_colors = anno_colors)  # 应用自定义颜色
  

• 关键注意事项

  • 样本名必须严格匹配，否则注释会出现错位
  • 如果需要将分类注释放在热图顶部（列方向），替换为annotation_col参数，但这里样本对应热图的行，所以用annotation_row
  • annotation_colors是可选参数，省略后pheatmap会自动分配分类颜色

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