test_that("The simple lnVR is computed correctly", {
  C_sd <- 1.5
  C_n <- 10
  A_sd <- 2.5
  A_n <- 10

  test_lnVR <- log(A_sd / C_sd) + (1 / (2 * (A_n - 1))) - (1 / (2 * (C_n - 1)))
  test_lnVRv <- (1 / (2 * (A_n - 1))) + (1 / (2 * (C_n - 1)))

  res <- .simple_lnVR(
    Ctrl_sd = C_sd, 
    Ctrl_n = C_n,
    A_sd = A_sd,
    A_n = A_n
  )

  expect_equal(res$simple_lnVR, test_lnVR, tolerance = 1e-6)
  expect_equal(res$simple_lnVRv, test_lnVRv, tolerance = 1e-6)
})


test_that("The main lnVR is computed correctly", {
  df <- data.frame(
    Ctrl_sd = 2.0,
    Ctrl_n = 12,
    A_sd = 2.8,
    A_n = 12,
    B_sd = 2.4,
    B_n = 12,
    AB_sd = 3.6,
    AB_n = 12
  )

  test_main_lnVR <- 0.37096
  test_main_lnVRv <- 0.045454

  res <- with(df, .main_lnVR(
    Ctrl_sd = Ctrl_sd,
    Ctrl_n = Ctrl_n,
    A_sd = A_sd,
    A_n = A_n,
    B_sd = B_sd,
    B_n = B_n,
    AB_sd = AB_sd,
    AB_n = AB_n
  ))

  expect_equal(res$main_lnVR, test_main_lnVR, tolerance = 1e-4)
  expect_equal(res$main_lnVRv, test_main_lnVRv, tolerance = 1e-4)
})


test_that("The interaction lnVR is computed correctly", {
  df <- data.frame(
    Ctrl_sd = 2.0,
    Ctrl_n = 12,
    A_sd = 2.8,
    A_n = 12,
    B_sd = 2.4,
    B_n = 12,
    AB_sd = 3.6,
    AB_n = 12
  )

  test_inter_lnVR <- 0.06899
  test_inter_lnVRv <- 0.18181

  res <- with(df, .interaction_lnVR(
    Ctrl_sd = Ctrl_sd,
    Ctrl_n = Ctrl_n,
    A_sd = A_sd,
    A_n = A_n,
    B_sd = B_sd,
    B_n = B_n,
    AB_sd = AB_sd,
    AB_n = AB_n
  ))

  expect_equal(res$inter_lnVR, test_inter_lnVR, tolerance = 1e-4)
  expect_equal(res$inter_lnVRv, test_inter_lnVRv, tolerance = 1e-4)
})