在之前的文章中,分享了箱线图的绘制方法:
进一步,分享一下带误差棒的箱线图的绘制方法。
先来看一下成品效果:
特别提示:本期内容『数据+代码』已上传资源群中,加群的朋友请自行下载。有需要的朋友可以在公众号后台回复关键词【全家桶】查看加群方式。
此部分主要是读取原始数据,并进行绘图前的预处理。
# 加载必要的库
library(readxl)
library(tidyr)
library(ggplot2)
# 修改工作路径
# setwd("D:/R")
# 读取数据
data <- read_excel("data_box.xlsx")
data = data[,c(2:10)]
# 将数据转换为数据框格式
data <- as.data.frame(data)
# 将数据转换为长格式
data_long <- data %>%
pivot_longer(cols = -c(Sample))
作图不配色就好比做菜不放盐,总让人感觉少些味道。
但颜色搭配比较考验个人审美,需要多加尝试。
这里直接使用TheBestColors配色工具中的SCI顶刊配色:
library(TheBestColors)
map <- Best100(1)
获取方式:公众号(阿昆的科研日常)后台回复100
利用‘geom_boxplot’和‘stat_boxplot’函数,绘制初始带误差棒的箱线图。
p <- ggplot(data_long, aes(x = name, y = value)) +
stat_boxplot(geom="errorbar",width=0.2)+
geom_boxplot(aes(fill = name),width=0.4)+
scale_fill_manual(values=map)+
labs(x = 'Sample', y = 'Weight')
print(p)
为了插图的美观,对坐标轴细节等进行美化,并按期刊所需分辨率、格式输出图片。
# 细节调整
p+scale_y_continuous(breaks = seq(0, 5, 1), # y轴刻度线位置
expand = c(0,0), # y坐标轴与绘图区域边缘的扩展量
limits = c(0, 5) ) + # y轴显示范围
theme_classic(base_size = 8) + # 设置基本字号为8pt
theme(
axis.line = element_blank(), # 删除原始轴线
panel.border = element_rect(size = 0.5,fill = 'transparent'), # 添加全框线
panel.grid.major.y = element_line(size = 0.2, linetype = 'solid', color = "gray"), # 开启y轴主网格线
axis.ticks = element_line(size = 0.5), #设置刻度线粗细
axis.ticks.length = unit(1.5, "mm"), # 设置刻度线的长度
legend.position = 'NULL',# 不显示图例
axis.text = element_text(size = 10), # 设置坐标轴标签字号和颜色
axis.title = element_text(size = 10), # 设置坐标轴标题字号
)
# 保存图形
ggsave("带误差棒的箱线图.png", width = 15, height = 12, units = "cm", dpi = 300)
以上。
如果你觉得我的分享对你有帮助的话,欢迎大家在这里点赞、在看、分享。当然,也欢迎大家在这里打赏。互动越多,更新越快哦~
声明:本公众号的所有原创内容,在未经允许的情况下,不得用于任何商业用途,违者必究。