tests/testthat/test-kthPair.R

In robcp: Robust Change-Point Tests

context("kth Pair")

test_that("Output of kthPair and medianDiff is of the correct format", 
{
  skip_on_os(os = "solaris")
  
  expect_error(kthPair(1:3, 1:3, 0))
  expect_error(kthPair(1:3, 1:3, 10))
  
  expect_error(kthPair(1:3, 1:3, 1, 3))
  expect_error(kthPair(1:3, 1:3, 0, 3))
  expect_error(kthPair(1:3, 1:3, 1, 0))
  
  ## medianDiff
})

test_that("kthPair returns the correct value",
{
  skip_on_os(os = "solaris")
  
  x <- rnorm(10)
  y <- runif(10)
  
  res1 <- sapply(1:100, function(k) kthPair(x, y, k))
  res2 <- sort(apply(expand.grid(x, y), 1, sum), decreasing = TRUE)
  
  expect_equal(res1, res2)
  
  res1 <- sapply(1:99, function(k) kthPair(x, y, k, k+1))
  res2 <- filter(res2, c(0.5, 0.5))[-100]
  
  expect_equal(res1, res2)
  
  x <- c(1, 2, 3)
  kthPair(x, -x, 1)
  expect_equal(x, c(1, 2, 3))
})

test_that("medianDiff returns the correct value",
{
  skip_on_os(os = "solaris")
  
  x <- sample(1:100, replace = TRUE)
  
  expect_equal(medianDiff(x, 1), median(x - 1))
  expect_equal(medianDiff(1, x), median(1 - x))
  
  y <- rnorm(100)
  
  expect_equal(medianDiff(x, y), median(apply(expand.grid(x, -y), 1, sum)))
  
  x <- runif(99)
  
  expect_equal(medianDiff(x, 1), median(x - 1))
  expect_equal(medianDiff(1, x), median(1 - x))
  
  y <- rnorm(99)
  
  expect_equal(medianDiff(x, y), median(apply(expand.grid(x, -y), 1, sum)))
  
  
  x <- runif(3)
  y <- rnorm(47)
  
  res1 <- medianDiff(x, y)
  res2 <- median(apply(expand.grid(x, -y), 1, sum))
  
  expect_equal(res1, res2)
  
  res1 <- medianDiff(y, x)
  res2 <- median(apply(expand.grid(y, -x), 1, sum))
  
  expect_equal(res1, res2)
})