R语言ENTREZID ID 转换 gene SYMBOL

时间:2020-06-23 22:51:22   收藏:0   阅读:656
library(data.table)
library(org.Hs.eg.db)
library(clusterProfiler)

setwd("example")
 
dt <- read.table("VanAllen.self_subtract", sep = "\t", header = T)
dt$geneid <- rownames(dt)

# transform id  
map_dt <- bitr(dt$geneid, fromType = "ENTREZID",toType = c( "SYMBOL"),OrgDb = org.Hs.eg.db)
dt_merge <- merge(map_dt,dt, by.y = "geneid", by.x = "ENTREZID")
rownames(dt_merge) <- dt_merge$SYMBOL
dt_merge$SYMBOL <- NULL
dt_merge$ENTREZID <- NULL

# write out
write.table(dt_merge, file = "example42.txt", sep = "\t", col.names = T, row.names = T, quote = F)

 

原文:https://www.cnblogs.com/ywliao/p/13184572.html

评论(0
© 2014 bubuko.com 版权所有 - 联系我们:wmxa8@hotmail.com
打开技术之扣,分享程序人生!