Abstract

The basic assumption of the parametric control charts is that the underlying process follows the specific distribution. The appropriateness of parametric control charts is questionable when different industrial applications do not support this desired assumption. Recently, nonparametric control charts have been introduced to overcome this deficiency of parametric control charts. Nonparametric control charts have the same in-control run-length characteristics for all continuous distributions and are in-control robust. On the other hand, the ranked set sampling technique is preferred over the simple random sampling technique with control charts because it reduces variability and improves the control chart’s performance. So, this study aims to propose a nonparametric triple homogenously weighted moving average control chart under Wilcoxon signed-rank test with a ranked set sampling technique (regarded as ) to further enhance the process location monitoring. Monte Carlo simulations are used for computational purposes. The proposed control chart’s run-length performance is compared with competing control charts, like TEWMA-SR, TEWMA-, TEWMA-SN, TEWMA-, and control charts. The comparison revealed that the proposed control chart outperformed the other competing control charts, particularly for small to moderate shifts in process location. Finally, a real-life application is presented for quality practitioners to illustrate the effectiveness of the proposed control chart.

1. Introduction

The main concern in the manufacturing process is the product’s quality, and the variations are the primary source that influences the quality of the product. The variations are generally categorized into the natural cause of variations and unnatural causes of variations. Natural variations are harmless while unnatural variations seriously impact the quality of the ultimate product. The process remains in-control (IC) with natural variations, whereas the process is regarded as an out-of-control (OOC) process in the presence of unnatural variations. The existence of unnatural variations causes the shift in process parameters (location and/or dispersion). Control charts are used to monitor these shifts in the process parameters. Control charts are typically divided into memoryless and memory-type control charts based on their structural system. Shewhart and Van Nostrand [1] presented the first memoryless control chart, regarded as the Shewhart control chart, whereas Page [2] and Roberts [3] introduced the concept of memory-type control charts, like cumulative sum (CUSUM) and exponentially weighted moving average (EWMA) control charts, respectively.

Several extensions to the main structure of classical control charts have already been made in order to improve their ability to detect earlier shifts. For example, Alevizakos et al. [4] offered a triple EWMA (TEWMA) control chart that outperformed the EWMA and DEWMA control charts. Similarly, Chatterjee et al. [5] investigated the use of a new TEWMA control chart for detecting shifts in process dispersion. Likewise, Alevizakos et al. [6] and Alevizakos et al. [7] introduced the TEWMA-SR and the TEWMA-SN control charts, respectively, to enhance their shift detection ability in the process location. Also, Mahmood et al. [8] suggested a control chart for joint monitoring of process location and dispersion. Later, Alevizakos et al. [9] presented one and two-sided TEWMA control chart, observing the time between events using the gamma distribution. After that, Rasheed et al. [10] and Rasheed et al. [11] presented two control charts for efficiently monitoring the process mean.

Conventional parametric control charts are often used when an ongoing process follows a specified probability distribution. In many cases, the ongoing process may deviate from a particular distribution, or the distribution of the current process may be in question. Nonparametric (NP) control charts are a good alternative to parametric control charts. The IC run-length (RL) features of the NP control charts remain the same for all symmetric continuous distributions. The sign (SN), as well as Wilcoxon signed-rank (SR), are quite well NP approaches used in statistical process monitoring (SPM). Likewise, simple random sampling (SRS) and ranked set sampling (RSS) techniques are quite often used in SPM with both parametric and NP control charts to observe the ongoing process data. The SPM literature supports RSS because it reduces variability and enhances the efficiency of the associated control charts (Haq et al. [12] and Noor-ul-Amin and Tayyab [13]). Many researchers, including Tsai et al. [14], Abid et al. [15], Abbasi [16], Chakraborti and Graham [17], Rasheed et al. [18], and Abbasi et al. [19], also used these sampling techniques along with NP control charts.

Hunter [20] pointed out that the EWMA control chart gives more weight to recent observations and less weight to former observations. Abbas [21] introduced the homogeneously weighted moving average (HWMA) control chart to monitor shifts in process location effectively. It gives a specific weight to current observations and the remaining weight to all previous observations equally. The distribution of weights in such a way enhances the HWMA control chart’s capacity to detect early shifts. Similarly, Adeoti and Koleoso [22] developed a hybrid HWMA (HHWMA) control chart to monitor the process location’s shifts. Furthermore, Anwar et al. [23] suggested a double HWMA (DHWMA) control chart for observing shifts in process location. They reported that the DHWMA control chart detects shifts better than the HWMA control chart. Also, Riaz et al. [24] designed a triple HWMA (THWMA) control chart to identify the small shifts in the process location. The THWMA control chart outperforms the HWMA and DHWMA control charts for monitoring process location.

Takahasi and Wakimoto [25] suggested that an RSS-based mean estimator is more unbiased and efficient than an SRS-based estimator. Similarly, the THWMA control chart is more effective for monitoring process location parameters. A thorough literature review observed that no one had explored an NP THWMA control chart based on RSS to date. To fill this gap in this study, we are aimed to propose an NP THWMA-SR control chart under RSS () for monitoring shifts in process location for continuous and symmetric distribution. The proposed control chart’s characteristics, including average RL (), median RL (), and standard deviation of RL (), are measured using various distributions like normal, student’s t, contaminated normal (CN), Laplace, and logistic distributions. These characteristics of the proposed control chart is compared to the other competing control charts, like TEWMA-SR, TEWMA-, TEWMA-SN, TEWMA-, and control charts.

The remainder of the paper is ordered as follows: Section 2 offers the design structure of the competing and the proposed control chart. Section 3 defines the proposed control chart’s IC and out-of-control (OOC) performance. Section 4 covers a comparative study of the proposed control chart, whereas Section 5 provides a real-life application. Section 6 concludes with closing remarks.

2. Competing and Proposed Control Charts

This section provides the design structure of the competing and the proposed control chart. More details can be found in the subsequent subsections.

2.1. Control Chart

Various researchers like Kim and Kim [26] and Abid et al. [15] used the following RSS-based Wilcoxon signed-rank statistic to monitor the shift in process location.where symbolizes the process median and denote the number of samples, observations, and cycles used in the RSS approach, respectively. The mean and variance of statistic are and , respectively, where r denotes the number of replications and can be defined as . In this study, we assume for a valid assessment of the proposed control chart to the competitors, so , which is . The quantity is used to improve the efficiency of the control chart and can be defined as . The values of can be obtained by solving the following mathematical expression . The plotting statistic of the NP HWMA under RSS () control chart can be defined as follows:where is a smoothing parameter and meets the condition . Also, the term is the mean of of samples. The mean and variance of the statistic are andrespectively. The control limits of the control chart are as follows:

If or , the process is considered OOC; otherwise, it remains IC.

2.2. Control Chart

The plotting statistic of the control chart is defined as follows:

The simplified version of the plotting statistic of the proposed control chart is expressed as follows:

The is 0, whereas the IC mean and variance of the plotting statistic areand

The control limits of the proposed control chart are defined as follows:

Plot against and . If or the underlying process is OOC; else, it is IC.

2.3. Proposed Control Chart

This subsection describes the design structure of the proposed control chart to monitor shifts in process location. The plotting statistic of the control chart based on RSS is defined as follows:

The simplified version of (10) is

The mean and variance of the statistic are and

Based on this mean and variance, the control limits for the control chart are as follows:

When , the process is in IC state; otherwise, it is OOC.

3. Implementation of the Proposed Control Chart

The performance metrics of the proposed control chart for monitoring shifts in process location is given in this section. Subsection 3.1 offers the performance metrics of the proposed control chart. Similarly, the proposed control chart’s robustness, IC, and OOC performance are presented in Subsection 3.2.

3.1. Performance Metrics

The is commonly used to evaluate the control chart’s performance. The is the expected number of sample points before the first OOC signal from the control chart. In this study, we set , with sample sizes and 10. Monte Carlo simulations with 50,000 simulations in the R program are used to determine the characteristics. To examine the performance behaviour of the proposed control chart, various values of (0.10, 0.25, 0.50, and 0.75) and (0.025, 0.05, 0.075, 0.10, 0.25, 0.50, 0.75, 1.00, 1.50, 2.00, 2.50, 3.00, and 5.00) are used. The following steps are used to obtain the nominal (prespecified) value:(i)To select a random sample from the evaluated distribution.(ii)Specify the process parameters ().(iii)Determine the plotting statistics using equation (8).(iv)Find and from equation (9).(v)Plot the plotting statistic against and over .(vi)If or note this sample of statistic as an RL. For instance, at , if or record 195 as a first RL.(vii)Repeat steps (ii) to (vi) 50,000 times and record RLs.(viii)Using 50,000 simulations, calculate and record RLs.(ix)If ; otherwise, adjust the constant in step (ii) as needed, and then, repeat steps (ii) to (ix) to achieve .(x)To acquire values, draw a shifted sample from the considered distribution again and repeat steps (ii) to (ix).

3.2. Robustness, IC, and OOC Performances of the Control Chart

This subsection highlights the proposed control chart’s robustness, IC, and OOC behaviour when a process location is shifted. Normal and nonnormal continuous and symmetrical distributions are used to evaluate such properties. For this study, the distributions used were the standard normal distribution, i.e., , Laplace or double exponential, i.e., , student’s t distribution, i.e., , logistic distribution i.e., , contaminated normal (CN) distribution, which is the mixture of and (see Table 1). Tables 2 and 3 demonstrate the RL characteristics of the proposed control chart for location shift. For comparison purposes, all these distributions were reparametrized with zero mean/median and unit variance. For all symmetric continuous distributions, the IC RL characteristics of the NP control chart remain constant.

The same parameters are used for comparison purposes as reported in numerous relevant articles. The measures are used to compare the proposed and competing control charts. Table 2 and Figures 1 and 2 provide the following outcomes:(i)The proposed control chart’s IC distribution looks remarkably similar to all distributions examined in this study. For example, at , and the for all investigated distributions (see Tables 2 and 3).(ii)As the smoothing parameter reduces, the proposed control chart becomes more effective in detecting shifts. This illustrates that the proposed control chart is more sensitive to small smoothing parameters (see Figure 1).(iii)As the sample size increases, the proposed control chart’s performance increases.(iv)The Laplace distribution outperforms the other distributions regarding OOC RL performance (see Figure 2).(v)The proposed control chart’s values increase as increases at a certain size of the shift. For instance, under normal distribution at , the , whereas when , the (see Tables 2 and 3).(vi)The distribution of RL values is positively skewed, i.e., (see Tables 2 and 3).(vii)The values of the proposed control chart are smaller than competing control charts with different shift sizes in process location (see Figure 3).

4. Comparative Study

This section compares the proposed control chart with other competing control charts, including TEWMA-SR, TEWMA-, TEWMA-SN, TEWMA-, and .

4.1. Proposed versus TEWMA-SR Control Chart

The profile demonstrates that the proposed control chart outperforms the TEWMA-SR control chart. For example, under normal distribution, at , and , the values of the proposed control chart are , whereas the values of the TEWMA-SR control charts are (see Tables 2 and 4). Similarly, when we consider Laplace distribution for comparison, we observed the same behaviour of the proposed control chart. For instance, at , and the values of the proposed control chart are while the values of the TEWMA-SR control chart are (see Tables 2 and 4). Figure 4 depicts the efficacy of the proposed control chart over the TEWMA-SR control chart. These findings demonstrate that the proposed control chart performs better than the TEWMA-SR control chart for detecting shifts in process location.

4.2. Proposed versus TEWMA-SN Control Chart

The proposed control chart performs better than the TEWMA-SN control chart. The RL profiles of and TEWMA-SN control charts are shown in Tables 2 and 5 for a nominal of 370. It is preferable to use the control chart with the smaller value in a given shift because it detects the shift more quickly. The proposed control charts’ efficiency can be observed under the normal distribution. For example, at , the values for the proposed control charts are whereas the values for the TEWMA-SN control charts are (see Tables 2 and 5). Figure 4 shows the supremacy of the proposed control chart over the TEWMA-SN control chart.

4.3. Proposed Versus TEWMA- Control Chart

The proposed control chart’s ARL values are more sensitive at all combinations of and than the TEWMA- control chart. For instance, in case of logistic distribution, when , and 0.05, 0.10, 0.25, 0.50, and 0.75, the values of the proposed and TEWMA- control charts are 22.05, 10.14, 3.69, 1.81, and 1.26 and 182.82, 76.01, 18.40, 5.95, and 2.94, respectively (see Tables 2 and 6). The supremacy of the proposed control chart to the TEWMA- can also be seen in Figure 4. When compared to the normal distribution, the and TEWMA- control charts reduce 82.12% and 8.80%, respectively, when and , (see Tables 2 and 6). In all investigated distributions, the proposed control chart is more responsive to identifying shifts than the TEWMA- control chart.

4.4. Proposed versus TEWMA- Control Chart

The study reveals that the proposed control chart outperforms the TEWMA- control chart for all possible values of (see Tables 2 and 7). As an illustration, for CN distribution, at , the values of the and TEWMA- control charts are and , respectively (see Tables 2 and 7). Under normal distribution, at , a 2.5% increase in process location parameter reduces the by 23.31% for the TEWMA- control chart while the reduces by 56.31% (see Tables 2 and 7). Figure 4 also shows the superiority of the control chart over the TEWMA- control chart. The results show that the control chart is superior to the TEWMA- control chart in process location’s shift detection.

4.5. Proposed versus Control Chart

The control chart outperforms the control chart in detecting shifts in process location. For example, using t(8) distribution at , the values of the proposed and control charts are and , respectively (see Tables 2 and 8). Similarly, under Laplace distribution, when , the values of the proposed control chart are while the values for control chart are , respectively (see Tables 2 and 8). As illustrated in Figure 4, the control chart outperforms the control chart.

5. Real-Life Application

This section illustrates a real-life application of the nonisothermal continuous stirred tank reactor (CSTR) process to demonstrate the applicability of the proposed control chart. Adegoke et al. [27] recently considered this real-life data. The CSTR process has nine different variables, one of which we choose as the variable of interest (), which represents the output temperature, and this variable is used in real-life application with parameters and . It started with 1000 observations, 600 of which were made while the process was in an IC state. After the 24th observation, a shift in the process location is introduced following Anwar et al. [28] and Anwar et al. [29]. The estimation is carried out with the help of the mentioned parameters, and the control limits are obtained. In order to use the proposed control chart and the existing (TEWMA-) control chart in practise, the variable of interest () is used. We used the RSS approach to generate 40 paired observations of size  = 5 and from a normal distribution. The parameters of the proposed and TEWMA- control charts used for real-life analysis are and , respectively. Figure 5 indicates that the proposed control chart triggers the first OOC signal at sample number 29 while the TEWMA- control chart detects the first OOC point at sample number 35. Similarly, the proposed control chart detects overall 12 OOC points, whereas the TEWMA- control chart detects 06 OOC points (see Figure 5).

6. Conclusions and Recommendations

The basic assumption of the parametric control charts is that the underlying process follows the specific distribution. When the features do not follow a particular distribution, using nonparametric (NP) control charts to monitor the quality characteristics is recommended. The triple homogeneously weighted moving average (THWMA) is the advanced version of the triple exponentially weighted moving average (TEWMA) control chart for enhanced process location monitoring. Similarly, the ranked set sampling (RSS) technique is preferred over simple random sampling (SRS) when the sampling procedure is costly and complicated. To this end, this study uses the NP THWMA control chart and the RSS technique, which presents an NP THWMA Wilcoxon signed-rank control chart (denoted by ) for enhanced monitoring of process location shifts. The proposed control chart’s performance is evaluated in ARL, MDRL, and SDRL. The findings demonstrated that the proposed control chart outperforms competitor control charts such as TEWMA-SR, TEWMA-, TEWMA-SN, TEWMA-, and . Furthermore, a real-life application is presented to illustrate the proposed control chart’s application in practice. The proposed study can be extended to the more ranked set sampling techniques and multivariate scenarios.where .

Data Availability

The [nonisothermal continuous stirred tank reactor] data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.