tests/testthat/test-numeric_utils.R

In OncoBayes2: Bayesian Logistic Regression for Oncology Dose-Escalation Trials

context("numeric utils")

test_that("log_inv_logit_fast evaluates correctly", {
  l <- seq(-20, 20, length = 41)
  ob2 <- OncoBayes2:::log_inv_logit_fast(l)

  bfam <- binomial()
  inv_logit <- bfam$linkinv

  expect_equal(ob2, log(inv_logit(l)), tolerance = 1E-8)
})

test_that("log1m_exp_max0_fast evaluates correctly", {
  l <- seq(-40, 1, length = 41)
  ob2 <- OncoBayes2:::log1m_exp_max0_fast(l)

  suppressWarnings(ref <- log1p(-exp(l)))

  expect_equal(ob2, ref, tolerance = 1E-8)
})