在 R 中绘制时间序列数据的每日分布

Question

我有一个时间序列数据（日期列和值列）。我正在尝试每日分布图。

下图中是每周分布图，绘制了一周中几天的值。同样，我正在尝试绘制每日分布图，其中 x 轴为月，y 轴为值，该图有 10 行，其中每行为您提供日期 1、日期 2、日期 3 等等，直到日期 10（因为在一个子图中 30 天会很笨拙，所以我想将地块分为 3、1-10、11-20 和 21-31）

每周分发代码供参考：

#dummy data
start_date <- as.Date("2020-01-01")
end_date <- as.Date("2021-12-31")
date_seq <- seq(from = start_date, to = end_date, by = "day")
set.seed(123)
value <- round(runif(length(date_seq), min = 10000, max = 100000000), 0)
df <- data.frame(date = date_seq, value = value)

df$week_number <- as.numeric(format(as.Date(df$date), "%U")) + 1
df$weekday <- weekdays(as.Date(df$date))
df$year <- as.numeric(format(as.Date(df$date), "%Y"))
years <- unique(df$year)

# Create a list of ggplots, one for each year
plots <- lapply(years, function(y) {
  year_df <- df[df$year == y, ]
  ggplot(year_df, aes(x = week_number, y = value, color = weekday)) +
    geom_line() +
    scale_color_discrete(limits = c("Monday", "Tuesday", "Wednesday", "Thursday", "Friday", "Saturday", "Sunday")) +
    ggtitle(paste("Weekday Distribution", y)) +
    xlab("Week number") +
    ylab("Value") +
    theme(legend.key.size = unit(0.4, "cm")) +
    theme(plot.title = element_text(hjust = 0.5, vjust = 1.5))
library(cowplot)
plot_grid(plotlist = plots, ncol = 1)

所以最后，将有三个图（1 到 10 个日期、11 到 20 个日期和 21 到 31 个日期），每个图将包含 2 个子图（因为日期范围从 2020 年到 2021 年）。谁能帮我这个？

Answer 1

下面我将如何做到这一点。 lubridate 包是你的朋友。对于分组，使用削减。

结果是（在我看来）非常无用的混乱线条。但这不是我不认可这种可视化的唯一原因。我觉得这在某种程度上打败了时间序列的要点……其中一点是可视化数据的自相关性。人为地从每个月中仅分离出特定日期会极大地影响使用时间序列的这一特殊优势（也许是原因）。您不仅会丢失信息，还会使您自己的分析生活变得更加复杂。

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

df %>%
  mutate(day = mday(date), 
         day_group = cut(day, c(1,11,21, 31), incl = T),
         month = month(date, label = T, abbr = T))  %>%
ggplot(aes(x = month, y = value, color = day, group=interaction(day, day_group))) +
  geom_line() +
  theme(legend.key.size = unit(0.4, "cm"),
        plot.title = element_text(hjust = 0.5, vjust = 1.5), 
        axis.text.x = element_text(angle = 90)) +
  facet_wrap(year~day_group)