tests/testthat/test.hdi.R

In genphen: A tool for quantification of associations between genotypes and phenotypes in genome wide association studies (GWAS) with Bayesian inference and statistical learning

context("Test HDI")


test_that("getHdi output is numeric vector of length 2", {
  expect_is(class = "numeric",
            object = getHdi(vec = runif(n = 10^3, min = 0, max = 1), 
                            hdi.level = 0.99))
  expect_equal(object = length(getHdi(vec = runif(n = 10^3, min = 0, max = 1), 
                                      hdi.level = 0.99)), 
               expected = 2)
  expect_equal(object = getHdi(vec = runif(n = 10^3, min = 0, max = 1), 
                               hdi.level = 0.99), 
               expected = c(0, 1), c(0, 1)+0.05)
})