tests/testthat/test_attrCUSUM_getContl_zipois.R

In attrCUSUM: Tools for Attribute VSI CUSUM Control Chart

#library(attrCUSUM)
#library(testthat)

context("Tests related to getContl_zipois")

# NULL, NA, TRUE/FALSE, NaN, Inf/-Inf, Int, Real, Complex, Character, List, Expression
test_that("Arguments checking : getContl_zipois", {
  # Check : rho (numeric value in [0,1])
  # NULL, NA, TRUE/FALSE, NaN, Inf/-Inf
  expect_error(getContl_zipois(rho = NULL, lambda = 4, anss.target = 200,
                               refv = 5,
                               c.zero = 0, maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = NA, lambda = 4, anss.target = 200,
                               refv = 5,
                               c.zero = 0, maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = TRUE, lambda = 4, anss.target = 200,
                               refv = 5,
                               c.zero = 0, maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = NaN, lambda = 4, anss.target = 200,
                               refv = 5,
                               c.zero = 0, maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = Inf, lambda = 4, anss.target = 200,
                               refv = 5,
                               c.zero = 0, maxndec = 7L, maxnumsubI = 1000L))
  # Invalid numeric value
  expect_error(getContl_zipois(rho = -0.1, lambda = 4, anss.target = 200,
                               refv = 5,
                               c.zero = 0, maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 1.1, lambda = 4, anss.target = 200,
                               refv = 5,
                               c.zero = 0, maxndec = 7L, maxnumsubI = 1000L))
  # Complex, Character, List, Expression
  expect_error(getContl_zipois(rho = 0.5 + 1i, lambda = 4, anss.target = 200,
                               refv = 5,
                               c.zero = 0, maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = "0.5", lambda = 4, anss.target = 200,
                               refv = 5,
                               c.zero = 0, maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = list(0.5), lambda = 4, anss.target = 200,
                               refv = 5,
                               c.zero = 0, maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = expression(0.5), lambda = 4,
                               anss.target = 200, refv = 5,
                               c.zero = 0, maxndec = 7L, maxnumsubI = 1000L))

  # Check : lambda (positive numeric value)
  # NULL, NA, TRUE/FALSE, NaN, Inf/-Inf
  expect_error(getContl_zipois(rho = 0.5, lambda = NULL, anss.target = 200,
                               refv = 5,
                               c.zero = 0, maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = NA, anss.target = 200,
                               refv = 5,
                               c.zero = 0, maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = TRUE, anss.target = 200,
                               refv = 5,
                               c.zero = 0, maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = NaN, anss.target = 200,
                               refv = 5,
                               c.zero = 0, maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = Inf, anss.target = 200,
                               refv = 5,
                               c.zero = 0, maxndec = 7L, maxnumsubI = 1000L))
  # Invalid numeric value
  expect_error(getContl_zipois(rho = 0.5, lambda = 0, anss.target = 200,
                               refv = 5,
                               c.zero = 0, maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = -1, anss.target = 200,
                               refv = 5,
                               c.zero = 0, maxndec = 7L, maxnumsubI = 1000L))
  # Complex, Character, List, Expression
  expect_error(getContl_zipois(rho = 0.5, lambda = 1 + 1i, anss.target = 200,
                               refv = 5,
                               c.zero = 0, maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = "1", anss.target = 200,
                               refv = 5,
                               c.zero = 0, maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = list(1), anss.target = 200,
                               refv = 5,
                               c.zero = 0, maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = expression(1),
                               anss.target = 200, refv = 5,
                               c.zero = 0, maxndec = 7L, maxnumsubI = 1000L))

  # Check : anss.target (positive numeric value in (1,50000])
  # NULL, NA, TRUE/FALSE, NaN, Inf/-Inf
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = NULL,
                               refv = 5, c.zero = 0,
                               maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = NA,
                               refv = 5, c.zero = 0,
                               maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = TRUE,
                               refv = 5, c.zero = 0,
                               maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = NaN,
                               refv = 5, c.zero = 0,
                               maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = Inf,
                               refv = 5, c.zero = 0,
                               maxndec = 7L, maxnumsubI = 1000L))
  # Invalid numeric value
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 1,
                               refv = 5, c.zero = 0,
                               maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4,
                               anss.target = .Machine$double.xmax,
                               refv = 5, c.zero = 0, maxndec = 7L,
                               maxnumsubI = 1000L))
  # Complex, Character, List, Expression
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200 + 1i,
                               refv = 5, c.zero = 0,
                               maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = "200",
                               refv = 5, c.zero = 0,
                               maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = list(200),
                               refv = 5, c.zero = 0,
                               maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4,
                               anss.target = expression(200), refv = 5, c.zero = 0,
                               maxndec = 7L, maxnumsubI = 1000L))

  # Check : refv (positive numeric value)
  # NULL, NA, TRUE/FALSE, NaN, Inf/-Inf
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = NULL, c.zero = 0,
                               maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = NA, c.zero = 0,
                               maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = TRUE, c.zero = 0,
                               maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = NaN, c.zero = 0,
                               maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = Inf, c.zero = 0,
                               maxndec = 7L, maxnumsubI = 1000L))
  # Invalid numeric value
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 0, c.zero = 0,
                               maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = -1.5, c.zero = 0,
                               maxndec = 7L, maxnumsubI = 1000L))
  # Complex, Character, List, Expression
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5 + 1i, c.zero = 0,
                               maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = "5", c.zero = 0,
                               maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = list(5), c.zero = 0,
                               maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = expression(5), c.zero = 0,
                               maxndec = 7L, maxnumsubI = 1000L))

  # Check : c.zero (positive numeric value in [-'refv', 'contl'))
  # NULL, NA, TRUE/FALSE, NaN, Inf/-Inf
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = NULL,
                               maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = NA,
                               maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = TRUE,
                               maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = NaN,
                               maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = Inf,
                               maxndec = 7L, maxnumsubI = 1000L))
  # Invalid numeric value
  # c.zero < -'refv'
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = -5.1,
                               maxndec = 7L, maxnumsubI = 1000L))
  # Complex, Character, List, Expression
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0 + 1i,
                               maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = "0",
                               maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = list(0),
                               maxndec = 7L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = expression(0),
                               maxndec = 7L, maxnumsubI = 1000L))

  # Check : maxndec (non-negative integer value <= 7L)
  # NULL, NA, TRUE/FALSE, NaN, Inf/-Inf
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0,
                               maxndec = NULL, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0,
                               maxndec = NA, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0,
                               maxndec = TRUE, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0,
                               maxndec = NaN, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0,
                               maxndec = Inf, maxnumsubI = 1000L))

  # Invalid integer value
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0,
                               maxndec = -1L, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0,
                               maxndec = 8L, maxnumsubI = 1000L))
  # Real, Complex, Character, List, Expression
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0,
                               maxndec = 5.5, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0,
                               maxndec = 5L + 1i, maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0,
                               maxndec = "7L", maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0,
                               maxndec = list(7L), maxnumsubI = 1000L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0,
                               maxndec = expression(7L), maxnumsubI = 1000L))

  # Check : maxnumsubI (positive integer value in [100,6000])
  # NULL, NA, TRUE/FALSE, NaN, Inf/-Inf
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0,
                               maxndec = 7L, maxnumsubI = NULL))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0,
                               maxndec = 7L, maxnumsubI = NA))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0,
                               maxndec = 7L, maxnumsubI = TRUE))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0,
                               maxndec = 7L, maxnumsubI = NaN))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0,
                               maxndec = 7L, maxnumsubI = Inf))

  # Invalid integer value
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0,
                               maxndec = 7L, maxnumsubI = 99L))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0,
                               maxndec = 7L, maxnumsubI = 10000L))
  # Real, Complex, Character, List, Expression
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0,
                               maxndec = 7L, maxnumsubI = 5000.1))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0,
                               maxndec = 7L, maxnumsubI = 1000L + i))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0,
                               maxndec = 7L, maxnumsubI = "1000L"))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0,
                               maxndec = 7L, maxnumsubI = list(1000L)))
  expect_error(getContl_zipois(rho = 0.5, lambda = 4, anss.target = 200,
                               refv = 5, c.zero = 0,
                               maxndec = 7L, maxnumsubI = expression(1000L)))
})

test_that("Results : getContl_zipois", {
  # Test for results
  rho <- 0.9
  lambda <- 4
  anss.target <- 200L
  refv <- 5
  c.zero <- 0
  maxndec <- 7L
  maxnumsubI <- 1000L
  res <- getContl_zipois(rho = rho, lambda = lambda, anss.target = anss.target,
                         refv = refv, c.zero = c.zero, maxndec = maxndec,
                         maxnumsubI = maxnumsubI)
  expect_equal(res$refv.act, 5)
  expect_equal(res$c.zero.act, 0)
  expect_equal(unname(round(res$sol1, 2L)), c(3, 189.56))
  expect_equal(unname(round(res$sol2, 2L)), c(4, 438.19))

  rho <- 0.9
  lambda <- 4.1
  anss.target <- 100L
  refv <- 5 + (1 / 3)
  c.zero <- 1 + (1 / 3)
  maxndec <- 1L
  maxnumsubI <- 1000L
  res <- getContl_zipois(rho = 0.5, lambda = lambda, anss.target = anss.target,
                         refv = refv, c.zero = c.zero, maxndec = maxndec,
                         maxnumsubI = maxnumsubI)
  expect_equal(round(res$refv.act, 2L), 5.3)
  expect_equal(round(res$c.zero.act, 2L), 1.3)
  expect_equal(unname(round(res$sol1, 2L)), c(3.70, 69.05))
  expect_equal(unname(round(res$sol2, 2L)), c(3.80, 115.97))
})