基于R语言ggplot2实现多SPI指标分面绘图的技术求助

问题背景 #

我拥有包含year、SPI 3、SPI 6、SPI 9、SPI 12的数据集，目标是将4个SPI指标绘成分面可视化图以同步分析，但初始代码生成的图形不符合预期。

样本数据 #

# 示例数据集结构
sample_data <- data.frame(
  year = 2000:2020,
  `SPI 3` = rnorm(21),
  `SPI 6` = rnorm(21),
  `SPI 9` = rnorm(21),
  `SPI 12` = rnorm(21)
)

初始代码（存在问题） #

library(ggplot2)
ggplot(sample_data, aes(x = year)) +
  geom_line(aes(y = `SPI 3`)) +
  geom_line(aes(y = `SPI 6`)) +
  geom_line(aes(y = `SPI 9`)) +
  geom_line(aes(y = `SPI 12`)) +
  facet_wrap(~ year) # 分面逻辑错误，未按SPI指标分类

当前效果图 #

问题表现：所有SPI指标叠加在同一面板，或分面按年份拆分，无法实现4个指标独立对比的需求

目标效果图 #

预期效果：4个SPI指标分别作为独立面板，按整齐布局排列，每个面板仅展示对应指标随年份的变化趋势

修正代码（感谢@r2evans 协助） #

library(tidyr)
library(ggplot2)

# 将宽格式数据转换为长格式
long_spi_data <- pivot_longer(
  sample_data,
  cols = starts_with("SPI"), # 匹配所有SPI开头的列
  names_to = "SPI_Scale",    # 新列名：SPI时间尺度
  values_to = "SPI_Value"    # 新列名：SPI对应数值
)

# 绘制分面图
ggplot(long_spi_data, aes(x = year, y = SPI_Value)) +
  geom_line(color = "#2c3e50", linewidth = 0.8) +
  facet_wrap(~ SPI_Scale, ncol = 2) + # 按SPI尺度分面，设置2列布局
  labs(
    x = "年份",
    y = "SPI值",
    title = "不同时间尺度SPI年度变化趋势"
  ) +
  theme_minimal()

最终效果图 #

修正后效果：4个SPI指标各自独立成面板，布局规范，可直观同步对比不同时间尺度的SPI变化趋势

核心技术要点 #

• 数据格式适配：ggplot2的分面功能依赖长格式数据，宽格式（每个指标为一列）无法直接按指标分面，需用tidyr::pivot_longer()完成格式转换
• 分面变量设置：分面参数facet_wrap()的变量需对应转换后的分类列（SPI_Scale），而非原数据的年份列
• 布局控制：通过ncol/nrow参数可调整分面的行列数，优化可视化布局

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