跟着NC学绘图-ggplot2批量绘制配对箱线图

library(tidyverse)
library(patchwork)

inflam <- read_tsv("data.tsv") %>% 
  arrange(Group) %>% 
  mutate(id = rep(1:(n()/3),each=3)) %>% drop_na() %>% 
  mutate(id=as.character(id)) %>% 
  mutate(across(where(is.numeric),log2),
         id=as.numeric(id))

variables <- c("GM_CSF", "IFNy", "IL_1B", "IL_2", "IL_4", "IL_5", 
               "IL_6", "IL_8", "IL_10", "IL_12p70", "IL_13", 
               "IL_17A", "IL_23", "TNFa")

plot_variable <- function(variable, data) {
  ggplot(data, aes(x = Time, y = .data[[variable]], fill = Group, color = Group)) +
    geom_boxplot(outlier.shape = NA,staplewidth = 0.5) +
    geom_line(aes(group = id), size = 0.5, alpha = 0.2) +
    geom_jitter(shape = 21, size = 1, stroke = 0.5,
                position = position_jitterdodge(dodge.width = 0.8,
                                                jitter.width = 0.2)) +
    facet_wrap(~factor(Group, levels = c('IF-P', 'CR'))) +

    scale_fill_manual(name = NULL, values = c("magenta3", "darkturquoise")) +
    scale_color_manual(name = NULL, values = c("black", "black")) +
    labs(y = paste("Log2(", variable, " Concentration)", sep = ""), x = "Time (weeks)") +
    scale_x_discrete(labels = c("WK0" = "0", "WK4" = "4", "WK8" = "8")) +
    theme_bw() +
    theme(legend.position = "none", 
          strip.background = element_rect(colour = "black", fill = "black", size = 1),
          strip.text = element_text(colour = 'white', size = 6, face = "bold"),
          axis.text = element_text(size = 6, color = "black"),
          axis.title = element_text(size = 8, color = "black",face="bold"),
          panel.border = element_rect(color = "black", fill = NA, size = 1),
          plot.tag = element_text(size =12,face="bold",color="black")) +
    ylim(NA, max(data[[variable]], na.rm = TRUE) * 1.1)
}

plots <- lapply(variables[c(5, 7, 8, 11)], function(variable) {
  plot_variable(variable, inflam)
})

do.call("wrap_plots", c(plots, ncol = 4))+
  plot_annotation(tag_levels = 'a')

ggsave(filename = "./Fig2.png", 
       plot = do.call("wrap_plots", c(plots, ncol = 4))+
         plot_annotation(tag_levels = 'A'),
       height = 6, width = 18, units = "cm")