前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >TxDb/org.Hs.eg.db 使用方法

TxDb/org.Hs.eg.db 使用方法

原创
作者头像
用户9781992
发布2022-05-27 16:39:21
1.1K0
发布2022-05-27 16:39:21
举报
文章被收录于专栏:R学习R学习
代码语言:javascript
复制
library("TxDb.Hsapiens.UCSC.hg19.knownGene")
library("GenomicFeatures")
#org.Hs.eg.db数据对象里面包含着各大主流数据库的数据,一般人都比较熟悉的entrez ID 和ensembl 数据库的ID。
library("org.Hs.eg.db")



txdb <- TxDb.Hsapiens.UCSC.hg19.knownGene

gene=genes(txdb) #提取23056个基因信息,数据以Grangs格式显示。
exon=exons(txdb) #提取exons信息,共289969个exon。
trans=transcripts(txdb) #提取转录本信息,共82960个转录本。
cds=cds(txdb) #获取cds区域信息,提取到237533个cds信息。

Tran_gene=transcriptsBy(txdb,by="gene") #通过基因分组获取每个基因的转录本信息。分成了23459个元素的list。
Exon_gene=exonsBy(txdb,by="gene") #基于基因进行外显子区域信息获取。
Cds_gene=cdsBy(txdb,by="gene") #基于基因的CDS区域信息获取。

gene1=mcols(gene) #mcols 函数可以查看GRanges object右边的metadata column

原创声明:本文系作者授权腾讯云开发者社区发表,未经许可,不得转载。

如有侵权,请联系 cloudcommunity@tencent.com 删除。

原创声明:本文系作者授权腾讯云开发者社区发表,未经许可,不得转载。

如有侵权,请联系 cloudcommunity@tencent.com 删除。

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档