先放最终效果图

结果图
# bar_plot 2021/01/19
# 导入所需的包
library(ggplot2)
library(ggthemr)
library(ggsignif)
library(tidyverse)
library(dplyr)
library(ggpubr)
library(devEMF)
# 导入并处理数据,需要两张表,一个是原始汇总表格,另外一些是每一个测量数据的均值和标准差
data1 <- read.csv(file = "../all_enzymatic.txt", sep = "\t", header = T)
data2 <- dplyr::select(data1,1,3,4,5,6)
# SOD
data_SOD <- dplyr::select(data2,1,5)
SOD_mean <- data_SOD %>%
dplyr::group_by(group) %>%
dplyr::summarize(
count=n(),
mean = mean(SOD),
sd = sd(SOD)
)