如何在 ggplot2 中对齐不同年份的周数？

这也许是一种使用 .从本质上讲，将每个以休息时间作为周数绘制，并重新应用开始周日期作为标签，以便图对齐：patchwork

library(tidyverse)
library(ggthemes)

dates_weeks <- df |>
  mutate(date = mdy(paste0(date_in, year)),
         week = week(date)) 

min_week <- min(dates_weeks$week)
max_week <- max(dates_weeks$week)

year_breaks <- 2010:2013 |>
  set_names() |>
  map(\(year_n) {
    
    year_dat <- dates_weeks |> 
      filter(year == year_n) |> 
      select(date_in, week) |> 
      unique()
    
    week_map <- character(max_week - min_week + 1)
    names(week_map) <- seq(min_week, max_week)
    
    week_map[as.character(year_dat$week)] <- as.character(year_dat$date_in)
    
    
    week_map
    
  })

library(patchwork)

dates_weeks |>
  mutate(week = factor(week), Treatment = factor(Treatment)) |> 
  nest(data = -year) |>
  mutate(pl = map2(data, year, \(data, year) {
    # browser()
    
    ggplot(data, aes(x = week, y = total_count, fill = Treatment)) +
      geom_col() +
      scale_x_discrete(
        labels = year_breaks[[as.character(year)]],
        drop = FALSE
      ) +
      facet_wrap(year, strip.position = "right") +
      scale_fill_manual(values = c("IP" = "blue", "LD" = "red"), drop = FALSE) +
      theme_few(base_size = 10)
    
  })) |>
  pull(pl) |>
  wrap_plots(ncol = 1) +
  plot_layout(guides = "collect") &
  theme(
    legend.position = "top",
    axis.text.x = element_text(
      angle = 90,
      hjust = 1,
      vjust = 0.5,
      color = "black"
    ),
    axis.text.y = element_text(vjust = 0.5, color = "black"),
    axis.title = element_text(color = "black"),
    strip.text = element_text(face = "bold", color = "black")
  )

更新：

如果您需要根据日期编号对齐日期，那么我们需要获取周数并将其用于 .ggplot

我不知道有一种方法可以在分面中拥有自定义轴刻度标签，所以我只使用了轴上的周数并为实际日期添加了标签。不过，您可以创建多个绘图并将它们拼接在一起。

library(dplyr)
library(tidyr)
library(ggplot2)
library(ggthemes)
  
df %>% 
  mutate(across(where(is.factor), as.character),
         week_no = strftime(paste(year, date_in, sep = "/"), format = "%V")) %>% 
  filter(total_count != 0) %>%
  pivot_wider(id_cols = c(week_no, year, date_in), 
              names_from = c(Treatment), values_from = c(total_count)) %>% 
  full_join(expand.grid(week_no = unique(.$week_no), year = unique(.$year)), .) %>% 
  replace_na(., list(date_in = "", LD = 0, IP = 0)) %>% 
  pivot_longer(LD:IP, names_to = "Treatment", values_to = "total_count") %>% 
  mutate(date_in = ifelse(Treatment == "IP", date_in, "")) %>% 
 ggplot(., aes(x = week_no, y = total_count, fill = Treatment)) +
  geom_col() +
  geom_text(aes(y=0, label=date_in), vjust=-0.5, color = "darkgray") +
  facet_wrap(~ year, ncol = 1, scales = "free_y", strip.position = "right") + 
  scale_fill_manual(values = c("IP" = "navyblue", "LD" = "darkred")) +
  xlab("Week no.") + ylab("Total Count") +
  theme_few(base_size = 10, base_family = "Arial") +
  theme(axis.text.x = element_text(angle = 90, hjust = 1, 
                                   vjust = 0.5, color = "black"),
        axis.ticks.x = element_blank(),
        axis.text.y = element_text(vjust = 0.5, color = "black"),
        strip.text = element_text(face = "bold", color = "black"),
        legend.position = "top")

^{创建于 2023-09-25 使用 reprex v2.0.2}

原文答案：

您可以过滤掉空行，并仅显示具有 .free_xdate_intotal_count

library(ggplot2)
library(ggthemes)

ggplot(df[df$total_count != 0,], aes(x = date_in, y = total_count,
                                     fill = Treatment)) +
  geom_col() +
  facet_wrap(~ year, ncol = 1, scales = "free_y", strip.position = "right") + 
  scale_fill_manual(values = c("IP" = "blue", "LD" = "red")) +
  theme(axis.text.x = element_text(angle = 90, hjust = 1, vjust = 0.5, 
                                   color = "black"),
        axis.text.y = element_text(vjust = 0.5, color = "black"),
        axis.title = element_text(color = "black"),
        strip.text = element_text(face = "bold", color = "black"),
        legend.position = "top") +
  theme_few(base_size = 10, base_family = "Arial")

^{创建于 2023-09-23 使用 reprex v2.0.2}