xsewma.arl {spc} | R Documentation |
Compute ARLs of simultaneous EWMA control charts (mean and variance charts)
Description
Computation of the (zero-state) Average Run Length (ARL)
for different types of simultaneous EWMA control charts
(based on the sample mean and the sample variance S^2
)
monitoring normal mean and variance.
Usage
xsewma.arl(lx, cx, ls, csu, df, mu, sigma, hsx=0, Nx=40, csl=0,
hss=1, Ns=40, s2.on=TRUE, sided="upper", qm=30)
Arguments
lx |
smoothing parameter lambda of the two-sided mean EWMA chart. |
cx |
control limit of the two-sided mean EWMA control chart. |
ls |
smoothing parameter lambda of the variance EWMA chart. |
csu |
upper control limit of the variance EWMA control chart. |
df |
actual degrees of freedom, corresponds to subgroup size (for known mean it is equal to the subgroup size, for unknown mean it is equal to subgroup size minus one. |
mu |
true mean. |
sigma |
true standard deviation. |
hsx |
so-called headstart (enables fast initial response) of the mean chart – do not confuse with the true FIR feature considered in xewma.arl; will be updated. |
Nx |
dimension of the approximating matrix of the mean chart. |
csl |
lower control limit of the variance EWMA control chart; default value is 0;
not considered if |
hss |
headstart (enables fast initial response) of the variance chart. |
Ns |
dimension of the approximating matrix of the variance chart. |
s2.on |
distinguishes between |
sided |
distinguishes between one- and two-sided two-sided EWMA- |
qm |
number of quadrature nodes used for the collocation integrals. |
Details
xsewma.arl
determines the Average Run Length (ARL) by
an extension of Gan's (derived from ideas already published by Waldmann)
algorithm. The variance EWMA part is treated
similarly to the ARL calculation method
deployed for the single variance EWMA charts in Knoth (2005), that is, by means of
collocation (Chebyshev polynomials). For more details see Knoth (2007).
Value
Returns a single value which resembles the ARL.
Author(s)
Sven Knoth
References
K. Waldmann (1986), Bounds for the distribution of the run length of geometric moving average charts, J. R. Stat. Soc., Ser. C, Appl. Stat. 35, 151-158.
F. F. Gan (1995), Joint monitoring of process mean and variance using exponentially weighted moving average control charts, Technometrics 37, 446-453.
S. Knoth (2005),
Accurate ARL computation for EWMA-S^2
control charts,
Statistics and Computing 15, 341-352.
S. Knoth (2007), Accurate ARL calculation for EWMA control charts monitoring simultaneously normal mean and variance, Sequential Analysis 26, 251-264.
See Also
xewma.arl
and sewma.arl
for zero-state ARL computation of
single mean and variance EWMA control charts, respectively.
Examples
## Knoth (2007)
## collocation results in Table 1
## Monte Carlo with 10^9 replicates: 252.307 +/- 0.0078
# process parameters
mu <- 0
sigma <- 1
# subgroup size n=5, df=n-1
df <- 4
# lambda of mean chart
lx <- .134
# c_mu^* = .345476571 = cx/sqrt(n) * sqrt(lx/(2-lx)
cx <- .345476571*sqrt(df+1)/sqrt(lx/(2-lx))
# lambda of variance chart
ls <- .1
# c_sigma = .477977
csu <- 1 + .477977
# matrix dimensions for mean and variance part
Nx <- 25
Ns <- 25
# mode of variance chart
SIDED <- "upper"
arl <- xsewma.arl(lx, cx, ls, csu, df, mu, sigma, Nx=Nx, Ns=Ns, sided=SIDED)
arl