1. OUQ技术分享首页
  2. 生物技术
  3. 生物信息技术
  4. R

R:gene symbol与id转换

• R

利用R将gene symbol与id转换

安装所需包

if (!requireNamespace(“BiocManager”, quietly = TRUE))
install.packages(“BiocManager”)

BiocManager::install("org.Hs.eg.db")

R:gene symbol与id转换

运行代码

library("org.Hs.eg.db") #加载人源包，鼠org.Mm.eg.db
rt=read.table("symbol.txt",sep="\t",check.names=F,header=T) #读取文件
genes=as.vector(rt[,1])
entrezIDs <- mget(genes, org.Hs.egSYMBOL2EG, ifnotfound=NA) #找出基因对应的id
entrezIDs <- as.character(entrezIDs)
out=cbind(rt,entrezID=entrezIDs)
write.table(out,file="id.txt",sep="\t",quote=F,row.names=F) #输出结果

本站原创，如若转载，请注明出处：https://www.ouq.net/110.html

赞 (1)
打赏 微信打赏，为服务器增加50M流量 微信打赏，为服务器增加50M流量 支付宝打赏，为服务器增加50M流量 支付宝打赏，为服务器增加50M流量
上一篇 03/16/2020 00:41
下一篇 03/17/2020 00:37

相关推荐