ggplot2分面图实现：每个面板嵌入对照组参考图

没问题，我来帮你搞定这个需求！核心思路是让对照组数据不受分面变量的过滤，在每个分面面板中完整显示，下面是具体的实现步骤和代码：

关键思路 #

我们需要把对照组数据单独提取出来，然后在绘制患者组分面图的基础上，叠加对照组的图层——因为对照组数据的class是ctr，而分面变量是患者的class（a和b），所以每个面板都会自动加载完整的对照组数据作为参考。

完整代码实现 #

library(dplyr)
library(ggplot2)

# 你的示例数据
df <- data_frame(
 id = c(rep('ctr',40), rep('pat',80)),
 class = c(rep('ctr',40), rep(c('a','b'), each = 40)),
 rank = rep (1:20,6),
 mean = c(rep(seq(3,-3, length.out = 20),2), rep(seq(1,-4, length.out = 20),2), rep(seq(-2,-8, length.out = 20),2)),
 sd = rep(seq(1.2,0.8, length.out = 20), times = 6),
 exam = rep(c('blue','red'), each = 20, times = 3))

# 单独提取对照组数据，方便后续重复调用
df_ctr <- df %>% filter(id == 'ctr')

# 绘制分面图，每个class面板包含对照组参考
ggplot() +
  # 1. 绘制患者组的折线：按exam区分颜色，class区分线型
  geom_line(data = df %>% filter(id == 'pat'),
            aes(x = rank, y = mean, color = exam, linetype = class)) +
  # 2. 绘制患者组的置信带：用interaction(class, exam)确保分组正确
  geom_ribbon(data = df %>% filter(id == 'pat'),
              aes(x = rank, ymax = mean + sd, ymin = mean - sd,
                  fill = exam, group = interaction(class, exam)),
              alpha = 0.1) +
  # 3. 叠加对照组的参考折线：每个面板都会显示
  geom_line(data = df_ctr,
            aes(x = rank, y = mean, color = exam),
            linetype = 1, size = 1) +
  # 4. 叠加对照组的置信带：每个面板都会显示
  geom_ribbon(data = df_ctr,
              aes(x = rank, ymax = mean + sd, ymin = mean - sd, fill = exam),
              alpha = 0.1) +
  # 5. 设置分面：只针对患者的class（a和b），避免出现ctr面板
  facet_wrap(~class, data = df %>% filter(id == 'pat')) +
  # 6. 保持你原来的颜色和填充设置
  scale_colour_manual(values = c("#00b6eb","#eb0041")) +
  scale_fill_manual(values = c("#00b6eb","#eb0041")) +
  # 可选：优化主题
  theme_bw()

代码说明 #

• 单独提取对照组数据：df_ctr是独立于分面变量的数据集，所以在每个分面面板中都会完整渲染，不会被class过滤。
• 分面数据源指定：facet_wrap里的data = df %>% filter(id == 'pat')确保分面只生成a和b两个面板，不会出现对照组的ctr面板。
• 患者组置信带分组：用interaction(class, exam)作为分组变量，保证每个class+exam组合的置信带正确绘制，避免因为分面后class固定而出现分组错误。
• 对照组图层强化：给对照组折线加了size = 1，让参考线更突出，你可以根据需求调整。

如果需要同时按class和exam分面（比如每个面板是单个class+exam的组合，里面显示对应exam的对照组参考），只需要把分面代码改成facet_grid(class ~ exam)即可，其他部分不需要调整。

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