R:gene symbol与id转换

利用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/rgene-symbol%e4%b8%8eid%e8%bd%ac%e6%8d%a2.html

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

相关推荐