Hello!

1)I have created TxDb  using makeTxDbFromGFF.

2)Then I want to extract only transcript with length >100.

3)Then I want to export.gff3 to save a new gff3 file utilizing rtracklayer. 

How to do the second step with length parameter???

Thank you in advance,

Alex

Extract the transcripts from the TxDb and subset for the width you are interested in

> tx = transcripts(TxDb.Hsapiens.UCSC.hg38.knownGene)
> length(tx)
[1] 197782
> tx100 = tx[width(tx) > 100]
> length(tx100)
[1] 192874