## Not run:
# bamfile <- system.file("extdata", "SRR027894subRBM17.bam", package="biovizBase")
# library(TxDb.Hsapiens.UCSC.hg19.knownGene)
# txdb <- TxDb.Hsapiens.UCSC.hg19.knownGene
# data(genesymbol)
# exons <- exonsBy(txdb, by = "tx")
# exons.rbm17 <- subsetByOverlaps(exons, genesymbol["RBM17"])
# plotSpliceSum(bamfile, exons.rbm17)
# plotSpliceSum(bamfile, exons.rbm17, weighted = FALSE, offset = 0.01)
# plotSpliceSum(bamfile, txdb, which = genesymbol["RBM17"])
# plotSpliceSum(bamfile, txdb, which = genesymbol["RBM17"], offset = 0.01)
# plotSpliceSum(bamfile, txdb, which = genesymbol["RBM17"],
# show.label = TRUE,
# label.type = "count")
# ## End(Not run)
Run the code above in your browser using DataLab