symbol ID转换

rm(list = ls())
read.table('a.txt')
library(org.Hs.eg.db)
g2s=toTable(org.Hs.egSYMBOL)
g2e=toTable(org.Hs.egENSEMBL)
library(stringr)
a$ensembl_id=unlist(lapply(a$V1, function(x){strsplit(as.character(x),'[.]')[[1]][1]}))
b=merge(a,g2e,by='ensembl_id')
d=merge(b,g2s,by='gene_id')
View(d)
©著作权归作者所有,转载或内容合作请联系作者
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

推荐阅读更多精彩内容