本文是回到B站关注者的提问，他的问题如下
原始数据存储在一个excel文件里，这个excel文件里有三个子表格，每一个子表格的数据如下：
image.png
总的数据格式
image.png
现在的需要是做如下的图
image.png
接下来就介绍如何利用原始数据到最终的图的ggplot2的代码
首先是将3个子表格的数据整理到一张表格里
比如这里我新建了一个子表格sheet4,数据最终的格式如下
image.png
接下来R语言里操作
首先是读取数据
library(readxl)
df<-read_excel("prac.xlsx",
sheet = "sheet4")
df
将宽格式转换为长格式
library(tidyverse)
df %>%
pivot_longer(!var4) -> new_df
new_df
定义误差线函数
这里用到的是标准误
ebtop<-function(x){
return(mean(x)+sd(x)/sqrt(length(x)))
}
ebbottom<-function(x){
return(mean(x)-sd(x)/sqrt(length(x)))
}
ggplot2作图
library(ggplot2)
ggplot(data=new_df,aes(x=name,y=value,fill=var4))+
stat_summary(geom = "bar",fun = "mean",
position = position_dodge(0.9))+
stat_summary(geom = "errorbar",
fun.min = ebbottom,
fun.max = ebtop,
position = position_dodge(0.9),
width=0.2)+
scale_y_continuous(expand = expansion(mult = c(0,0.1)))+
theme_bw()+
theme(panel.grid = element_blank())+
scale_fill_manual(values = c("#e20612","#ffd401","#00b0eb"),
name="")+
labs(x="XXXXX",y="YYYYY")
image.png
调整不同分组之间的顺序
new_df$name<-factor(new_df$name,
levels = c("var2","var3","var1"))
ggplot(data=new_df,aes(x=name,y=value,fill=var4))+
stat_summary(geom = "bar",fun = "mean",
position = position_dodge(0.9))+
stat_summary(geom = "errorbar",
fun.min = ebbottom,
fun.max = ebtop,
position = position_dodge(0.9),
width=0.2)+
scale_y_continuous(expand = expansion(mult = c(0,0.1)))+
theme_bw()+
theme(panel.grid = element_blank())+
scale_fill_manual(values = c("#e20612","#ffd401","#00b0eb"),
name="")+
labs(x="XXXXX",y="YYYYY")
image.png
调整组内柱子之间的顺序
new_df$var4<-factor(new_df$var4,
levels = c("group2",
"group3",
"group1"))
ggplot(data=new_df,aes(x=name,y=value,fill=var4))+
stat_summary(geom = "bar",fun = "mean",
position = position_dodge(0.9))+
stat_summary(geom = "errorbar",
fun.min = ebbottom,
fun.max = ebtop,
position = position_dodge(0.9),
width=0.2)+
scale_y_continuous(expand = expansion(mult = c(0,0.1)))+
theme_bw()+
theme(panel.grid = element_blank())+
scale_fill_manual(values = c("#e20612","#ffd401","#00b0eb"),
name="")+
labs(x="XXXXX",y="YYYYY")
image.png
这里新学到一个知识点是，柱子默认是不贴底的，如果要贴底使用函数scale_y_continuous(expand = expansion(mult = c(0,0.1)))
mult对应的两个值一个是控制下面，一个是控制上面，贴底就设置为0 就可以了
公众号推文由获取数据和代码的方式
欢迎大家关注我的公众号
小明的数据分析笔记本
小明的数据分析笔记本 公众号 主要分享：1、R语言和python做数据分析和数据可视化的简单小例子；2、园艺植物相关转录组学、基因组学、群体遗传学文献阅读笔记；3、生物信息学入门学习资料及自己的学习笔记！
