Generate plots to visualy assess the quality of select colonies
or sequencing cycles
Usage
## S4 method for signature 'RolexaRun'
CombinedPlot(run=Rolexa.env, int, seq, scores, colonies = 1:4, par = list())
CombinedPlot(run,...)
## S4 method for signature 'SolexaIntensity'
ChannelHistogram(int, cycles = c(1,18,36),
threemodes = FALSE, par = list())
ChannelHistogram(int,...)
Arguments
run
a RolexaRun object defining the run parameters
int
a SolexaIntensity
object
seq
a ShortRead
object
scores
a matrix of base quality scores (one column per base,
one row per sequence)
cycles
the list of cycles to plot
colonies
the list of rows to select for plotting
threemodes
fit and plot a mixture of 3 gaussians (2 by default)
par
parameters for the plotting function
...
additional arguments, ignored
Details
CombinedPlot creates one plot for each selected colony
with the sequence along the x axis, the four intensities
plotted as barplots above each base and the quality scores as a line
plot below the sequence.
ChannelHistogram plots histograms and signal-noise
thresholds for each of the four intensity channels on selected
cycles. Fits to 2 or 3 gaussians are overlaid on the histograms.
Author(s)
Jacques Rougemont, Arnaud Amzallag, Christian Iseli, Laurent Farinelli, Ioannis Xenarios, Felix Naef
References
Probabilistic base calling of Solexa sequencing data,
BMC Bioinformatics 2008, 9:431
R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(Rolexa)
Error in library(Rolexa) : there is no package called 'Rolexa'
Execution halted