R语言混合传统绘图与ggplot2添加曼哈顿图标签失败求助

你遇到的核心问题是基础R绘图系统和ggplot2是完全独立的两套框架——你用plot()和points()画的是基础图形设备上的图，而geom_label_repel()是ggplot2专属的图层函数，根本无法直接叠加到基础图上，这就是标签完全不显示的原因。

下面给你两种可行的解决方案，优先推荐转用ggplot2，这是适配ggrepel的最优方式：

方案1：完全转用ggplot2绘制（推荐） #

我们把你原来的基础绘图逻辑重写成ggplot2代码，同时集成ggrepel的自动标签功能，完美匹配你的需求：

1. 数据预处理（和你的原始逻辑对齐） #

先处理染色体的位置偏移，让不同染色体在x轴上连续排列：

library(ggplot2)
library(ggrepel)

# 读取主数据
d <- read.table("a.txt", header=T, fill=T, sep=" ")
dmht <- data.frame(
  chrom = d[,9],
  txStart = d[,11],
  neg_log10_p = -log10(d[,5])
)
# 按染色体和位置排序
dmht <- dmht[order(dmht$chrom, dmht$txStart), ]

# 读取需要标注的基因数据
ds <- read.table("selected2label.txt", header=T, fill=T, sep=" ")
dsel <- data.frame(
  chrom = ds[,9],
  txStart = ds[,11],
  neg_log10_p = -log10(ds[,5]),
  GENE_ID = ds[,1]
)

# 定义染色体顺序和显示标签
chrs <- c('chr1','chr2','chr3','chr4','chr5','chr6','chr7','chr8','chr9','chr10','chr11','chr12','chr13','chr14','chr15','chr16','chr17','chr18','chr19','chrX','chrY')
chrLabels <- c('1', '2', '3', '4', '5', '6', '7', '8', '9', '10', '11', '12', '13', '14', '15', '16', '17', '18', '19', 'X', 'Y')

# 计算每个染色体的起始偏移（让x轴连续无重叠）
chr_len <- sapply(chrs, function(chr) {
  max(dmht$txStart[dmht$chrom == chr]) - min(dmht$txStart[dmht$chrom == chr]) + 1
})
chr_start <- c(0, cumsum(chr_len)[-length(chr_len)])
names(chr_start) <- chrs

# 给主数据和标注数据添加x轴的实际位置
dmht$x_pos <- dmht$txStart + chr_start[dmht$chrom]
dsel$x_pos <- dsel$txStart + chr_start[dsel$chrom]

# 计算x轴刻度位置（每个染色体的中点）
chr_mid <- chr_start + chr_len/2

2. 绘制曼哈顿图+自动防重叠标签 #

# 设置染色体交替颜色
chr_colors <- rep(c("blue", "green", "cyan"), length.out = length(chrs))
names(chr_colors) <- chrs

png("result/test/mhtTest.png", width=1600)
ggplot(dmht, aes(x = x_pos, y = neg_log10_p)) +
  # 绘制每个染色体的散点
  geom_point(aes(color = chrom), size = 1) +
  # 绑定染色体颜色与标签
  scale_color_manual(values = chr_colors, breaks = chrs, labels = chrLabels) +
  # 添加自动避免重叠的基因标签
  geom_label_repel(
    data = dsel,
    aes(label = GENE_ID, x = x_pos, y = neg_log10_p),
    size = 5,
    box.padding = unit(0.35, "lines"),
    point.padding = unit(0.5, "lines"),
    show.legend = FALSE
  ) +
  # 自定义x轴为染色体标签（替代原始的连续位置）
  scale_x_continuous(
    breaks = chr_mid,
    labels = chrLabels,
    expand = c(0, 0)
  ) +
  # 设置轴标签和标题
  labs(x = "Chromosome", y = "-log10(PValue)", title = "") +
  # 调整主题，让x轴标签垂直（对应你原来的las=3）
  theme_bw() +
  theme(
    axis.text.x = element_text(angle = 90, vjust = 0.5, hjust = 1),
    panel.grid.major.x = element_blank()
  )
dev.off()

方案2：坚持用基础绘图（不推荐） #

如果实在不想切换到ggplot2，你可以用基础绘图兼容的标签防重叠工具，比如textplot包的textrepel()函数，替换你原来的text()调用：

# 先安装依赖包（如果未安装）
install.packages("textplot")
library(textplot)

# 在你原来的标注循环里，替换text()为：
textrepel(
  x = as.numeric(dchr[,2])+chrStart[i],
  y = as.numeric(dchr[,3]),
  labels = dchr[,4],
  cex = 1.2,
  col = "red"
)

关键提醒 #

• 绝对不要混用基础绘图和ggplot2的函数，它们的底层渲染逻辑完全不同，无法互相叠加。
• ggplot2+ggrepel是绘制带自动标签曼哈顿图的行业标准方案，后续调整样式（比如添加阈值线、修改配色）也会更灵活。

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