Enhanced Cumulative Sum Charts for Monitoring Process Dispersion

The cumulative sum (CUSUM) control chart is widely used in industry for the detection of small and moderate shifts in process location and dispersion. For efficient monitoring of process variability, we present several CUSUM control charts for monitoring changes in standard deviation of a normal process. The newly developed control charts based on well-structured sampling techniques - extreme ranked set sampling, extreme double ranked set sampling and double extreme ranked set sampling, have significantly enhanced CUSUM chart ability to detect a wide range of shifts in process variability. The relative performances of the proposed CUSUM scale charts are evaluated in terms of the average run length (ARL) and standard deviation of run length, for point shift in variability. Moreover, for overall performance, we implore the use of the average ratio ARL and average extra quadratic loss. A comparison of the proposed CUSUM control charts with the classical CUSUM R chart, the classical CUSUM S chart, the fast initial response (FIR) CUSUM R chart, the FIR CUSUM S chart, the ranked set sampling (RSS) based CUSUM R chart and the RSS based CUSUM S chart, among others, are presented. An illustrative example using real dataset is given to demonstrate the practicability of the application of the proposed schemes.


Introduction
Cumulative sum (CUSUM) control charts have received a great deal of attention in modern industries to monitor unusual variations in manufacturing and service processes. There is a CUSUM location chart for monitoring changes in process mean and CUSUM scale or dispersion chart for monitoring changes in process variability. Many authors have contributed to the theory of the CUSUM charts with emphasis on location charts. However, in reality, monitoring changes in process variability is equally important [1]. An increase in process variance could lead to increase in the number of defective items, and a decrease in process variability made more items closer to the target value resulting in an improved process capability [2]. Scale or dispersion charts are also important in the interpretation of location charts, as the process standard deviation is assumed constant [2].
A number of researchers with the pioneering work of Page in 1963 have studied the use of the CUSUM control chart to monitor process variability. Page [3] studied a CUSUM chart based on subgroup range to monitor increases in the process dispersion. Tuprah [4] evaluated the performance of a CUSUM chart based on subgroup standard deviation. Chang and Gan [5] studied CUSUM chart based of logarithmic transformation of subgroup variance. Acosta-Mejia et al. [2] examined and compared the performance of several CUSUM charts for monitoring process dispersion. Some recent contributions to the advancement of CUSUM chart for monitoring process dispersion include [6,7] and references therein.
Ranked set sampling (RSS) is a statistical technique for data collection that is currently gaining popularity in the quality control chart. McIntyre [8] introduced the scheme while Takahashi and Wakimoto [9] developed its mathematical theory. RSS has been used to enhance the performance of Shewhart Àx charts for monitoring the process mean [10,11,12]. Recently, the RSS methodology has found its way into the memory control charts such as the exponentially weighted moving average (EWMA) and the CUSUM [13,14,15] among others. It is well known that these memory charts are more efficient in the detection of small and moderate process shifts than the traditional Shewhart-type charts.
The RSS scheme employs judgment ordering and utilizes extra information from each specific unit to obtain actual samples. Hence, it provides true representative samples about the target population when compared to the traditional simple random sampling (SRS) for the same number of sample observations. In literature, most of the RSS applications in control charting focus on the location chart for monitoring the process mean. Monitoring the process variance based on ranked data have received less attention and is the focus of this paper. In this article, new CUSUM control charts based on extreme RSS, extreme double RSS and double extreme RSS are proposed to further enhance the overall performance of a CUSUM scale chart at detecting a wide range of shifts in process variability.
The rest of the article is organized as follows: Section 2 introduces the RSS method and some of its extreme variations. In Section 3, we present the design structure for the proposed CUSUM chart using the extreme RSS scheme and variants. In Section 4, we discuss the computations of average run length (ARL) and average extra quadratic loss (AEQL) for the proposed charts via Monte Carlo simulation. A comparison of several control charts for monitoring changes in the process dispersion is given in Section 5. Some graphic representations are also given to aid comparison. In Section 6, we present a practical application example of the proposed control schemes using a real dataset. Finally, some concluding remarks are given in Section 7.

RSS Properties and Variants
The method of RSS introduced by McIntyre involves the drawing of n simple random samples each of subgroup size n units from the target population. Rank the units in each set by visual comparison, some less-expensive method, or using auxiliary variables, with respect to a variable of interest. The smallest ranked unit is selected from the first set; the second smallest ranked unit is selected from the second set and so on until the largest ranked unit is selected from the last set. The cycle may be repeated r times to obtain a ranked set sample of nr units. It is important to point out that although RSS requires the ranking of n by n sample units, only n observations are actually measured accurately, and hence, makes the comparison with SRS of same size n meaningful [16].
Let X i1 , X i2 , X i3 ,. . ., X in , for i = 1, 2,. . ., n be n independent simple random samples each of size n from a continuous distribution with probability density function f(x) and cumulative distribution function F(x). The simple random sample estimates of the population range, R, mean, μ, and variance, σ 2 , are respectively given by the sample range R srs , sample mean ÀX srs and sample variance S 2 srs . Let X rss = {X (i:n)j , i = 1, 2,. . ., n, j = 1, 2, . . ., r} denotes the ith order statistic in ith sample of size n in jth cycle under perfect ranking. Then the sample range based on ranked set sample of subgroup size n in jth cycle is defined as are control chart constants that are subgroup size n dependants. We have provided the estimates for the values of d Ã 2 and d Ã 3 for sample sizes of two to ten, via Monte Carlo simulation of 100,000 iterations, with r = 100 based of different sampling schemes (cf. Table 1).
For the estimate of population variance based on ranked set samples, Stokes [17] gave a natural estimator define by where X À rss ¼ 1=nr ð Þ P n i¼1 P r j¼1 X ði:nÞj is an unbiased estimator of the population mean, [9]. The expected value of S 2 rssðsÞ is E S 2 where μ (i:n) is the mean of ith order statistic defined by However, unlike S 2 srs which is an unbiased estimator of population variance, S 2 rssðsÞ is biased. To obtain an estimator, which is unbiased, a modified version of S 2 rssðsÞ [16,18], is expressed as where u rss ¼ 1=n ð Þ P n i ¼ 1 u ði:nÞ is a correction constant that depends on sample size n and υ (i: n) = [μ (i:n) -μ]/σ. We also define E S rss ð Þ ¼ c Ãrss 4 s and s S rss ð Þ ¼ c Ãrss 5 s, where c Ãrss 4 and c Ãrss 5 are the S chart constants that depends on subgroup size n. We have as well provided the estimates via Monte Carlo simulation of 10 7 iterations with single cycle, for the values of c Ã 4 and c Ã 5 based on sample sizes of n = 2 to 10 using different sampling strategies in Table 2.
In RSS scheme, the ranking of units in a subgroup with respect to the variable of interest may be perfect (without ranking errors) or imperfect (with ranking errors or using auxiliary variables). In reality, imperfect ranking may be unavoidable and the accuracy of ranking depends on the linear relationship between the variable of interest (X) and its auxiliary variable (Y). Suppose that (X, Y) has a bivariate normal distribution and assume the regression of the main variable of interest X on the auxiliary variable Y is linear. Then following Stokes [19], where μ x , μ y , σ x and σ y are the population means and standard deviations of X and Y, respectively, and ρ xy is the correlation between X and Y variables. The ε is the error term with mean 0 and variance s 2 . Let X [i:n]j denotes the ith judgment order statistic of the variable of interest X based on auxiliary variable Y, then the above equation can be written as where Y [i:n]j is the ith order statistic based on perfect ranking in jth cycle of subgroup size n.
In practice, ranking of quality characteristics of interest of subgroup sizes greater than four may be difficult, [8]. A very important variation of RSS, which is easier to execute and has the potential of reducing errors associated with ranking, most especially in symmetric distributions, is the extreme RSS (ERSS) introduced by [20]. Some contributions to the theory of ERSS can be found in [21,22,23,24]. The procedure of the scheme is summarized as follows: • Randomly select n random samples each of size n units from the target population.
• Rank the units within each set by visual inspection or some less-expensive method, with respect to the quality characteristics of interest.
• If the subgroup size is even, select the lowest rank unit from the first n/2 samples in a subgroup and the highest rank unit from the second n/2 samples.
• For odd subgroup size, measure the median observation from one sample, and select the lowest rank from the first (n-1)/2 samples and the highest rank from the other (n-1)/ 2 samples.
• The steps may be repeated r times to obtain a sample of nr units of ERSS data. There are also double extensions of the ERSS schemes namely, extreme double RSS (EDRSS) and double ERSS (DERSS) which are more efficient than the RSS and ERSS. The EDRSS is a two-stage procedure in which RSS is selected on the first stage and ERSS on the second stage. The DERSS, on the other hand, is a two-stage ERSS. That is, select ERSS on both the first and second stages [25]. These extreme variations of RSS are very effective in the estimation of population standard deviation σ than the RSS and its median variations [26], if the underlying distribution is normal.

New CUSUM Charts for Process Dispersion
Since the first study of Page [3] in which he applied a CUSUM on a subgroup sample ranges to monitor process variability, many researchers have proposed different enhancements to the CUSUM dispersion chart-for example, see [2,5,6,7]. Most of these studies, however, are based on the assumption that samples are drawn from a production line using the traditional SRS. Recent studies by Muttlak and Al-Sabah [22], Abujiya and Muttlak [11], Al-Sabah [27], Jafari and Mirkamali [28], Abujiya et al. [15], Mehmood et al. [12], Abujiya et al. [29] and Mehmood et al. [30] among others, have shown that the use of a well-structured sampling technique such as RSS in control charting had significantly improved the performance of the charts over the random sampling.
To increase the sensitivity of CUSUM scale charts at detecting wide range of variability, we studied CUSUM charts based on R and S statistics for monitoring the dispersion of a normal process using ERSS, EDRSS and DERSS. For simplicity, these will be referred to as scheme I, II and III, respectively. Throughout this article, we assume that the in-control mean μ 0 and standard deviation σ 0 are known or can be estimated from the preliminary in-control process samples. Also, let σ 1 6 ¼ σ 0 be the standard deviation that needs to be detected assuming that the process mean does not change.

The CUSUM-R chart
In this subsection, we give the design structure of CUSUM of the range R for monitoring the process depression based on the extreme variations of RSS. Let R erssj , R edrssj and R derssj denote the sample ranges of the jth cycle based on schemes I, II and III respectively. Then the two one -sided standardized CUSUM statistics for detecting increases and decreases in process variability using scheme I, can be expressed as (cf. [2,6]): where k erss is called reference value defined by k erss ¼ d Ãerss ÞE R erssj and τ = σ 1 / σ 0 is the shift in process standard deviation. Here, E(R erssj ) can be estimated by ÀR erss ¼ 1=r ð Þ P r j ¼ 1 R erssj . If either C þ erssj or C À erssj exceeds the predetermined control limit h for the first time, then an out-of-control signal is given and j is the run length of the distribution. For a zero head start, we set the starting values C þ erss 0 and C À erss 0 equal to zero. The advantage of standardizing CUSUM is that it makes parameters of the chart independent of σ [2]. Similarly, for schemes II and III, we considered Although sample range is the simplest measure of variation, it is not always as efficient as the standard deviation, which is the natural measure of variation, [1]. Therefore, we present the CUSUM of S chart in the next section.

The CUSUM-S chart
Let X erssj denotes the extreme ranked set samples of subgroup size n in jth cycle. Then the estimator of the population standard deviation can be expressed as where ÀX erssj ¼ 1=n ð Þ P n i¼1 X erssj and υ erss is the bias correction constant similar to υ rss of section 2. The standardized two one-sided CUSUM of standard deviation S for monitoring the variation of a process based on scheme I is given by In a similar manner, we define the standardized CUSUM statistics for scheme II as In addition, for scheme III, we have Þ. An alarm is given when any of the CUSUM statistics in Eqs (13) to (18) exceeds a control limit for the first time. We now evaluate the performance of the proposed CUSUM charts for monitoring process variability.

Performance Measure of the New Charts
The statistical performances of control charts are generally evaluated in terms of the ARL, the mean of the run length (RL) distribution. In other words, ARL represents the average number of observations plotted on a control chart until the chart gives an out-of-control signal or produces a false alarm. This performance measure may be computed by integral equation approach originally introduced by Page [31], Markov Chain approach of Brook and Evans [32], or a Monte Carlo simulation approach adopted by Hawkins [33], Chang and Gan [5], Mehmood et al. [12], Haq et al. [14] among others. In this article, we used the Monte Carlo simulation approach through an algorithm developed in FORTRAN.
In addition to the computation of ARL, we have also computed the standard deviation of RL (SDRL) distribution often used by researchers as a supplementary indicator of the effectiveness of a control chart. Usually, control charts with smaller out-of-control ARL (ARL 1 ) or SDRL (SDRL 1 ) value at more points are considered more efficient than the others. But Reynolds and Stoumbos [34] pointed out that, in terms of ARL 1 and SDRL 1 , no chart will give better performance than the others for all shift values τ. Since our objective is to enhance the performance CUSUM control charts at detecting a wide range of process variability across all shifts rather than at a particular point, following Wu et al. [35] and Wu et al. [36], we used the AEQL.
The AEQLmeasures the overall performance of a control chart across a range of shifts rather than at a particular shift [36]. It is defined by: where ARL(τ) is the ARLvalues of a particular control chart at standard deviation shift τ. τ max and τ min are the upper and lower bounds of shifts in process standard deviation, respectively. The AEQL is based on a loss function, which takes all factors that affect the quality cost into account [36]. The AEQL of a CUSUM control chart can be computed by integration over a range of shifts in standard deviation. The smaller the AEQL value of a control chart, the better the overall statistical performance of the chart. If we assume that the observations are normally distributed with mean μ 0 and standard deviation σ 0 and that the process is initially in-control and remains in-control until it goes out-ofcontrol with a shift in process standard deviation from σ 0 to σ 1 = τσ 0 . The shift in standard deviation can be expressed as τ = σ 1 /σ 0 . If τ = 1, the process is in-control, while an out-of-control process is indicated by τ > 1 and τ < 1, representing increase and decease in process variability respectively. Furthermore, we assume that the in-control process has mean μ 0 = 0 and standard deviation σ 0 = 1 without loss of generality. All the proposed charts were designed to have incontrol ARL 0 values of 200 with subgroup size of n = 5. Using the right combinations of the chart parameters, h and k, each chart is design to detect 20%, 30%, 40% and 50% increases and decreases in process standard deviation. We have evaluated the run-length properties using simulations of 100,000 iterations for each shift τ in process standard deviation, and the results are summarized in Tables 3-6. Also presented in Tables 3-6, is the overall performance measures in terms of AEQL.
Based on the results in Tables 3-6, we summarize our findings for the proposed CUSUM scale control charts as follows: 1. There are no significant differences between the in-control ARL 0 values of the proposed charts and the corresponding SDRL 0 values (cf. Tables 3-6).
2. As the shifts in standard deviation increase (τ > 1), Tables 3 and 4, or decrease (τ < 1), Tables 5 and 6, both the out-of-control ARL 1 and SDRL 1 decrease rapidly, with h inversely proportional to k (for increases in dispersion) and directly proportional to k (for decreases in dispersion).
3. All the proposed charts have greatly enhanced the detection ability of the CUSUM scale chart to a wide range of increases and decrease in the dispersion of a process (cf. Tables 3-6).
4. The proposed schemes are generally more effective at detecting decreases in the process variability (cf. Tables 5 and 6) than for monitoring increases (cf. Tables 3 and 4) in standard deviation, in terms of low, ARL 1 and SDRL 1 .
5. The out-of-control ARL 1 and SDRL 1 of proposed CUSUM of R and S charts for monitoring increases in the dispersion are, however, smaller than those of decreases in a very small region. Less than fifteen and twenty-five percentage increases for R and S charts, respectively.
6. The ARL performance for increases in standard deviation indicates that the proposed R and S charts are particularly very effective at detecting small changes in variation when the  Tables 3 and 4) and vice versa for decreases in variation.
7. The overall performance of the proposed charts, in terms of AEQL, indicate significant gains in efficiency of the proposed CUSUM charts, particularly at a detection ability of 40% increase or 30% decrease in standard deviation except for increases in variation based on scheme I (cf. Tables 3-6).
8. The proposed CUSUM of S charts appears to outperform all the corresponding proposed CUSUM of R charts in terms of ARL 1 , SDRL 1 and AEQL, particularly in the detection of increases in variability (cf. Tables 3-6).
9. The significant improvement archived, in terms of ARL 1 , SDRL 1 and AEQL, when monitoring decreases in process variability using CUSUM of S chart (cf. Table 5) are practically comparable to CUSUM of R chart (cf. Table 6).
10. The proposed CUSUM of R and S charts based on scheme III appeared to be more effective than the corresponding CUSUM of R and S charts based on schemes I and II (cf. Tables 3-6).
11. Of all the proposed charts, the CUSUM of S charts based on scheme III dominates all other schemes at detecting various changes in process standard deviation for all combinations of h and k. The proposed scheme I is the least performing (cf. Tables 3-6). 12. In general, all the proposed charts not only maintain the advantages of CUSUM to detect small shifts but have also greatly increased the sensitivity of the charts in detecting very large process shifts in standard deviation (cf. Tables 3-6).

Comparisons of the Control Charts
In this section, we compare the performance of the proposed CUSUM scale charts with several control charts for monitoring process dispersion. These include: the classical CUSUM of the range R, the classical CUSUM of the standard deviation S, the FIR CUSUM of R chart, the FIR CUSUM of S chart, the RSS CUSUM of R chart and the RSS CUSUM of S chart. The discussions here are based on subgroup size of n = 5, and the cases based on other sample sizes have similar conclusions, [5]. Assuming that the underlying distribution is normal and using an incontrol ARL 0 of 500, the control charts are designed to detect 30% increase and decrease in process standard deviation. The performances of these control charts in monitoring changes in the dispersion are compared in terms of their ARL values. For a better comparison, we also implore the overall performance measures-the AEQL, the Average Ratio of ARL (ARARL) and the Performance Comparison Index (PCI). The ARARLmeasures the overall performance of a control chart across a range of shifts rather than at a particular shift [35]. It is defined by: where ARL(τ) and ARL(τ) benchmark are respectively, the ARL values of a particular control chart and a benchmark control chart at τ. The benchmark chart is the one with the smallest ARL 1 . τ max and τ min are the upper and lower bounds of shifts process standard deviation, respectively. The ARARLof a CUSUM control chart can be computed by numerical methods. Similar to ARARL, the PCI is the ratio between the AEQL of a particular control chart and the AEQLof the benchmark control chart. However, unlike the ARARL, it does not involve integration over a range of shifts. It is defined as: where the AEQL benchmark is the AEQL value of the best-performing control chart [37]. A control chart with minimal ARARL gives smaller AEQL and PCI values and is subsequently more effective in detecting changes in the process standard deviation. The ARL, ARARL, AEQL and

Comparison with classical CUSUM R chart
The classical CUSUM of the range R chart introduced by Page [3] outperforms the traditional Shewhart R chart as well as the warning line scheme for the detection of changes in process variability [2,3]. The chart is based on SRS and not only has the ability to give the position of the departure from the target value but also the amount of departures. The ARL performances of the classical one-sided CUSUM R chart are presented in column two of Tables 7 and 8. Comparison of the control charts indicates that the proposed schemes have smaller ARL 1 than the classical CUSUM of R chart. This is equally supported by the ARARL, AEQL and PCI values. All the proposed schemes have greater sensitivity of detecting changes in process standard deviation than the classical CUSUM of R chart (cf . Tables 7 and 8).

Comparison with classical CUSUM S chart
Acosta et al. [2] considered the classical CUSUM of the standard deviation S chart. This control chart provides better sensitivity for detecting small and moderate changes in the variation of a process than the Shewhart S chart [3]. The ARL performance of the classical CUSUM S chart in column three of Tables 7 and 8 also indicates the superiority of the chart over the classical CUSUM R chart but not as effective as the proposed schemes. In other words, all the proposed charts have smaller ARL 1 than the classical CUSUM S chart. Furthermore, the values of ARARL, AEQL and PCI in Tables 7 and 8 also reveal the dominance of the proposed schemes over the classical CUSUM S chart.

Comparison with RSS CUSUM R chart
Al-Sabah [27] presented the RSS based CUSUM chart for location to further enhance the detection ability of a standard CUSUM chart for monitoring changes in process mean level. As an enhancement to the CUSUM scale chart, we studied the performance of the CUSUM of R chart under RSS. The ARL values of RSS based CUSUM of R chart is presented in column four of Tables 7 and 8. Results show that this scheme has smaller ARL 1 compared to the classical CUSUM of R and S charts but is less effective than the proposed schemes. From the overall point of view, all the proposed schemes perform significantly better than the RSS CUSUM R chart (cf . Tables 7 and 8). For example, the RSS based CUSUM R chart for monitoring increases in standard deviation is inferior to its CUSUM R chart counterpart based scheme I by 29%, scheme II by 47% and scheme III by over 100% in terms of ARARL and PCI (cf. Table 7).

Comparison with RSS CUSUM S chart
In the same line with the RSS based CUSUM R chart, we equally studied the performance of the CUSUM of S chart using ranked set samples, and the results obtained for detection of  Tables 7 and 8 also reveal the superiority of the proposed schemes over the RSS CUSUM S control chart.

Comparison with FIR CUSUM R chart
Lucas and Crosier [38] suggested and analyzed the effect of the fast initial response (FIR) feature on CUSUM chart. The feature sets the initial value of the upper and lower CUSUM statistic to a specified positive value, usually h/2. Here, h is the predetermined control limit of a CUSUM chart. The FIR feature generally minimizes the out-of-control ARL 1 values and as an enhancement, the feature is implemented with the CUSUM of R chart. We present the results for monitoring increases and decrease in variability in column six of Tables 7 and 8, respectively. The ARL 1 values of the FIR CUSUM R chart shows that significant improvement is achieved. The chart outperforms the RSS CUSUM S chart in all shift points except in a very small region when percentage increase or decrease is less than 30 and 15 respectively. The FIR CUSUM R chart has larger ARL 1 values than the proposed schemes and is substantially less effective at detecting ranges of shifts in process standard deviation than the proposed charts (cf .  Tables 7 and 8).

Comparison with FIR CUSUM S chart
Poetrodjojo et al. [39] studied the FIR CUSUM of standard deviation S chart using the statistics V i ¼ ffiffiffi n p S srsi s 1 =s 2 0 . They used h/2 as a head start when controlling increases in process standard deviation. Following Acosta et al. [2], we adopted the statistics S i = S srsi /σ 0 for the development of FIR CUSUM of S chart and the results are given in column seven of Tables 7 and 8. This scheme is certainly more sensitive to departures from target standard deviation than the earlier CUSUM scale charts but not as effective as the proposed schemes, in terms of ARL 1 . Furthermore, the overall performance of the proposed schemes, in terms of ARARL, AEQL and PCI also indicate the superiority of the proposed charts over the FIR CUSUM of S chart except for the CUSUM R chart of scheme I which are comparable. For example, the ARARL and PCI values of upper one-sided CUSUM chart reveal that the FIR CUSUM of S chart is less effective than the corresponding CUSUM of S chart based on scheme I by an average 10%, scheme II by an average of 31% and scheme III by average of 62% (cf. Table 7).

The issue of imperfectness in ranking
In practice, ranking of main variable of interest may not always be perfect since the accuracy of the ranking depends on operator judgment. The objective of this sub-section is to investigate the performance of the proposed schemes when ranking of units in each subgroup is not perfect or when ranking of units is carried out with help of an auxiliary variable. We replace the perfect ranking statistics X (i:n)j in design structure of the proposed CUSUM scale charts (cf. Section 3) with the imperfect ranking statistics X [i:n]j , Eq 5. It is well known that perfect ranking and random sampling are special cases of imperfect ranking with correlation coefficients of ρ xy = 1 and ρ xy = 0 respectively, [19]. The ARL values of imperfect samples using schemes I, II and III are obtained via Monte Carlo simulation based on data of size n = 5 generated from a bivariate standard normal population using ρ xy = 0.25, 0.5 and 0.75 with an in-control ARL 0 = 200. The computed ARLs, designed to detect 20% increases in standard deviation, are compared to some existing schemes-the CUSUM of lnS 2 , the CUSUM of the variance S 2 , the EWMA of lnS 2 and the EWMA of the variance S 2 , using graphical displays of ARL curves (cf . Figs 1-3).
The displayed ARL curves of imperfect cases clearly show that the proposed schemes are not doing badly even in the presence of ranking errors or not very strong correlation between the main quality characteristic of interest and its auxiliary variable. The proposed schemes II and III outperform all the control charts investigated in the comparative studies in the detection of upward shifts in standard deviation (cf . Figs 2 and 3). Scheme I appeared not to be doing very well when ρ xy < 0.75 but still dominates other control schemes when ρ xy < 0.75 (cf. Fig 1).

Practical Applications of the New Schemes
To illustrate the practical application of the proposed CUSUM schemes, we used a real dataset based on the problem of filling bottles on a Pepsi-Cola production line. The original dataset,  Enhanced CUSUM Charts for Monitoring Process Dispersion based on an example in Muttlak and Al-Sabah [10], is from a production line of the Pepsi-Cola production company, Al-Khobar, Saudi Arabia, when the process was in-control. Using the resampling approach of Takahasi and Wakimoto [9], thirty random samples each of subgroup size n = 5 were collected based on schemes I, II and III, and they reasonably satisfied the normality assumption. Fig 4 displays the standard deviations of the thirty data points for SRS, scheme I, II and III.  To measure effectively how the proposed schemes respond to the changes in process standard deviation, some disturbances were introduced to the data by increasing each of the last ten data points by thirty percent to make the process out-of-control. In other words, we assume that the process in-control until the twentieth sample value. Following Hawkins and Zamba  To see the effect of imperfect ranking using real dataset, a second set of data is collected using schemes I, II and III with errors in ranking of the bottles with respect to the volume of liquid. In perfect ranking procedure, we placed five randomly selected Pepsi-Cola bottles next to each other, and then ranked them visually with respect to the level of liquid. While in imperfect ranking, a random set of five bottles next to each other in a moving production line is identified and ranked visually as the bottles are moving on the production line. Like the first dataset, thirty random imperfect samples each of size n = 5 were obtained using the resampling method of Takahasi and Wakimoto [9], from the original imperfect RSS and ERSS dataset in an example in Muttlak and Al-Sabah [10]. We follow the design structure for the first dataset to develop CUSUM S charts based on imperfect ranking dataset and the results are displayed in Figs 8-10.
Based on the real dataset used in this application, we observe that the proposed schemes have less variation and are more effective at detecting shifts in process standard deviation than the classical CUSUM S chart. Although all the schemes indicate that the process drifts out-ofcontrol from the twenty-first observation, the classical CUSUM S chart only gives an out-ofcontrol signal on the twenty-seventh sample, allowing some tiny changes escape detection (cf. Figs 5-7). The perfect cases of propose schemes I, II and III give out-of-control signals at twenty-fifth, twenty-third and twenty-first sample points, respectively. This means that the proposed schemes I, II and III are more effective in detecting tiny changes than the classical scheme with the proposed scheme III giving the correct signal. The imperfect cases given in Figs 8-10, also indicates that the proposed charts are performing better than the classical scheme even with the ranking errors. Conclusions based on CUSUM R chart using the real dataset were found to be similar.

Conclusions
This article proposed some enhanced CUSUM of R and S control charts based on ERSS, EDRSS and DERSS sampling techniques for monitoring process dispersion. Using Monte Carlo simulations, the performance measures in terms of ARL, SDRL, AEQL, ARARL and PCI, of the proposed CUSUM of R and S schemes were computed. The performances of the proposed schemes were compared with classical CUSUM R, classical CUSUM S, FIR CUSUM R, FIR CUSUM S, RSS CUSUM R and RSS CUSUM S control charts when standard deviation shifts from its target value. All the proposed schemes were found to be more efficient in detecting a wide range of shifts in the process dispersion than other schemes.
We also examined the performance of the proposed schemes when ranking of the main characteristic of interest is based on auxiliary variable and compared them to CUSUM lnS 2 , CUSUM S 2 , EWMA lnS 2 and EWMA S 2 . It is observed that even in the presence of ranking errors, the proposed schemes still dominate other procedures. Of all the proposed, the CUSUM S chart based on scheme III have the best overall performance and CUSUM R chart based on scheme I is the least efficient performance scheme. We have provided real data application of the proposed schemes to demonstrate the practicability of the procedures. In addition, provided are the constant values for R and S charts of normal order statistics based on different sampling schemes can be used as design aid for the charts. The recommendation is to use the proposed scheme III since it performs better than any other control chart in detecting both increases and decreases in process standard deviation. The scope of this work may be extended to the EWMA charts and multivariate control chart structures.