R语言 转换metascape GO结果中的gene name为ensembl ID




# 安装并加载必要的包
if (!requireNamespace("AnnotationDbi", quietly = TRUE)) {
    install.packages("AnnotationDbi")
}
if (!requireNamespace("org.Hs.eg.db", quietly = TRUE)) {
    install.packages("org.Hs.eg.db")
}
library(AnnotationDbi)
library(org.Hs.eg.db)
 
# 定义函数,将基因名转换为对应的Ensembl ID
convert_gene_name_to_ensembl <- function(gene_name_vector) {
    # 使用 egSYMBOLtoENTREZID 函数转换基因名到Entrez ID
    entrez_id_vector <- org.Hs.egSYMBOLtoENTREZID(gene_name_vector)
    
    # 使用 keys 函数从org.Hs.eg.db数据库中获取Ensembl ID
    ensembl_id_vector <- org.Hs.egENZYMEIDtoENTREZID(entrez_id_vector, "ENSEMBL")
    
    # 返回Ensembl ID向量
    return(ensembl_id_vector)
}
 
# 示例使用
gene_names <- c("TP53", "BRCA1", "BRCA2")
ensembl_ids <- convert_gene_name_to_ensembl(gene_names)
print(ensembl_ids)

这段代码定义了一个函数convert_gene_name_to_ensembl,它接受一个基因名的向量作为输入,并返回对应的Ensembl ID的向量。代码使用了org.Hs.eg.db包中的egSYMBOLtoENTREZIDegENZYMEIDtoENTREZID函数来完成从基因名到Ensembl ID的转换。最后,代码示例展示了如何使用这个函数来转换一个包含几个基因名的向量。

最后修改于:2024年08月17日 20:06

评论已关闭

推荐阅读

DDPG 模型解析,附Pytorch完整代码
2024年11月24日
DQN 模型解析,附Pytorch完整代码
2024年11月24日
AIGC实战——Transformer模型
2024年12月01日
Socket TCP 和 UDP 编程基础(Python)
2024年11月30日
python , tcp , udp
如何使用 ChatGPT 进行学术润色?你需要这些指令
2024年12月01日
AI
最新 Python 调用 OpenAi 详细教程实现问答、图像合成、图像理解、语音合成、语音识别(详细教程)
2024年11月24日
ChatGPT 和 DALL·E 2 配合生成故事绘本
2024年12月01日
omegaconf,一个超强的 Python 库!
2024年11月24日
【视觉AIGC识别】误差特征、人脸伪造检测、其他类型假图检测
2024年12月01日
[超级详细]如何在深度学习训练模型过程中使用 GPU 加速
2024年11月29日
Python 物理引擎pymunk最完整教程
2024年11月27日
MediaPipe 人体姿态与手指关键点检测教程
2024年11月27日
深入了解 Taipy:Python 打造 Web 应用的全面教程
2024年11月26日
基于Transformer的时间序列预测模型
2024年11月25日
Python在金融大数据分析中的AI应用(股价分析、量化交易)实战
2024年11月25日
AIGC Gradio系列学习教程之Components
2024年12月01日
Python3 `asyncio` — 异步 I/O,事件循环和并发工具
2024年11月30日
llama-factory SFT系列教程:大模型在自定义数据集 LoRA 训练与部署
2024年12月01日
Python 多线程和多进程用法
2024年11月24日
Python socket详解,全网最全教程
2024年11月27日