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)
代码库8-获得hg38启动子并添加gene symbol
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
- 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
推荐阅读更多精彩内容
- 代码移动到 HG 管理后 .gitignore文件不能用了然后创建.hgignore文件 忽略不了下边的这个每次都...
- 作者:童蒙编辑:angelica 由于基因一般有别名、或者因为excel的转换、人为错误等其他原因,使得基因名容易...
- 由于公司的代码库并不是使用 git 而是使用 Mercurial,网上常用的 jenkins 和 hg 配置的东西...
- Swift代码库之storyboard添加tabbar 创建一个项目,swift_tabbar_storyboar...