histogram: Plot histogram of STR sequences

In ShawHahnLab/microsat: Computational, High-throughput Individual Identification through Microsatellite Profiling

Plot histogram of STR sequences

Description

Given processed STR sequences and optionally the filtered version, plot a histogram of counts per sequence length.

Usage

histogram(
  seq_data,
  sample_data = NULL,
  main = "Number of Reads by Sequence Length",
  xlim = range(seq_data$Length),
  cutoff_fraction = NULL,
  render = TRUE
)

Arguments

seq_data	data frame of dereplicated sequences as created by analyze_seqs.
sample_data	data frame of filtered and categorized sequences as created by analyze_sample.
main	title of the plot.
xlim	numeric range for x-axis.
cutoff_fraction	numeric threshold for the fraction of locus-matching counts needed to call an allele. Used to draw a horizontal line if sample_data is given.
render	Should the plot be drawn to the display device?

Value

list of data frames for the sets of counts-versus-length bars drawn in the plot, split by category.

ShawHahnLab/microsat documentation built on Aug. 25, 2023, 11:16 p.m.