# NOT RUN {
##-------------------------------------##
## A) Ensembl
##-------------------------------------##
ensfile <- system.file("extdata",
"hs.ensembl.62.small.RData", package="refGenome")
ens <- loadGenome(ensfile)
enex <- refExons(ens)
# }
# NOT RUN {
saveGenome(enex,"enex.RData", useBasedir=FALSE)
er <- loadGenome("enex.RData")
# }
# NOT RUN {
# Unify exons
uex <- unifyRanges(enex)
##-------------------------------------##
## B) UCSC
##-------------------------------------##
ucfile <- system.file("extdata", "hs.ucsc.small.RData", package="refGenome")
uc <- loadGenome(ucfile)
ucex <- refExons(uc)
# }
# NOT RUN {
saveGenome(ucex, "ucex.RData", useBasedir=FALSE)
ur <- loadGenome("ucex.RData")
# }
Run the code above in your browser using DataLab