在线时间:8:00-16:00
迪恩网络APP
随时随地掌握行业动态
扫描二维码
关注迪恩网络微信公众号
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)
|
请发表评论