Calculates a peak cutoff value given an FDR, assuming a Poisson noise distribution estimated from the frequency of singleton and doubleton islands.

1 | ```
peakCutoff(cov, fdr.cutoff = 0.001, k = 2:20)
``` |

`cov` |
The coverage object, e.g.,
an |

`fdr.cutoff` |
The maximum-allowed FDR for calculating the cutoff. |

`k` |
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 |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.