plotSpliceSum: Plot Splice Summary from RNA-seq data

Description

Plot splice summary by simply counting overlaped junction read in weighted way or not.

Usage

## For character,GRangesList
"plotSpliceSum"(data, model, ..., weighted = TRUE)
## For character,TxDb
"plotSpliceSum"(data, model, which, ..., weighted = TRUE)

Arguments

data

A character specifying the bam file path of RNA-seq data.

model

A GRangesList which represting different isoforms, or a TxDb object. In the second case, users need to pass "which" argument which is a GRanges object to specify the region. And we get connonical model internally.

which

A GRanges object specifying the region you want to get model from the TxDb object.

weighted

If TRUE, weighted by simply add 1/cases matched to each model and if FALSE, simply add 1 to every case.

...

Extra arugments passed to qplot function. such as, offset which control the height of chevron.

Value

A ggplot object.

Details

Internally we use biovizBase:::spliceSummary for simple counting, but we encourage users to use their own robust way to make slicing summary and store it as GRangesList, then plot the summary by qplot function.

Examples

Run this code

## 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)