Entering edit mode
8.3 years ago
bsmith030465
▴
240
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!
This is a user guide:
http://bioconductor.org/packages/devel/bioc/vignettes/Gviz/inst/doc/Gviz.pdf
You need a transformation between UCSC gene_ids and ensemble gene_ids, am I right?
There is a post about it:
Ensembl And Ucsc Gtfs To Match
That's what I'm trying to find out. How was the dataset (used in Gviz) created.
thanks.