pick.best.index: pick best index

In tdhock/bams: Breakpoint annotation model smoothing

Description

Minimizer for local models, described in article section 2.3 "Picking the optimal model"

Usage

pick.best.index(err)

Arguments

err	Vector of errors to minimize.

Value

Integer index of the minimal error.

Author(s)

Toby Dylan Hocking

Examples

stopifnot(pick.best.index(rep(0,100))==50)

err <- rep(1,100)
err[5] <- 0
stopifnot(pick.best.index(err)==5)

## should pick the middle
err <- rep(1,100)
err[40:60] <- 0
stopifnot(pick.best.index(err)==50)

## should pick the biggest
err <- rep(1,100)
err[1:60] <- 0
stopifnot(pick.best.index(err)==60)

## should pick the smallest
err <- rep(1,100)
err[50:100] <- 0
stopifnot(pick.best.index(err)==50)

tdhock/bams documentation built on May 31, 2019, 7:34 a.m.