之前的文章中,分享了R语言折线图的绘制模板:
带标记点的折线图:
进一步,分享一下带置信区间的折线散点图的绘制模板。
先来看一下成品效果:
特别提示:本期内容『数据+代码』已上传资源群中,加群的朋友请自行下载。有需要的朋友可以在公众号后台回复关键词【全家桶】查看加群方式。
1. 数据准备
此部分主要是读取原始数据,并进行绘图前的预处理。
# 加载库
library(readxl)
library(tidyr)
library(ggplot2)
library(ggpubr)
# 修改工作路径
# setwd("D:/R")
# 读取数据
data <- read_excel("data_CI.xlsx")
# 将数据转换为数据框格式
data <- as.data.frame(data)
2. 颜色定义
作图不配色就好比做菜不放盐,总让人感觉少些味道。
但颜色搭配比较考验个人审美,需要多加尝试。
这里直接使用TheBestColors配色工具中的SCI顶刊配色:
library(TheBestColors)
map <- Best100(16)
CC <- map[c(1,4,3)]
3.带置信区间的折线散点图绘制
利用‘stat_smooth’、‘stat_regline’和‘geom_point’函数,绘制初始带置信区间的折线散点图。
# 带置信区间的折线散点图绘制
formula = y ~ poly(x, 2, raw = TRUE) # 定义拟合公式
p <- ggplot(data, aes(GSmax, GPPsat, color=Data)) +
geom_point(size = 3) + # 散点
stat_smooth(method = "lm", formula = formula,size = 1.2) + # 拟合线与置信区间
stat_regline_equation(
aes(group=Data, label = paste(..eq.label.., ..adj.rr.label.., sep = "~~~~")),
formula = formula, size = 3,
label.x.npc = 0.4,
label.y.npc = 0.23,
show.legend = FALSE
)+ # 公式
scale_color_manual(values = CC)+ # 散点、折线颜色
labs(x = 'Day', y = 'Weight change(%)') # 坐标轴标题
print(p)
4. 细节优化
为了插图的美观,对坐标轴细节等进行美化,并按期刊所需分辨率、格式输出图片。
# 细节调整
p+theme_classic(base_size = 8) + # 设置基本字号为8pt
theme(
axis.line = element_blank(), # 删除原始轴线
panel.border = element_rect(size = 0.5,fill = 'transparent'), # 添加全框线
axis.ticks = element_line(size = 0.5), #设置刻度线粗细与颜色
axis.ticks.length = unit(1.5, "mm"), # 设置刻度线的长度
legend.position = 'top', # 设置图例位置框
axis.text = element_text(size = 10), # 设置坐标轴标签字号和颜色
legend.text = element_text(size = 10), # 设置图例文本字号
axis.title = element_text(size = 10), # 设置坐标轴标题字号
)
# 图片导出
ggsave("带置信区间的折线散点图.png", width = 12, height = 12, units = "cm", dpi = 600)
以上。