tests/testthat/test-quantiles.R

In brr: Bayesian Inference on the Ratio of Two Poisson Rates

context("Quantiles")

test_that("Quantiles PGIB and PGB2", {
  p <- 0.5
  q <- qPGIB(p, 10, 2, 12, 2.1) 
  expect_true(pPGIB(q, 10, 2, 12, 2.1) > p && pPGIB(q-1, 10, 2, 12, 2.1) < p)
  q <- qPGB2(p, 10, 2, 12, 2.1) 
  expect_true(pPGB2(q, 10, 2, 12, 2.1) > p && pPGB2(q-1, 10, 2, 12, 2.1) < p)
})

test_that("Quantiles dprior_x", {
  p0 <- 0.5
  q <- qprior_x(p0, 2, 3, 4, 5, 10)
  p <- pprior_x(q, 2, 3, 4, 5, 10)
  pprevious <- pprior_x(q-1, 2, 3, 4, 5, 10)
  expect_true(p > p0 && pprevious < p0)
})