tests/testthat/test-racusum_betabinomial_crit_sim.R

In vlad: Variable Life Adjusted Display and Other Risk-Adjusted Quality Control Charts

context("racusum_betabinomial_crit_sim")

set.seed(1234)
L0 <- 100
RQ <- 1
RA <- 2
maxS <- 71
g0 <- -3.6798
g1 <- 0.0768
shape1 <- 1
shape2 <- 3
tol <- 0.3
m <- 1e3

test_that("Different input values for RA", {
  RAtest <- list(-1, 0, "0", NA)
  lapply(RAtest, function(x) {
    expect_error(do.call(x, racusum_betabinomial_crit_sim, L0=L0, shape1=shape1, shape2=shape2, coeff=c(g0, g1),
                                             rs = maxS, RQ = RQ, m=m, RA = x))})
})

test_that("Different simulation algorithms, detecting deterioration", {
  skip_on_cran()
  skip_if(SKIP == TRUE, "skip this test now")
  expect_equal(racusum_betabinomial_crit_sim(L0=L0, shape1=shape1, shape2=shape2, coeff=c(g0, g1), RA = RA,
                                             rs = maxS, RQ = RQ, verbose=TRUE, m=m), 1.3491, tolerance=tol)
  expect_equal(racusum_betabinomial_crit_sim(L0=L0, shape1=shape1, shape2=shape2, coeff=c(g0, g1), RA = RA,
                                             rs = maxS, RQ = RQ, verbose=FALSE, m=m), 1.3491, tolerance=tol)
})