大橙子网站建设,新征程启航
为企业提供网站建设、域名注册、服务器等服务
小编给大家分享一下如何利用R语言的ggplot2包绘制KEGG气泡图,希望大家阅读完这篇文章之后都有所收获,下面让我们一起去探讨吧!
专注于为中小企业提供网站建设、网站制作服务,电脑端+手机端+微信端的三站合一,更高效的管理,为中小企业温江免费做网站提供优质的服务。我们立足成都,凝聚了一批互联网行业人才,有力地推动了上千余家企业的稳健成长,帮助中小企业通过网站建设实现规模扩充和转变。
一 载入数据集和R包
library(ggplot2)pathway = read.csv("KEGG.csv",header=TRUE,check.names = FALSE)head(pathway)
不同软件得到的KEGG结果的列名称可能不一致,但是这几列几乎都会有。
二 绘制KEGG气泡图
2.1初始化数据并绘制散点图
ggplot(pathway,aes(Pvalue,PATHWAY)) +geom_point()
可在以下几个方面进行优化:
A:标题,横纵坐标轴;
B:按照通路上基因的多少定义点的大小;
C:根据P值定义点的颜色;
2.2 修改点的大小
#按照Gene个数定义点的大小ggplot(pathway,aes(Pvalue,PATHWAY)) +geom_point() +geom_point(aes(size=Gene))
2.3 修改点的颜色
#定义连续型的配色ggplot(pathway,aes(Pvalue,PATHWAY))+geom_point(aes(size=Gene,color=-1*log10(Qvalue)))+scale_color_gradient(low="green",high = "red")
三 汇总展示
ggplot(pathway,aes(Pvalue,PATHWAY))+geom_point(aes(size=Gene,color=-1*log10(Qvalue)))+scale_color_gradient(low="green",high = "red")+ #labs(color=expression(-log[10](Qvalue)),size="Gene", ##expression函数定义函数样式 []添加下标,^添加上标x="Pvalue", ##自定义标轴y="Pathway name",title="Pathway enrichment")+ ##自定义坐标轴theme_bw() #去掉背景
看完了这篇文章,相信你对“如何利用R语言的ggplot2包绘制KEGG气泡图”有了一定的了解,如果想了解更多相关知识,欢迎关注创新互联行业资讯频道,感谢各位的阅读!