R语言ggplot2分组箱线图添加误差线的简单小例子

用户7010445

发布于 2021-10-13 11:46:02

4.9K0

发布于 2021-10-13 11:46:02

文章被收录于专栏：小明的数据分析笔记本

ggplot2 更改图例的顺序

箱线图添加误差线这个自己老是记不住，每次作图都得现查，今天的推文记录一下实现代码，方便自己以后查看。

首先是示例数据集

示例数据集还是使用鸢尾花数据集

部分如下

image.png

收下是读取数据

df<-read.csv("iris.csv")
head(df)

宽格式转换为长格式

reshape2::melt(df,id.vars="Species") -> dfa

基本箱线图

ggplot(data = dfa,
       aes(x=Species,y=value,fill=Species))+
  geom_boxplot()

image.png

添加误差线

这里使用到的是stat_boxplot()函数

ggplot(data = dfa,
       aes(x=Species,y=value,fill=Species))+
  geom_boxplot()+
  stat_boxplot(geom = "errorbar",
               width=0.3)

image.png

这样多了一个垂直线，不好看，我们把误差线的图层放到最下层，就是把代码写到boxplot的前面，然后加一些基本的美化

ggplot(data = dfa,
       aes(x=Species,y=value,fill=Species))+
  
  stat_boxplot(geom = "errorbar",
               width=0.3)+
  geom_boxplot(width=0.5)+
  scale_fill_material_d()+
  theme_bw()

image.png

这里还有一个小知识点是更改图例的顺序，现在图例从上到下依次是 setosa versicolor virginica,如果要反过来可以加一行代码guides(fill=guide_legend(reverse = T))

ggplot(data = dfa,
       aes(x=Species,y=value,fill=Species))+
  
  stat_boxplot(geom = "errorbar",
               width=0.3)+
  geom_boxplot(width=0.5)+
  scale_fill_material_d()+
  theme_bw()+
  guides(fill=guide_legend(reverse = T))

image.png

分组箱线图

ggplot(data = dfa,
       aes(x=Species,y=value,fill=variable))+
  
  stat_boxplot(geom = "errorbar",
               width=0.3,
               position = position_dodge(0.5))+
  geom_boxplot(width=0.5,
               position = position_dodge(0.5))+
  scale_fill_material_d()+
  theme_bw()+
  guides(fill=guide_legend(reverse = T))

image.png

分组箱线图需要我们在stat_boxplot()函数和geom_boxplot()函数里分别加上position = position_dodge(0.5)参数，目的是使误差线和箱子的位置对应上

如果不加这两个参数，效果如下

image.png

这里还有一个疑问是箱线图中间的线好像是中位数，如何把这个线更改为平均值呢？暂时没有想明白

今天推文的示例数据和代码可以直接留言20210929获取

欢迎大家关注我的公众号

小明的数据分析笔记本

小明的数据分析笔记本公众号主要分享：1、R语言和python做数据分析和数据可视化的简单小例子；2、园艺植物相关转录组学、基因组学、群体遗传学文献阅读笔记；3、生物信息学入门学习资料及自己的学习笔记！

本文参与腾讯云自媒体同步曝光计划，分享自微信公众号。

原始发表：2021-09-29，如有侵权请联系 cloudcommunity@tencent.com 删除

数据分析

本文分享自小明的数据分析笔记本微信公众号，前往查看

如有侵权，请联系 cloudcommunity@tencent.com 删除。

本文参与腾讯云自媒体同步曝光计划，欢迎热爱写作的你一起参与！

数据分析

登录后参与评论

0 条评论

热度