使用 R 中 ggplot2 的 facet_wrap 功能对多个数据集进行箱线图比较？

Question

我有一个data.frame，如下所示。

library(tidyverse)
library(lubridate)

set.seed(129)

DF <- data.frame(Date = seq(as.Date("2001-01-01"), to = as.Date("2005-12-31"), by = "month"),
                 O1 = runif(60,1,5), R1 = runif(60,1,5), C1 = runif(60,1,5),
                 O2 = runif(60,1,5), R2 = runif(60,1,5), C2 = runif(60,1,5),
                 O3 = runif(60,1,5), R3 = runif(60,1,5), C3 = runif(60,1,5))

目标：

我想生成 3 个facetboxplot，其中facet 一个是比较O1, R1, and C1 的数据集。在facet 2 上，我希望看到boxplot 的O2, R2, and C2，同样也希望看到第三个facet。

示例：

我正在寻找plot 喜欢的附件。示例情节是2 facet，而我正在寻找3 facet。

Answer 1

更新：

DF1 <- DF %>% 
  as_tibble() %>% 
  mutate(Date= ymd(Date)) %>% 
  mutate(month = month(Date),
         year = year(Date)) %>% 
  pivot_longer(
    cols = -c(Date, month, year)
  ) %>% 
  mutate(facet = case_when(name %in% c("O1", "R1", "C1") ~ "facet1",
                           name %in% c("O2", "R2", "C2") ~ "facet2",
                           name %in% c("O3", "R3", "C3") ~ "facet3")) 

  ggplot(DF1, aes(x=factor(month), y=value)) +
  geom_boxplot(aes(fill=name)) +
  facet_wrap(.~facet, nrow = 3) 

第一个答案： 像这样？

library(tidyverse)
DF %>% 
  pivot_longer(
    cols = -Date
  ) %>% 
  mutate(facet = case_when(name %in% c("O1", "R1", "C1") ~ "facet1",
                           name %in% c("O2", "R2", "C2") ~ "facet2",
                           name %in% c("O3", "R3", "C3") ~ "facet3")) %>% 
  ggplot(aes(name, value)) +
  geom_boxplot() +
  facet_wrap(.~facet, nrow = 3)

Answer 2

试试：

DF %>% 
  pivot_longer(-Date) %>% 
  mutate(month = factor(month.abb[month(Date)], month.abb),
         groups = readr::parse_number(name)) %>% 
  group_by(groups) %>% 
  mutate(facet_groups = paste0(unique(name), collapse = ","),
         name = fct_reorder(name, groups)) %>% 
  ggplot(aes(x = month, y = value, fill = name)) + 
  geom_boxplot() +
  facet_wrap(facet_groups ~ .,
             ncol = 1) +
  labs(y = "Monthly Precipitation",
       x = element_blank(),
       fill = element_blank())