library(TxDb.Hsapiens.UCSC.hg38.knownGene)
library(org.Hs.eg.db)
txdb <- TxDb.Hsapiens.UCSC.hg38.knownGene
gene_gr <- genes(txdb)
promoter_gr <- promoters(gene_gr)
promoter_gr$gene_symbol <- mapIds(x = org.Hs.eg.db,
keys = promoter_gr$gene_id,
keytype = "ENTREZID", #需要转换的类型
column = "SYMBOL") #需要转换为的类型
head(promoter_gr)
网友评论