Learn R Programming

R453Plus1Toolbox (version 1.22.0)

homopolymerHist: Create A Histogram Of The Homopolymer Stretches

Description

This function creates a histogram for the different lengths of the homopolymer stretches with one bar for each nucleotide in the flow. With the height giving the number of this homopolymer stretch in the flowgram.

Usage

homopolymerHist(x, range=c(0, length(flowgram(x))), xlab="Homopolymer length", ylab="Number of homopolymers", col=c(A="black", C="red", G="blue", T="green"), ...)

Arguments

x
An object of class SFFRead.
range
Two positions between which the flows should be plotted.
xlab
The X axis label.
ylab
The Y axis label.
col
The colors in which the four nucleotids should be plotted.
...
Arguments to be passed to methods, such as graphical parameters (see ‘par’).