投稿
  1. 极客之音首页
  2. 前端开发

R语言|ggplot2| 绘制KEGG气泡图

小半 • • 前端开发 • 阅读 387

如果你不相信努力和时光，那么成果就会是第一个选择辜负你的。不要去否定你自己的过去，也不要用你的过去牵扯你现在的努力和对未来的展望。不是因为拥有希望你才去努力，而是去努力了，你才有可能看到希望的光芒。R语言|ggplot2| 绘制KEGG气泡图，希望对大家有帮助，欢迎收藏，转发！站点地址：www.bmabk.com，来源：原文

在 RStudio 中使用 BiocManager 安装包

install.packages("BiocManager")
install.packages("ggplot2")

安装对应版本的 Rtools

详细指导说明
安装连接

安装完成后，使用命令安装 KEGGREST 以及所需要的相关包

BiocManager::install("KEGGREST")
BiocManager::install("fmcsR")
~~devtools::install_git("https://github.com/cran/RbioRXN.git")~~

包加载

library(KEGGREST)
~~library(RbioRXN)~~

查看KEGG数据库包含的数据

listDatabases()

获取单个数据集中的数据

pathway<- keggList("pathway")
head(pathway)

对单个数据库进行组织的选择

org <-keggList("pathway","hsa")
head(org)

如下图所示：
R语言|ggplot2| 绘制KEGG气泡图

从上面可以看出keggList不仅可以提取单个数据集还可以获取对应物种的信息。在这里我们发现同样的通路编码ID却不一样，map+num泛指KEGG中的所有通路；has+num指的是人类物种的通路信息。

获取所有的代谢反应和化合物数据

keggAll = get.kegg.all()
save(keggAll,file="C:/data/metabolism/database/KEGG/keggAll.Rdata")

提取数据

reaction=keggAll$reaction
write.csv(reaction," reaction.csv")
compound=keggAll$compound
write.csv(compound," compound.csv")

至此我们就可以将KEGG中的数据提取到本地进行接下来的分析处理。
提取方法参考
数据提取参考

绘制气泡图

参考

版权声明：本文内容由互联网用户自发贡献，该文观点仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请发送邮件至 举报，一经查实，本站将立刻删除。

文章由极客之音整理，本文链接：https://www.bmabk.com/index.php/post/203763.html

赞 (0)
小半的头像小半
0 0

相关推荐

极客之音——专业性很强的中文编程技术网站，欢迎收藏到浏览器，订阅我们！