如何在ggplot2中将X轴年份每3年分组并显示无数据区间？

完全可以实现，核心思路是先生成覆盖所有目标区间的完整分组列表，再将原数据匹配到这些分组中，最后基于完整分组绘图，就能显示无数据的区间。以下是具体步骤和代码示例：

1. 预处理数据：生成完整分组并匹配 #

假设你的数据框为df，包含year（年份）和diversity（物种多样性）两列：

library(dplyr)

# 1. 确定数据的年份范围
min_year <- min(df$year)
max_year <- max(df$year)

# 2. 生成所有3年一组的起始年份（确保覆盖从最小到最大年份的完整区间）
start_years <- seq(
  from = floor((min_year - 1)/3)*3 + 1,
  to = ceiling(max_year/3)*3,
  by = 3
)

# 3. 创建带区间标签的完整分组数据框
group_df <- tibble(
  year_group = factor(paste(start_years, start_years + 2, sep = "-"),
                     levels = paste(start_years, start_years + 2, sep = "-")),
  start = start_years,
  end = start_years + 2
)

# 4. 将原数据匹配到分组中，无数据的组会自动填充NA
full_df <- df %>%
  mutate(year_group = cut(
    year,
    breaks = c(group_df$start, max(group_df$end) + 1),
    labels = group_df$year_group,
    include.lowest = TRUE
  )) %>%
  right_join(group_df, by = "year_group")

2. 用ggplot2绘图 #

基于预处理后的full_df绘图，所有3年区间都会显示在X轴，无数据的区间对应位置会呈现空白（以柱状图为例）：

library(ggplot2)

ggplot(full_df, aes(x = year_group, y = diversity)) +
  geom_col(fill = "#2c3e50", na.rm = FALSE) + # na.rm=FALSE保留无数据的位置
  labs(x = "年份区间", y = "物种多样性") +
  theme(axis.text.x = element_text(angle = 45, hjust = 1)) # 旋转X轴标签避免重叠

关键说明 #

• 通过right_join确保所有预设的3年分组都被保留，无数据的分组对应的diversity为NA
• 绘图时设置na.rm = FALSE，让无数据的分组在X轴上仍占据位置，直观体现无数据区间

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