# NOT RUN {
##-------------------------------------##
## A) Ensembl
##-------------------------------------##
ensfile <- system.file("extdata",
"hs.ensembl.62.small.RData", package="refGenome")
ens <- loadGenome(ensfile)
gt <- getGeneTable(ens)
gene_id <- as.character(gt$gene_id[1])
gm <- geneModel(ens, gene_id)
tr <- getTranscript(gm, 1)
tr
getExonData(tr)
getCdsData(tr)
##-------------------------------------##
## B) UCSC
##-------------------------------------##
ucfile <- system.file("extdata", "hs.ucsc.small.RData", package="refGenome")
uc <- loadGenome(ucfile)
gt <- getGeneTable(uc)
gene_id <- as.character(gt$gene_id[1])
gm <- geneModel(uc, gene_id)
tr <- getTranscript(gm, 1)
tr
# }
Run the code above in your browser using DataLab