ggplot2绘制气泡图+分类+趋势线

最新推荐文章于 2024-05-29 16:09:04 发布

小Li先生

最新推荐文章于 2024-05-29 16:09:04 发布

阅读量4.8k

点赞数 4

分类专栏：科研绘图文章标签：数据可视化 r语言气泡图趋势线和置信区间

本文链接：https://blog.csdn.net/LeaningR/article/details/112560787

版权

科研绘图专栏收录该内容

15 篇文章 2 订阅

订阅专栏

今天有小伙伴让我给他做个气泡图，我让他给我发了个例子，作图基本分为两部分，数据整理和画图，数据整体一般注意变量类型、排列、分类等，画图细节比较多，先出图，后调试（坐标轴、图例、背景、图标、标签等），我的观点是先出图，查看大体样子，然后根据需求再选择调试与否，节省时间。因为调一个图真的太花时间。

ggplot2绘制气泡图+分类+趋势线

1、模板

左边这个是气泡图+标签+分类+趋势线，
右边这个是点图+分类+趋势线+置信区间+标签
在这里插入图片描述

2、数据列

三个变量，IN.rate和C.rate为分类变量，因此需要宽转长

> data.demon1 <- data.demon[c(5,8,12)];head(data.demon1)##提取需要数据
  INevent.rate Cevent.rate endpoint
1         1.46        1.41     0.78
2         1.26        1.48     1.37
3         2.23        2.26     1.39
4         1.17        1.46     1.40
5         1.10        1.60     1.60
6         2.30        2.91     1.71

> gather(data = data.demon1,key = "IN.C",value = "event.rate",-3)
   endpoint         IN.C event.rate
1      0.78 INevent.rate       1.46
2      1.37 INevent.rate       1.26
3      1.39 INevent.rate       2.23
4      1.40 INevent.rate       1.17
5      1.60 INevent.rate       1.10
6      1.71 INevent.rate       2.30
7      1.89 INevent.rate       1.86
8      1.99 INevent.rate       1.16
9      2.09 INevent.rate       1.72

3、画图

注意：
1）ggrepel包geom_text_repel()标注
2）

##加载包
library(ggplot2)
library(ggrepel)
library(dplyr)
library(tidyr)
##转化提取数据
data.demon2 <- gather(data = data.demon1,key = "IN.C",value = "event.rate",-3)
str(data.demon2)##查看变量类型
##画图
p.demon <- ggplot(data = data.demon2,aes(x=endpoint,y=event.rate))+##数据列
  geom_point(aes(size=event.rate,color=IN.C),alpha=0.7,show.legend = TRUE)+##点图大小、颜色、透明度、图例
  geom_smooth(aes(x=endpoint,y=event.rate,color=IN.C),method = lm,linetype=1,se=FALSE,span=1)+##趋势线、颜色、方法、线型、置信区间
  guides(color=guide_legend(title=NULL))+##去除图例标题
  scale_size(range = c(1, 10),guide=FALSE)+##气泡大小区间，图例标题去除
  scale_color_manual(values=c("CornflowerBlue","Gold"),
                     breaks=c("INevent.rate","Cevent.rate"),
                     labels=c("1","2"))+##气泡颜色、图例名称修改)+##气泡颜色
  labs(x='Endpoint',y='Increase ratio')+##横纵坐标
  geom_text_repel(aes(label=event.rate))+##气泡标注
  annotate("text",x=1.0,y=2.8,label="atop(Y==3*X+5,R^2==0.9)",parse=T)+##添加公式
  #theme(legend.title = element_blank())+
  theme_bw();p.demon##设置背景theme
ggsave(p.demon,file="p_demon.png", height=7, width=10)##数据保存默认dpi300