test_that("Dependence measures work accurately", {
  dep1 <- dep_bitl(xi1 = 1.5, xi2 = 1.2, delta = 0.6)
  expect_s3_class(dep1, "dep_bitl")
  expect_equal(dep1$kendall_tau, (2/9) * 0.6, tolerance = 1e-10)
  expect_equal(dep1$spearman_rho, 0.6 / 3, tolerance = 1e-10)
  expect_equal(dep1$blomqvist_beta, 0.6 / 4, tolerance = 1e-10)
  expect_match(dep1$orthant_dependence, "NOD")

  dep2 <- dep_bitl(xi1 = 1.5, xi2 = 1.2, delta = -0.4)
  expect_match(dep2$orthant_dependence, "POD")

  dep3 <- dep_bitl(xi1 = 1.5, xi2 = 1.2, delta = 0)
  expect_match(dep3$orthant_dependence, "Independent")

  dep_loc <- dep_bitl(xi1 = 1.5, xi2 = 1.2, delta = 0.5, x = c(0.2, 0.5), y = c(0.3, 0.6))
  expect_false(is.null(dep_loc$local_dependence))
  expect_true(all(dep_loc$local_dependence$holland_wang > 0))
})