R语言pheatmap基因显示不全及样本注释报错求助

问题1：heatmap函数显示全部差异基因的设置方法 #

当使用heatmap()绘制426个差异基因的热图时，标签被截断或仅显示部分基因，核心原因是默认参数下标签字体过大、图形尺寸不足，导致重叠或被截断。可通过以下步骤解决：

1. 确保矩阵仅包含目标差异基因
   先从标准化计数矩阵中提取筛选出的426个差异基因，再完成z-score转换：

# 假设sigs数据框包含差异基因的gene_id列，norm_counts是标准化计数矩阵
diff_gene_z <- t(scale(t(norm_counts[sigs$gene_id, ])))

2. 调整图形设备尺寸与热图参数
   打开足够大的图形设备，同时缩小行标签字体、调整边距，强制显示所有基因名：

# 生成png格式热图，设置足够宽高
png("diff_gene_heatmap.png", width = 1200, height = 1800)
heatmap(diff_gene_z,
        cexRow = 0.4,  # 缩小行标签字体，可根据基因数量调整（比如0.3-0.6）
        labRow = rownames(diff_gene_z),  # 强制显示所有基因名
        margins = c(10, 16),  # 增大左右边距，避免标签被截断
        col = colorRampPalette(c("blue", "white", "red"))(100))  # 可选：自定义颜色梯度
dev.off()

• 若仍有部分标签不显示，可进一步减小cexRow的值，或增加height参数；
• 若使用RStudio绘图，可手动放大绘图窗口后再导出高清图片。

问题2：pheatmap添加样本条件注释的正确方法 #

报错Error in check.length("fill") : 'gpar' element 'fill' must not be length 0的核心原因是注释数据框的行名与热图矩阵行名不匹配，或注释列无有效颜色映射。以下是正确实现步骤：

1. 预处理注释数据框
   确保targets数据框的行名与样本距离矩阵的行名（样本ID）完全一致：

# 设置targets的行名为样本ID（Sample列的值）
rownames(targets) <- targets$Sample
# 提取仅包含Condition的注释数据框
sample_annot <- targets[, "Condition", drop = FALSE]

2. 验证行名匹配
   确保样本距离矩阵的行名与注释数据框行名完全一致，避免匹配失败：

# 假设sample_dist_matrix是样本距离矩阵（由dist(t(norm_counts))转换而来）
stopifnot(all(rownames(sample_dist_matrix) == rownames(sample_annot)))

3. 绘制带注释的样本距离热图
   可选自定义注释颜色，确保每个Condition有对应颜色，避免默认映射出错：

# 定义注释颜色映射（根据你的Condition类别调整）
annot_colors <- list(
  Condition = c(
    Control = "#1f78b4",
    Treatment = "#e31a1c"
  )
)

# 绘制热图（样本距离矩阵需转为矩阵格式）
pheatmap(as.matrix(sample_dist_matrix),
         annotation_row = sample_annot,
         annotation_col = sample_annot,  # 行和列都是样本，可同时添加列注释
         annotation_colors = annot_colors,
         main = "Sample Distance Heatmap",
         treeheight_row = 0, treeheight_col = 0)  # 可选：关闭聚类树（样本距离热图常不需要）

• 如果不自定义annotation_colors，pheatmap会自动生成颜色，但需确保Condition列是因子或字符型，无空值；
• 若仍报错，检查sample_dist_matrix的行名是否有重复或拼写错误，以及sample_annot的Condition列是否存在空值。

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