Figure 2B source data.xlsx

the codelibrary(ggplot2)library(readxl)library(dplyr)library(tidyr)library(gridExtra)library(ggpubr)data_long <- data %>% pivot_longer(cols = starts_with("value"), names_to = "measurement", values_to = "value")comparison_list <- list(c("Control", "Heat stress"))plot_list <- lapply(unique(data_long$sample), function(sample) { df <- data_long %>% filter(sample == !!sample) ggplot(df, aes(x = group, y = value, color = group, fill = group)) + geom_jitter(width = 0.15) + geom_boxplot(width = 0.4, alpha = 0.2) + scale_color_manual(values = c("#6F99AD", "#BF5960")) + scale_fill_manual(values = c("#6F99AD", "#BF5960")) + theme_test() + stat_compare_means(comparisons = comparison_list, method = "t.test", label = "p.signif") + labs(title = sample, y = sample, x = "") + scale_y_continuous(expand = expansion(c(0.1, 0.1))) + theme(axis.text = element_text(color = "black", size = 10), axis.line = element_blank(), legend.position = "none")})grid.arrange(grobs = plot_list, nrow = 2, ncol = 4)