R可视化——分组散点图的绘制以及拟合曲线和回归方程的添加！

 

点击上方

“科研后花园”

关注我们

代码如下

1、加载绘图包

rm(list=ls())#clear Global Environment
setwd('D:/桌面/分组散点图')#设置工作路径

#加载R包
library(ggplot2) # Create Elegant Data Visualisations Using the Grammar of Graphics
library(ggpmisc) # Miscellaneous Extensions to 'ggplot2'
library(RColorBrewer) # ColorBrewer Palettes
library(grid) # The Grid Graphics Package
library(scales) # Scale Functions for Visualization

2、读取数据——以Chiplot绘图平台数据为例

df <- read.table(file="data.txt",sep="\t",header=T,check.names=FALSE)

3、自定义颜色

#自定义颜色
col<-c("#be0027", "#cf8d2e")
#构建背景色
color1 <- colorRampPalette(brewer.pal(11,"PiYG"))(30)
color2 <- colorRampPalette(brewer.pal(11,"PuOr"))(30)

4、绘图

p <- ggplot(df,aes(x,y,fill=group))+
  geom_point(shape=21,size=3,alpha=0.5)+
  #添加回归曲线并添加置信区间
  geom_smooth(method = "lm",aes(color=group), se=T, 
              formula = y ~ x,
              linetype=1,alpha=0.5)+
  #添加回归方程
  stat_poly_eq(formula = y ~ x, 
               aes(color=group,label = paste(after_stat(eq.label),
                                 ..rr.label..,sep = "~~~")), parse = TRUE) +
  scale_fill_manual(values = col)+
  scale_color_manual(values = col)+
  theme_bw()+
  theme(panel.grid=element_blank(),
        axis.text=element_text(color='#333c41',size=12),
        legend.text = element_text(color='#333c41',size=12),
        legend.title = element_blank())+
  labs(x=NULL,y=NULL)
p

5、添加背景

grid.raster(alpha(color1, 0.2), 
            width = unit(1, "npc"), 
            height = unit(1,"npc"),
            interpolate = T)

6、分面并添加背景

p+facet_grid(~group, scales = "free")+
  theme(legend.position = "none")
#添加背景
grid.raster(alpha(color2, 0.2), 
            width = unit(1, "npc"), 
            height = unit(1,"npc"),
            interpolate = T)

温馨提示

如果你喜欢本文，请分享到朋友圈，想要获得更多信息，请关注我。