tests/testthat/test_quantileFunctions.R

In TauStar: Efficient Computation and Testing of the Bergsma-Dassios Sign Covariance

library(TauStar)
context("Testing the quantile functions.")

test_that("qHoeffInd function returns correct values", {
  set.seed(9234)
  xVals <- rHoeffInd(5)
  for (i in 1:length(xVals)) {
    p <- pHoeffInd(xVals[i])
    expect_equal(qHoeffInd(p, 10^-3), xVals[i], tolerance = 10^-3)
  }

  expect_error(qHoeffInd(-.01))
  expect_error(qHoeffInd(1.1))
})

test_that("qDisHoeffInd function returns correct values", {
  set.seed(9234)
  sims <- 10
  for (i in 1:sims) {
    p <- runif(sample(1:5, 1))
    p <- p / sum(p)
    q <- runif(sample(1:5, 1))
    q <- q / sum(q)
    x <- rDisHoeffInd(1, p, q)
    expect_equal(qDisHoeffInd(pDisHoeffInd(x, p, q), p, q, 10^-4),
      x,
      tolerance = 10^-3
    )
  }

  expect_error(qDisHoeffInd(-.01, 1, 1))
  expect_error(qDisHoeffInd(1.1, 1, 1))
})

test_that("qMixHoeffInd function returns correct values", {
  set.seed(9234)
  sims <- 10
  for (i in 1:sims) {
    p <- runif(sample(1:5, 1))
    p <- p / sum(p)
    x <- rMixHoeffInd(1, p)
    expect_equal(qMixHoeffInd(pMixHoeffInd(x, p), p, 10^-4),
      x,
      tolerance = 10^-3
    )
  }

  expect_error(qDisHoeffInd(-.01, 1, 1))
  expect_error(qDisHoeffInd(1.1, 1, 1))
})