如何用ggplot2的facet_wrap实现独立变量的多面板IMT关联图？

需求背景 #

需要用ggplot2绘制以年龄（H2_lft）为X轴、平均IMT值（MeanIMT_alg）为Y轴的图形，分析吸烟（H2_Roken）、性别（H1_geslacht）、种族（H1_EtnTotaal）这三个心血管风险因素与IMT随年龄变化的关联。原本将所有变量的拟合线放在同一图中过于杂乱，希望生成三个独立面板：每个面板对应一个风险因素的分组拟合线（比如一个面板展示种族分组的拟合线，一个展示吸烟状态的拟合线，一个展示性别的拟合线），而非变量交叉组合的面板，但现有代码生成了交叉组合的面板，无法满足需求。

核心解决方案：转换为长格式数据 #

ggplot2对长格式数据的支持更好，需要先将宽格式数据转换为长格式，把三个独立的分组变量合并为统一的“风险因素”和“分组类别”列，再基于长数据绘图。

完整代码示例 #

首先加载所需包（未安装则先执行安装步骤）：

# 首次运行安装所需包
install.packages(c("ggplot2", "tidyr"))

# 加载包
library(ggplot2)
library(tidyr)

转换数据格式并绘图：

# 将宽数据转换为长格式
Dataset_long <- Dataset %>%
  pivot_longer(
    cols = c(H1_EtnTotaal, H2_Roken, H1_geslacht), # 指定要转换的三个分组变量
    names_to = "risk_factor", # 新列：存储风险因素名称（种族/吸烟/性别）
    values_to = "group_level" # 新列：存储对应分组的具体类别
  )

# 绘制多面板图
t <- ggplot(data = Dataset_long, aes(x = H2_lft, y = MeanIMT_alg)) +
  # 按分组类别绘制拟合线
  geom_smooth(se = FALSE, aes(group = group_level, colour = group_level)) +
  # 添加全数据的基准拟合线（黑色虚线）
  stat_smooth(method = lm, se = FALSE, colour = "black", linetype = "dashed") +
  # 按风险因素分面板，每个面板对应一个变量的分组
  facet_wrap(~risk_factor, scales = "free_y") +
  theme_classic() +
  labs(
    x = "Age (years)", 
    y = "Mean IMT (mm)", 
    title = "IMT vs Age by Cardiovascular Risk Factors",
    caption = "Figure 2: Mean IMT by Age and Risk Factors",
    color = "Group"
  )

# 显示图形
print(t)

代码说明 #

1. 数据转换：pivot_longer把原来的三个分组列合并成两列：risk_factor记录当前是哪个风险因素，group_level记录该因素下的具体分组值（比如种族的不同类别、吸烟的是/否等）。
2. 绘图逻辑：
   • geom_smooth通过group和colour映射group_level，实现每个分组单独绘制拟合线；
   • stat_smooth添加全数据的基准拟合线，方便对比各组与整体趋势的差异；
   • facet_wrap(~risk_factor)直接按风险因素拆分面板，每个面板只展示对应变量的分组拟合线，不会产生交叉组合；
   • scales = "free_y"允许每个面板的Y轴范围独立，更贴合各组数据的分布。

原代码错误原因 #

之前的facet_wrap(~H1_EtnTotaal + ~H2_Roken + ~H1_geslacht)是将三个变量的所有交叉组合作为面板（比如“种族A+吸烟+男性”会成为一个独立面板），这会生成大量冗余的交叉面板，完全不符合需求。

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