find_peak: Find peaks

In grafab/gsrc: Genome Structure Rearrangement Calling in Genomes with High Synteny

Description

Find peaks

Usage

find_peak(dat, npeaks = NULL, breaks = round(length(dat) * 0.1),
  check = FALSE, method = "density")

Arguments

dat	Vector, containing theta values for one sample and chromosome.
npeaks	Integer, Number of peaks to be detected.
breaks	Integer, Number of breaks for the histogram.
check	Logical, if TRUE it is checked if the central peak is approximately in the middle of the two other peaks. If the data has not heterozygotes, the third peak might be very close to one of the homozygous peaks. In that case the peak is set to the middle of the other two peaks.
method	Character, "mixture", "density" or "histogram". Defines which method is used to transform the data into a distribution.

grafab/gsrc documentation built on May 17, 2019, 8:18 a.m.