gsrc (version 1.1)

find_peak: Find peaks

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.