tests/testthat/test-logitnavg.R

In kaybrodersen/micp: Bayesian mixed-effects inference on classification performance

test_that("logitnavgcdf rejects bad input", {
  expect_error(logitnavgcdf("foo", 0, 1, 0, 1),
               "not a numeric or integer vector")
  expect_error(logitnsumcdf(matrix(c(0, 0.1, 0.2, 0.3), nrow = 2), 0, 1, 0, 1),
               "not an atomic vector")
})

test_that("logitnavgcdf returns NA for NA input", {
  expect_equal(logitnavgcdf(NA_real_, 0, 1, 0, 1), NA_real_)
  expect_equal(logitnavgcdf(NA, 0, 1, 0, 1), NA_real_)
  expect_equal(logitnavgcdf(c(NA, NA), 0, 1, 0, 1), c(NA_real_, NA_real_))
})

test_that("logitnavgcdf returns expected result", {
  expect_equal(logitnavgcdf(c(0.8, 0.9), 0.1, 1.3, -0.7, 1.1),
               logitnsumcdf(c(1.6, 1.8), 0.1, 1.3, -0.7, 1.1))
  expect_equal(logitnavgcdf(0, 0, 1, 0, 1), 0)
  expect_equal(round(logitnavgcdf(0.5, 0, 1, 0, 1), 2), 0.5)
  expect_equal(logitnavgcdf(1, 0, 1, 0, 1), 1)
})

test_that("logitnavginv rejects bad input", {
  expect_error(logitnavginv("foo", 0, 1, 0, 1),
               "not a numeric or integer vector")
})

test_that("logitnavginv returns NA for NA input", {
  expect_equal(logitnavginv(NA_real_, 0, 1, 0, 1), NA_real_)
  expect_equal(logitnavginv(NA, 0, 1, 0, 1), NA_real_)
  expect_equal(logitnavginv(c(NA, NA), 0, 1, 0, 1), c(NA_real_, NA_real_))
})

test_that("logitnavginv is the inverse of logitnavgcdf", {
  q <- logitnavgcdf(0.6, 2.2, 1.3, -1.7, 1.1)
  expect_equal(logitnavginv(q, 2.2, 1.3, -1.7, 1.1), 0.6, tolerance = 1e-4)
})

test_that("logitnavginv works on vector input", {
  q1 <- logitnavgcdf(0.6, 2.2, 1.3, -1.7, 1.1)
  q2 <- logitnavgcdf(0.7, 2.3, 1.4, -1.8, 1.2)
  expect_equal(logitnavginv(c(q1, q2), c(2.2, 2.3), c(1.3, 1.4),
                            c(-1.7, -1.8), c(1.1, 1.2)), c(0.6, 0.7),
               tolerance = 1e-4)
})

test_that("logitnavgmean result is approximately the sample mean", {
  set.seed(1)
  x1 <- safesigm(rnorm(10000, 2, 1))
  x2 <- safesigm(rnorm(10000, 3, 2))
  sample_mean <- mean((x1 + x2) / 2)
  theoretical_mean <- logitnavgmean(2, 1, 3, 2)
  expect_equal(theoretical_mean, sample_mean, tolerance = 1e-3)
})

test_that("logitnavgmean works on vector input", {
  expect_equal(logitnavgmean(c(1, 1.1), c(2, 2.2), c(3, 3.3), c(4, 4.4)),
               c(logitnavgmean(1, 2, 3, 4), logitnavgmean(1.1, 2.2, 3.3, 4.4)))
})