inst/tinytest/test_C_pm.R

In pci: A Collection of Process Capability Index Functions

expect_error(C_pm(), 'argument "mu" is missing, with no default')
expect_error(C_pm(2.5), 'argument "sigma" is missing, with no default')
expect_error(C_pm(2.5, 1.), 'argument "target" is missing, with no default')
expect_error(C_pm(2.5, 1., 2.5), 'argument "lsl" is missing, with no default')
expect_error(C_pm(2.5, 1., 2.5, 2.),
             'argument "usl" is missing, with no default')

expect_error(C_pm(2.5, 1., 2.5, 2., 3.),
             'argument "dl" is missing, with no default')

expect_error(C_pm(2.5, 1., 2.5, 2., 3., 6., 7.), 'unused argument \\(7\\)')
expect_error(
  C_pm("2.5", 1., 2.5, 2., 3., 6.),
  '"mu" must be a base-R numeric vector and contain no Inf values'
)

expect_error(
  C_pm(2.5, "1.", 2.5, 2., 3., 6.),
  '"sigma" must be a base-R numeric vector and contain no Inf values'
)

expect_error(
  C_pm(2.5, 1., "2.5", 2., 3., 6.),
  '"target" must be a base-R numeric vector and contain no Inf values'
)

expect_error(
  C_pm(2.5, 1., 2.5, "2.", 3., 6.),
  '"lsl" must be a base-R numeric vector and contain no Inf values'
)

expect_error(
  C_pm(2.5, 1., 2.5, 2., "3.", 6.),
  '"usl" must be a base-R numeric vector and contain no Inf values'
)

expect_error(
  C_pm(2.5, 1., 2.5, 2., 3., "6."),
  '"dl" must be a base-R numeric vector and contain no Inf values'
)

expect_error(C_pm(2.5, 1., 2.5, 2:3, 3:5, 6.),
             'All vector arguments with length >1 must be of equal length')

expect_error(C_pm(2.5, -1., 2.5, 2., 3., 6.), '"sigma >= 0" is false')
expect_error(C_pm(2.5, 1., 2.5, 3., 2., 6.), '"usl > lsl" is false')
expect_error(C_pm(2.5, 1., 2.5, 2., 2., 6.), '"usl > lsl" is false')
expect_error(C_pm(2.5, 1., 2., 2., 3., 6.), '"target > lsl" is false')
expect_error(C_pm(2.5, 1., 3., 2., 3., 6.), '"target < usl" is false')
expect_error(C_pm(2.5, 1., 2.5, 2., 3., 0.), '"dl > 0" is false')

# Test int yields double.
expect_identical(C_pm(2.5, integer(0L), 2.5, 2., 3., 6.), double(0L))

# Test NA over NaN.
expect_identical(C_pm(2.5, 1/6, 2.5, 2., NaN, NA_real_), NA_real_)
expect_identical(C_pm(2.5, 1/6, 2.5, 2., NA_real_, NaN), NA_real_)

expect_identical(C_pm(2.5, 1/6, 2.5, 2., 3., 6.), 1.) # known output
# Test output with NA's.
expect_identical(C_pm(2.5, c(NA, 1/6), 2.5, 2., 3., 6.), c(NA, 1.))
# Test output with NaN's.
expect_identical(C_pm(2.5, c(NaN, 1/6), 2.5, 2., 3., 6.), c(NaN, 1.))

# sigma == 0 & mu == target yields NaN.
expect_identical(C_pm(2.5, 0., 2.5, 2., 3., 6.), NaN)
expect_identical(C_pm(2.5, c(1/6, 0.), 2.5, 2., 3., 6.), c(1., NaN))


# Assert that output won't carry names attribute.
expect_identical(C_pm(c(foo = 2.5), 1/6, 2.5, 2., 3., 6.), 1.)
expect_identical(C_pm(2.5, c(foo = 1/6), 2.5, 2., 3., 6.), 1.)
expect_identical(C_pm(2.5, 1/6, c(foo = 2.5), 2., 3., 6.), 1.)
expect_identical(C_pm(2.5, 1/6, 2.5, c(foo = 2.), 3., 6.), 1.)
expect_identical(C_pm(2.5, 1/6, 2.5, 2., c(foo = 3.), 6.), 1.)
expect_identical(C_pm(2.5, 1/6, 2.5, 2., 3., c(foo = 6.)), 1.)


# Test special cases.
expect_identical(C_pm(2.5, 1/12, 2.5, 2., 3., 6.), C_p(1/12, 2., 3., 6.))


# Test pci_info$expr_r.
e = str2expression(pci_info["C_pm", "expr_r"])
args1 = list(mu = 2.5, sigma = 1/6, target = 2.5, usl = 3.,
             lsl = 2., dl = 6.)

args2 = list(mu = 2.5, sigma = .314, target = 2.5, usl = 3.,
             lsl = 2., dl = 6.)

expect_identical(do.call(C_pm, args1), eval(e, args1))
expect_identical(do.call(C_pm, args2), eval(e, args2))
rm(e, args1, args2)