tests/testthat/test-distributions.R

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

context("Distributions")

test_that("dPGIB returns positive values", {
  expect_that(all(dPGIB(0:2,3,4,2,2.5)>0), is_true())
})

test_that("Check moment GB2", {
  a <- 3; c <- 4; d <- 5; tau <- 10
  expect_that((moment_GB2(1, a, c, d, tau) - summary_gamma(a, 1)$mean*summary_beta2(d, c, 1/tau)$mean) < .Machine$double.eps, 
              is_true())
})

test_that("Check moments GIB and PGIB", {
  a <- 3; alpha <- 4; beta <- 5; rho <- 10
  expect_equal(summary_GIB(a,alpha,beta,rho)$mean, summary_PGIB(a,alpha,beta,rho)$mean)
  expect_equal(summary_GIB(a,alpha,beta,rho)$mean+summary_GIB(a,alpha,beta,rho)$sd^2, summary_PGIB(a,alpha,beta,rho)$sd^2)
})