Research Article  Open Access
Performance of Synthetic Double Sampling Chart with Estimated Parameters Based on Expected Average Run Length
Abstract
A synthetic double sampling (SDS) chart is commonly evaluated based on the assumption that process parameters (namely, mean and standard deviation) are known. However, the process parameters are usually unknown and must be estimated from an incontrol PhaseI dataset. This will lead to deterioration in the performance of a control chart. The average run length (ARL) has been implemented as the common performance measure in process monitoring of the SDS chart. Computation of ARL requires practitioners to determine shift size in advance. However, this requirement is too restricted as practitioners may not have the experience to specify the shift size in advance. Thus, the expected average run length (EARL) is introduced to assess the performance of the SDS chart when the shift size is random. In this paper, the SDS chart, with known and estimated process parameters, was evaluated based on EARL and compared with the performance measure, ARL.
1. Introduction
The quality of products and services is an essential factor in the world of business [1]. In order to enhance the quality of products and services, statistical process control (SPC) is used to monitor and attain the process of manufacturing and services. Among the SPC techniques, control charts are one of the most prominent techniques for detecting shifts in a process. The first control chart was proposed by Dr. Walter A. Shewhart, and it was named the Shewhart chart. The Shewhart chart is frequently used to detect large process mean shifts due to its simplicity [2]. However, the main limitation of the Shewhart chart is that it is insensitive in detecting moderate and small process mean shifts.
In recent years, Khoo et al. [3] suggested synthetic double sampling (SDS) chart, which combines double sampling (DS) subchart and conforming run length (CRL) subchart. From the findings, the SDS chart is efficient in detecting moderate and small process mean shifts compared to the synthetic chart and double sampling chart. The implementation of the SDS chart is based on the assumption that the process parameters are known. Nevertheless, the process parameters are generally unknown in practice. This requires for an estimation of the process parameters from the incontrol PhaseI samples.
Saleh and Mahmoud [4] claimed that when the estimated process parameters are used in place of the known process parameters, the performance of the control chart is affected due to the existence of variability in the estimation. Woodall and Montgomery [5] recognised the importance of examining the effects of parameter estimation on the performance of a particular control chart. Therefore, the effects of parameter estimation need to be considered when designing a control chart. Among others, [6–8] have examined the impacts of PhaseI parameter estimation on the performance of a control chart.
The performance of a control chart is crucial in determining the appropriate control chart to be implemented in a process. A common performance measure in process monitoring is average run length (ARL). ARL is the number of samples (on average) plotted on a control chart before it signals an outofcontrol [9]. By using ARL as a performance measure, the chart’s user needs to determine the process shift size.
You et al. [10] investigated the ARL performance of the SDS chart when process parameters were estimated and this has motivated the current research work. In most practical situations, it is usual that the next shift size is unknown in advance [11, 12]. To circumvent this problem, the expected average run length (EARL) is proposed to examine the performance of the SDS chart when the process shift size is unknown and random. Hence, the performance of the SDS chart, with known and estimated process parameters, is investigated using EARL. In addition, the proposed alternative performance measure, i.e., EARL of the SDS chart, will be compared with the corresponding chart using ARL.
The rest of this paper is structured as follows: Section 2 presents the operation and steps to implement the SDS chart. Moreover, the run length properties of the SDS chart with known and estimated process parameters are also given in Section 2. Section 3 illustrates the performance comparison of the SDS chart, based on EARL and ARL, for known process parameters and that of the corresponding chart with estimated process parameters. Finally, concluding remarks are drawn in the last section.
2. Materials and Methods
Khoo et al. [3] developed the SDS chart, which comprises the DS subchart and a CRL subchart. The CRL subchart is an attribute chart, with one lower limit, i.e., . Figure 1 illustrates the operation of the SDS chart with known process parameters.
(a) First sample
(b) Combined samples
Step 1. Set the charting parameters , , , , , and .
Step 2. At sampling time , take the first sample of size and the sample mean, , is calculated.
Step 3. Compute the standardised statistic corresponding to the first sample.
Step 4. If is in , the th sampling time is conforming and the control flow returns to Step 2.
Step 5. If is in , the th sampling time is nonconforming and the control flow goes to Step 9.
Step 6. If is in , take a second sample with and compute the sample mean, .
Step 7. Calculate the sample mean of the combined samples and the standardised statistic corresponding to the combined samples.
Step 8. If Z_{i} is in , and the th sampling time is conforming, then return to Step 2. Otherwise, the sampling time is regarded as nonconforming and the control flow goes to Step 9.
Step 9. Count the number of inspected sampling times between the present and last nonconforming sampling times inclusive of the present nonconforming sampling time, and denote it as the CRL value.
Step 10. If CRL > , the process is incontrol and the control flow returns to Step 2. Otherwise, the process is outofcontrol and immediate actions are required to eliminate the assignable cause(s). Then, return to Step 2.
Without loss of generality, the incontrol mean, , and incontrol standard deviation, , are assumed as known. Let be the probability of deciding that a sampling time is nonconforming in the DS subchart. Note that and can be expressed as follows [3]:where and are the cumulative distribution function (cdf) and probability density function (pdf) for a standard normal random variable, respectively. Here, , , and
Finally, for the SDS chart when process parameters are known, the ARL is equal toMoreover, when the exact shift size is unknown, it is essential to consider the EARL for an overall range of shifts , where and indicate the lower and upper bounds of the mean shift, respectively. The EARL of the SDS chart with known process parameters iswhere ARL can be obtained from (2) and is the pdf of the shift size . Since the actual shape of is usually unknown, it is assumed that the shifts in process mean happen with equal probability; i.e., is uniformly distributed with . Therefore, (3) reduces toIn reality, and are unknown and need to be estimated from m PhaseI samples, each of size n, i.e., , for i = 1, 2, …, m. The estimators of and are [10]respectively. As the values of and are both unknown and need to be estimated using and , the standardised statistic for the first sample and combined samples at the sampling time of the DS subchart with estimated process parameters becomesrespectively.
Here, and represent the sample sizes of the SDS chart with the estimated process parameters. Similarly, , and correspond to the limits of the SDS chart when the process parameters are estimated.
Then, the probability for the DS subchart with the estimated process parameters to identify a nonconforming sampling time is given as ; i.e., [10],withHere, and are random variables denoted as respectively.
As , it can be deduced that the pdf of isFor the random variable , it is known that , i.e., the gamma distribution with parameters and . Using this property, the pdf of is as follows: where is the pdf of the gamma distribution with parameters and Note that, for complete and detailed derivation, reader can refer to You et al. [10].
Thus, the ARL of the SDS chart with estimated process parameters isConsequently, when the process parameters are estimated from the incontrol PhaseI samples, the computation of the EARL iswhere EARL can be obtained from (3) by replacing and with and , respectively.
3. Results and Discussion
In practice, the exact shift size of a process is unknown. In this situation, if the corresponding optimal charting parameters are employed based on a particular shift size, the performance of the control chart will be significantly different if different shift occurred in the process. Therefore, it is essential to evaluate the performance of the SDS chart using alternative performance measure, i.e., EARL. In this paper, the optimal charting parameters of the SDS chart were computed using a nonlinear minimisation problem, i.e., optimal statistical design that minimises the outofcontrol EARL . The programmes are written in the ScicosLab software version 4.4.2 (http://www.scicoslab.org).
The optimal charting parameters and the corresponding with different combinations of sample size, , the number of PhaseI samples, with and with the incontrol EARL, i.e., = 370.4, are presented in Table 1. Here, denotes the known process parameters case, while denotes the estimated process parameters case. The performance of the SDS chart for both the known and estimated process parameters cases was calculated using the optimal charting parameters in columns 4–9, which were obtained by minimising when the process parameters are known.

From Table 1, for the same , , and , the value of decreases with the increases in . This is due to the fact that as the more PhaseI samples are taken, the performance of the estimated process parameters SDS chart approaches to the corresponding chart with known process parameters; i.e., value decreases to indicate better performance. For instance, when , , and , the optimal charting parameters yield the lowest , when the process parameters are known. With these optimal charting parameters, for , respectively. It is noticeable that the for the estimated process parameters case is deviated from the known process parameters. However, the value approaches to the value which corresponds to the when the number of PhaseI samples increased. These findings show that more than 80 PhaseI samples are required to reduce the effects of process parameters estimation when estimating the process parameters from the incontrol PhaseI samples.
To illustrate the implementation of the proposed optimal charting parameters, Table 2 presents the optimal charting parameters for the known process parameters SDS chart based on minimising outofcontrol ARL and the corresponding for the same combinations of in Table 1. Here, the incontrol ARL, i.e., is intended. Note that are considered here to accommodate the that are considered in Table 1; i.e., and are included in = (0.2, 1.0) and , respectively. In Table 1, when , , and , the is obtained using the optimal charting parameters . Here, by using the same optimal charting parameters for , it yields when using the ScicolsLab program. It is observed that the value is almost the same to those in Table 2 when and , although the optimal charting parameters based on minimising are different, i.e., (see Table 2). This indicates that the optimal charting parameters obtained based on minimising can be employed as long as , i.e., when the practitioners do not have knowledge to determine the exact process shift size in advance.

4. Conclusions
In the production and manufacturing industries, it is a typical situation where quality practitioners are undecided about the process shift size to be implemented. The findings showed that the performance criterion EARL is capable of tackling the random shift size situation. Furthermore, the results also revealed that the performance of the SDS chart was adversely affected by process parameters estimation. This was proven when more than 80 PhaseI samples were needed for the chart with the estimated process parameters to behave similarly like the one with known process parameters. Therefore, future research works can consider the optimal charting parameters by minimising for the SDS chart when the process parameters are estimated.
Data Availability
The ScicosLab programs to compute the average run length and expected average run length of the synthetic double sampling chart can be requested from the author.
Conflicts of Interest
The author declares that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research is supported by the Universiti Kebangsaan Malaysia, Geran Galakan Penyelidik Muda, GGPM2017062.
References
 S. Shang, Q. Zhou, M. Liu, and Y. Shao, “Sample size calculation for controlling false discovery proportion,” Journal of Probability and Statistics, vol. 2012, Article ID 817948, 13 pages, 2012. View at: Publisher Site  Google Scholar
 D. C. Montgomery, A Modern Introduction, John Wiley & Sons, New York, NY, USA, 7th edition, 2013.
 M. B. C. Khoo, H. C. Lee, Z. Wu, C.H. Chen, and P. Castagliola, “A synthetic double sampling control chart for the process mean,” IIE Transactions, vol. 43, no. 1, pp. 23–38, 2010. View at: Publisher Site  Google Scholar
 N. A. Saleh and M. A. Mahmoud, “Accounting for phase I sampling variability in the performance of the {MEWMA} control chart with estimated parameters,” Communications in Statistics—Simulation and Computation, vol. 46, no. 6, pp. 4333–4347, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 W. H. Woodall and D. C. Montgomery, “Some current directions in the theory and application of statistical process monitoring,” Journal of Quality Technology, vol. 46, no. 1, pp. 78–94, 2014. View at: Publisher Site  Google Scholar
 S. Psarakis, A. K. Vyniou, and P. Castagliola, “Some recent developments on the effects of parameter estimation on control charts,” Quality and Reliability Engineering International, vol. 30, no. 8, pp. 1113–1129, 2014. View at: Publisher Site  Google Scholar
 P. Castagliola, P. E. Maravelakis, and F. O. Figueiredo, “The EWMA median chart with estimated parameters,” Institute of Industrial Engineers (IIE). IIE Transactions, vol. 48, no. 1, pp. 66–74, 2016. View at: Publisher Site  Google Scholar
 S. Du, X. Yao, and D. Huang, “Engineering modelbased Bayesian monitoring of rampup phase of multistage manufacturing process,” International Journal of Production Research, vol. 53, no. 15, pp. 4594–4613, 2015. View at: Publisher Site  Google Scholar
 S. Chakraborti, “Run length distribution and percentiles: The shewhart X̄ chart with unknown parameters,” Quality Engineering, vol. 19, no. 2, pp. 119–127, 2007. View at: Publisher Site  Google Scholar
 H. W. You, M. B. C. Khoo, M. H. Lee, and P. Castagliola, “Synthetic double sampling X chart with estimated process parameters,” Quality Technology and Quantitative Management, vol. 12, no. 4, pp. 579–604, 2016. View at: Publisher Site  Google Scholar
 Y. Ou, Z. Wu, K. M. Lee, and S. Chen, “An optimal design algorithm of the SPRT chart for minimizing weighted ATS,” International Journal of Production Economics, vol. 139, no. 2, pp. 564–574, 2012. View at: Publisher Site  Google Scholar
 W. L. Teoh, J. K. Chong, M. B. C. Khoo, P. Castagliola, and W. C. Yeong, “Optimal Designs of the Variable Sample Size XfiChart Based on Median Run Length and Expected Median Run Length,” Quality and Reliability Engineering International, vol. 33, no. 1, pp. 121–134, 2017. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Huay Woon You. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.