R ggplot2绘制含跨年日期值哑铃图线段异常如何解决

直接使用同一年的日期映射跨年洄游的起止点时，跨年记录的结束儒略日数值小于起始儒略日，geom_segment 会按照x轴坐标值从小到大自动连线，最终错误覆盖全年无洄游的时段。

提前预处理线段数据集，对跨年/非跨年分组做区分处理：

• 非跨年洄游分组：保留单段连线，直接连接起始、结束儒略日
• 跨年洄游分组：拆分为两段连线，第一段从起始儒略日连到当年年末（第365天），第二段从次年年初（第1天）连到结束儒略日
• 绘图时调用预处理完成的线段数据集，避免坐标自动排序导致的连线错误

完整修正代码 #

数据预处理 #

替换原代码中拆分debut、fin表之后的处理逻辑：

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

# 整理为每个物种-生活史阶段为一行的宽表，存储起止儒略日
segment_df <- table_fin %>%
  pivot_wider(names_from = deb_fin, values_from = valeur) %>%
  rowwise() %>%
  mutate(
    # 标记是否为跨年洄游
    is_cross = fin < debut,
    # 第一段线的起止坐标
    seg1_start = debut,
    seg1_end = ifelse(is_cross, 365, fin),
    # 第二段线的起止坐标，非跨年分组设为NA
    seg2_start = ifelse(is_cross, 1, NA),
    seg2_end = ifelse(is_cross, fin, NA)
  ) %>%
  ungroup() %>%
  # 转为长表，适配ggplot的geom_segment映射要求
  select(sp_stade, seg1_start:seg2_end) %>%
  pivot_longer(
    cols = -sp_stade,
    names_to = c("seg_id", "pos"),
    names_sep = "_"
  ) %>%
  pivot_wider(names_from = pos, values_from = value) %>%
  filter(!is.na(start))

绘图代码 #

替换原代码中的geom_segment层：

p <- ggplot(table_fin) +
  geom_segment(
    data = segment_df,
    aes(
      x = as.Date(start, origin = as.Date("2024-01-01")),
      xend = as.Date(end, origin = as.Date("2024-01-01")),
      y = sp_stade,
      yend = sp_stade
    ),
    color = "#aeb6bf",
    size = 4.5,
    alpha = 0.5
  ) +
  scale_x_date(date_breaks = "2 weeks", date_labels = "%d/%m") +
  geom_point(
    aes(
      x = as.Date(valeur, origin = as.Date("2024-01-01")),
      y = sp_stade,
      color = deb_fin
    ),
    size = 4
  )
p

注意事项 #

• 转日期类型时统一使用平年作为基准origin即可，不需要单独处理闰年，儒略日的相对位置不会出现偏差
• 修正后所有跨年分组的线段会在年末年初自然断开，不会出现横跨全年的错误连线，完全匹配银鳗、产卵后阶段、LPR/LPM等类群的跨年洄游规律

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