Description Usage Arguments Value Examples
plotBinCounts
visualizes counts or enrichment for a particular region across bins, datasets, and read groups.
Many graphical parameters can be adjusted to suit the needs of the user, see below.
1 2 3 4 5 |
SNPhood.o |
Object of class |
regions |
Numeric or Character or NULL. Default NULL. Regions that should be plotted, either specified as integer (such as 1, value must be
between 1 and the total number of regions as defined in the object) or their annotation (name must appear in the region names as obtained via the function |
readGroups |
Character or NULL. Default NULL. Read groups that should be plotted, specified by their name as obtained by the function |
datasets |
Numeric or Character or NULL. Default NULL. Datasets that should be used for plotting, either specified as integer (such as 1, value must be
between 1 and the total number of datasets as defined in the object) or their annotation (name must appear in the dataset names as obtained via the function |
readGroupColors |
Character or NULL. Default NULL. Colors of the read groups that appear in the plot.
If set to NULL, colors will be set automatically. The length of the vector must equal the total number of read groups
that are defined in the |
ylim |
Numeric(2). Default NULL. Range of the y-axis, as specified by a minimum and a maximum value. See ?ylim for details. |
addGenotype |
Logical(1). Default TRUE. Should the genotype distribution for each read group at the original user position be displayed in the legend in addition? See the Vignette for more details how this distribution is determined. |
plotGenotypeRatio |
Logical(1). Default FALSE. Should the ratio of the genotypes be plotted instead of the count or enrichment values? Only applicable if the number of read groups to be plotted is 2 and if one region is plotted. Setting this parameter to TRUE may result in ratios across bins that are interrupted due to zero counts (and a resulting division by zero, which can therefore not be displayed). Also, ratios cannot be plotted if the genotype for the selected regions could not be determined due to the lack of reads overlapping with the particular region (see the Vignette for details). |
addTitle |
Logical(1). Default TRUE. Should the plot contain a title that summarizes the genomic region that is visualized? |
colorPalette |
Character(1). Default "Set1". Name of the palette from the |
printPlot |
Logical(1). Default TRUE. Should the plots be printed? Only relevant if |
fileToPlot |
Character(1) or |
maxWidthLabels |
Numeric(1). Default NULL. Maximum width of the legend labels in number of characters. If the width of the legend labels are longer, they are shortened. Set to NULL to not shorten labels. |
verbose |
Logical(1). Default FALSE. Should the verbose mode (i.e., diagnostic messages during execution of the script) be enabled? |
the generated ggplot2 plot(s) as list for further processing. May contain multiple plots, depending on the function. The plot(s) can then be plotted individually or modified arbitrarily as the user wants. For example, if multiple plots are returned and the plots have been saved in a variable called plots.l, simply type plots.l[[1]] to view the first plot.
1 2 3 4 5 6 7 8 9 10 11 12 13 | data(SNPhood.o, package="SNPhood")
# Plot the first region, all parameters with their default values
plot = plotBinCounts(SNPhood.o)
# Plot the second region for the first dataset, using specific colors for the read groups.
plot = plotBinCounts(SNPhood.o, regions = 2, dataset = 1, readGroupColors = c("red","blue","gray"))
# Plot the first region for the first dataset and the genotype ratio. Save the plot in a variable
plot = plotBinCounts(SNPhood.o, regions = 1, readGroups = c("maternal", "paternal"), dataset = 1, plotGenotypeRatio = TRUE)
#' # Plot all regions for the first dataset and aggregate. Save the plot in a variable
plot = plotBinCounts(SNPhood.o, regions = NULL, readGroups = c("maternal", "paternal"), dataset = 1)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.