在R中绘制含大日期间隔的时间序列数据：消除间隙并加分隔虚线

问题核心 #

需要绘制包含2023年9月、12月，2024年4月、7月各6天的时间序列数据，要求：

• 视觉上不显示各时间段之间的日期间隙
• 用虚线分隔不同时间段
  当前代码因类型冲突（将日期转为因子后又用日期时间轴）导致效果不符合预期。

解决方案 #

步骤1：预处理数据，生成无间隙连续x轴变量 #

先给数据添加用于无间隙显示的连续位置变量，同时标记时间段分组：

library(dplyr)
library(lubridate)

# 预处理数据：生成全局连续x轴位置和时间段标记
AllSite_RedwoodResidual <- AllSite_RedwoodResidual %>%
  # 按年月划分时间段
  mutate(period = floor_date(Time, "month")) %>%
  # 每个时间段内生成相对顺序（连续数值）
  group_by(period) %>%
  mutate(relative_time = row_number()) %>%
  ungroup() %>%
  # 计算每个时间段的结束位置
  group_by(period) %>%
  mutate(period_end = max(relative_time)) %>%
  ungroup() %>%
  # 累加前面时间段的观测数，得到全局连续x轴位置
  mutate(global_x = relative_time + c(0, cumsum(diff(unique(period_end))))[match(period, unique(period))])

步骤2：绘制符合需求的图表 #

使用预处理后的global_x作为x轴，自动添加分隔虚线，调整x轴标签：

library(ggplot2)
library(scales)

ggplot(AllSite_RedwoodResidual, aes(x = global_x, y = AvgRedwood, color = ID)) +
  geom_line(aes(group = ID), na.rm = TRUE) +
  geom_point(aes(fill = ID), size = 2, shape = 21, color = "black", stroke = 0.1) +
  geom_errorbar(aes(ymin = AvgRedwood - AvgResidual, ymax = AvgRedwood + AvgResidual), 
                width = 0.2, color = "grey13", alpha = 0.3) +
  # 自动添加时间段分隔虚线
  annotate("segment", 
           x = unique(AllSite_RedwoodResidual$period_end)[-length(unique(AllSite_RedwoodResidual$period_end))] + 
             cumsum(c(0, diff(unique(AllSite_RedwoodResidual$period_end))))[-length(diff(unique(AllSite_RedwoodResidual$period_end)))],
           xend = .data$x, y = 0, yend = Inf,
           linetype = "dotted", color = "black", size = 0.8) +
  theme_bw() +
  theme(
    plot.title = element_text(hjust = 0.5),
    legend.title = element_text(hjust = 0.5, size = 16),
    legend.text = element_text(size = 12),
    axis.text.x = element_text(size = 12, angle = 45, hjust = 1, color = "black"),
    axis.text.y = element_text(size = 12, color = "black"),
    axis.title.x = element_text(size = 14, margin = margin(t = 10)), 
    axis.title.y = element_text(size = 14, margin = margin(r = 10))
  ) +
  scale_color_manual(
    values = c("BurnRed" = "red", "TimberRed" = "darkgreen", "ControlRed" = "dodgerblue3"),
    labels = c("Prescribed Fire", "Timber Harvest", "Control")
  ) +
  scale_fill_manual(
    values = c("BurnRed" = "red", "TimberRed" = "darkgreen", "ControlRed" = "dodgerblue3"),
    labels = c("Prescribed Fire", "Timber Harvest", "Control")
  ) +
  # 自定义x轴标签：显示小时+年月
  scale_x_continuous(
    breaks = AllSite_RedwoodResidual$global_x,
    labels = paste(format(AllSite_RedwoodResidual$Time, "%H:%M"), "\n", format(AllSite_RedwoodResidual$period, "%Y-%m"))
  ) +
  labs(x = "Date Range", y = "Average Redwood")

关键说明 #

• global_x将离散的时间段转换为连续无间隙的x轴数值，解决了原代码中因子与日期时间轴的类型冲突
• 分隔虚线位置根据时间段结束位置自动计算，无需手动指定x值
• x轴标签同时显示小时和年月，清晰区分不同时间段

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