Hi,

I was trying to make a geneModel object, as used in the Gviz package:

> data(geneModels)

head(geneModels,2)

chromosome start end width strand feature gene exon

1 chr7 26591441 26591829 389 + lincRNA ENSG00000233760 ENSE00001693369

2 chr7 26591458 26591829 372 + lincRNA ENSG00000233760 ENSE00001596777

transcript symbol

1 ENST00000420912 AC004947.2

2 ENST00000457000 AC004947.2

How can I get my data in this format using TxDb.Hsapiens.UCSC.hg38.knownGene, or some other appropriate package?

thanks!