geom_tile绘图无法排序实现对角线布局问题求助

核心问题：染色体的正确排序 #

默认字符串排序会把chr10排在chr2前，导致热力图对角线错位。解决的关键是给染色体列设置自定义有序因子，而非单纯做字符串排序。

1. 处理染色体排序 #

假设你的bands数据包含chrom列（染色体标识），先生成正确的排序列表：

library(gtools)

# 提取并按自然顺序排序染色体
sorted_chroms <- unique(mixedsort(bands$chrom))

# 将chrom列转为有序因子
bands$chrom <- factor(bands$chrom, levels = sorted_chroms)

2. 生成染色体组合并过滤三角区域 #

用expand.grid生成两两组合，通过因子的索引值过滤出上/下三角：

# 生成所有染色体对组合
chrom_pairs <- expand.grid(
  chrom1 = sorted_chroms,
  chrom2 = sorted_chroms,
  stringsAsFactors = FALSE
)

# 转为有序因子，保持排序
chrom_pairs$chrom1 <- factor(chrom_pairs$chrom1, levels = sorted_chroms)
chrom_pairs$chrom2 <- factor(chrom_pairs$chrom2, levels = sorted_chroms)

# 替换为你的实际热力值数据
set.seed(123)
chrom_pairs$value <- rnorm(nrow(chrom_pairs))

# 过滤上三角（含对角线）：chrom1的层级索引 ≤ chrom2的索引
upper_tri <- chrom_pairs[as.integer(chrom_pairs$chrom1) <= as.integer(chrom_pairs$chrom2), ]

# 过滤下三角（含对角线）则用 ≥
lower_tri <- chrom_pairs[as.integer(chrom_pairs$chrom1) >= as.integer(chrom_pairs$chrom2), ]

3. 绘制三角热力图 #

因为已经设置了因子顺序，ggplot会自动按正确的染色体顺序排列：

library(ggplot2)

# 上三角热力图示例
ggplot(upper_tri, aes(x = chrom1, y = chrom2, fill = value)) +
  geom_tile(color = "white") +
  scale_x_discrete(limits = sorted_chroms) +
  scale_y_discrete(limits = sorted_chroms) +
  scale_fill_gradient2(low = "blue", mid = "white", high = "red") +
  theme_minimal() +
  theme(axis.text.x = element_text(angle = 45, hjust = 1))

额外提示 #

如果你的bands是染色体上的区间数据而非整条染色体，需先按chrom+start排序后，再生成区间的两两组合，后续过滤三角区域的逻辑不变。

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