用org.HS.eg.db进行SYMBOL\GENENAME\ENTREZID的转换

######################entrezID转换##################################
#############教程:B站不吃糖的康康《ensembl_ID转换》################

library(org.Hs.eg.db)

keytypes(org.Hs.eg.db)#####查看org.Hs.eg.db中可以被选择/使用的类型
setwd("C:/Users/HUAWEI/Desktop/Rdata/bss/GSEA")
gsea <- read.csv("gsea.csv",header=T)
X <- gsea$ensembleID######第一列就是ensembleID,把他赋值给一个X
get <- c("SYMBOL","GENENAME","ENTREZID")
result <- select(org.Hs.eg.db,
                 keys = X,
                 columns = get,
                 keytype = "ENSEMBL")
                 
                 ###select选择一个函数(org函数)
                 ####keys=x是表示数据来源(来源于x,x是刚刚赋值的gsea$ensembleID)
                 ###column是需要获得的类型(需要获得get,get <- c("SYMBOL","GENENAME""ENTREZID")
                 ###keytype是输入的类型(ENSEMBL)
                 #####【'select()' returned 1:many mapping between keys and columns】不是一一映射的,是有重合的。
result <- na.omit(result)###去除result里面的NA值(有些在org包里匹配不到的)
table(duplicated(result$ENTREZID))#####有重复,原因是好几个探针覆盖同一个基因。可以看一下多少重复。

让我们谢谢b站~~~~!

感谢bilibili大学。

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值