如何在循环绘制的ggplot图中添加对应细胞因子的lloq/uloq水平线？

问题背景 #

已有细胞因子数据集，通过for循环绘制13种不同analyte的ggplot散点+均值图，需要为每个analyte的图添加对应其自身的最低定量限（lloq）和最高定量限（uloq）水平线，但尝试代码时要么报错Error in check_required_aesthetics():! geom_hline requires the following missing aesthetics: yintercept，要么会把所有analyte的lloq/uloq都画到同一张图里。

错误原因 #

之前的代码中geom_hline调用方式存在两个问题：

1. 未指定yintercept参数，导致缺失必要的美学映射；
2. 未对loq数据集进行子集化，直接使用完整数据集导致所有analyte的水平线被同时绘制。

正确解决方案 #

在循环中，先为当前analyte筛选对应的loq数据，再分别添加lloq和uloq水平线，同时可设置不同样式区分两条线：

# 获取唯一的analyte列表
analyte_list <- unique(cytokines$analyte)

for (i in analyte_list) {
  # 筛选当前analyte的细胞因子数据
  current_cyto <- subset(cytokines, analyte == i)
  # 筛选当前analyte的loq/uloq数据
  current_loq <- subset(loq, analyte == i)
  
  plot <- ggplot(data = current_cyto) +
    aes(x = treat, y = result, color = groep) +
    geom_point(shape = 1, position = position_jitterdodge(dodge.width = 0.5, jitter.width = 0.1)) +
    stat_summary(fun = mean, geom = "crossbar", width = 0.3, 
                 mapping = aes(group = groep), position=position_dodge(0.5)) +
    # 添加LLOQ水平线，红色虚线
    geom_hline(yintercept = current_loq$lloq, color = "red", linetype = "dashed") +
    # 添加ULOQ水平线，蓝色虚线
    geom_hline(yintercept = current_loq$uloq, color = "blue", linetype = "dashed") +
    scale_y_log10() +
    theme(legend.position="bottom") + 
    ggtitle(i) +
    labs(y = "pg/ml", x = "Time")  
  
  ggsave(plot, file=paste0("plot_", i, ".png"))
}

代码说明 #

• 提前筛选当前analyte对应的loq数据，避免重复子集化操作；
• 直接用yintercept = current_loq$lloq指定水平线位置，无需在aes()中映射（每个图仅需单条lloq/uloq线）；
• 为两条线设置不同颜色和线型，便于区分LLOQ与ULOQ。

优化方案（无需循环） #

如果不想用循环，可先合并cytokines和loq数据集，通过facet_wrap一次性生成所有图：

# 按analyte合并两个数据集
merged_data <- merge(cytokines, loq, by = "analyte")

ggplot(merged_data) +
  aes(x = treat, y = result, color = groep) +
  geom_point(shape = 1, position = position_jitterdodge(dodge.width = 0.5, jitter.width = 0.1)) +
  stat_summary(fun = mean, geom = "crossbar", width = 0.3, 
               mapping = aes(group = groep), position=position_dodge(0.5)) +
  geom_hline(aes(yintercept = lloq), color = "red", linetype = "dashed") +
  geom_hline(aes(yintercept = uloq), color = "blue", linetype = "dashed") +
  scale_y_log10() +
  theme(legend.position="bottom") + 
  labs(y = "pg/ml", x = "Time") +
  facet_wrap(~analyte, scales = "free_y") # 每个图的y轴独立缩放

ggsave("all_cytokine_plots.png", width = 12, height = 8)

这种方式更简洁，自动为每个analyte匹配对应的水平线，scales = "free_y"还能让每个图的y轴范围适配自身数据。

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