美文网首页RNAseq
计算基因长度

计算基因长度

作者: 牛顿大尉 | 来源:发表于2020-12-10 11:18 被阅读0次

用在转录组RPKM计算上
1载入GTF文件

library(GenomicFeatures)
txdb <- makeTxDbFromGFF("hg38.gtf",format="gtf")

2通过exonsBy获取每个gene上的所有外显子的起始位点和终止位点,然后用reduce去除掉重叠冗余的部分,最后计算长度

exons_gene <- exonsBy(txdb, by = "gene")
exons_gene_lens <- lapply(exons_gene,function(x){sum(width(reduce(x)))})

相关文章

网友评论

    本文标题:计算基因长度

    本文链接:https://www.haomeiwen.com/subject/covlgktx.html