如何在R的ggplot2棒棒糖图中为每条竖线底部添加SNP标签

问题需求 #

在ggplot2绘制的棒棒糖图中，需为每条对应SNP的竖线**底部（图外正下方）**添加唯一的SNP编号标签，避免重复，且标签需精准对应每条竖线的正下方（类似示例图前两条线黑框位置）。

解决方案思路 #

1. 去重标签数据：原始数据中每个SNP对应两个表型（Diabetes/Cardiac），直接添加文本会出现重复，需先提取唯一的SNP与BP_2组合，确保每条竖线仅对应一个标签。
2. 扩展绘图边界：通过调整y轴范围，为图外的标签预留显示空间。
3. 精准定位标签：将标签放在每条竖线的正下方（y轴范围外），使用geom_text而非geom_text_repel避免位置偏移，保证对齐。

完整代码实现 #

1. 导入并预处理数据 #

# 导入测试数据
data <- read.table(text = "SNP        BP_2         PP4    Phenotype   Present 
rs5478100  117.558703    0.5    Diabetes    1
rs6372078  117.563687     0    Diabetes     
rs5866152  117.564875   0.06   Diabetes     1
rs2263839  117.56644    0.05   Diabetes     1
rs6499109  117.568766   0.06   Diabetes     1
rs7863647  117.569046   0.55   Diabetes     1
rs12787678 117.569046   0.2    Diabetes     1
rs5478100  117.558703   0      Cardiac
rs6372078  117.563687   0.4    Cardiac      1
rs5866152  117.56487    0.2    Cardiac      1
rs2263839  117.56644    0.04   Cardiac      1
rs6499109  117.568766   0.19   Cardiac      1
rs7863647  117.569046   0.11   Cardiac      1
rs12787678 117.579457   0.37   Cardiac      1", header = TRUE, stringsAsFactors = FALSE)

# 提取唯一的SNP-BP_2组合，用于生成标签
label_data <- unique(data[, c("SNP", "BP_2")])

2. 绘制棒棒糖图并添加标签 #

library(ggplot2)

dot_plot <- ggplot(data, aes(BP_2, Phenotype)) +
  # 绘制对应SNP的竖线
  geom_line(aes(group = BP_2)) +
  # 绘制表型对应的数据点
  geom_point(aes(size = Present, colour = Phenotype)) +
  # 设置x轴刻度在顶部
  scale_x_continuous(breaks = seq(117.558000, 117.580000, by = 0.002000), position = 'top') +
  # 扩展y轴下边界，为图外标签预留空间
  scale_y_discrete(expand = expansion(add = c(1.2, 0))) +
  # 主题样式调整
  theme_bw() + 
  theme(panel.border = element_blank(), 
        panel.grid.major = element_blank(),
        panel.grid.minor = element_blank(), 
        axis.line.x = element_line(colour = "black"),
        axis.line.y = element_line(color = "black"),
        legend.position = "none") +
  labs(title="", x="", y = "") +
  # 添加SNP标签，定位到竖线正下方的图外位置
  geom_text(data = label_data, aes(x = BP_2, y = -0.3, label = SNP), 
            size = 3, color = "black", hjust = 0.5)

dot_plot

关键说明 #

• 去重标签数据：unique(data[, c("SNP", "BP_2")])确保每个SNP仅生成一个标签，彻底解决重复问题。
• y轴扩展：scale_y_discrete(expand = expansion(add = c(1.2, 0)))为图下方预留足够空间，保证标签能正常显示。
• 标签定位：y = -0.3将标签设置在最下方表型（Cardiac）的外侧，hjust = 0.5保证标签与对应竖线居中对齐。

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