tests/local/test-bnecfit.R

In AIMS/bayesnec: A Bayesian No-Effect- Concentration (NEC) Algorithm

nec_ <- pull_out(manec_example, "nec4param")
ecx_ <- pull_out(manec_example, "ecx4param")

test_that("Adding only works if either if bnecfit", {
  expect_error(1 + nec_)
  expect_error(NA + nec_)
  expect_error(NULL + nec_)
  expect_error("a" + nec_)
  expect_s3_class(ecx_ + nec_, c("bayesmanecfit", "bnecfit")) |>
    expect_message()
  expect_s3_class(manec_example + nec_, c("bayesmanecfit", "bnecfit")) |>
    expect_message() |>
    expect_warning() |>
    expect_warning() |>
    expect_warning() |>
    expect_warning()
  expect_s3_class(manec_example + manec_example,
                  c("bayesmanecfit", "bnecfit")) |>
    expect_message() |>
    expect_warning() |>
    expect_warning() |>
    expect_warning() |>
    expect_warning()
  expect_s3_class(manec_example + nec_ + ecx_,
                  c("bayesmanecfit", "bnecfit")) |>
    expect_message() |>
    expect_message() |>
    expect_warning() |>
    expect_warning() |>
    expect_warning() |>
    expect_warning() |>
    expect_warning() |>
    expect_warning() |>
    expect_warning() |>
    expect_warning()
})

test_that("Concat only works if datasets are identical", {
  expect_message(c(mod_1, mod_2)) # should work
  expect_message(c(mod_1, mod_2, mod_3)) # should work
  expect_error(c(mod_1, mod_2, mod_4)) # should fail
  expect_error(c(mod_1, mod_2, mod_5)) # should fail
  expect_error(c(mod_1, mod_2, mod_6)) # should fail
})

test_that("Concatenating only works if either if bnecfit", {
  expect_true(is.list(c(1, nec_)))
  expect_true(is.list(c(NA, nec_)))
  expect_true(is.list(c(NULL, nec_)))
  expect_true(is.list(c("a", nec_)))
  expect_s3_class(c(ecx_, nec_), c("bayesmanecfit", "bnecfit")) |>
    expect_message()
  expect_s3_class(c(manec_example, nec_), c("bayesmanecfit", "bnecfit")) |>
    expect_message() |>
    expect_warning() |>
    expect_warning() |>
    expect_warning() |>
    expect_warning()
  expect_s3_class(c(manec_example, manec_example),
                  c("bayesmanecfit", "bnecfit")) |>
    expect_message() |>
    expect_warning() |>
    expect_warning() |>
    expect_warning() |>
    expect_warning()
  expect_s3_class(c(manec_example, nec_, ecx_),
                  c("bayesmanecfit", "bnecfit")) |>
    expect_message() |>
    expect_warning() |>
    expect_warning() |>
    expect_warning() |>
    expect_warning()
})

test_that("Update works with regular fitting arguments", {
  # no recompilation
  expect_error(update(manec_example, chains = 1, iter = 50))
})

test_that("Different distribution triggers error or message", {
  expect_error(update(manec_example, newdata = nec_data, recompile = TRUE,
                      chains = 1, iter = 50, family = Beta(link = "identity")))
  expect_error(update(manec_example, newdata = nec_data, recompile = TRUE,
                      chains = 1, iter = 50))
  expect_error(update(manec_example, recompile = TRUE, chains = 1, iter = 50,
                      family = Beta(link = "identity")))
  expect_error(update(manec_example, newdata = nec_data, recompile = TRUE,
                      chains = 2, iter = 50, family = Beta(link = "identity")),
               "You either input new")
  expect_error(update(manec_example, newdata = nec_data, recompile = TRUE,
                      chains = 2, iter = 50), "You either input new")
  expect_error(update(manec_example, recompile = TRUE, chains = 2, iter = 50,
                      family = Beta(link = "identity")),
               "You either input new")
})

test_that("Update works with regular fitting arguments", {
  expect_s3_class(update(nec_, chains = 2, iter = 50, recompile = TRUE,
                         refresh = 0, verbose = FALSE), "bnecfit") |>
    suppressWarnings() |>
    suppressMessages()
})

test_that("Different distribution triggers error or message", {
  expect_error(update(manec_example, newdata = nec_data, recompile = TRUE,
                      chains = 2, iter = 50, family = Beta(link = "identity"),
                      force_fit = TRUE, refresh = 0, verbose = FALSE)) |>
    suppressWarnings() |>
    suppressMessages()
})