The one-sided and two-sided Shewhart w-of-w standard and improved runs-rules monitoring schemes to monitor the mean of normally distributed observations from independent and identically distributed (iid) samples are investigated from an overall performance perspective, i.e., the expected weighted run-length (EWRL), for every possible positive integer value of w. The main objective of this work is to use the Markov chain methodology to formulate a theoretical unified approach of designing and evaluating Shewhart w-of-w standard and improved runs-rules for one-sided and two-sided X- schemes in both the zero-state and steady-state modes. Consequently, the main findings of this paper are as follows: (i) the zero-state and steady-state ARL and initial probability vectors of some of the one-sided and two-sided Shewhart w-of-w standard and improved runs-rules schemes are theoretically similar in design; however, their empirical performances are different and (ii) unlike previous studies that use ARL only, we base our recommendations using the zero-state and steady-state EWRL metrics and we observe that the steady-state improved runs-rules schemes tend to yield better performance than the other considered competing schemes, separately, for one-sided and two-sided schemes. Finally, the zero-state and steady-state unified approach run-length equations derived here can easily be used to evaluate other monitoring schemes based on a variety of parametric and nonparametric distributions.
University of South Africa1. Introduction
Balakrishnan and Koutras [1] define a run as an uninterrupted sequence of the same elements bordered at each end by other types of elements. Supplementary runs-rules have been used since the 1950s to improve the performance of the basic Shewhart control charts; see some detailed discussions of some of these earlier works in [2–6]. Some of the commonly cited and recent research works on runs-rules are done in [7–21]. For a literature review on the parametric runs-rules charts that cover articles up to 2006, a reader is referred to Koutras et al. [22], whereas, for the full discussion of nonparametric runs-rules charts until 2017, see the book by Chakraborti and Graham [23]. While runs-rules have been mostly applied in statistical process control and monitoring to improve the detection rate of the basic Shewhart charts, more recently, these have also been used to further increase the detection rate of the exponentially weighted moving average (EWMA) and cumulative sum (CUSUM) schemes; see [24–28].
To differentiate between common and special causes of variation, control charts are the most used tools of statistical process control and monitoring to achieve this goal. That is, when a process has only common causes of variation present, a control chart will indicate that the process is in statistical control, or in short, in-control (IC); however, when a process has special causes of variation present, it is said to be in a state of out-of-control (OOC). Assume that {Yij:i≥1;j=1,2,…,n} is a sequence of samples from iid N(μ0,σ02) distribution. Let Y-i denote the plotting statistic calculated from {Yij} at sampling point i. In a production process, say, samples are usually taken at each sampling point to be inspected and then each of these samples is classified as either conforming or nonconforming depending on where the sample plots on the control charting regions are shown in Figure 1.
Charting zones of the one-sided and two-sided standard runs-rules (SRR) (i.e., left panel) and improved runs-rules (IRR) (i.e., on the right panel) monitoring schemes.
Consider Figure 1, given that μ0 and σ02 are the specified IC mean and variance (process parameters), respectively; let LCL and UCL denote the lower control limit and the upper control limit of some monitoring scheme with limits given by (1)UCLLCL=μY-,0±kσY-,0.In addition to the limits in (1) (however, with different charting constant (i.e., k value)), let LWL and UWL denote the lower and upper warning limits of the Y- monitoring scheme given by(2)UCLLCL=μY-,0±k1σY-,0,UWLLWL=μY-,0±k2σY-,0,where μY-,0 and μY-,0 are the specified IC mean and variance of the plotting statistic Y-, respectively. Note that the UCL and LCL in (1) and in (2) are not equal because the control limits constants have the following relation: k<k1. This is so that the resulting control limits yield the constraint that the actual average run-length (ARL) must equal the nominal IC ARL (denoted by ARL0); otherwise, if k>k1, the additional warning limits in (2) will lower the ARL0.
The charting regions in the left panel of Figure 1 correspond to the one-sided and two-sided standard runs-rules (i.e., standalone w out of the last w consecutive plotting statistics rule; denoted by SRR); see, for instance, [7, 8]. A one-sided upper (lower) SRR scheme issues an OOC signal when there are w consecutive plotting statistics that fall in Zone A (Zone C), respectively. A non-side-sensitive (denoted by NSS) two-sided SRR scheme issues an OOC signal when there are w consecutive plotting statistics that fall in Zone A or Zone C. A side-sensitive (denoted by SS) two-sided SRR scheme issues an OOC signal when there are w consecutive plotting statistics that fall in Zone A (Zone C), respectively.
However, the charting regions in the right panel of Figure 1 correspond to the one-sided and two-sided improved runs-rules (i.e., combination of the 1-of-1 and SRR; denoted by IRR); see, for instance, Khoo and Ariffin [10]. A one-sided upper (lower) IRR scheme issues an OOC signal when a single sampling point plots in Zone 1 (Zone 5) or w out of the last w consecutive plotting statistics fall in Zone 2 (Zone 4), respectively. A NSS two-sided SRR scheme issues an OOC signal when a single sampling point plots in Zone 1 or Zone 5 or w out of the last w consecutive plotting statistics fall in Zone 2 or Zone 4. A SS two-sided IRR scheme issues an OOC signal when a single sampling point plots in Zone 1 (Zone 5) or w out of the last w consecutive plotting statistics fall in Zone 2 (Zone 4), respectively.
The zero-state and steady-state mode of analysis are used to characterize the short-term and long-term run-length properties of a monitoring scheme. It should be noted that Champ [6] discussed the steady-state one-sided IRR schemes performance and derived some of its ARL expressions. Next, Balakrishnan and Koutras [1] showed that the zero-state NSS SRR scheme’s run-length distribution is the same as the geometric distribution of order w. Shmueli and Cohen [29] derived some closed-form ARL expressions of the SS two-sided SRR scheme. Acosta-Mejia [30] conducted an empirical zero-state ARL performance of the two-sided SS SRR and IRR schemes. More recently, Lim and Cho [31] conducted an extensive investigation into the empirical performance and derived the steady-state closed-form ARL expressions for the SS two-sided IRR scheme.
The main objective of this paper is to unify these publications (i.e., Champ [6], Balakrishnan and Koutras [1], Shmueli and Cohen [29], Khoo and Ariffin [10], Acosta-Mejia [30], and Lim and Cho [31]) and formulate a unified approach to evaluate one-sided and two-sided w-of-wX- SRR and IRR schemes for any possible integer value of w, separately, for the zero-state and the steady-state contexts. More specifically, in this paper, we show the following:
There is an ideal manner to define the w-of-w scheme’s transition probability matrices (TPM) so that it can easily be formulated for any possible integer value of w for both one- and two-sided schemes.
The design structure of the TPM and other run-length distribution properties of the upper/lower one-sided w-of-w SRR and IRR schemes are actually similar to those of the two-sided NSS w-of-w SRR and IRR schemes with different probability elements.
Derive initial probabilities and ARL vectors, so that we formulate the zero-state and steady-state closed-form ARL expressions for the one-sided and two-sided w-of-w SRR and IRR schemes for any possible integer value of w.
In the papers by [6, 10, 30, 31] that go into detail about w-of-wX- runs-rules monitoring schemes, it is not easy to figure out how one should select a specific best value of w to use as the ARL metric is based on a specific size shift which must be determined in advance. To bypass this problem, in this paper, we propose the use of overall performance measures to examine the performance of one- and two-sided SRR and IRR schemes for a range of small, medium, and large shift sizes. The overall performance measures are better measures than ARL when the quality practitioner does not know beforehand the magnitude of the target shift size, that is, when the shift size is random.
The rest of the paper is structured as follows: In Section 2, we illustrate the difference between the design structure of the one-sided and the two-sided (NSS and SS) SRR and IRR schemes’ TPMs. In Section 3, we describe some run-length properties as well as the overall performance metrics. In Sections 4 and 5, an empirical discussion of the one-sided and two-sided runs-rules schemes is done, respectively. An example is shown in Section 6 illustrating how the monitoring schemes discussed here are implemented in real life. In Section 7, some concluding remarks are given. Finally, in the Appendix, we derive the closed-form expressions of the expected run-length distribution for the one-sided and two-sided SRR and IRR schemes in a different approach from those that exist currently in the literature as separately documented in [1, 6, 29, 31]. Moreover, expressions of the false alarm rate (FAR) are derived in the Appendix for the one- and two-sided SRR and IRR monitoring schemes discussed here.
2. The Design of the w-of-w SRR and IRR Control Chart
Given the charting zones in Figure 1, consider Zone A. The probability of a charting statistics falling in Zone A may be calculated as follows: pA=PY-∈ZoneA=PY-≥UCL=1-Φ(k-δn), where Φ(∙) is the cumulative distribution function (cdf) of the standard normal distribution. Similarly, for the other charting zones, pΛ=PY-∈ZoneΛ, where Λ={A,B,C,D,E,O,1,2,…,9}, with Zone O = Zone A ∪ Zone C, Zone 8 = Zone 1 ∪ Zone 5, and Zone 9 = Zone 2 ∪ Zone 4.
The main requirement of the Markov chain procedure is the TPM of a w-of-w scheme of interest. To construct the TPM, we need to discretize the charting regions of each SRR and IRR monitoring scheme as done in Figure 1. The charting regions corresponding to each SRR scheme are as follows:
One-sided: Upper {A,D} and Lower {C,E}.
Two-sided: NSS {B, O, with O≡A∪C} and SS {A,B,C}.
However, the charting regions corresponding to each IRR scheme are as follows:
One-sided: Upper {1,2,6} and Lower {4,5,7}.
Two-sided: NSS {3,8,9} and SS {1,2,3,4,5}, with 8 ≡ 1∪5 and 9 ≡ 2∪4.
In Tables 1 and 2, we illustrate how the TPM is constructed for each SRR and IRR monitoring scheme when w = 3. That is, in Table 1, we give all the compound patterns, denoted by “OOC”, which depicts consecutive elements plotting on distinct zones in Figure 1 that result in OOC signaling events. The steps involved in constructing the TPMs of the SRR and IRR schemes are as follows for any w:
Step (i): Outline the absorbing states that lead to an OOC signal and denote these as OOC.
Step (ii): Define the conforming zone that represent the IC state, denoted by ϕ, where(3)ϕ=η1,for one-sided and two-sided NSS schemesηw,for two-sided SS schemes.
Step (iii): Decompose the absorbing states in Step (i) into their corresponding transient states and denote these as ηi.
Step (iv): Define the state space, denoted by Ω, which is an amalgamation of Steps (i) to (iii).
Decomposition of the state space of the 3-of-3 standard and improved runs-rules (SRR and IRR) schemes.
Type
OOC
ϕ
η
Ω
SRR
Upper
{AAA}
η1={D}
η2={A}, η3={AA}
{ϕ;η2,η3;OOC}
Lower
{CCC}
η1={E}
η2={C}, η3={CC}
{ϕ;η2,η3;OOC}
NSS
{OOO}
η1={B}
η2={O}, η3={OO}
{ϕ;η2,η3;OOC}
SS
{AAA},{CCC}
η3={B}
η1={AA}, η2={A}, η4={C}, η5={CC}
{η1,η2;ϕ;η4,η5;OOC}
IRR
Upper
{1},{222}
η1={6}
η2={2}, η3={22}
{ϕ;η2,η3;OOC}
Lower
{5},{444}
η1={7}
η2={4}, η3={44}
{ϕ;η2,η3;OOC}
NSS
{8},{999}
η1={3}
η2={9}, η3={99}
{ϕ;η2,η3;OOC}
SS
{1},{5},{222},{444}
η3={3}
η1={22}, η2={2}, η4={4}, η5={44}
{η1,η2;ϕ;η4,η5;OOC}
The transition probability matrices of the 3-of-3 standard and improved runs-rules (SRR and IRR) charts.
Therefore, following the latter description, the state spaces for each of the schemes are shown in Table 1 and these are used to construct each of the TPMs in Table 2 for the one-sided (upper and lower) and two-sided SRR and IRR monitoring schemes when w = 3.
We see from Table 2 that the TPM consists of absorbing and transient states defined within Ω, and its structure is such that, for any positive integer τ, it is given by an (τ+1)×(τ+1) matrix, P:(4)P=Q∣r---0′∣1where the τ×1 vector r satisfies r=1-Q1 with 1=(11…1)′, 0=(00…0)′, Q is the τ×τ essential TPM consisting of transient states, where, for both the w-of-w SRR and IRR monitoring schemes, we have (5)τ=w,for one-sided and two-sided NSS schemes2w-1,for two-sided SS schemes.The construction and properties of TPMs of the one-sided and two-sided SRR and IRR schemes for any possible integer value of w are thoroughly discussed in the Appendix.
3. Run-Length Characteristics of the w-of-w SRR and IRR Control Chart3.1. Some Run-Length Characteristics
Let N denote the run-length of some w-of-w control chart. Then N is the number of sample points plotted on the control chart until it gives an OOC signal for the first time. In this paper, we compute the expected run-length of the chart using the Markov chain technique best explained in Fu and Lou [32]; this is further discussed in the Appendix. The most used quantity to measure the performance of a monitoring scheme is the E(N), and we denote this here as ARL given by (6)ARL=ξ′∙Rwhere ξ is the τ×1 initial probability vector (see Section 3.2) that depends on whether a zero-state or a steady-state analysis is of interest and,(7)R=I-Q-1∙1,where I is the τ×τ identity matrix. The R closed-form expressions are formulated in the Appendix for each of the considered schemes.
3.2. Initial Probabilities Vectors
ξ=q is the vector of initial probabilities associated with the zero-state mode and it has a one in the component associated with the state in which the chart begins (i.e., state ϕ) and each of the other components of the vector is equal to zero; this is further shown in the Appendix.
ξ=s is the vector of initial probabilities associated with the steady-state mode and its elements are nonzero. There are a number of methods used to compute s, and in this paper, we focus on three of these steady-state probability vector (SSPV) methods which are denoted here by SSPV1, SSPV2, and SSPV3 (each of these is computed while the process is IC; i.e., δ = 0).
(i) SSPV1 Method. The SSPV1 method (by Crosier [33]) entails computing P∗, by altering P in (4) so that the control statistic is reset to the “initial state” whenever it goes into an “OOC state”. That is, the last row of the TPM is changed such that the value of one is moved to the respective initial state (i.e., state ϕ) instead of the OOC state. That is, (4) becomes P∗=Qrej′0, where ej is the jth unit vector. Note that ej corresponds to e1 for the one-sided (upper or lower) and the two-sided NSS SRR and IRR schemes. However, ej corresponds to ew for the two-sided SS SRR and IRR schemes. Consequently, we then use P∗ to find the (τ+1)×1 probability vector π such that the following equation is satisfied: π=P∗′π subject to 1′π=1. Finally, the SSPV1 method yields (8)s=1′z-1∙zwhere z is the τ×1 vector obtained from π by deleting the (τ+1)th component associated with the absorbing state.
(ii) SSPV2 Method. The SSPV2 method (by Champ [6]) is given by (8), however, with zτ×1=G-Q′-1ej and the matrix G is given by G=ej∙1′+I.
(iii) SSPV3 Method. The SSPV3 method (used by [31, 34–36], etc.) is obtained by dividing each element of Q by its corresponding row sum, so that we may have an ergodic altered version of the essential TPM called the conditional essential TPM, which is denoted by QC. Consequently, the SSPV3 method is a vector such that s∙QC=s subject to ∑i=1τsi=1.
The SSPV1, SSPV2, and SSPV3 are each formulated in the Appendix for each of the runs-rules schemes discussed in this paper. Calculations in this paper were done using SSPV2 method.
3.3. Overall Performance Measures
A number of authors have argued that if a control chart is designed based on one specific size of a shift, it would perform poorly when the actual size of a mean shift is significantly different from the assumed size; see [36–41]. Hence, they recommend that control charts should be designed in terms of the overall performance rather than a specific shift size performance. The expected weighted run-length (EWRL) is a quality loss function that describes the relationship between the shift size and the quality impact of a control chart, overall; and this is given by (9)EWRL=Ewδ×ARLδ=∫δminδmaxwδ×ARLδ×fδdδ,where δ follows some probability distribution function with a density function fδ and a range [δmin, δmax], where δmin and δmax are the lower and upper bound of the range of δ, and wδ is a weight function associated with δ. Note that the EWRL is a generalized quality loss function and by assigning different weight functions, it yields the following different common quality loss function metrics: (10)(i)Extra quadratic loss EQL if wδ=δ2;(ii)Expected ARLEARL if wδ=1.
Note that the logic behind the EQL weight function is that the larger the shift size, the greater the quality loss, whereas the EARL assigns the same weight on each ARL value, irrespective of the shift size.
Here we compute both the zero-state and steady-state EQL and EARL to investigate whether different EWRL functions have a similar or different effect on the choice of the optimal value of w for each of the w-of-w SRR and IRR schemes. Moreover, we consider only the case where fδ follows a Uniform (0, 1) distribution, which in a way implies that the objective function (i.e., (9)) that needs to be minimized can equivalently be written as(11)EWRL=Ewδ×ARLδ=1δmax-δmin∑δ=δminδmaxwδ×ARLδΔδ.
Throughout this paper, we use the increment in the shift, i.e., Δδ, equal to 0.1. Finally, for any competing schemes, the best scheme will be the one that yields the smallest EWRL value.
4. Performance of the One-Sided w-of-w SRR and IRR Monitoring Schemes
A monitoring scheme is designed such that when the process is IC, the ARL0 is set at some desirable level (or equivalently, the significance level is set at some standard value). For instance, a significance level of sizes 0.005, 0.0027, 0.0020, and 0.0010 implies that the ARL0 = 200, 370.4, 500, and 1000, respectively. Due to writing space constraint, only the performance relating to ARL0 = 370.4 is illustrated and for the other ARL0 values, a similar conclusion follows.
We conduct the analysis of the OOC performance by separately looking at two run-length characteristics, i.e., the ARL and EWRL.
Based on the ARL: note that, for w > 7, there is no k > 0 such that the actual IC ARL is equal to 370.4 in Table 3. Next, we use (A.4) and (A.5) to compute the zero-state and steady-state ARLs which are shown in Table 3. In zero-state, each one-sided w-of-w scheme converges to a lower bound ARL value equal to w for any large shift value; that is, a one-sided w-of-w scheme can only signal after exactly w sampling points. Note though, in steady-state, the lower bound is slightly less than the value of w. For small shifts, i.e., δ < 1, the higher the value of w, the better, as this yields smaller OOC ZSARL and SSARL values. However, for large shifts, increasing w is not advisable due to the lower bound just explained. For δ > 1.5, the basic one-sided X- chart tends to be more competitive, as it outperforms the one-sided w-of-w schemes with higher values of w. Due to a lack of a single monitoring scheme outperforming the rest, for all shift values, separately in zero-state and steady-state modes, it is not easy to choose the optimal value of w.
Based on the EWRL: using (A.12) and (A.13) in the Appendix, we calculate the zero-state and steady-state EARL and EQL given in Table 4. In both states, as δmax decreases, the optimal w increases; see the boldfaced values that yield a minimum EWRL for a given range of w values. We see that, in each state, the EARL of the w-of-w schemes is better than that of the one-sided X- chart. On the contrary, only w = 2, 3, 4 in both states yield EQLs less than that of the X- chart when δmax = 3.
The zero- and steady-state ARL for the one-sided w-of-w SRR X- charts when ARL0 = 370.4.
Zero-state
Steady-state
w
1
2
3
4
5
6
7
2
3
4
5
6
7
k
2.7822
1.6134
1.0509
0.6927
0.4348
0.2358
0.0752
1.6141
1.0520
0.6942
0.4367
0.2381
0.0778
δ
0
370.4
370.4
370.4
370.4
370.4
370.4
370.4
370.4
370.4
370.4
370.4
370.4
370.4
0.1
273.4
251.4
240.7
234.5
230.4
227.6
225.7
251.3
240.4
233.9
229.6
226.6
224.4
0.2
203.7
173.9
160.7
153.5
149.0
146.1
144.3
173.7
160.2
152.7
148.0
144.8
142.7
0.3
153.2
122.5
110.2
103.8
100.1
97.9
96.6
122.3
109.6
103.0
99.0
96.5
94.9
0.4
116.2
87.9
77.5
72.5
69.7
68.2
67.5
87.7
76.9
71.6
68.6
66.9
65.8
0.5
89.0
64.3
55.9
52.1
50.3
49.4
49.1
64.0
55.4
51.4
49.2
48.1
47.5
0.6
68.7
47.8
41.3
38.6
37.5
37.1
37.1
47.6
40.8
37.9
36.5
35.8
35.6
0.7
53.6
36.2
31.3
29.4
28.8
28.7
29.0
36.0
30.8
28.8
27.9
27.6
27.6
0.8
42.1
27.9
24.2
23.0
22.7
22.9
23.4
27.7
23.8
22.4
21.9
21.9
22.1
0.9
33.4
21.9
19.2
18.4
18.5
18.8
19.4
21.7
18.8
17.9
17.7
17.8
18.2
1
26.8
17.4
15.5
15.1
15.3
15.8
16.5
17.3
15.1
14.6
14.6
14.9
15.3
1.1
21.6
14.1
12.8
12.7
13.0
13.6
14.3
14.0
12.4
12.2
12.3
12.7
13.2
1.2
17.6
11.6
10.7
10.8
11.3
12.0
12.7
11.5
10.4
10.3
10.6
11.1
11.7
1.3
14.5
9.7
9.1
9.4
10.0
10.7
11.5
9.6
8.9
9.0
9.3
9.9
10.5
1.4
12.0
8.2
7.9
8.3
9.0
9.7
10.5
8.1
7.7
7.9
8.3
8.9
9.5
1.5
10.0
7.0
7.0
7.5
8.1
8.9
9.8
6.9
6.7
7.1
7.6
8.2
8.8
1.6
8.4
6.1
6.2
6.8
7.5
8.3
9.2
6.0
6.0
6.4
6.9
7.6
8.2
1.7
7.2
5.4
5.6
6.2
7.0
7.8
8.7
5.3
5.4
5.9
6.5
7.1
7.8
1.8
6.1
4.8
5.1
5.8
6.6
7.5
8.4
4.7
4.9
5.4
6.1
6.7
7.4
1.9
5.3
4.3
4.7
5.5
6.3
7.2
8.1
4.2
4.5
5.1
5.8
6.4
7.2
2
4.6
3.9
4.4
5.2
6.0
6.9
7.8
3.8
4.2
4.8
5.5
6.2
6.9
2.1
4.0
3.6
4.2
4.9
5.8
6.7
7.6
3.5
4.0
4.6
5.3
6.0
6.8
2.2
3.6
3.3
3.9
4.8
5.6
6.6
7.5
3.2
3.8
4.4
5.1
5.9
6.6
2.3
3.2
3.1
3.8
4.6
5.5
6.4
7.4
3.0
3.6
4.3
5.0
5.8
6.5
2.4
2.8
2.9
3.6
4.5
5.4
6.3
7.3
2.8
3.4
4.2
4.9
5.7
6.4
2.5
2.6
2.7
3.5
4.4
5.3
6.3
7.2
2.7
3.3
4.1
4.8
5.6
6.4
2.6
2.3
2.6
3.4
4.3
5.2
6.2
7.2
2.6
3.2
4.0
4.8
5.5
6.3
2.7
2.1
2.5
3.3
4.2
5.2
6.1
7.1
2.5
3.2
3.9
4.7
5.5
6.3
2.8
2.0
2.4
3.3
4.2
5.1
6.1
7.1
2.4
3.1
3.9
4.7
5.4
6.2
2.9
1.8
2.3
3.2
4.1
5.1
6.1
7.1
2.3
3.0
3.8
4.6
5.4
6.2
3
1.7
2.3
3.2
4.1
5.1
6.1
7.0
2.2
3.0
3.8
4.6
5.4
6.2
The zero-state and steady-state EARL and EQL (with δmin=0 and δmax=3, 2, 1) for the one-sided w-of-w, w∈{2,3,…,7}, SRR X- charts when ARL0 = 370.4.
Zero-state
Steady-state
δmax
X-
2
3
4
5
6
7
2
3
4
5
6
7
EARL
3
521.3
441.5
418.6
411.2
410.3
412.8
417.5
440.3
415.6
406.6
403.6
404.1
406.5
2
768.9
648.4
610.2
594.8
588.8
587.8
590.0
646.9
606.7
589.4
581.2
578.1
577.9
1
1430.5
1221.6
1146.9
1111.3
1092.7
1082.9
1079.0
1219.7
1142.2
1104.6
1083.4
1071.3
1064.5
EQL
3
193.4
158.3
170.4
192.3
217.9
245.8
274.5
156.3
164.0
181.5
201.6
223.3
245.9
2
209.5
149.2
141.3
145.1
152.9
162.9
174.3
147.6
137.6
139.0
144.3
151.4
159.8
1
197.3
138.6
121.9
115.5
113.4
113.3
114.6
137.9
120.2
113.2
110.2
109.2
109.4
Based on this example, it is apparent that the different EWRL functions do lead to different recommended values of w. Hence, the choice between any EWRL function needs to depend on each user as per weight function structure preference in (10) and the magnitude of shifts of interest. That is, we recommend w found using EARL approach when all the shifts are equally important (i.e., the quality practitioner is interested in all magnitudes of shifts) and recommend w found using EQL approach when the magnitude of the shift is more important (i.e., the quality practitioner is interested in shifts according to their magnitude). Thus, moving forward, we separately present both the results of the EARL and EQL so that we may see the resulting optimal values in each case.
Next, the w-of-w IRR schemes have two design parameters for each w (i.e., k1 and k2). We proceed as similarly done in Tables 3 and 4 (i.e., use (A.12) and (A.13) and then compute the EARL (see Table 5) and EQL (see Table 6)) for w=2,3,…, 11, k1=2.8,2.9,…,4.0 when δmin = 0 and δmax = 3.
The zero- and steady-state EARL (with δmin=0 and δmax=3) for the one-sided w-of-w IRR X- charts when ARL0 = 370.4.
w
Zero-state
Steady-state
2
3
4
5
6
7
8
9
10
11
2
3
4
5
6
7
8
9
10
11
k1
2.8
500.8
487.6
479.9
475.0
471.8
469.6
468.0
467.0
466.3
466.0
500.7
487.3
479.4
474.3
470.8
468.4
466.7
465.4
464.6
464.0
2.9
463.5
442.3
432.1
426.4
423.3
421.8
421.1
421.0
463.1
441.3
430.4
424.2
420.6
418.5
417.3
416.8
3.0
449.8
427.1
416.8
411.8
409.3
408.4
408.6
449.1
425.7
414.5
408.6
405.5
403.8
403.3
3.1
443.1
419.9
410.0
405.4
403.5
403.2
404.0
442.4
418.2
407.2
401.7
398.9
397.7
397.7
3.2
439.9
416.4
406.7
402.6
401.3
401.5
439.0
414.5
403.7
398.3
396.0
395.2
3.3
438.5
414.8
405.4
401.6
400.8
401.5
437.4
412.7
402.0
397.0
394.9
394.5
3.4
437.9
414.3
405.1
401.6
401.1
402.4
436.8
412.0
401.4
396.7
394.9
394.9
3.5
438.0
414.3
405.3
402.1
402.1
403.7
436.8
411.8
401.5
396.9
395.4
395.7
3.6
438.3
414.6
405.8
403.0
403.2
405.1
437.1
412.1
401.8
397.5
396.3
396.7
3.7
438.7
415.2
406.5
403.9
404.3
406.6
437.5
412.5
402.3
398.2
397.2
398.0
3.8
439.2
415.6
407.1
404.8
405.6
408.1
437.9
413.0
403.0
399.0
398.2
399.1
3.9
439.6
416.1
407.8
405.6
406.6
409.5
438.3
413.5
403.4
399.6
399.0
400.2
4.0
440.0
416.5
408.5
406.4
407.7
410.8
438.6
413.8
404.0
400.3
399.9
401.1
The zero- and steady-state EQL (with δmin=0 and δmax=3) for the one-sided w-of-w IRR X- charts when ARL0 = 370.4.
w
Zero-state
Steady-state
2
3
4
5
6
7
8
9
10
11
2
3
4
5
6
7
8
9
10
11
k1
2.8
173.3
163.2
160.1
159.6
160.0
160.9
162.0
163.2
164.4
165.5
173.1
162.5
159.0
158.1
158.3
158.9
159.8
160.8
161.9
162.9
2.9
151.0
143.9
144.5
147.1
150.2
153.4
156.5
159.4
150.3
142.2
141.8
143.6
146.1
148.8
151.4
154.0
3.0
145.7
141.2
144.0
148.5
153.2
157.8
162.2
144.6
138.7
140.2
143.6
147.5
151.4
155.1
3.1
144.1
141.7
146.3
152.4
158.6
164.5
170.1
142.9
138.7
141.7
146.4
151.5
156.5
161.3
3.2
144.2
143.4
149.6
157.3
164.9
172.2
142.7
140.0
144.3
150.3
156.6
162.7
3.3
145.1
145.7
153.4
162.6
171.7
180.4
143.4
141.9
147.5
154.7
162.1
169.4
3.4
146.3
148.2
157.4
168.0
178.6
188.7
144.5
144.1
150.8
159.2
167.8
176.2
3.5
147.6
150.8
161.3
173.3
185.4
197.0
145.7
146.3
154.2
163.7
173.5
183.1
3.6
149.0
153.2
165.0
178.5
192.1
205.2
147.0
148.5
157.4
168.1
179.1
189.8
3.7
150.3
155.6
168.6
183.4
198.5
213.1
148.3
150.6
160.5
172.3
184.4
196.3
3.8
151.6
157.7
171.9
188.0
204.5
220.6
149.5
152.6
163.4
176.2
189.5
202.5
3.9
152.7
159.7
174.9
192.3
210.1
227.6
150.5
154.3
166.0
179.8
194.1
208.3
4.0
153.7
161.4
177.6
196.1
215.2
234.2
151.5
155.9
168.4
183.1
198.4
213.6
For instance, when w = 3 and given that we use k1 = 2.8,2.9,…,5.0, the couple (k1, k2) = (3.5, 1.0711) yields the lowest steady-state EARL equal to 411.8 over the given range of k1 values when δmin = 0 and δmax = 3. Continuing in a similar manner as done in Tables 5 and 6, for w = 2,3,…,11, k1 = 2.8, 2.9, …, 4.0, and δmax = 3, 2, 1, then at each possible value of w, the minimum EWRL values (i.e., boldfaced values in Tables 5 and 6) are plotted in Figure 2. In each graph in Figure 2, the EWRL value of the one-sided X- chart (when δmin = 0 and δmax = 3, 2, 1) is overlaid as a reference line, respectively.
The zero-state and steady-state EARL and EQL values (with δmin=0) of the one-sided SRR and IRR schemes when ARL0 = 370.4.
δmax=3
δmax=2
δmax=1
Firstly, as expected, in each state and for each w (wherever both the one-sided w-of-w SRR and IRR exist), the IRR scheme has a better overall performance. Secondly, using either the EARL or EQL, we see that, for each w, the zero-state one-sided SRR scheme has the worst performance (or least improvement from the basic one-sided X- chart) whereas the steady-state one-sided IRR scheme has the best performance. Thirdly, in terms of EARL, the one-sided w-of-w SRR and IRR schemes always outperform the X- chart; however, in terms of the EQL, the zero-state and steady-state one-sided SRR schemes are outperformed by the X- chart once w≥ 5 when δmax = 3. Fourthly, in general, for small shifts (i.e., δmax = 1), there seem to be a small difference in the performance of the different runs-rules schemes using either EARL or EQL; however, for large shifts (i.e., δmax = 3), there seem to be a noticeable significant difference among the one-sided runs-rules schemes, especially when using the EQL metric due to the weight structure in (10). Finally, we observe that each of these EWRL functions (i.e., (A.12)–(A.13)) tends to decrease and then at some point, the curve increases (i.e., concave up function); hence these minimum turning points represent the value of w that yields the lowest EWRL for that particular SRR or IRR scheme.
Thus, based on the EARL, we recommend the use of steady-state mode one-sided IRR scheme with w = 7 for all shift types; however, based on the EQL, we recommend the steady-state one-sided IRR scheme with w = 7, 4, and 3, for small, moderate, and large shifts, respectively. The steady-state mode performance is slightly better than the corresponding zero-state mode; hence, we recommend the steady-state mode to evaluate the performance of the w-of-w monitoring schemes.
5. Performance of the Two-Sided w-of-w SRR and IRR Monitoring Schemes
Similar to the calculations done in Tables 3, 4, 5, and 6 and Figure 2 in the one-sided case, in Figures 3 and 4, we show the zero-state and steady-state EWRL for the two-sided SS case using (A.14) and (A.15). For the two-sided NSS and SS IRR schemes, at each w value, using the corresponding optimal design parameters (k1, k2) and (A.12) to (A.15), the zero- and steady-state EARL and EQL that satisfy mink1∈3.1,3.2,…5.0EARL and mink1∈3.1,3.2,…5.0EQL, for a given δmax = 1, 2, 3 (with δmin = 0) are the ones that are plotted in Figures 3 and 4, respectively.
The zero-state and steady-state EARL and EQL values (with δmin=0) of the two-sided non-side-sensitive SRR and IRR schemes when ARL0 = 370.4.
δmax=3
δmax=2
δmax=1
The zero-state and steady-state EARL and EQL values (with δmin=0) of the two-sided side-sensitive SRR and IRR schemes when ARL0 = 370.4.
δmax=3
δmax=2
δmax=1
The two-sided NSS SRR and IRR schemes in Figure 3 are not recommended, at all, because
increasing w leads to deteriorating overall performance because the EARL and EQL increase as w increases;
for values of w > 3, the two-sided NSS SRR schemes are outperformed by the basic X- chart.
Unlike the two-sided NSS schemes, the two-sided SS SRR and IRR schemes given in Figure 4 have a similar general behavior as those discussed in Figure 2. Thus, following a similar argument as in Figure 2, based on the EARL, we recommend the use of the two-sided SS IRR scheme with w = 8 for all shift types; however, based on the EQL, we recommend the use of the two-sided SS IRR scheme with w = 8, 4, and 3 for small, moderate, and large shift sizes, respectively.
6. Application Example
To illustrate the use and the application of the one-sided and two-sided w-of-w SRR and IRR schemes, we consider a well-known dataset from Montgomery [42] on the inside diameters of piston rings manufactured by a forging process. This data set contains 25 retrospective or Phase I samples, each of size 5, that were collected when the process was thought to be IC. These data are considered to be the Phase I reference data for which a goodness of fit test for normality is not rejected. This data set also contains 15 prospective (Phase II) samples each of 5 observations (i.e., n = 5). Note that when the distribution parameters of a particular process are unknown, it is generally accepted that there are two phases of application for a monitoring scheme, namely, Phase I (for estimation of distribution parameters) and Phase II (continuous monitoring using the parameters estimated in Phase I); see the book by Chakraborti and Graham [23] for further discussion on these phases of application. Using Phase I techniques, with an IC data, we estimate that the mean and standard deviation of the piston rings data are equal to 74.0011 and 0.0048, respectively.
Consequently, the limits in (1) and (2) are given in Table 7 for the upper one-sided SRR and IRR schemes and the side-sensitive two-sided SRR and IRR schemes with w = 4 that yield an IC ARL equal to 370.4.
The design parameters and limits of the piston ring data.
Using the limits in Table 7, we construct the corresponding monitoring schemes in Figure 5 for the upper one-sided and side-sensitive two-sided schemes. In Phase I, all the monitoring schemes depict processes that have some suspect samples but none are OOC according to SRR and IRR guidelines. We observe that, in Phase II, for this specific dataset, the upper one-sided and two-sided schemes issue an OOC signal for the first time at time points 40 and 38 (or, on time points 15 and 13 on Phase II) for the SRR and IRR, respectively, showing the improvement that is brought by the IRR design over the SRR design.
Upper one-sided and two-sided SRR and IRR schemes to monitor the mean of piston ring size.
7. Concluding Remarks
In this paper, we revisited the design of the w-of-w standard and improved runs-rules schemes for one-sided and two-sided charts based on the mean of the normal distribution from iid samples. Then, we implemented a unified approach in designing these schemes and unlike the existing studies which are based on the ARL only (see [6, 10, 30, 31]), we base our recommendations on the overall performance, using specifically, the extra quadratic loss and the expected average run-length. Using these overall performance measures, we show that the one-sided and the two-sided side-sensitive steady-state improved runs-rules schemes have a much better performance than the other competing one-sided and two-sided schemes considered here, respectively. Moreover, we showed that the two-sided non-sided-sensitive standard and improved runs-rules schemes should never be used as they yield a uniformly deteriorating overall performance as w increases.
Furthermore, for ease of calculating expected run-length characteristics, in the Appendix, we derived some closed-form expressions (in a slightly different manner as currently available in the literature) that can easily be used to obtain the zero-state and steady-state average run-length values of the one-sided and two-sided standard and improved runs-rules schemes. These closed-form expressions are valuable because any user with or without prior knowledge of Markov chain or simulation or possessing any advanced statistical software can easily use a pocket calculator to compute the performance measurements of the schemes considered here.
The empirical results in this article are based on the assumption of iid normally distributed samples; hence the study is limited to these assumptions. If the assumptions are violated, the above results may have to be reexamined. The empirical run-length properties of the proposed schemes under the violation of these assumptions are under investigation and will be reported in a separate article. Note though the theoretical expressions derived here may easily be extended for other one-sided or two-sided schemes with symmetric control limits like (i) nonparametric schemes (sign, signed-rank, precedence, minimum, median, etc.), (ii) attributes schemes (number nonconforming, nonconformities per unit item, etc.), and (iii) parametric schemes (t-distribution chart, variable sample size, and interval runs-rules chart, etc.).
Appendix
In this Appendix, we follow a similar line of argument in deriving general TPMs and closed-form expressions as that done for the synthetic charts as well as 2-of-(H+1) runs-rules (with H a positive integer greater than 0) in Shongwe and Graham [43–45] using Markov chain design approach. In Section A, we discuss the theoretical run-length properties of the w-of-w upper or lower one-sided as well as the two-sided NSS monitoring schemes based on the SRR and IRR designs. The theoretical run-length properties of the two-sided SS monitoring schemes based on the SRR and IRR designs are discussed in Section A. Finally, in Section A, the expressions used to calculate the empirical values of the EWRL in Section 4 are shown.
A. A.1. Run-Length Properties of the One-Sided and Two-Sided NSS w-of-w SRR and IRR Schemes
By following a similar procedure as that done in Tables 1 and 2, it follows that, for any w, the TPM in (4) for the upper or lower one-sided as well as the two-sided NSS w-of-w SRR and IRR schemes is given by(A.1)ϕη2η3η4⋯ηw-1ηwOOCϕp∗p#00⋯00pλη2p∗0p#0⋯00pλη3p∗00p#⋯00pλ⋮⋮⋮⋮⋮⋱⋮⋮⋮ηw-2p∗000⋯p#0pληw-1p∗000⋯0p#pληwp∗000⋯00p#+pλOOC0000⋯001with each of the probability elements as given in Table 8.
The elements of the TPM of the upper or lower one-sided and the two-sided non-side-sensitive SRR and IRR w-of-w monitoring schemes.
Type
p∗
p#
pλ
SRR
Upper one-sided
pDδ
pAδ
0
Lower one-sided
pEδ
pCδ
0
Two-sided non-side-sensitive
pBδ
pOδ
0
IRR
Upper one-sided
p6δ
p2δ
p1δ
Lower one-sided
p7δ
p4δ
p5δ
Two-sided non-side-sensitive
p3δ
p9δ
p8δ
Based on the notation introduced in Table 8, the false alarm rate (FAR) of the upper or lower one-sided as well as the two-sided NSS SRR and IRR monitoring schemes at sampling time point t is given by (A.2)FAR=pλwhen t=1,2,…,w-1pλ+p#wwhen t≥w.Using (7), it follows that the ARL vector, R, of the upper or lower one-sided as well as the two-sided NSS SRR and IRR monitoring schemes is given by(A.3)R=ς1δς2δς3δ⋮ςw-2δςw-1δςwδ=11-p∗∑i=0w-1p#i∑i=0w-1p#i∑i=0w-2p#i∑i=0w-3p#i⋮∑i=02p#i∑i=01p#i∑i=00p#i.The zero-state, w×1 initial probability vector, q=(100…0), i.e., 1 in the first position and zero elsewhere, then it follows that, using (6), the corresponding zero-state ARL (ZSARL) of the upper or lower one-sided as well as the two-sided NSS SRR and IRR monitoring schemes is given by(A.4)ZSARL=∑i=0w-1p#i1-p∗∑i=0w-1p#i.Next, as discussed in Section 3.2, the SSPV1, SSPV2, and SSPV3 of the upper or lower one-sided as well as the two-sided non-side-sensitive SRR and IRR monitoring schemes are given in Table 9, where y#=p#/(p∗+p#).
The steady-state initial probability vectors of the one-sided and two-sided NSS w-of-w SRR and IRR schemes.
s
SSPV1: SRR and IRR
SSPV3: IRR only
SSPV2: SRR and IRR
SSPV3: SRR only
s1s2s3⋮sw-2sw-1sw
1∑i=0w-1p#i1p#p#2⋮p#w-3p#w-2p#w-1
1∑i=0w-1y#i1y#y#2⋮y#w-3y#w-2y#w-1
Therefore, the corresponding steady-state ARL (SSARL) of the upper or lower one-sided as well as the two-sided NSS SRR and IRR w-of-w monitoring schemes is given by(A.5)SSARL=∑j=1wsjςjδwith ςjδ as given in (A.3) and sj as shown in Table 9.
A.2. Run-Length Properties of the Side-Sensitive Two-Sided w-of-w SRR and IRR Schemes
By following a similar procedure as that done in Tables 1 and 2, it follows that, for any w, the TPM in (4) for the SS two-sided w-of-w SRR and IRR schemes is given by(A.6)η1η2⋯ηw-4ηw-3ηw-2ηw-1ϕηw+1ηw+2ηw+3ηw+4⋯η2w-3η2w-2η2w-1OOCη1q∗qβqα+qλη2qαq∗qβqλη3qαq∗qβqλ⋮⋱⋮⋮⋮ηw-3qαq∗qβqληw-2qαq∗qβqληw-1qαq∗qβqλϕqαq∗qβqληw+1qαq∗qβqληw+2qαq∗qβqληw+3qαq∗qβqλ⋮⋮⋮⋱⋮η2w-4qαq∗qβqλη2w-3qαq∗qβqλη2w-2qαq∗qβqλη2w-1qαq∗qβ+qλOOC1with each of the probability elements as given in Table 10.
The elements of the TPM of the two-sided side-sensitive SRR and IRR w-of-w monitoring schemes.
Type
qα
q∗
qβ
qλ
Two-sided side-sensitive SRR
pAδ
pBδ
pCδ
0
Two-sided side-sensitive IRR
p2δ
p3δ
p4δ
p1δ+p5δ
Based on the notation introduced in Table 10, the FAR of the two-sided side-sensitive SRR and IRR monitoring schemes at sampling time point t is given by (A.7)FAR=qλwhen t=1,2,…,w-1qλ+qαw+qβwwhen t≥w.Using (7), it follows that the ARL vector, R, of the two-sided side-sensitive SRR and IRR monitoring schemes is given by (A.9) with (A.8)G=1-q∗-q∗∑i=1w-1qαi-q∗∑i=1w-1qβi-q∗∑i=1w-1qαi∑i=1w-1qβi-∑i=1w-1qαi∑i=1w-1qβi.(A.9)R=ς1δς2δς3δ⋮ςw-2δςw-1δςwδςw+1δςw+2δ⋮ς2w-3δς2w-2δς2w-1δ=1G∑i=00qαi∑i=0w-1qβi∑i=01qαi∑i=0w-1qβi∑i=02qαi∑i=0w-1qβi⋮∑i=0w-3qαi∑i=0w-1qβi∑i=0w-2qαi∑i=0w-1qβi∑i=0w-1qαi∑i=0w-1qβi∑i=0w-2qβi∑i=0w-1qαi∑i=0w-3qβi∑i=0w-1qαi⋮∑i=02qβi∑i=0w-1qαi∑i=01qβi∑i=0w-1qαi∑i=00qβi∑i=0w-1qαiThe zero-state, (2w-1)×1 initial probability vector, q=(00⋯010⋯00), i.e., 1 in the wth position and zero elsewhere, then it follows that, using (6), the corresponding ZSARL of the two-sided SS w-of-w SRR and IRR schemes is given by(A.10)ZSARLδ=∑i=0w-1qαi∑i=0w-1qβiG,with G defined in (A.8).
The steady-state initial probability vectors, s, of the two-sided SS w-of-w SRR and IRR schemes are given in Table 11, where q=qα=qβ (as s is calculated when the process is IC, i.e., δ = 0), c1=q∗/(q∗+q), c2=q∗/(q∗+2q), c3=q/(q∗+2q), F0=∑i=0w-1qi, F1=2∑i=0w-2qi+2c1qw-2+2q∗∑i=0w-3qi1-q∗-1 and F2=2∑i=0w-2c3i+2c1c3w-2+2c2∑i=0w-3c3i1-c2-1.
The steady-state initial probability vectors of the two-sided SS w-of-w SRR and IRR schemes.
Therefore, the corresponding steady-state ARL of the two-sided SS w-of-w SRR and IRR schemes is given by(A.11)SSARLδ=swςwδ+∑j=1w-1sjςjδ+ς2w-jδwith ςjδ as given in (A.9) and sj as shown in Table 11.
A.3. Zero-State and Steady-State EWRL Expressions
Based on the expressions derived in Appendices A.1 and A.2, in this section, we give the expressions of the EARL and EQL used to calculate the zero-state and steady-state empirical values in Section 4 for the one-sided and two-sided SRR and IRR monitoring schemes.
For the upper or lower one-sided as well as the two-sided NSS SRR and IRR w-of-w monitoring schemes, using the notation in Table 8, the zero-state EWRL is given by(A.12)1δmax-δmin∫δminδmaxwδ×∑i=0w-1p#i1-p∗∑i=0w-1p#idδwhereas the steady-state EWRL is given by(A.13)1δmax-δmin∫δminδmaxwδ×∑j=1wsjςjδdδwith ςjδ as given in (A.3), sj as shown in Table 9, and wδ defined in (10).
For the two-sided SS SRR and IRR w-of-w monitoring schemes, using the notation in Table 10, the zero-state EWRL is given by(A.14)1δmax-δmin∫δminδmaxwδ×∑i=0w-1qαi∑i=0w-1qβiGdδwith G defined in (A.8), whereas the steady-state EWRL is given by(A.15)1δmax-δmin∫δminδmaxwδ×swςwδ+∑j=1w-1sjςjδ+ς2w-jδdδwith ςjδ as given in (A.9), sj as shown in Table 11, and wδ defined in (10).
Data Availability
SAS 9.3 programs are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors would like to thank University of South Africa for financial support.
BalakrishnanN.KoutrasM. V.2002New York, NY, USAJohn Wiley & Sons, Inc.MR1882476WeilerH.The Use of Runs to Control the Mean in Quality Control1953482648168252-s2.0-034775631910.1080/01621459.1953.10501203PageE. S.Control charts with warning lines1955421-2243257Western Electric Handbook1956Indianapolis, Ind, USAChampC. W.WoodallW. H.Exact results for Shewhart control charts with supplementary runs rule198729439339910.1080/00401706.1987.104882662-s2.0-0023451218ChampC. W.Steady-state run length analysis of a shewhart quality control chart with supplementary runs rules19922137657772-s2.0-004085503210.1080/03610929208830813Zbl0800.62635DermanC.RossS. M.1997San Diego, Calif, USAAcademic PressKleinM.Two alternatives to the Shewhart X- control chart200032442743110.1080/00224065.2000.119800282-s2.0-0038723914ZhangS.WuZ.Designs of control charts with supplementary runs rules200549176972-s2.0-2314445466210.1016/j.cie.2005.02.002KhooM. B. C.AriffinK. N. B.Two improved runs rules for the Shewhart X- control chart200618217317810.1080/089821106005675172-s2.0-33645140625ChakrabortiS.EryilmazS.A nonparametric Shewhart-type signed-rank control chart based on runs200736233535610.1080/03610910601158427MR2370905AntzoulakosD. L.RakitzisA. C.The modified r out of m control chart200837239640810.1080/03610910701501906MR2422896AntzoulakosD. L.RakitzisA. C.The revised m-of-k runs rule2008201758110.1080/089821107016364012-s2.0-37549045606ZhangY.CastagliolaP.Run rules X- charts when process parameters are unknown201017438139910.1142/S02185393100038582-s2.0-78249285976KritzingerP.HumanS. W.ChakrabortiS.Improved Shewhart-type runs-rules nonparametric sign charts201443224723474810.1080/03610926.2012.729637MR3276972Zbl06414142KhilareS. K.ShirkeD. T.Fraction nonconforming control charts with m-of-m runs-rules2015785-81305131410.1007/s00170-014-6735-12-s2.0-84940007858Malela-MajikaJ.-C.ChakrabortiS.GrahamM. A.Distribution-free precedence control charts with improved runs-rules201632442343910.1002/asmb.2159MR35396582-s2.0-84992293253Malela-MajikaJ.-C.ChakrabortiS.GrahamM. A.Distribution-free Phase II Mann–Whitney control charts with runs-rules201686172373510.1007/s00170-015-8083-12-s2.0-84951801028RakitzisA. C.On the performance of modified runs rules charts with estimated parameters20174621360138010.1080/03610918.2014.1002846MR3606887KumarN.ChakrabortiS.RakitzisA. C.Improved Shewhart-type charts for monitoring time between events201749327829610.1080/00224065.2017.119179952-s2.0-85031087463MehmoodR.LeeM. H.HussainS.RiazM.On efficient construction and evaluation of runs rules-based control chart for known and unknown parameters under different distributions201810.1002/qre.2423KoutrasM. V.BersimisS.MaravelakisP. E.Statistical process control using Shewhart control charts with supplementary runs rules20079220722410.1007/s11009-007-9016-8MR2404264Zbl1171.623662-s2.0-34248161040ChakrabortiS.GrahamM. A.2019WileyAbbasN.RiazM.DoesR. J. M. M.Enhancing the performance of EWMA charts20112768218332-s2.0-8005293645710.1002/qre.1175RiazM.AbbasN.DoesR. J. M. M.Improving the performance of CUSUM charts20112744154242-s2.0-7995747964410.1002/qre.1124KhooM. B. C.CastagliolaP.LiewJ. Y.TeohW. L.MaravelakisP. E.A study on EWMA charts with runs-rules - the Markov chain approach201645154156418010.1080/03610926.2014.917187MR3505437ArshadW.AbbasN.RiazM.HussainZ.Simultaneous Use of Runs Rules and Auxiliary Information With Exponentially Weighted Moving Average Control Charts20173323233362-s2.0-8501183548110.1002/qre.2007MaravelakisP. E.CastagliolaP.KhooM. B. C.Run length properties of run rules EWMA chart using integral equations201710.1080/16843703.2017.1372853ShmueliG.CohenA.Run-length distribution for control charts with runs and scans rules200332247549510.1081/STA-120018196MR1958524Zbl1024.62049Acosta-MejiaC. A.Two sets of runs rules for the X- chart200719212913610.1080/175134707012636412-s2.0-34047220204LimT.-J.ChoM.Design of control charts with m-of-m runs rules20092581085110110.1002/qre.10232-s2.0-73449099980FuJ. C.LouW. Y. W.2003SingaporeWorld Scientific Publishing10.1142/4669MR2000533CrosierR. B.A new two-sided cumulative sum quality control scheme198628218719410.1080/00401706.1986.104881262-s2.0-0022758777DavisR. B.WoodallW. H.Evaluating and improving the synthetic control chart20023422002082-s2.0-164251982010.1080/00224065.2002.11980146HaqA.KhooM. B. C.A new synthetic control chart for monitoring process mean using auxiliary information201686153068309210.1080/00949655.2016.1150477MR3523162MachadoM. A. G.CostaA. F. B.A side-sensitive synthetic chart combined with an X- chart201452113404341610.1080/00207543.2013.8792212-s2.0-84898013561WuZ.YangM.JiangW.KhooM. B. C.Optimization designs of the combined Shewhart-CUSUM control charts200853249650610.1016/j.csda.2008.08.032MR2649104RyuJ.-H.WanH.KimS.Optimal design of a cusum chart for a mean shift of unknown size20104233113262-s2.0-7995172364210.1080/00224065.2010.11917826TranK. P.CastagliolaP.BalakrishnanN.On the Performance of Shewhart median Chart in the Presence of Measurement Errors2017335101910292-s2.0-8499568538210.1002/qre.2087YouH. W.Run length distribution of synthetic double sampling chart201712241426814272YouH. W.Performance of the synthetic double sampling chart with estimated parameters based on expected average run length201820186758361010.1155/2018/7583610MR3812833MontgomeryD. C.20137thSingaporeJohn Wiley & Sons, LtdShongweS. C.GrahamM. A.Synthetic and Runs-Rules Charts Combined With an X- Chart: Theoretical Discussion201733173510.1002/qre.19872-s2.0-84966355946ShongweS. C.GrahamM. A.Some theoretical comments regarding the run-length properties of the synthetic and runs-rules X- monitoring schemes – Part 1: Zero-state201710.1080/16843703.2017.1389141MR3740814ShongweS. C.GrahamM. A.Some theoretical comments regarding the run-length properties of the synthetic and runs-rules X- monitoring schemes – Part 2: Steady-state201710.1080/16843703.2017.1389142MR3740814