tests/testthat/test-bootstrap.R

In BSW: Fitting a Log-Binomial Model Using the Bekhit–Schöpe–Wagenpfeil (BSW) Algorithm

test_that("Bootstrap process runs successfully", {

  # Create example data
  set.seed(123)
  x <- rnorm(100, 50, 10)
  y <- rbinom(100, 1, exp(-4 + x * 0.04))

  # Fit a model
  fit <- bsw(formula = y ~ x, data = data.frame(y = y, x = x))

  # Run the bootstrap procedure
  result <- bootbsw(fit, ci_level = 0.95, R = 1000L)

  # Extract the bootstrap object
  bootobj <- result$Bootstrap_Object

  # Check that the original estimates were computed
  expect_true(is.numeric(bootobj$t0))

  # Check that 1000 bootstrap replicates were performed
  expect_equal(length(bootobj$t[,1]), 1000)

  # Check that at least 95% of bootstrap replicates are complete
  expect_gt(sum(complete.cases(bootobj$t)), 950)
})