inst/tests/test-edges.R

In lianos/biosignals: Smoothing and peak/edge finding in 1D data (signal vectors).

context("edge detection")

## This example is taken from the supplement of the FACADE paper:
## http://nar.oxfordjournals.org/content/38/15/e157.full
test_that("idealized step edge detected", {
  ## Idealized step
  x <- c(rep(0, 15), rep(20, 15), -1, rep(0, 20))
  x <- x + rnorm(length(x)) ## with some noise

  edges <- detectEdges(x, bandwidth=15)
  
  ## Give some slack on what the start/end calls are given the randomness
  ## in the noise added to x, the real start/end should be 16 and 31
  expect_true(any(c(15:17) %in% edges$start))
  expect_true(any(c(30:32) %in% edges$end))
})