温馨提示×

温馨提示×

您好,登录后才能下订单哦!

密码登录×
登录注册×
其他方式登录
点击 登录注册 即表示同意《亿速云用户服务条款》

r语言中如何使用ggplot2绘制气泡图

发布时间:2022-03-18 15:11:20 来源:亿速云 阅读:1103 作者:小新 栏目:开发技术

这篇文章主要为大家展示了“r语言中如何使用ggplot2绘制气泡图”,内容简而易懂,条理清晰,希望能够帮助大家解决疑惑,下面让小编带领大家一起研究并学习一下“r语言中如何使用ggplot2绘制气泡图”这篇文章吧。

Pathway富集结果常用气泡图来显示,其基本代码如下:

案例数据:

> dat
     ko_id                                Kegg_pathway Rich_factor       Pvalue DEGs
1  ko00195                              Photosynthesis    3.543363 1.000000e-15    5
2  ko00196           Photosynthesis - antenna proteins    5.575221 1.500000e-10    6
3  ko03008           Ribosome biogenesis in eukaryotes    2.464696 4.920000e-06    7
4  ko00710 Carbon fixation in photosynthetic organisms    2.567654 4.510000e-05    9
5  ko01200                           Carbon metabolism    1.738592 7.350000e-05   11
6  ko03030                             DNA replication    2.453097 1.772599e-03   14
7  ko00030                   Pentose phosphate pathway    2.397575 1.772599e-03   12
8  ko00860        Porphyrin and chlorophyll metabolism    2.574238 1.772599e-03   14
9  ko00062                       Fatty acid elongation    2.857551 2.316650e-03   18
10 ko00630     Glyoxylate and dicarboxylate metabolism    2.129425 5.299956e-03   19
11 ko00230                           Purine metabolism    1.671398 5.299956e-03    4

绘图常利用geom_point几何对象函数实现气泡图绘制过程,同时实现坐标轴翻转

注意显著性大小-log10(Pvalue)映射到的点颜色上、差异基因数映射到点大小上。

kegg_point=ggplot(dat,aes(x=Kegg_pathway,y=Rich_factor))+
  geom_point(aes(color=-log10(Pvalue),size=DEGs),alpha=0.8)+
  coord_flip()+
  scale_color_gradient(low = "green",high = "red")
kegg_point


不过需要注意的是在部分情况下,对应的DEGs数量较小,自动匹配的点的大小较小,在图片上难以显示明显,可以利用scale_size_continuous(range)参数控制,

譬如:

kegg_point=kegg_point+scale_size_continuous(range=c(4,10))
kegg_point

size对应的range参数可以有效控制点大小的范围。

以上是“r语言中如何使用ggplot2绘制气泡图”这篇文章的所有内容,感谢各位的阅读!相信大家都有了一定的了解,希望分享的内容对大家有所帮助,如果还想学习更多知识,欢迎关注亿速云行业资讯频道!

向AI问一下细节

免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:is@yisu.com进行举报,并提供相关证据,一经查实,将立刻删除涉嫌侵权内容。

AI