ReactomePA数据库下载及后续格式处理

原创

于 2024-05-16 16:43:32 发布 · 1.1k 阅读

10 ·

CC 4.0 BY-SA版权

文章标签：

#python #笔记 #经验分享

Reactome分析需要用到ReactomePA包处理

由于硬件能力有限，只能借助多个脚本实现（主打一个那个方便上手用哪个），总之自用能解决问题，及时记录一下，可能后面还能用的到的

最终想要得到的结果类型如下：

文件1：pathway_TERM2NAME.txt

文件2：pathway_TERM2SYMBOL.txt

第一步ReactomePA数据库下载：

rm(list = ls())
# 安装 ReactomePA
BiocManager::install("ReactomePA")
library(ReactomePA)
library(reactome.db)
ls("package:reactome.db")
keytypes(reactome.db)
library(dplyr)

PathwayId_Genes <- as.list(reactomePATHID2EXTID) %>% .[grep("HSA",names(.))]
PathwayId_PathwayName <- as.list(reactomePATHID2NAME) %>% .[grep("HSA",names(.))]
PathwayName_PathwayId <- as.list(reactomePATHNAME2ID) %>% .[grep("HSA",.)]
Genes_PathwayId <- as.list(reactomeEXTID2PATHID) %>% .[grep("HSA",.)]
name