getAve_binomial {attrCUSUM} | R Documentation |
Computation of ANSSs, ATSs and other related informations on (zero-inflated) binomial CUSUM chart.
getAve_binom(size, prob, refv, contl, c.zero = 0, warnl = 0, ds = 1, dl = NULL, di = NULL, maxndec = 7L, maxnumsubI = 6000L) getAve_zibinom(rho, size, prob, refv, contl, c.zero = 0, warnl = 0, ds = 1, dl = NULL, di = NULL, maxndec = 7L, maxnumsubI = 6000L)
rho |
A length-one numeric vector of zero-inflated parameter in [0,1). |
size |
A length-one non-negative integer-valued vector of number of trials. |
prob |
A length-one numeric vector of probability of success in [0,1]. |
refv |
A reference value of CUSUM statistic. |
contl |
A control limit in CUSUM control scheme. |
c.zero |
An initial state of CUSUM statistic (default is 0). |
warnl |
A warning limit in VSI control scheme (default is 0). |
ds |
The shorter sampling interval in VSI CUSUM control scheme (default is 1). |
dl |
The longer sampling interval in VSI CUSUM control scheme. |
di |
The sampling interval at initial state in VSI CUSUM control scheme. |
maxndec |
The maximum number of decimal places
of |
maxnumsubI |
The maximum number of sub-intervals in [100,6000] to be used in applying the Markov chain approach. |
getAve_binom
and getAve_zibinom
are
wrapper functions of getAve
for (zero-inflated) binomial CUSUM chart.
A list from getAve
.
See getAve
for more information.
Rakitzis et al. (2016). CUSUM Control Charts for the Monitoring of Zero-inflated Binomial Processes, Quality and Reliability Engineering International, 32, 465-483.
# Example 1: Reproduction of results from Rakitzis et al. (2016) rho <- 0.9 size <- 100 prob0 <- 0.01 refv <- 0.26 contl <- 3.86 deltas <- seq(1, 2, 0.1) ANSSs <- numeric(length(deltas)) for(i in seq(deltas)) { prob1 <- deltas[i] * prob0 ANSSs[i] <- getAve_zibinom(rho = rho, size = size, prob = prob1, refv = refv, contl = contl)$ANSS } names(ANSSs) <- deltas ANSSs <- round(ANSSs, 2) ANSSs # Example 2: ANSS profiles (h in seq(10L)) rho <- 0.9 size <- 100 prob0 <- 0.01 refv <- 0.26 h <- seq(10L) ANSSs <- numeric(10) for(i in seq(ANSSs)) { ANSSs[i] <- getAve_zibinom(rho = rho, size = size, prob = prob0, refv = refv, contl = h[i], ds = 1)$ANSS } ANSSs