查看原文
其他

R绘图模板——散点+回归线+误差线+显著性!!!

王志山 科研后花园 2023-09-08
 

点击上方

“科研后花园”

关注我们

图片灵感来源:

代码如下:1、设置工作环境:
rm(list = ls())setwd("D:/桌面/散点+回归线+误差线+显著性")
#加载R包library(ggplot2) # Create Elegant Data Visualisations Using the Grammar of Graphicslibrary(ggpubr) # 'ggplot2' Based Publication Ready Plotslibrary(ggpmisc) # Miscellaneous Extensions to 'ggplot2'2、加载数据(随机编写,无实际意义):
df <- read.table("data.txt", header = 1, check.names = F)

3、数据均值计算及整理:
#分别求取size,X和Y的均值X <- aggregate(X ~ group, df, function(x) c(mean = mean(x), sd = sd(x)))Y <- aggregate(Y ~ group, df, function(x) c(mean = mean(x), sd = sd(x)))size <- aggregate(size ~ group, df, function(x) c(mean = mean(x), sd = sd(x)))# 保留两位小数X$X <- round(X$X, 2)Y$Y <- round(Y$Y , 2)size$size <- round(size$size , 2)#整理绘图数据data <- data.frame( group=X$group, X=X$X[,1], X_error=X$X[,2], Y=Y$Y[,1], Y_error=Y$Y[,2], size=size$size[,1])

4、绘图:
#绘图ggplot(data, aes(X, Y))+ #水平方向误差线 geom_errorbar(aes(xmin=X-X_error,xmax=X+X_error), linewidth = 0.8, width = 0)+ #垂直方向误差线 geom_errorbar(aes(ymin=Y-Y_error,ymax=Y+Y_error), linewidth = 0.8, width = 0)+ #散点 geom_point(aes(color = group, size = size))+ #回归线 geom_smooth(method = "lm", color="red", fill = "red", se=T, formula = y ~ x, linetype=1, alpha=0.2)+ #回归方程及R2、p值 # stat_cor(method = "pearson",label.x = 17, label.y = 27, size=4.5)+ stat_poly_eq(formula = y ~ x, aes(label = paste(eq.label,after_stat(adj.rr.label),..p.value.label..,sep = "~~~~")), parse = TRUE, size=4.5, label.x = "left")+ #点大小范围 scale_size_continuous(range = c(3, 12))+ #图例 guides(color=guide_legend(override.aes = list(size=5,alpha=1)), size = "none")+ #标题 labs(x = "X_lab", y = "Y_lab", color = NULL)+ #主题 theme_bw()+ theme(panel.grid=element_blank(), axis.text=element_text(color='black',size=12), axis.title = element_text(color='black',size=14), legend.position = c(0.92,0.35))+ #颜色 scale_color_manual(values = c("#ff0000","#ffed00","#ff0092","#00b2a9", "#00c7f2","#dc5034","#a626aa","green","blue"))

PS: 以上内容是小编个人学习代码笔记分享,仅供参考学习,欢迎大家一起交流学习。


完整代码(附带注释)及原始数据获取:

夸克网盘(手机端下载夸克APP进行搜索提取):

链接:https://pan.quark.cn/s/f691cd1630a0

提取码:Ngw8

更多推荐

R绘图模板——散点+分组+size+显著性!!!

R可视化——基于MicrobiotaProcess包进行物种差异分析!

R绘图模板——气泡图+注释+配对连线!!!

基于MicrobiotaProcess包进行PCoA分析!

R绘图模板——世界地图的绘制及拓展!!!

跟着Nature学绘图——柱状图+散点+配对连线+显著性!!!

R绘图模板——中国地图+散点+柱状图!!!

R绘图模板——热图和柱状堆积图组合图!!!

R绘图模板源代码获取!

不会用代码做RDA分析?没关系,Canoco软件帮你实现!

扩增子测序数据分析还不会?小编整理的全套R语言代码助您轻松解决问题!

您可能也对以下帖子感兴趣

文章有问题?点此查看未经处理的缓存