如何在ggplot中仅标记同时满足GSEA与PSEA阈值的点？

在R语言的ggplot绘图场景中，需要仅标记同时满足以下两个阈值条件的点：

• GSEA（-LOG10 FDR q-value）>0.6
• PSEA（-LOG10 Kuiper q-value）>1.0

当前代码中使用了|（逻辑或）运算符，导致满足任一条件的点都被标记，不符合需求；同时还需要优化标签位置（位于点左侧）、解决标签重叠和绘图边缘截断的问题。当前完整代码如下：

ggplot(df, aes(x=PSEALOG, y=GSEALOG, color=GSEAFDRqvalue)) +
  geom_point()+
  scale_colour_gradient(low = "blue", high = "red", na.value = NA) +
  theme(text=element_text(size=10)) +
  xlab("PSEA (-LOG10 Kuiper q-value)") + ylab("GSEA (-LOG10 FDR q-value)" ) +
  geom_text(data=subset(df, GSEALOG > 0.6 | PSEALOG > 1.0),
            aes(PSEALOG, GSEALOG, label=NAME))

1. 修正双阈值筛选逻辑 #

将subset函数中的逻辑运算符从|（逻辑或）改为&（逻辑与），即可仅筛选出同时满足两个阈值的点：

geom_text(data=subset(df, GSEALOG > 0.6 & PSEALOG > 1.0),
          aes(PSEALOG, GSEALOG, label=NAME))

2. 优化标签显示效果 #

让标签位于点左侧 #

在geom_text中添加hjust=1（文本右对齐，相对点的位置就会靠左），搭配nudge_x参数微调偏移距离，让标签和点保持合适间距：

geom_text(data=subset(df, GSEALOG > 0.6 & PSEALOG > 1.0),
          aes(PSEALOG, GSEALOG, label=NAME),
          hjust=1, nudge_x=-0.1)

解决标签重叠问题 #

推荐使用ggrepel包的geom_text_repel替代原生geom_text，它会自动调整标签位置避免重叠：

1. 先安装并加载包：

install.packages("ggrepel")
library(ggrepel)

2. 替换为geom_text_repel并配置参数：

geom_text_repel(data=subset(df, GSEALOG > 0.6 & PSEALOG > 1.0),
                aes(PSEALOG, GSEALOG, label=NAME),
                hjust=1,
                nudge_x=-0.1,
                force=0.5) # 可根据实际情况调整标签排斥力度

避免标签被边缘截断 #

可以通过扩展绘图范围解决，有两种方式：

• 方式一：在geom_text_repel中设置xlim和ylim参数，手动扩展轴范围：

geom_text_repel(data=subset(df, GSEALOG > 0.6 & PSEALOG > 1.0),
                aes(PSEALOG, GSEALOG, label=NAME),
                hjust=1,
                nudge_x=-0.1,
                xlim=c(min(df$PSEALOG)-0.5, max(df$PSEALOG)),
                ylim=c(min(df$GSEALOG)-0.2, max(df$GSEALOG)))

• 方式二：在ggplot主图层中添加expand_limits扩展整体绘图范围：

ggplot(df, aes(x=PSEALOG, y=GSEALOG, color=GSEAFDRqvalue)) +
  geom_point()+
  scale_colour_gradient(low = "blue", high = "red", na.value = NA) +
  theme(text=element_text(size=10)) +
  xlab("PSEA (-LOG10 Kuiper q-value)") + ylab("GSEA (-LOG10 FDR q-value)" ) +
  expand_limits(x = c(min(df$PSEALOG)-0.5, max(df$PSEALOG)),
                y = c(min(df$GSEALOG)-0.2, max(df$GSEALOG))) +
  geom_text_repel(data=subset(df, GSEALOG > 0.6 & PSEALOG > 1.0),
                  aes(PSEALOG, GSEALOG, label=NAME),
                  hjust=1, nudge_x=-0.1)

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