R中时间序列坐标轴调整：按年拆分绘图与添加月份标签

Hey there! Let's tackle your two main issues: splitting your time series into yearly plots, and adding proper month labels even with irregular sampling dates.

First, Prep Your Data #

First off, let's add a year column to your data frame—this makes grouping and plotting way easier. We'll use the lubridate package for clean date handling:

library(lubridate)
library(ggplot2)
library(dplyr)

# Add a year column to your tibble
DWPhyto <- DWPhyto %>%
  mutate(year = year(date))

Option 1: Faceted Plots (Most Efficient) #

Instead of manually splitting plots (and dealing with xlim headaches), use facet_wrap to create a grid of yearly plots. The scales="free_x" argument ensures each subplot's x-axis automatically adapts to that year's date range (no more missing axes!).

To add month labels, use scale_x_date to set monthly breaks and format the labels to show month names:

ggplot(DWPhyto, aes(x = date, y = data)) +
  geom_line(color = "#2c3e50") +
  # Split into yearly plots, independent x-axes
  facet_wrap(~year, scales = "free_x", ncol = 2) +
  # Format x-axis to show monthly labels
  scale_x_date(
    breaks = "1 month",  # Add a break every month
    labels = scales::date_format("%b"),  # Show 3-letter month abbreviations (e.g., Jan, Feb)
    name = "Month"
  ) +
  labs(
    y = "Data Value",
    title = "Yearly Breakdown of Phyto Data (2008-2017)",
    subtitle = "Irregular Sampling Dates"
  ) +
  theme_minimal() +
  theme(
    axis.text.x = element_text(angle = 45, hjust = 1),  # Rotate labels for readability
    strip.background = element_rect(fill = "#f8f9fa"),
    strip.text = element_text(face = "bold")
  )

Why your xlim approach caused missing axes: When you set xlim manually, if any of your yearly dates fall outside the range you specified (even by a day), ggplot will drop the axis entirely. Using scales="free_x" avoids this by auto-adjusting each plot's x-range to its own data.

Option 2: Save Individual Yearly Plots #

If you need separate files for each year instead of a faceted grid, use purrr::walk to loop through each year, generate a plot, and save it:

library(purrr)

# Get all unique years from your data
unique_years <- unique(DWPhyto$year)

# Loop through each year to create and save plots
walk(unique_years, function(current_year) {
  # Filter data to the current year
  yearly_data <- filter(DWPhyto, year == current_year)
  
  # Create the plot
  yearly_plot <- ggplot(yearly_data, aes(x = date, y = data)) +
    geom_line(color = "#2c3e50") +
    scale_x_date(
      breaks = "1 month",
      labels = scales::date_format("%b"),
      name = "Month"
    ) +
    labs(
      y = "Data Value",
      title = paste("Phyto Data -", current_year)
    ) +
    theme_minimal() +
    theme(axis.text.x = element_text(angle = 45, hjust = 1))
  
  # Save the plot to your working directory
  ggsave(
    filename = paste0("phyto_data_", current_year, ".png"),
    plot = yearly_plot,
    width = 8,
    height = 5,
    dpi = 300
  )
})

This will generate a separate PNG file for each year, each with properly formatted month labels on the x-axis.

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