Calculate a peak cutoff
Calculates a peak cutoff value given an FDR, assuming a Poisson noise distribution estimated from the frequency of singleton and doubleton islands.
peakCutoff(cov, fdr.cutoff = 0.001, k = 2:20)
The coverage object, e.g.,
The maximum-allowed FDR for calculating the cutoff.
The coverage levels at which to estimate an FDR value. The maximal
value that is less than
A numeric value to use for calling peaks
D. Sarkar and M. Lawrence
islandDepthPlot for the graphical equivalent; the
vignette for a bit more explanation.
1 2 3