如何在R的ggplot森林图中添加定位准确的标签与标题

一、Y轴"Analyte"标签精准对齐 #

默认ylab()很难让标签和Assay水平线完美对齐，换用annotation_custom()手动定位更靠谱：

1. 加载grid包，用textGrob()生成旋转后的标签
2. 隐藏默认y轴标题，将自定义标签放在Assay水平线正上方

library(grid)

# 替换原有ylab("Analyte")，添加这段到ggplot代码中
annotation_custom(textGrob("Analyte", rot = 90, gp = gpar(fontsize = 12)),
                  xmin = -0.15, xmax = -0.15, # 根据你的绘图范围调整x值，确保在y轴外侧
                  ymin = mean(unique(df$Assay)), ymax = mean(unique(df$Assay))) +
theme(axis.title.y = element_blank())

ymin/ymax用所有Assay值的均值，保证标签在y轴垂直中心位置，和所有分析物的水平线对齐。

二、面板上方"p-value"标签添加与对齐 #

给每个面板的QM_P_trend上方加标签，用数据驱动的geom_text()最灵活，不会错位：

1. 从原数据框生成每个面板的标签位置数据，定位在该面板p-value正上方
2. 将标签加入绘图

library(dplyr)

# 生成标签位置数据框，替换"分组变量"为你实际的分面列名
label_df <- df %>% 
  group_by(分组变量) %>% 
  summarize(
    x = max(QM_P_trend) + 0.1, # 调整x值到p-value右侧，避免重叠
    y = max(Assay) + 0.3 # 调整y值到该面板最高分析物的上方
  )

# 在ggplot代码中添加geom_text
geom_text(data = label_df, aes(x = x, y = y, label = "*p*-value"), 
          fontface = "italic", hjust = 0.5, vjust = 0)

如果你的p-value固定在面板某个位置，直接把x设为对应数值即可，不用取max。

三、完整修改示例（结合原有代码） #

假设你的原有代码框架如下：

library(ggplot2)
df <- read.csv("你的数据文件.csv")

ggplot(df, aes(x = OR, y = Assay, color = 分组变量)) +
  geom_point() +
  geom_errorbarh(aes(xmin = OR_low, xmax = OR_high)) +
  facet_wrap(~分组变量) +
  theme_bw() +
  ylab("Analyte")

修改后代码：

library(ggplot2)
library(grid)
library(dplyr)

df <- read.csv("你的数据文件.csv")

# 生成p-value标签位置数据
label_df <- df %>% 
  group_by(分组变量) %>% 
  summarize(
    x = max(QM_P_trend) + 0.1,
    y = max(Assay) + 0.3
  )

ggplot(df, aes(x = OR, y = Assay, color = 分组变量)) +
  geom_point() +
  geom_errorbarh(aes(xmin = OR_low, xmax = OR_high)) +
  # 添加p-value标签
  geom_text(data = label_df, aes(x = x, y = y, label = "*p*-value"), 
            fontface = "italic", hjust = 0.5, vjust = 0) +
  # 添加Analyte标签
  annotation_custom(textGrob("Analyte", rot = 90, gp = gpar(fontsize = 12)),
                    xmin = -0.15, xmax = -0.15,
                    ymin = mean(unique(df$Assay)), ymax = mean(unique(df$Assay))) +
  facet_wrap(~分组变量) +
  theme_bw() +
  # 隐藏默认y轴标题，调整刻度标签避免重叠
  theme(axis.title.y = element_blank(),
        axis.text.y = element_text(margin = margin(r = 10)))

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