plotWin: Plot the number of reads vs the proportion of '+' stranded...

In UofABioinformaticsHub/strandCheckR: Calculate strandness information of a bam file

Description

Plot the number of reads vs the proportion of '+' stranded reads of all windows from the input data frame.

Usage

plotWin(windows, split = c(10, 100, 1000), threshold = c(0.6, 0.7, 0.8,
  0.9), save = FALSE, file = "win.pdf", groupBy = NULL,
  useCoverage = FALSE, ...)

Arguments

windows	data frame containing the strand information of the sliding windows. Windows should be obtained using the function getStrandFromBamFile to ensure the correct data structure.
split	an integer vector that specifies how you want to partition the windows based on coverage. By default split = c(10,100,1000), partition windows into 4 groups based on these values.
threshold	a numeric vector between 0.5 & 1 that specifies which threshold lines to draw on the plot. The positive windows above the threshold line (or negative windows below the threshold line) will be kept when using filterDNA.
save	if TRUE, then the plot will be save into the file given by file parameter
file	the file name to save to plot
groupBy	the column that will be used to split the data (which will be used in the facets method of ggplot2).
useCoverage	if TRUE then plot the coverage strand information, otherwise plot the number of reads strand information. FALSE by default
...	used to pass parameters to facet_wrap during plotting

Details

This function will plot the proportion of '+' stranded reads for each window, against the number of reads in each window. The threshold lines indicate the hypothetical boundary where windows will contain reads to kept or discarded using the filtering methods of filterDNA. Any plot can be easily modified using standard ggplot2 syntax (see Examples)

Value

The plot will be returned as a standard ggplot2 object

See Also

getStrandFromBamFile, plotHist

Examples

bamfilein = system.file('extdata','s2.sorted.bam',package = 'strandCheckR')
windows <- getStrandFromBamFile(file = bamfilein,sequences = '10')
plotWin(windows)

# Change point colour using ggplot2
library(ggplot2)
plotWin(windows) + 
scale_colour_manual(values = rgb(seq(0, 1, length.out = 4), 0, 0))

UofABioinformaticsHub/strandCheckR documentation built on Aug. 15, 2021, 9:08 a.m.