ggplot中Facet Wrap绘图异常：多日期SPEI时间序列数据问题

问题根源 #

你的问题大概率出在数据重塑逻辑错误或NA值未处理：

• 不同SPEI起始日期不同，重塑后会存在大量NA值，ggplot绘制面积图时NA会导致断层、空白；
• 若用melt时未拆分站点与SPEI类型，分面变量variable可能只对应SPEI类型，导致所有站点数据混在一起，可视化混乱。

具体解决步骤 #

1. 正确重塑宽格式数据（推荐用pivot_longer替代melt）

假设你的原始CSV是宽格式，列结构为date, site1_SPEI3, site1_SPEI6, site1_SPEI12, site2_SPEI3,...，先拆分站点与SPEI类型，并过滤NA值：

library(tidyverse)

# 读取数据并转换日期格式
df <- read.csv("your_data.csv", stringsAsFactors = FALSE)
df$date <- as.Date(df$date)  # 确保日期为Date类型

# 重塑为长格式，拆分站点和SPEI类型
long_df <- df %>%
  pivot_longer(
    cols = -date,  # 排除日期列，其余列转为长格式
    names_to = c("site", "spei_type"), 
    names_sep = "_",  # 按下划线拆分列名（需匹配你的实际列名分隔符）
    values_to = "spei_value"
  ) %>%
  drop_na(spei_value)  # 移除无有效SPEI值的行（解决不同起始日期导致的NA问题）

2. 绘制正确的分面面积图

根据你的需求选择分面方式：

方式1：每个站点单独分面，展示3种SPEI的堆叠/叠加面积图

ggplot(long_df, aes(x = date, y = spei_value, fill = spei_type)) +
  geom_area(alpha = 0.7, position = "identity") +  # identity为叠加，stack为堆叠
  facet_wrap(~site, ncol = 2) +  # 按站点分面，2列排列
  labs(x = "日期", y = "SPEI值", fill = "SPEI类型") +
  theme_bw() +
  scale_x_date(date_labels = "%Y", date_breaks = "2 years")  # 自定义日期轴显示

方式2：每个站点的每种SPEI单独分面

ggplot(long_df, aes(x = date, y = spei_value)) +
  geom_area(fill = "steelblue", alpha = 0.7) +
  facet_wrap(site ~ spei_type, ncol = 3) +  # 行=站点，列=SPEI类型
  labs(x = "日期", y = "SPEI值") +
  theme(axis.text.x = element_text(angle = 45, hjust = 1))  # 旋转x轴标签避免重叠

3. 异常排查

• 若分面后数据混乱：检查site和spei_type列的取值是否正确，比如列名拆分是否匹配你的实际命名规则（如果列名是SPEI3_site1，则names_sep和names_to顺序要调换）；
• 若面积图有断层：确认已用drop_na()移除NA值，或在geom_area中添加na.rm = TRUE；
• 若日期轴显示异常：确保date列是Date类型，可通过class(long_df$date)验证。

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