从gff3文件提取基因长度

library(GenomicFeatures)
txdb <- makeTxDbFromGFF("../genome/me.gff3",format="auto")
exons_gene <- exonsBy(txdb, by = "gene")
exons_gene_lens <- lapply(exons_gene,function(x){sum(width(reduce(x)))})
length=t(as.data.frame(exons_gene_lens))
write.table(length,'../genome/gene_length.txt',col.names=F,row.names=T,quote=F,sep='\t')

---

💌lixiang117423@foxmail.com
💌lixiang117423@gmail.com