如何用ggplot+facet_wrap绘制含Starch参考线的多面板折线图？

1. 数据预处理 #

要让每个面板都显示Starch参考线，核心是把Starch的观测数据和每个目标底物（Arabinan、Glucomannan、Xylan、Xyloglucan）做关联，确保每个面板的数据集里都包含Starch的行。

假设你的数据集名为df，先加载必要的包，再处理数据：

library(ggplot2)
library(dplyr)
library(tidyr)

# 定义需要单独做面板的目标底物
target_substrates <- c("Arabinan", "Glucomannan", "Xylan", "Xyloglucan")

# 提取Starch的所有观测
starch_data <- df %>% filter(Substrate == "Starch")

# 组合数据：把目标底物的观测和对应分组绑定，再加入每个分组对应的Starch数据
plot_data <- df %>%
  filter(Substrate %in% target_substrates) %>%
  mutate(facet_group = Substrate) %>%
  bind_rows(
    starch_data %>%
      crossing(facet_group = target_substrates)
  )

2. 绘制带参考线的分面板图 #

用ggplot结合facet_wrap，通过颜色和线型区分目标底物与Starch参考线，让参考线更醒目：

ggplot(plot_data, aes(x = Time, y = OD, color = Substrate, linetype = Substrate)) +
  geom_line(linewidth = 1) +
  # 自定义颜色：Starch用黑色，目标底物用不同彩色
  scale_color_manual(values = c(
    "Starch" = "black",
    "Arabinan" = "#E69F00",
    "Glucomannan" = "#56B4E9",
    "Xylan" = "#009E73",
    "Xyloglucan" = "#F0E442"
  )) +
  # 自定义线型：Starch用虚线，目标底物用实线
  scale_linetype_manual(values = c(
    "Starch" = "dashed",
    "Arabinan" = "solid",
    "Glucomannan" = "solid",
    "Xylan" = "solid",
    "Xyloglucan" = "solid"
  )) +
  # 按分组生成4个面板
  facet_wrap(~facet_group) +
  # 设置坐标轴和图例标签
  labs(x = "时间", y = "OD值", color = "底物", linetype = "底物") +
  # 用简洁的主题
  theme_bw()

关键说明 #

• 数据处理时，crossing函数会让Starch数据和每个目标底物分组一一配对，确保每个面板都能加载到Starch的折线数据。
• 通过color和linetype的双重映射，能清晰区分参考线与当前面板的目标底物曲线，避免混淆。

内容的提问来源于stack exchange，提问作者aetobatus-narinari