pois.ewma.arl: Compute ARLs of Poisson EWMA control charts

In spc: Statistical Process Control -- Calculation of ARL and Other Control Chart Performance Measures

Compute ARLs of Poisson EWMA control charts

Description

Computation of the (zero-state) Average Run Length (ARL) at given mean mu.

Usage

pois.ewma.arl(lambda, AL, AU, mu0, z0, mu, sided="two", rando=FALSE, gL=0, gU=0,
mcdesign="transfer", N=101)

Arguments

lambda	smoothing parameter of the EWMA p control chart.
AL, AU	factors to build the lower and upper control limit, respectively, of the Poisson EWMA control chart.
mu0	in-control mean.
z0	so-called headstart (give fast initial response).
mu	actual mean.
sided	distinguishes between one- and two-sided EWMA control chart by choosing "upper", "lower", and "two", and "zwei", respectively.
rando	Switch between the standard limit treatment, FALSE, and an additional randomisation (to allow ‘perfect’ ARL calibration) by setting TRUE. If randomisation is used, then set the corresponding probailities, gL and gU, appropriately.
gL, gU	If the EWMA statistic is at the limit (approximately), then an alarm is triggered with probability gL and gU for the lower and upper limit, respectively.
mcdesign	choose either "classic" which follows Borror, Champ and Rigdon (1998), or the more sophisticated "transfer" which improves the accuracy heavily.
N	number of states of the approximating Markov chain; is equal to the dimension of the resulting linear equation system.

Details

The monitored data follow a Poisson distribution with mu. The ARL values of the resulting EWMA control chart are determined by Markov chain approximation. We follow the algorithm given in Borror, Champ and Rigdon (1998). However, by setting mcdesign="transfer" (now the default) from Morais and Knoth (2020), the accuracy is considerably improved.

Value

Return single value which resembles the ARL.

Author(s)

Sven Knoth

References

C. M. Borror, C. W. Champ and S. E. Rigdon (1998) Poisson EWMA control charts, Journal of Quality Technonlogy 30(4), 352-361.

M. C. Morais and S. Knoth (2020) Improving the ARL profile and the accuracy of its calculation for Poisson EWMA charts, Quality and Reliability Engineering International 36(3), 876-889.

See Also

later.

Examples

## Borror, Champ and Rigdon (1998), Table 2, PEWMA column
mu0 <- 20
lambda <- 0.27
A <- 3.319
mu1  <- c(2*(3:15), 35)
ARL1 <- rep(NA, length(mu1))
for ( i in 1:length(mu1) )
  ARL1[i] <- pois.ewma.arl(lambda, A, A, mu0, mu0, mu1[i], mcdesign="classic")
print(cbind(mu1, round(ARL1, digits=1)))

## the same numbers with improved accuracy
ARL2 <- rep(NA, length(mu1))
for ( i in 1:length(mu1) )
  ARL2[i] <- pois.ewma.arl(lambda, A, A, mu0, mu0, mu1[i], mcdesign="transfer")
print(cbind(mu1, round(ARL2, digits=1)))

## Morais and Knoth (2020), Table 2, lambda = 0.27 column
lambda <- 0.27
AL <- 3.0870
AU <- 3.4870
gL <- 0.001029
gU <- 0.000765
mu0 <- 20
mu1  <- c(16, 18, 19.99, mu0, 20.01, 22, 24)
ARL3 <- rep(NA, length(mu1))
for ( i in 1:length(mu1) )
  ARL3[i] <- pois.ewma.arl(lambda,AL,AU,mu0,mu0,mu1[i],rando=TRUE,gL=gL,gU=gU, N=101)
print(cbind(mu1, round(ARL3, digits=1)))

spc documentation built on Oct. 24, 2022, 5:07 p.m.