使用 R 的 facet_wrap 功能逐年比较观察数据和模拟数据？

Posted 2023-02-16

【中文标题】使用 R 的 facet_wrap 功能逐年比较观察数据和模拟数据？

【英文标题】：Year by year comparison of observed and simulated data using facet_wrap functionality of R?

【发布时间】：2021-10-04 11:13:55

【问题描述】：

我有以下data.frame

library(lubridate)
library(tidyverse)

set.seed(123)

DF <- data.frame(Date = seq(as.Date("2001-01-01"), to = as.Date("2005-12-31"), by = "day"),
                 Sim = runif(1826,1,5), Obs = runif(1826,3,7)) %>% 
  separate(Date, into = c("Year", "Month","Day")) %>% 
  mutate(Sim_Vol = Sim*86400, Obs_Vol = Obs*86400 ) %>% 
  filter(between(Month, 3,5))

我想逐年使用ggplot 的facet_wrap 功能为每个Year 找到Sim_Vol 和Obs_Vol 中的total sum，然后是bar plot 和Yearly total data比较。 附上一个示例图以进一步澄清

【参考方案1】：

这样的？

DF %>%
  group_by(Year) %>%
  summarise(across(c(Sim_Vol, Obs_Vol), sum)) %>%
  pivot_longer(cols = ends_with("Vol")) %>%
  ggplot(aes(x = name)) +
  geom_bar(aes(weight = value)) +
  facet_wrap(~Year)