tests/testthat/test_local.minmax.R

In helen-zhu/ConsensusPeaks: A Method for Peak Merging And Distribution Fitting Across Biological Replicates

test_that("local.minmax works on non consecutive numbers", {
  x <- c(1,2,1,4,1)
  x.minmax <- local.minmax(x)
  expect_equal(x.minmax$min.ind, c(1,3,5))
  expect_equal(x.minmax$max.ind, c(2,4))
})


test_that("local.minmax works on consecutive starting numbers", {
  x <- c(1, 1,2,1,4,1)
  x.minmax <- local.minmax(x)
  expect_equal(x.minmax$min.ind, c(2,4,6))
  expect_equal(x.minmax$max.ind, c(3,5))

  x <- c(1,1,1, 1,2,1,4,1)
  x.minmax <- local.minmax(x)
  expect_equal(x.minmax$min.ind, c(4,6,8))
  expect_equal(x.minmax$max.ind, c(5,7))
})

test_that("local.minmax works on consecutive ending numbers", {
  x <- c(1,2,1,4,1, 1)
  x.minmax <- local.minmax(x)
  expect_equal(x.minmax$min.ind, c(1,3,5))
  expect_equal(x.minmax$max.ind, c(2,4))

  x <- c(1,2,1,4,1, 1,1,1)
  x.minmax <- local.minmax(x)
  expect_equal(x.minmax$min.ind, c(1,3,5))
  expect_equal(x.minmax$max.ind, c(2,4))
})

test_that("local.minmax works on consecutive starting and ending numbers", {
  x <- c(1, 1,2,1,4,1, 1)
  x.minmax <- local.minmax(x)
  expect_equal(x.minmax$min.ind, c(2,4,6))
  expect_equal(x.minmax$max.ind, c(3,5))

  x <- c(1,1,1,1, 1,2,1,4,1, 1,1,1)
  x.minmax <- local.minmax(x)
  expect_equal(x.minmax$min.ind, c(5,7,9))
  expect_equal(x.minmax$max.ind, c(6,8))
})

test_that("local.minmax works on even consecutive numbers in middle of sequence", {
  x <- c(1,2,2,1,1,4,4, 1)
  x.minmax <- local.minmax(x)
  expect_equal(x.minmax$min.ind, c(1,4,8))
  expect_equal(x.minmax$max.ind, c(2,6))
})


# NOTE: This may not be the behavior that we want but is one way of dealing with the consecutive sequences
test_that("local.minmax works on odd consecutive numbers in middle of sequence", {
  x <- c(1,2,2,2,1,1,1,4,4,4, 1)
  x.minmax <- local.minmax(x)
  expect_equal(x.minmax$min.ind, c(1,5,11))
  expect_equal(x.minmax$max.ind, c(2,8))

  x <- c(4,2,2,2,1,1,1,4,4,4,5)
  x.minmax <- local.minmax(x)
  expect_equal(x.minmax$min.ind, c(2,5, 10))
  expect_equal(x.minmax$max.ind, c(1,4,8,11))
})