Next Article in Journal
A Single-Averaged Model for the Solar Radiation Pressure Applied to Space Debris Mitigation Using a Solar Sail
Next Article in Special Issue
Two-Sided Beneficial Value-Added Service Investment and Pricing Strategies in Asymmetric/Symmetric Investment Scenarios
Previous Article in Journal
Hyperstructure Theory Applied to BF-Algebras
Previous Article in Special Issue
Daily Semiparametric GARCH Model Estimation Using Intraday High-Frequency Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Enhanced Dispersion Monitoring Structures Based on Modified Successive Sampling: Application to Fertilizer Production Process

1
Department of Economics and Statistics, Dr Hasan Murad School of Management (HSM), University of Management and Technology, Lahore 54770, Pakistan
2
Department of Statistics, Virtual University of Pakistan, Lahore 54000, Pakistan
3
Industrial and Systems Engineering Department, King Fahd University of Petroleum and Minerals, Dhahran 31261, Saudi Arabia
4
Centre for Smart Mobility and Logistics, King Fahd University of Petroleum and Minerals, Dhahran 31261, Saudi Arabia
5
Department of Mathematics, King Fahd University of Petroleum and Minerals, Dhahran 31261, Saudi Arabia
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(5), 1108; https://doi.org/10.3390/sym15051108
Submission received: 13 April 2023 / Revised: 12 May 2023 / Accepted: 16 May 2023 / Published: 18 May 2023
(This article belongs to the Special Issue Mathematical Models and Methods in Various Sciences)

Abstract

:
In this era of Industry 4.0, efficient and affordable monitoring solutions are needed for the surveillance of manufacturing/service operations. In general, memory-type control charts outperform memoryless control charts when it comes to determining the changes in location and dispersion parameters of symmetrically distributed processes. Before monitoring the process location, it is essential to monitor the process dispersion, since the latter presumes that the process variance remains stable. In practice, the modified successive sampling (MSS) mechanism is preferred over simple random sampling for its cost-effectiveness and efficiency. This study was designed in order to propose moving average and double moving average control charts based on the MSS mechanism for monitoring the dispersion parameter. The performance of the proposed charts is evaluated using run-length measures, and a comparison is made with an existing control chart based on MSS and repetitive sampling. Furthermore, the application of the designed moving and double moving average charts is demonstrated using a case study related to fertilizer production. It is observed that the proposed double moving average control chart performs better than the other control charts designed under the MSS and repetitive sampling schemes.

1. Introduction

In statistical process control (SPC), a control chart is the most important and frequently used tool to monitor process parameters (such as location and/or dispersion). Quality practitioners mostly prefer using control charts to identify sustainable variations in the process parameters. In a process, the main application of control charts is the visual detection of unusual variations for which educative action is needed to move the process back into the in-control (IC) state [1]. Shewhart [2] invented the first control chart, named the Shewhart chart, which was designed based on a current sample. Therefore, Shewhart-type (memoryless) control charts are more efficient for detecting large shifts in the process and are less efficient for detecting shifts of small magnitude. On the other hand, memory-type control charts, such as cumulative sum (CUSUM) by Page [3], exponentially weighted moving average (EWMA) by Roberts [4], moving average (MA) by Wong et al. [5] and double moving average (DMA) by [6], perform more efficiently for the timely detection of small shifts in the process. The reason for this is that their structures are based on current as well as past observations of the process [7].
Many quality practitioners have studied MA and DMA control charts in detail. MA and combined MA-Shewhart design procedures were designed by Wong et al. [5]. To monitor both process parameters (mean and variance), Khoo and Yap [8] designed an MA chart for the joint monitoring of both upwards and downwards shift. For the monitoring of process variability, Adeoti et al. [9] designed a double moving average-S chart that detects small to moderate shifts. Khoo and Wong [10] presented a double exponentially weighted moving average (DEWMA) chart to monitor the location parameter of a process. For monitoring defective products, Areepong [11] explicated formulas for the average run length (ARL) of an MA chart, and Phant et al. [12] gave an expression for the ARL of a DMA chart under the integer-valued autoregression of order one (INAR(1)) processes. Adeoti et al. [9] and Alevizakos et al. [6] studied a double moving average (DMA) chart to enhance the detection ability of a moving average (MA) chart under normal distribution and assuming a simple random sampling (SRS) mechanism.
In the literature, many studies have presented more accurate and efficient estimators, some of which use sampling schemes other than SRS. Their findings reveal that such sampling mechanisms are more reliable, time- and cost-effective as compared to SRS. Salazar and Sinha [13] designed memoryless control charts using ranked set sampling (RSS). Al-Nasser and Al-Rawwash [14] developed Shewhart control charts based on a robust RSS mechanism. Abujiya et al. [15] suggested an EWMA structure based on median RSS, while Munir and Haq [16] presented a CUSUM structure using an ordered and double ordered RSS scheme. For location monitoring, Nawaz et al. [17], Nawaz and Han [18] and Hussain et al. [19] studied different median- and mean-based Shewhart, CUSUM, EWMA and homogeneously weighted moving average (HWMA) charts under ranked set sampling (RSS) and neoteric ranked set sampling (NRSS) schemes. Reynolds Jr and Lou [20] evaluated a GLR control chart for the monitoring of process mean. The generalized control chart is based on the likelihood ratio of normal distribution under SRS. The average time to signal measure is used to evaluate the performance of the proposed charts. A comparison with the Shewhart, CUSUM and the combined Shewhart-CUSUM charts is provided. Furthermore, Sheriff et al. [21] performed a comparative study on PCA-based GLR control charts for process monitoring. Their method is based on the likelihood ratio of principal component analysis for high-dimensional data. The proposed chart is used to monitor location and/or dispersion shifts. Riaz et al. [22] designed a new HWMA control chart for the monitoring of the dispersion parameter using wind farm data. Anwar et al. [23] and Anwar et al. [24] proposed mixed memory-type control charts for the simultaneous monitoring of more than one parameters using auxiliary information. A double generally weighted moving average (GWMA) control chart for monitoring dispersion parameters was designed by Alevizakos et al. [25]. Akhtar et al. [26] evaluated the EWMA control chart by using log-normal distributions with estimated parameters to monitor process variability. A GWMA maximum chart for joint monitoring was designed by Chatterjee et al. [27]. Similarly, Chatterjee et al. [28] introduced an efficient control charting structure named double GWMA (DGWMA) for the monitoring of location and dispersion parameters. Khan et al. [29] developed a hybrid EWMA control chart based on a ranked set sampling scheme under a Bayesian approach using hard-bake process data. For the monitoring of variability in a process, an exponentially weighted moving average control chart based on generalized fast initial response was designed by Ajibade et al. [30]. Recently, a new weighted adaptive CUSUM chart using a bivariate normal process to monitor generalized variance was proposed by Haq and Abbasi [31].
All these aforementioned studies under the Shewhart setup and their analyses are conducted under the assumption that all n sample observations are available and no information from the previous sample(s) is used for making decisions. The CUSUM and EWMA structures overcome this issue partially by defining the plotting statistics as functions of current and previous samples. There are many practical situations in which it is desirable to conduct control charting analysis with a sample size n but only n c observations are available at each time point, where c is a positive integer. For such situations, a cost-effective and time-saving sampling procedure named modified successive sampling (MSS) was first introduced by Yaqub et al. [32]. Yaqub et al. [32] and Abbas et al. [33] proposed Shewhart charts based on MSSS schemes for the detection of location and dispersion parameters, respectively. The modified form of the successive sampling scheme is more effective than other existing sampling techniques and is specially designed to monitor process parameters (location or dispersion). For the real-time surveillance of the location parameter, memory-type control charts based on MSS were designed by Hyder et al. [34].
The situation of a reduced sample size can also occur in the case of dispersion monitoring. To the best of our knowledge, there are no studies in the literature that provide a solution for such situations under memory-type structures. Filling this research gap, the current study is focused on designing efficient memory-type control charts (MA- S M S S ( S ) 2 , DMA- S M S S ( S ) 2 ) based on the MSS scheme for the monitoring of symmetrical process variability. The performance of the designed control charts is compared with the existing Shewhart chart for process dispersion monitoring under various MSS schemes.
The rest of this article is organized as follows. In Section 2, a summary of MSS schemes is given, along with the construction of the designed charts based on these schemes. Section 3 presents the performance of proposed charts, which is evaluated through simulations. A comparative analysis is presented in Section 4, Section 5 illustrates the implementation of the designed charts using real-life data and the concluding remarks of the article and future recommendations are given in Section 6.

2. Methodology

This section briefly defines the structure and several schemes of MSS. Additionally, for the efficient monitoring framework of the scale parameter, we construct MA- S M S S ( S ) 2 and DMA- S M S S ( S ) 2 charts based on MSS in this section.

2.1. Modified Successive Sampling

The timely inspection of defective items may improve the cost-efficiency of any industrial process. To improve cost-efficiency, Jessen [35] developed the successive sampling technique for various inventory problems. However, in most surveys, SRS is suggested for a single occasion. However, an item’s quality is regularly examined in industrial practices. Therefore, to obtain reliable estimates, successive sampling plays an important role in such repeated assessments. The structure of successive sampling is that the initial sample is selected on the 1st occasion and the next sample (having some sample points from the previous sample) is drawn on the 2nd occasion, and so on. This sampling procedure will be carried out until desired sample size is achieved. Patterson [36], Rao and Graham [37], Choudhary et al. [38], Yaqub et al. [32], Abbas et al. [33] and Hyder et al. [34] described some modifications of successive sampling.
Abbas et al. [33] discussed the structure of modified successive sampling (MSS) for memoryless control charts. The procedure of MSS is outlined as follows:
  • Firstly, select a sample  y 1 , 1 , y 1 , 2 , , y 1 , n of size n from a symmetric distribution by using simple random sampling (SRS).
  • New observations  n c are drawn by using SRS for the 2nd sample y 2 , 1 , y 2 , 2 , , y 2 , n c and the rest of the c observations are some selected quantiles of the previous sample as [ y 2 , 1 , y 2 , 2 , , y 2 , n c , y 2 , n c + 1 , y 2 , n c + 2 , , y 2 , n ], where y 2 , n c + 1 = f 1 y 1 , 1 , y 1 , 2 , , y 1 , n , y 2 , n c + 2 = f 2 y 1 , 1 , y 1 , 2 , , y 1 , n and so on, up to y 2 , n = f c y 1 , 1 , y 1 , 2 , . , y 1 , n .
  • Using the SRS scheme again, take  n c new observations as the 3rd sample and the leftover c observations are selected by using some quantile points of the 2nd sample. Similarly, this procedure will continue for the complete production process run.
The M S S is symbolized as M S S n , c , f 1 , f 2 , f c . , where n and c denote the size of the sample and the number of values chosen from the preceding sample and f p is a function of the preceding sample   p = 1 , 2 , , c , respectively. For the monitoring of the dispersion parameter, following Abbas et al. [33], we have considered ten different MSS schemes, as given below:
(a)
For c = 2 , the n 2 new observations are selected by using the SRS scheme, and the other two observations are specific quantile points ( Q 1 , Q 2 ) of the previous sample. Notationally, this is defined as M S S n , 2 , Q 1 , Q 2 . In this study, different choices of quantile pairs such as ( Q 0.25 , Q 0.75 ) , Q 0.30 , Q 0.70 , Q 0.35 , Q 0.65 , Q 0.40 , Q 0.60 and ( Q 0.45 , Q 0.55 )   a r e   u s e d . By using these quantile pairs, the various M S S ( S ) schemes are considered, which are described below:
  • M S S n , 2 , Q 0.25 , Q 0.75 , is the first M S S ( 1 ) scheme, where S = 1 .
  • M S S n , 2 , Q 0.30 , Q 0.70 , is the second M S S ( 2 ) scheme, where S = 2 .
  • M S S n , 2 , Q 0.35 , Q 0.65 , is the third M S S ( 3 ) scheme, where S = 3 .
  • M S S n , 2 , Q 0.40 , Q 0.60 , is the fourth M S S ( 4 ) scheme, where S = 4 .
  • M S S n , 2 , Q 0.45 , Q 0.55 , is the fifth M S S ( 5 ) scheme, where S = 5 .
(b)
When  c = 3 , n 3 new observations are drawn by using the SRS scheme, and the remaining three observations are specific quantiles points ( Q 1 , Q 2 , Q 3 ) of the previous sample. Symbolically, This is defined as M S S n , 3 , Q 1 , Q 2 , Q 3 . In this study, different choices of quantile points such as ( Q 0.25 , Q 0.50 , Q 0.75 ) , Q 0.30 , Q 0.50 , Q 0.70 , Q 0.35 , Q 0.50 , Q 0.65 , Q 0.40 , Q 0.50 , Q 0.60 and ( Q 0.45 , Q 0.50 , Q 0.55 )   a r e   u s e d . By using these quantile pairs, the other M S S ( S ) schemes are described as follows:
  • M S S n , 3 , Q 0.25 , Q 0.50 , Q 0.75 , is the sixth M S S ( 6 ) scheme, where S = 6 .
  • M S S n , 3 , Q 0.30 , Q 0.50 , Q 0.70 , is the seventh M S S ( 7 ) scheme, where S = 7 .
  • M S S n , 3 , Q 0.35 , Q 0.50 , Q 0.65 , is the eighth M S S ( 8 ) scheme, where S = 8 .
  • M S S n , 3 , Q 0.40 , Q 0.50 , Q 0.60 , is the ninth M S S ( 9 ) scheme, where S = 9 .
  • M S S n , 3 , Q 0.45 , Q 0.50 , Q 0.55 , is the tenth M S S ( 10 ) scheme, where S = 10 .
It should be noted that the subscript (S) in M S S ( S ) describes the different pairing schemes of modified successive sampling (MSS).

2.2. Existing Shewhart- S M S S ( S ) 2 Control Chart

For the detection of large shifts in the dispersion parameter, Abbas et al. [33] introduced a Shewhart-S2 chart based on modified successive sampling. The plotting statistic of the Shewhart- S M S S ( S ) 2 control chart is as follows:
S i 2 = j = 1 n ( X i j X ) 2 n 1
The lower and upper control limits of the Shewhart- S M S S ( S ) 2 chart are defined as follows:
L C L M S S = S s 2 L M S S M S E S 2
U C L M S S = S s 2 + L M S S M S E S 2
where S s 2 is the mean of S 2 , M S E S 2 is the mean square error of S 2 and L M S S is the charting coefficient of the existing Shewhart chart based on MSS.

2.3. Proposed MA- S M S S ( S ) 2 Control Chart

Initially, the moving average (MA) control chart structure was proposed by Wong et al. [5]. For the monitoring of the dispersion parameter, the statistic of the designed chart (MA- S M S S ( S ) 2 ) under different MSS schemes is described as follows:
M A i ( S ) = j = 1 i S j ( S ) 2 i , f o r   i < w , j = i w + 1 i S j ( S ) 2 w , f o r   i w .
where S j ( S ) 2 indicates the variance of the j t h sample with the specific MSS schemes shown by the subscript (S), i is the sampling time and w denotes the span of the moving average. The mean of the M A i S statistic is E M A S = S s 2 , and the variance of M A i S is given as
V a r M A ( S ) = M S E S 2 n i , f o r   i < w , M S E S 2 n w , f o r   i w
The control limits of the MA- S M S S ( S ) 2 chart are calculated as follows:
L C L M A ( S ) = 0
U C L M A ( S ) = S s 2 + L M S S ( S ) V a r ( M A S )
where L M S S S is a control charting coefficient observed against the specified values of A R L 0 and sample size ( n ) . If MA- S M S S ( S ) 2 statistic is plotted outside the control limits ( L C L M S S S / U C L M S S S ) , then the process is declared to be out of control (OOC).

2.4. Proposed DMA- S M S S ( S ) 2 Control Chart

Alevizakos et al. [6] proposed a memory-type double moving average (DMA) control chart for efficient monitoring. The statistic of DMA- S M S S ( S ) 2 under MSS schemes is defined as follows:
D M A i ( S ) = j = 1 i M A j ( S ) i , f o r   i < w , j = i w + 1 i M A j ( S ) w , f o r   i w .
The expected value of the D M A ( S ) statistic is obtained as
E D M A S = S s 2
and the variance of the D M A ( S ) statistic is evaluated as follows:
for i < w ,
V a r D M A ( S ) = j = 1 i a j 2 M S E S 2 n i 2 ,
for w i < 2 w 1 ,
V a r D M A S = M S E S 2 n w 2 j 1 = i w + 1 w 1 1 j 1 + i w + 1 j 11 < j 12 w 1 2 j 12 + j 1 = i w + 1 w 1 j 2 = w i 2 j 1 j 2 + w j 1 w + i w + 1 w + w j 21 < j 22 i 2 j 21 j 22 + w w 2 .
and for i 2 w 1 ,
V a r D M A ( S ) = M S E S 2 n w 2 1 + i w + 1 j 1 < j 2 i 2 j 1 j 2 + w w 2 .
The mean and variances of the DMA statistic are derived by Alevizakos et al. [6]. The control limits of the DMA- S M S S ( S ) 2 control chart are defined as follows:
L C L D M A S = 0
U C L D M A ( S ) = S s 2 + K M S S ( S ) V a r ( D M A S )
The charting constant is denoted as K M S S ( S ) . If the statistic D M A i ( S ) lies between these (lower and upper) control limits, the process is called IC; otherwise, it is OOC.

3. Performance Evaluations of Designed Charts

In this section, performance evaluations and comparative analysis of the designed charts with the existing Shewhart control charts proposed by Abbas et al. [33] have been conducted. All the comparisons are carried out under the assumption of known parameters. However, the effect of parameter estimation on some moving average control charts is studied in detail by Jones et al. [39] and Noorossana et al. [40]. The performance of the designed MA- S M S S ( S ) 2 and DMA- S M S S ( S ) 2 charts are evaluated by using run length (RL) metrics. The RL metrics include various performance measuring criteria such as average RL (ARL), standard deviation RL (SDRL) and median RL (MDRL). The average RL (ARL) is the average of the number of plotting statistics before detecting an OOC signal [41,42]. A R L 0 is defined as the ARL of the IC process, while A R L 1 is the ARL of an OOC process. In practice, the control limits of their corresponding charts are set against the pre-specified value of A R L 0 , and then the performance is evaluated based on the A R L 1 values [43,44]. It should be noted here that there are several short production runs with low-volume manufacturing environments where the quality of characteristics is monitored [45]. The shifts in such processes, if any, occur at or close to the start. The current study is targeted towards such processes where early change point detection is desired for shifts that occur close to the beginning of the process. Accordingly, the algorithm of zero-state ARL computation is used where the possible change point is at i = 1 .
Assuming normal distribution (symmetrical distribution) with μ = 0 , σ = 1 , the average run length of an IC process is pre-fixed at 370, i.e., A R L 0 = 370 , the moving span is set at w = 2 and δ = 1 (this means that the process variability has no shift). For a fixed A R L 0 = 370 , the values of charting coefficients ( L M S S ( S ) and K M S S ( S ) ) of the designed MA- S M S S ( S ) 2 and DMA- S M S S ( S ) 2 charts, S s 2 and M S E S 2 under different schemes of MSS are presented in Table 1 for the sample size choices n = 5   a n d   7 .
The control chart with the fewest values of A R L 1 is said to be more efficient and have higher detection ability than the others. For this instance, the complete profiles of RL with different shifts δ in process variability, sample sizes n and moving span w are presented under various MSS schemes in Table 2 and Table 3.

Algorithm for Control Charting Constants

In this section, under various MSS schemes, the following basic procedure has been carried out to find the suitable values of the charting coefficients ( L M S S ( S ) and K M S S ( S ) ) of the designed MA- S M S S ( S ) 2 and DMA- S M S S ( S ) 2 charts, respectively.
Step 1: The first step of this procedure is to select the values of the parameters ( S s 2 , M S E S 2 , w , n ) for the designed chart with pre-fixed A R L 0 = 370 . This article also considers the various choices of the above-mentioned design parameters (i) sample size ( n = 5   a n d   7 ) and (ii) moving average span w = 2 , 3 , 5 , 10   a n d   15 for both the proposed MA- S M S S ( S ) 2 and DMA- S M S S ( S ) 2 charts.
Step 2: A normally distributed random sample of size n (symmetrical sample) is generated using μ = 0 , σ = δ × 1 , and the remaining samples with size n will be generated through the MSS technique. It should be noted that δ refers to the size of the shift; hence, to calculate A R L 0 , it is considered to be equal to 1; otherwise, for the A R L 1 study, it is set to be more than 1.
Step 3: The variance given in Equation (1) is computed for each sample, which is further used to compute the MA- S M S S ( S ) 2 and DMA- S M S S ( S ) 2 plotting statistics using Equations (4) and (8).
Step 4: The arbitrary values of the charting coefficients ( L M S S ( S ) and K M S S ( S ) ) of the designed MA- S M S S ( S ) 2 and DMA- S M S S ( S ) 2 charts are set, and the control limits of the designed charts are computed.
Step 5: The proposed charts’ statistics are plotted against their control limits to obtain a single RL value by the number of samples. Afterwards, the process is declared to be OOC.
The R language software is used to run steps 1–5 100,000 times and obtain the A R L 0 value. If the obtained A R L 0 is not equal to the pre-specified value of 370, then the L M S S ( S ) and K M S S ( S ) values are adjusted and this procedure is repeated to find the desired value of A R L 0 . Once the proper choice of L M S S ( S ) and K M S S ( S ) are obtained, then the shifts in variance and obtained A R L 1 values are added.

4. Comparative Analysis

In practice, if the variation in the quality characteristics is reduced, this leads to an efficient production process. However, if the variation in the quality characteristics increases, the production process deteriorates. Hence, most quality experts are interested in diagnosing the effect of process deterioration. Therefore, this study is also designed to evaluate the performance of the designed charts (MA- S M S S ( S ) 2 and DMA- S M S S ( S ) 2 ) by using incremental shifts ( δ = 1 , 1.1 , 1.2 , 1.3 , 1.4 , 1.5 , 1.6 , 1.7 , 1.8 , 1.9 , 2 ) in the dispersion parameter. The shift is defined as δ = σ 1 / σ 0 for the dispersion parameter, where σ 0   a n d   σ 1 are the standard deviation of IC and OOC processes, respectively. Under different MSS schemes, a comparative analysis has been made between these designed charts along with the existing Shewhart- S M S S ( S ) 2 chart (cf. Table 2 and Table 3). The A R L 1 , M D R L 1 and S D R L 1 are used as performance measures to examine the performances of the above-mentioned control charts.
Under different sampling schemes of MSS at constant c = 2 , the results of the Shewhart- S M S S ( S ) 2 , MA- S M S S ( S ) 2 and DMA- S M S S ( S ) 2 charts are shown in Table 2. The primary findings are listed below:
It is observed that the DMA- S M S S ( S ) 2 chart outperformed the Shewhart- S M S S ( S ) 2 and MA- S M S S ( S ) 2 charts on the same amount of shift in the dispersion parameter. For example, when n = 5 , a 30% increase in the change in the dispersion parameter reduces the A R L 1 values of the Shewhart- S M S S ( 1 ) 2 , MA- S M S S ( 1 ) 2 and DMA- S M S S ( 1 ) 2 charts to 72.18, 65.18 and 64.29, respectively. However, when n = 7 and δ = 1.5 , t h e   A R L 1 values 31.10, 21.91 and 21.64 are reported for the Shewhart- S M S S ( 5 ) 2 , MA- S M S S ( 5 ) 2 and DMA- S M S S ( 5 ) 2 charts, respectively.
The findings also reveal that increasing the sample size leads to better performance of the charts. For example, the DMA- S M S S ( S ) 2 chart outperformed the others; therefore, at δ = 1.4 , t h e   A R L 1 value of the DMA- S M S S ( 2 ) 2 chart is reported as 44.36 when n = 5 , while for the same settings with n = 7 , it is reported as 31.60.
The findings of the Shewhart- S M S S ( S ) 2 , MA- S M S S ( S ) 2 and DMA- S M S S ( S ) 2 charts under the MSS scheme with c = 3 are shown in Table 3. The primary findings are listed below:
It is noted that the DMA- S M S S ( S ) 2 chart outperformed the Shewhart- S M S S ( S ) 2 and MA- S M S S ( S ) 2 charts at a fixed amount of shift in the dispersion parameter. For example, when n = 5 , a 60% increase in the dispersion parameter reduces the A R L 1 values of the Shewhart- S M S S ( 7 ) 2 , MA- S M S S ( 7 ) 2 and DMA- S M S S ( 7 ) 2 charts to 34.21, 22.74 and 17.00, respectively. However, when n = 7 and δ = 1.7 , t h e   A R L 1 values 28.73, 20.29 and 10.00 are reported for the Shewhart- S M S S ( 9 ) 2 , MA- S M S S ( 9 ) 2 and DMA- S M S S ( 9 ) 2 charts, respectively.
Moreover, the results also show that the charts’ performance improved as the sample size increased. For example, the DMA- S M S S ( S ) 2 chart outperformed the others; therefore, at δ = 1.2 , t h e A R L 1 value of the DMA- S M S S ( 8 ) 2 chart is reported as 115.87 when n = 5 , while for the same settings with n = 7 , it is reported as 88.79.
Overall, from Table 2 and Table 3, it is also revealed that based on different MSS schemes, the DMA- S M S S ( S ) 2 chart performed better than the Shewhart- S M S S ( S ) 2 and MA- S M S S ( S ) 2 charts with respect to A R L 1 at a pre-fixed A R L 0 = 370 .
Furthermore, to check the effect of moving average span on the detection ability of the designed MA- S M S S ( S ) 2 and DMA- S M S S ( S ) 2 charts, sensitivity analysis was performed based on various choices of moving average span ( w = 2 , 3 , 5 , 10   a n d   15 ) . The results of the sensitivity analysis for the MA- S M S S ( S ) 2 and DMA- S M S S ( S ) 2 charts with c = 2   a n d   3 are presented in Figure 1, Figure 2, Figure 3 and Figure 4. At the pre-specified A R L 0 = 370 , it is noted that the performance ability of the MA- S M S S ( S ) 2 and DMA- S M S S ( S ) 2 charts increased with the increase in the moving average span. For example, it can be easily observed in Figure 1, Figure 2, Figure 3 and Figure 4 that at w = 15 , both the designed charts have the lowest curves of A R L 1 for all the above-mentioned MSS schemes as compared to all the choices of w for c = 2   a n d   3 .
Furthermore, a comparative study of the best proposed MA- S M S S ( 1 ) 2 and DMA- S M S S ( 1 ) 2 charts with the existing S 2 chart was designed based on the repetitive sampling S r e p e t i t i v e 2 proposed by Aslam et al. [46]. For the pre-fixed A R L 0 = 370 and c = 2 , the run length curves of the MA- S M S S ( 1 ) 2 , DMA- S M S S ( 1 ) 2 and S r e p e t i t i v e 2 are plotted in Figure 5. It is observed that for small amounts of shift (i.e., 1.1 to 1.5), the MA- S M S S ( 1 ) 2 and DMA- S M S S ( 1 ) 2 charts have comparatively better detection ability, while for large shifts (i.e., more than 1.5), the S r e p e t i t i v e 2 chart showed relatively lower ARL values. The dominance of the S r e p e t i t i v e 2 chart in the presence of large shifts is natural because of the Shewhart structure.

5. Case Study

This section contains real-life applications of the designed control charts to monitor process dispersion in a ZA fertilizer production-based dataset. Mashuri et al. [47] also used this dataset to illustrate the better performance of proposed multivariate control charts in terms of monitoring the process variability. ZA fertilizer is frequently used to increase the concentration of sulfur and nitrogen in soil. It is also called ammonium sulfate fertilizer. The production process of ZA fertilizer is based on six stages (cf. Figure 6): carbonation, reaction and gas scrub, filtration, neutralization, evaporation and crystal, and drying and cooling stages. Air and carbon dioxide are emitted in this production process.
In this case study, the concentration of C O 2 ( g / L ) is the variable of interest in the fertilizer production process. This variable has 526 values consisting of the concentration of C O 2 ( g / L ) in the body and in the glassy coating on ceramics called glaze. These observations related to the body are considered in-control (239 values) data, while observations related to glaze are considered out-of-control (287 values) data. To assess the normality of both datasets, probability plots were prepared, which are shown in Figure 7. It should be noted that the in-control data are normally distributed, with a mean of 0.3975 and 0.2115 standard deviation. However, the out-of-control data are not normally distributed; their mean is 0.5459 with 0.424 standard deviation. By adopting the mechanism of MSS, 174 samples of size 5 were drawn from this dataset, which has two subgroups (79 in control and 95 out of control). To evaluate the performance of both the designed charts, the in-control subgroup of this dataset was utilized to calculate the proposed MA- S M S S ( 1 ) 2 and DMA- S M S S ( 1 ) 2 charts’ control limits at the pre-defined A R L 0 = 370 .
Figure 8 presents a graphical view of both the proposed charts. On the X-axis, the subgroups (sample numbers) are plotted, and on the y-axis, the plotting statistics of both the proposed charts are plotted against their respective control limits. The pink shaded area represents in-control points while white colored is referring to the out-of-control region. It shows that the MA- S M S S ( 1 ) 2 control chart detects an out-of-control signal after inspecting 92 samples, while the DMA- S M S S ( 1 ) 2 control chart detects out-of-control points after observing 90 samples. Both control charts have alternatively in-control and out-of-control points in out-of-control regions; therefore, a specific pattern can be seen in both the MA- S M S S ( 1 ) 2 and DMA- S M S S ( 1 ) 2 charts (cf. Figure 8). From the above simulation study, it is noted that the DMA- S M S S ( S ) 2 chart has outperformed the Shewhart- S M S S ( S ) 2 and MA- S M S S ( S ) 2 charts with a fixed amount of shift in the dispersion parameter. Similarly, by implementing the real-life example, we observe that the detection ability of the DMA- S M S S ( 1 ) 2 control chart is high, and it is more sensitive than the MA- S M S S ( 1 ) 2 control chart for the monitoring of process variability.
Furthermore, the proposed charts have wide real-life applications in different major fields of life. For example, we can utilize this control charting structure in the chemical industry, medical industry, glass field, ceramic industry, engineering, fertilizer production, hard-bake process and wind farm data (cf. Riaz et al. [22], Anwar et al. [24] and Khan et al. [29]).

6. Conclusions

Usually, memory-type control charts are used to detect small to moderate shifts in the location or dispersion process parameters. The process is said to be OOC depending on the amount of shift detected in the location or dispersion parameters. Usually, it is preferable to detect a shift in the dispersion parameter of the process before the detection of a shift in the location parameter. In this study, we have designed memory-type MA- S M S S ( S ) 2 and DMA- S M S S ( S ) 2 control charts for the monitoring of the process dispersion parameter using a modified successive sampling technique. Additionally, a comparative analysis with the existing Shewhart- S M S S ( S ) 2 chart has been presented using different performance measures (ARL, MDRL and SDRL). From the results, we have observed that the value of A R L 1 decreases with the increasing value of dispersion shift and sample size at a fixed value of A R L 0 = 370 (cf. Table 2 and Table 3). Moreover, the findings lead to the conclusion that the designed DMA- S M S S ( S ) 2 control chart performs better than Shewhart- S M S S ( S ) 2 and MA- S M S S ( S ) 2 under all the aforementioned MSS schemes. We have also noted that increasing values of the moving average span ( w ) show a declining trend in A R L 1 curves under all MSS schemes (cf. Figure 1, Figure 2, Figure 3 and Figure 4).
Generally, the structure of mixed EWMA and CUSUM charts are more efficient compared to traditional MA, DMA, EWMA and CUSUM charts because they are more sensitive to the detection of small shifts in any process (cf. [48]). Therefore, one could extend this study by using different choices of MSS schemes and sample sizes to design a new mixed structure of the EWMA and CUSUM control charts. Moreover, the current study is designed based on the known parameter case (K-case), while the unknown parameter case (U-case) is recommended for future study. Furthermore, the MSS structure is proposed under normal distribution (symmetric distribution); however, one could propose MSS structure based on asymmetrical distribution and further extend the performance ability of the proposed control charts under the new MSS structure based on asymmetrical distributions.

Author Contributions

Conceptualization, M.H. and T.M.; methodology, M.H., T.M. and S.M.M.R.; software, M.H. and T.M.; validation, M.H., S.M.M.R., T.M. and N.A.; formal analysis, M.H.; writing—original draft preparation, M.H., N.A., T.M. and S.M.M.R.; writing—review and editing, M.H., N.A., T.M. and S.M.M.R.; visualization, M.H.; supervision, T.M. and S.M.M.R.; funding acquisition, T.M. All authors have read and agreed to the published version of the manuscript.

Funding

This publication is based on work supported by King Fahd University of Petroleum and Minerals. Tahir Mahmood and Nasir Abbas at KFUPM acknowledge the Interdisciplinary Research Center for Smart Mobility and Logistics for the support received under grant no. INML2205.

Data Availability Statement

Data sharing is not applicable to this article as no datasets were generated or analyzed during the current study.

Acknowledgments

We are thankful to King Fahd University of Petroleum and Minerals for providing us with research facilities for conducting this research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Arslan, M.; Anwar, S.; Gunaime, N.M.; Shahab, S.; Lone, S.A.; Rasheed, Z. An Improved Charting Scheme to Monitor the Process Mean Using Two Supplementary Variables. Symmetry 2023, 15, 482. [Google Scholar] [CrossRef]
  2. Shewhart, W.A. Quality control charts. Bell Syst. Tech. J. 1926, 5, 593–603. [Google Scholar] [CrossRef]
  3. Page, E.S. Continuous inspection schemes. Biometrika 1954, 41, 100–115. [Google Scholar] [CrossRef]
  4. Roberts, S. Control chart tests based on geometric moving averages. Technometrics 1959, 42, 97–101. [Google Scholar] [CrossRef]
  5. Wong, H.; Gan, F.; Chang, T. Designs of moving average control chart. J. Stat. Comput. Simul. 2004, 74, 47–62. [Google Scholar] [CrossRef]
  6. Alevizakos, V.; Chatterjee, K.; Koukouvinos, C.; Lappa, A. A double moving average control chart: Discussion. Commun. Stat.-Simul. Comput. 2020, 51, 6043–6057. [Google Scholar] [CrossRef]
  7. Hu, X.; Sun, G.; Xie, F.; Tang, A. Monitoring the Ratio of Two Normal Variables Based on Triple Exponentially Weighted Moving Average Control Charts with Fixed and Variable Sampling Intervals. Symmetry 2022, 14, 1236. [Google Scholar] [CrossRef]
  8. Khoo, M.B.; Yap, P. Joint monitoring of process mean and variability with a single moving average control chart. Qual. Eng. 2004, 17, 51–65. [Google Scholar] [CrossRef]
  9. Adeoti, O.A.; Akomolafe, A.A.; Adebola, F.B. Monitoring process variability using double moving average control chart. Ind. Eng. Manag. Syst. 2019, 18, 210–221. [Google Scholar] [CrossRef]
  10. Khoo, M.B.; Wong, V. A double moving average control chart. Commun. Stat.—Simul. Comput.® 2008, 37, 1696–1708. [Google Scholar] [CrossRef]
  11. Areepong, Y. Explicit formulas of average run length for a moving average control chart for monitoring the number of defective products. Int. J. Pure Appl. Math. 2012, 80, 331–343. [Google Scholar]
  12. Phant, S.; Sukparungsee, S.; Areepong, Y. Explicit formula for average run length of double moving control chart for INAR (1) processes. Preprint 2016. [Google Scholar] [CrossRef]
  13. Salazar, R.; Sinha, A. Control chart X based on ranked set sampling. Comun. Tecica 1997, 1, 1–97. [Google Scholar]
  14. Al-Nasser, A.D.; Al-Rawwash, M. A control chart based on ranked data. J. Appl. Sci. 2007, 7, 1936–1941. [Google Scholar]
  15. Abujiya, M.a.R.; Lee, M.H.; Riaz, M. Improving the performance of exponentially weighted moving average control charts. Qual. Reliab. Eng. Int. 2014, 30, 571–590. [Google Scholar] [CrossRef]
  16. Munir, W.; Haq, A. New cumulative sum control charts for monitoring process variability. J. Stat. Comput. Simul. 2017, 87, 2882–2899. [Google Scholar] [CrossRef]
  17. Nawaz, T.; Raza, M.A.; Han, D. A new approach to design efficient univariate control charts to monitor the process mean. Qual. Reliab. Eng. Int. 2018, 34, 1732–1751. [Google Scholar] [CrossRef]
  18. Nawaz, T.; Han, D. Monitoring the process location by using new ranked set sampling-based memory control charts. Qual. Technol. Quant. Manag. 2020, 17, 255–284. [Google Scholar] [CrossRef]
  19. Hussain, S.; Mahmood, T.; Riaz, M.; Nazir, H.Z. A new approach to design median control charts for location monitoring. Commun. Stat.-Simul. Comput. 2022, 51, 3553–3577. [Google Scholar] [CrossRef]
  20. Reynolds, M.R., Jr.; Lou, J. An evaluation of a GLR control chart for monitoring the process mean. J. Qual. Technol. 2010, 42, 287–310. [Google Scholar] [CrossRef]
  21. Sheriff, M.Z.; Karim, M.N.; Nounou, H.N.; Nounou, M.N. Process monitoring using PCA-based GLR methods: A comparative study. J. Comput. Sci. 2018, 27, 227–246. [Google Scholar] [CrossRef]
  22. Riaz, M.; Abbasi, S.A.; Abid, M.; Hamzat, A.K. A new HWMA dispersion control chart with an application to wind farm data. Mathematics 2020, 8, 2136. [Google Scholar] [CrossRef]
  23. Anwar, S.M.; Aslam, M.; Riaz, M.; Zaman, B. On mixed memory control charts based on auxiliary information for efficient process monitoring. Qual. Reliab. Eng. Int. 2020, 36, 1949–1968. [Google Scholar] [CrossRef]
  24. Anwar, S.M.; Aslam, M.; Zaman, B.; Riaz, M. Mixed memory control chart based on auxiliary information for simultaneously monitoring of process parameters: An application in glass field. Comput. Ind. Eng. 2021, 156, 107284. [Google Scholar] [CrossRef]
  25. Alevizakos, V.; Chatterjee, K.; Koukouvinos, C.; Lappa, A. A double generally weighted moving average control chart for monitoring the process variability. J. Appl. Stat. 2023. [Google Scholar] [CrossRef]
  26. Akhtar, N.; Abid, M.; Amir, M.W.; Abbas, Z.; Nazir, H.Z.; Raza, Z.; Riaz, M. Design and analysis of exponentially weighted moving average control charts for monitoring the variability of log-normal processes with estimated parameters. Qual. Reliab. Eng. Int. 2022, 38, 1590–1611. [Google Scholar] [CrossRef]
  27. Chatterjee, K.; Koukouvinos, C.; Lappa, A.; Roupa, P. A joint monitoring of the process mean and variance with a generally weighted moving average maximum control chart. Commun. Stat.-Simul. Comput. 2023. [Google Scholar] [CrossRef]
  28. Chatterjee, K.; Koukouvinos, C.; Lappa, A. Monitoring process mean and dispersion with one double generally weighted moving average control chart. J. Appl. Stat. 2023, 50, 19–42. [Google Scholar] [CrossRef]
  29. Khan, I.; Khan, D.M.; Noor-ul-Amin, M.; Khalil, U.; Alshanbari, H.M.; Ahmad, Z. Hybrid EWMA Control Chart under Bayesian Approach Using Ranked Set Sampling Schemes with Applications to Hard-Bake Process. Appl. Sci. 2023, 13, 2837. [Google Scholar] [CrossRef]
  30. Ajibade, G.A.; Ajadi, J.O.; Kuboye, O.J.; Alih, E. Generalized new exponentially weighted moving average control charts (NEWMA) for monitoring process dispersion. Int. J. Qual. Reliab. Manag. 2023. [Google Scholar] [CrossRef]
  31. Haq, A.; Abbasi, A.A. New weighted adaptive CUSUM charts for monitoring the generalized variance of a bivariate normal process. J. Stat. Comput. Simul. 2023, 93, 604–633. [Google Scholar] [CrossRef]
  32. Yaqub, M.; Abbas, N.; Riaz, M.; Shabbir, J. On modified successive sampling based control charting schemes. Qual. Reliab. Eng. Int. 2016, 32, 2491–2497. [Google Scholar] [CrossRef]
  33. Abbas, N.; Riaz, M.; Mahmood, T. An improved S2 control chart for cost and efficiency optimization. IEEE Access 2017, 5, 19486–19493. [Google Scholar] [CrossRef]
  34. Hyder, M.; Mahmood, T.; Butt, M.M.; Raza, S.M.M.; Abbas, N. On the location-based memory type control charts under modified successive sampling scheme. Qual. Reliab. Eng. Int. 2021, 38, 2200–2217. [Google Scholar] [CrossRef]
  35. Jessen, R.J. Statistical Investigation of a Sample Survey for Obtaining Farm Facts; Jessen, R.J., Ed.; Lowa State University: Ames, IA, USA, 1942. [Google Scholar]
  36. Patterson, H. Sampling on successive occasions with partial replacement of units. J. R. Stat. Society. Ser. B. 1950, 12, 241–255. [Google Scholar] [CrossRef]
  37. Rao, J.; Graham, J.E. Rotation designs for sampling on repeated occasions. J. Am. Stat. Assoc. 1964, 59, 492–509. [Google Scholar] [CrossRef]
  38. Choudhary, R.; Bathla, H.; Sud, U. On non-response in sampling on two occasions. J. Indian Soc. Agric. Stat. 2004, 58, 331–343. [Google Scholar]
  39. Jones, L.A.; Champ, C.W.; Rigdon, S.E. The performance of exponentially weighted moving average charts with estimated parameters. Technometrics 2001, 43, 156–167. [Google Scholar] [CrossRef]
  40. Noorossana, R.; Fathizadan, S.; Nayebpour, M.R. EWMA control chart performance with estimated parameters under non-normality. Qual. Reliab. Eng. Int. 2016, 32, 1637–1654. [Google Scholar] [CrossRef]
  41. Mahmood, T.; Erem, A. A bivariate exponentially weighted moving average control chart based on exceedance statistics. Comput. Ind. Eng. 2023, 175, 108910. [Google Scholar] [CrossRef]
  42. Omar, M.H.; Arafat, S.Y.; Hossain, M.P.; Riaz, M. Inverse Maxwell Distribution and Statistical Process Control: An Efficient Approach for Monitoring Positively Skewed Process. Symmetry 2021, 13, 189. [Google Scholar] [CrossRef]
  43. Erem, A.; Mahmood, T. A bivariate CUSUM control chart based on exceedance statistics. Qual. Reliab. Eng. Int. 2023, 39, 1172–1191. [Google Scholar] [CrossRef]
  44. Iqbal, A.; Mahmood, T.; Ali, Z.; Riaz, M. On Enhanced GLM-Based Monitoring: An Application to Additive Manufacturing Process. Symmetry 2022, 14, 122. [Google Scholar] [CrossRef]
  45. Berardinelli, C. Short-Run Statistical Process Control Techniques. Available online: https://www.isixsigma.com/control-charts/short-run-statistical-process-control-techniques/ (accessed on 6 April 2023).
  46. Aslam, M.; Khan, N.; Jun, C.-H. A new S2 control chart using repetitive sampling. J. Appl. Stat. 2015, 42, 2485–2496. [Google Scholar] [CrossRef]
  47. Mashuri, M.; Haryono, H.; Aksioma, D.F.; Wibawati, W.; Ahsan, M.; Khusna, H. Tr (R2) control charts based on kernel density estimation for monitoring multivariate variability process. Cogent Eng. 2019, 6, 1665949. [Google Scholar] [CrossRef]
  48. Zaman, B.; Abbas, N.; Riaz, M.; Lee, M.H. Mixed CUSUM-EWMA chart for monitoring process dispersion. Int. J. Adv. Manuf. Technol. 2016, 86, 3025–3039. [Google Scholar] [CrossRef]
Figure 1. The effect of moving average span parameter (w) on the performance of the MA- S M S S ( S ) 2 control chart under MSS schemes for constant c = 2; (a) M S S ( 1 ) , (b) M S S ( 2 ) , (c) M S S ( 3 ) , (d) MSS(4), (e) M S S ( 5 ) .
Figure 1. The effect of moving average span parameter (w) on the performance of the MA- S M S S ( S ) 2 control chart under MSS schemes for constant c = 2; (a) M S S ( 1 ) , (b) M S S ( 2 ) , (c) M S S ( 3 ) , (d) MSS(4), (e) M S S ( 5 ) .
Symmetry 15 01108 g001
Figure 2. The effect of moving average span parameter (w) on the performance of the MA- S M S S ( S ) 2 control chart under MSS schemes for constant c = 3; (a) M S S ( 6 ) , (b) M S S ( 7 ) , (c) M S S ( 8 ) , (d) M S S ( 9 )  , (e) M S S ( 10 ) .
Figure 2. The effect of moving average span parameter (w) on the performance of the MA- S M S S ( S ) 2 control chart under MSS schemes for constant c = 3; (a) M S S ( 6 ) , (b) M S S ( 7 ) , (c) M S S ( 8 ) , (d) M S S ( 9 )  , (e) M S S ( 10 ) .
Symmetry 15 01108 g002
Figure 3. The effect of moving average span parameter (w) on the performance of the DMA- S M S S ( S ) 2 control chart under MSS schemes for constant c = 2; (a) M S S ( 1 ) , (b) M S S ( 2 ) , (c) M S S ( 3 ) , (d) M S S ( 4 ) , (e) M S S ( 5 ) .
Figure 3. The effect of moving average span parameter (w) on the performance of the DMA- S M S S ( S ) 2 control chart under MSS schemes for constant c = 2; (a) M S S ( 1 ) , (b) M S S ( 2 ) , (c) M S S ( 3 ) , (d) M S S ( 4 ) , (e) M S S ( 5 ) .
Symmetry 15 01108 g003
Figure 4. The influence of moving average span parameter (w) on the performance of the DMA- S M S S ( S ) 2 control chart under MSS schemes for constant c = 3; (a) M S S ( 6 ) , (b) M S S ( 7 ) , (c) M S S ( 8 ) , (d) M S S ( 9 ) , (e) M S S ( 10 ) .
Figure 4. The influence of moving average span parameter (w) on the performance of the DMA- S M S S ( S ) 2 control chart under MSS schemes for constant c = 3; (a) M S S ( 6 ) , (b) M S S ( 7 ) , (c) M S S ( 8 ) , (d) M S S ( 9 ) , (e) M S S ( 10 ) .
Symmetry 15 01108 g004
Figure 5. At a fixed constant c = 2 , a comparison of the best proposed MA- S M S S ( 1 ) 2 and DMA- S M S S ( 1 ) 2 charts with the S 2 chart based on repetitive sampling S r e p e t i t i v e 2 .
Figure 5. At a fixed constant c = 2 , a comparison of the best proposed MA- S M S S ( 1 ) 2 and DMA- S M S S ( 1 ) 2 charts with the S 2 chart based on repetitive sampling S r e p e t i t i v e 2 .
Symmetry 15 01108 g005
Figure 6. Stages of ZA fertilizer production process.
Figure 6. Stages of ZA fertilizer production process.
Symmetry 15 01108 g006
Figure 7. The probability plots of the concentration of C O 2 ( g / L ) with respect to body and glaze.
Figure 7. The probability plots of the concentration of C O 2 ( g / L ) with respect to body and glaze.
Symmetry 15 01108 g007
Figure 8. Demonstration for (a) MA- S M S S ( S ) 2 and (b) DMA- S M S S ( S ) 2 control charts using ZA fertilizer production dataset.
Figure 8. Demonstration for (a) MA- S M S S ( S ) 2 and (b) DMA- S M S S ( S ) 2 control charts using ZA fertilizer production dataset.
Symmetry 15 01108 g008
Table 1. Control charting coefficients L M S S ( S ) and K M S S ( S ) for the proposed control charts under different schemes of MSS (w = 2).
Table 1. Control charting coefficients L M S S ( S ) and K M S S ( S ) for the proposed control charts under different schemes of MSS (w = 2).
Scheme n = 5 n = 7
S S 2 M S E S 2 L M S S ( S ) K M S S ( S ) S S 2 M S E S 2 L M S S ( S ) K M S S ( S )
M S S ( 1 ) 1.080.3710.1817.140.970.2610.7217.74
M S S ( 2 ) 0.860.349.3915.640.860.269.99416.48
M S S ( 3 ) 0.730.388.3713.840.80.289.3615.35
M S S ( 4 ) 0.670.407.7612.80.760.298.9514.58
M S S ( 5 ) 0.640.427.4912.30.740.308.7814.27
M S S ( 6 ) 0.760.39.616.10.780.259.816.32
M S S ( 7 ) 0.550.407.1211.850.680.298.514.04
M S S ( 8 ) 0.460.486.0810.030.630.327.8112.8
M S S ( 9 ) 0.430.515.759.530.600.347.512.23
M S S ( 10 ) 0.410.525.579.170.590.357.3812.02
Table 2. Run length profile of the proposed charts under MSS at pre-specified A R L 0 = 370 , w = 2 and c = 2 .
Table 2. Run length profile of the proposed charts under MSS at pre-specified A R L 0 = 370 , w = 2 and c = 2 .
Schemes δ Shewhart - S M S S ( S ) 2 MA - S M S S ( S ) 2 DMA - S M S S ( S ) 2
n = 5 n = 7 n = 5 n = 7 n = 5 n = 7
ARLMDRLSDRLARLMDRLSDRLARLMDRLSDRLARLMDRLSDRLARLMDRLSDRLARLMDRLSDRL
M S S ( 1 ) 1369.56309.5369.08369.21308.1371.56371.03254.00375.92369.91260.00368.71370.61261.00367.68369.37251.00363.10
1.1193.53162194.28177.70149177.88181.63126.00180.77165.03113.00165.99180.83129.00184.70156.43108.00156.00
1.2113.2894.5113.4197.628297.61102.8372.00102.6785.1159.0085.59102.6473.0098.5183.5758.0082.99
1.372.1860.572.0959.705060.0965.1844.0062.9949.3634.0049.4264.2946.0064.2547.6134.0045.74
1.449.4541.549.4839.1632.539.3144.3030.0042.5831.7222.0029.8243.7931.0044.2531.6722.0031.76
1.535.5329.535.4627.242327.2631.8723.0030.2722.3516.0020.7031.8322.0031.5022.2715.0022.18
1.626.9722.526.9020.181720.3724.2717.0022.7016.5812.0014.8424.0917.0023.8216.3811.0016.44
1.721.2317.521.1315.3612.515.3318.7314.0017.7412.889.0011.4618.2013.0018.1012.549.0012.25
1.817.0414.516.9612.231012.1515.3811.0013.9510.588.009.1315.2511.0014.7810.187.009.93
1.914.0911.514.009.948.59.8113.3910.0011.948.957.007.5412.839.0012.168.096.007.83
211.889.511.788.326.58.2211.158.009.967.556.006.1310.848.0010.207.105.006.79
M S S ( 2 ) 1370.31310371.98370.72310.1375.07369.67249.00368.53369.49252.00369.52369.95250.00373.62369.56253.00357.75
1.1200.19167202.15183.68153.5185.50181.62124.00185.11161.31111.00163.25180.49127.00186.63160.13113.00159.44
1.2119.8899.5121.58102.4786103.13102.9671.00104.1085.6259.0086.67101.6872.00103.1784.0859.0084.42
1.377.3364.578.6862.625263.4666.0645.0067.8049.2134.0049.9564.6045.0063.2848.9535.0047.53
1.453.304454.1941.5934.542.5644.4631.0043.1932.3723.0031.2344.3830.0046.0231.6022.0031.74
1.538.633239.7129.2424.529.9432.2623.0031.2622.3816.0021.2731.7622.0032.7322.1215.0022.46
1.629.2924.530.1121.331821.8524.6018.0023.2516.7012.0015.4923.9517.0024.3216.1311.0016.35
1.722.691923.3216.3313.516.6119.6414.0018.8312.859.0011.5018.2713.0018.7512.449.0012.73
1.818.4915.519.1812.9310.513.2715.5311.0014.2710.538.009.2615.1610.0015.789.667.009.74
1.915.1512.515.6910.51910.7713.149.0012.188.886.007.7012.558.0012.948.065.008.22
212.7010.513.198.6578.8211.328.0010.237.545.006.2810.477.0010.666.554.006.56
M S S ( 3 ) 1370.12313.5375.77369.67306.5367.15370.35250.50375.53369.46256.00372.18370.45261.00372.14369.14260.00367.74
1.1206.36172.5208.55184.49154185.84186.20128.00189.26163.16113.00167.37184.92134.00193.63160.93112.00157.70
1.2124.14103.5126.43103.3286105.21105.8673.00108.6985.5359.0087.62103.2977.00108.0083.5358.0082.80
1.381.446883.0163.945365.2566.6044.0070.4649.6133.0051.8365.2348.0070.3849.3235.0047.16
1.456.464758.1742.173543.1144.6430.0046.7132.4323.0030.6843.5133.0046.8731.8122.0032.94
1.541.1034.542.6529.712530.5933.1523.0032.2122.6616.0021.7532.0821.0034.0521.9315.0023.10
1.631.102632.4922.0618.522.7625.6218.0025.1216.6412.0015.3023.9316.0025.3416.0711.0016.91
1.724.322025.5816.7013.517.3819.8314.0019.0013.159.0011.8418.7012.0019.9112.318.0013.00
1.819.491620.6113.2010.513.8416.1111.0015.2810.497.009.3815.1410.0016.269.656.0010.09
1.916.0013.516.9910.63911.1213.359.0012.448.906.007.7212.408.0013.487.845.008.28
213.401114.318.8979.2811.598.0010.727.485.006.3010.126.0010.886.604.006.91
M S S ( 4 ) 1370.03313.5379.63370.33313.5378.75369.81252.50371.42369.00256.00369.53370.51254.00376.92370.17257.00374.43
1.1207.91173.5211.00189.03157.5191.70185.47126.00192.29163.40111.00170.20185.00131.00191.72162.91114.50165.71
1.2127.41106.5129.93107.0289108.71105.9872.00109.4186.7860.0086.12105.2577.00110.2486.1059.0089.59
1.384.137086.4266.265567.7666.8345.0070.2951.1636.0049.7865.3549.0069.4449.7833.0052.16
1.458.4148.560.5844.0236.545.1645.1433.0045.7533.1824.0031.3244.4630.0047.3131.7121.0033.36
1.542.5135.544.1030.8625.531.9334.4724.0033.8922.9716.0022.0131.9421.0034.6822.1915.0023.34
1.632.4126.534.1322.661923.5726.1518.0025.7516.8912.0015.8123.7216.0025.6816.0711.0016.88
1.725.322126.9317.351418.3320.4914.0019.9713.019.0011.9018.5912.0020.3012.218.0013.18
1.820.2116.521.4913.681114.4216.4812.0015.6510.798.009.8614.9210.0016.319.656.0010.29
1.916.721417.9411.109.511.7213.7610.0013.089.016.008.0112.268.0013.697.855.008.36
214.0211.515.269.097.59.6211.628.0011.107.545.006.3910.216.0011.236.444.006.84
M S S ( 5 ) 1370.87311.5373.74370.04316377.20369.75250.00378.34371.98258.00375.82369.95252.00367.78369.82261.00358.28
1.1208.89174.5211.68188.64157.5190.19187.49129.00192.44167.65115.00172.97186.50129.00188.01161.28111.00162.45
1.2128.23107131.36107.9090109.52106.9773.00111.6486.9959.0091.39105.1675.00108.4485.0860.0082.69
1.385.327188.5666.7855.568.1967.2345.0071.0950.5435.0048.9966.9048.0067.6650.1234.0051.97
1.459.264961.3644.493746.0445.1930.0048.6832.9823.0032.1744.2134.0047.1332.1622.0033.95
1.543.0735.545.1431.102632.1732.1521.0035.3621.9116.0021.8231.5924.0034.0921.6415.0023.81
1.632.602734.4823.101924.0125.9918.0025.7616.7612.0015.5224.0116.0026.7516.1211.0017.26
1.725.612127.4817.5514.518.4220.5314.0019.8513.259.0012.2618.6012.0020.6012.398.0013.36
1.820.5416.522.1813.871114.5916.6412.0015.9610.678.009.6814.999.0016.929.636.0010.59
1.916.7913.518.2711.229.511.9013.8110.0013.348.916.007.9612.278.0013.727.955.008.55
214.0611.515.449.227.59.8111.438.0010.757.595.006.579.956.0011.236.454.006.80
Table 3. Run length profile of the proposed charts under MSS at pre-specified A R L 0 = 370 , w = 2 and c = 3 .
Table 3. Run length profile of the proposed charts under MSS at pre-specified A R L 0 = 370 , w = 2 and c = 3 .
Schemes δ Shewhart - S M S S ( S ) 2 MA - S M S S ( S ) 2 DMA - S M S S ( S ) 2
n = 5 n = 7 n = 5 n = 7 n = 5 n = 7
ARLMDRLSDRLARLMDRLSDRLARLMDRLSDRLARLMDRLSDRLARLMDRLSDRLARLMDRLSDRL
M S S ( 6 ) 1370.83311379.94370.13313377.80370.40257.00371.47369.86250.00370.83370.31254.00375.76370.99256.00368.59
1.1208.43173213.70187.60156191.13195.24134.00201.60167.49115.00171.74194.09133.00194.52165.25124.00176.65
1.2126.35105130.14106.3989108.87116.9979.00123.0788.4466.0091.71115.6380.00115.9190.6061.0093.61
1.382.6268.585.9865.7554.568.3274.9750.0078.5855.5638.0055.7174.6751.0077.0055.0537.0057.81
1.458.0247.561.3143.293645.2550.9934.0054.3335.9625.0035.3750.6836.0052.9435.7524.0038.52
1.542.323545.4830.702532.4237.2225.0039.9025.2117.0024.7637.0926.0037.3524.1816.0025.77
1.631.922634.5722.4718.523.8028.7719.0029.3718.8213.0018.3827.7220.0028.6217.9212.0019.15
1.725.012027.4217.151418.3523.0616.0023.2615.0511.0014.3721.6514.0023.8813.399.0014.45
1.820.2216.522.2213.511114.6418.5012.0018.8611.708.0011.1517.2911.0019.0810.787.0011.59
1.916.6313.518.4410.928.511.9615.2610.0015.809.777.008.8814.629.0016.328.565.009.37
213.881115.639.097.59.8813.099.0013.188.196.007.4112.298.0013.817.194.007.65
M S S ( 7 ) 1369.52307376.50370.10313379.64370.76256.00388.48369.39248.00372.77370.38255.00373.82369.72264.00371.97
1.1211.02175.5219.56190.71158195.74198.16135.00211.65168.35111.00172.04197.03135.00203.24167.05121.00166.65
1.2131.91109139.26108.7590112.75115.5977.00126.9492.5960.0094.35114.0981.00126.9891.5764.0091.41
1.387.4671.593.8068.3256.571.9974.2347.0083.8057.3036.0057.0473.8552.0078.4455.9739.0056.26
1.461.645066.8945.4937.548.1349.7131.0056.7137.9622.0037.6448.3736.0056.0435.6024.0035.41
1.545.1336.550.2032.1726.534.7236.5822.0042.1727.3214.0026.1835.9825.0039.9025.1117.0025.47
1.634.2127.538.5823.621925.9127.7419.0030.2918.8113.0018.7827.2417.0031.5617.0611.0019.33
1.726.682130.6617.9814.519.7722.5314.0024.4214.3410.0014.1420.6312.0024.9613.138.0014.82
1.821.611725.2114.121115.8318.9412.0020.2111.538.0011.2016.429.0020.1010.096.0011.45
1.917.6513.520.8511.36912.8115.4010.0016.459.556.009.0513.727.0017.058.305.009.44
214.7811.517.649.40710.6913.118.0014.148.035.007.6011.306.0014.256.834.007.94
M S S ( 8 ) 1370.11311.5385.43369.13306.5372.65369.21251.00383.21369.71248.00381.95370.38260.00376.17370.92255.00378.80
1.1217.49180.5228.07190.47158195.58197.04130.00216.21168.85112.00179.77196.47136.00207.60167.63121.00175.42
1.2135.84111.5145.65110.4792114.50116.0675.00131.4089.4759.0095.15115.8781.00125.9588.7965.0092.54
1.391.237499.8469.395773.4773.8446.0084.7854.4436.0059.8172.2752.0082.8853.1038.0057.29
1.464.465271.8146.3838.549.5348.8029.0058.1634.7922.0039.3947.1135.0058.3133.7225.0037.26
1.547.2837.553.5832.682735.6336.2321.0043.9823.4214.0026.6935.5425.0042.8022.6017.0026.13
1.635.712841.3924.3319.526.8726.5515.0032.3217.9113.0018.9825.7319.0032.8217.2110.0019.93
1.728.142233.3718.4514.520.7320.5811.0025.7214.6210.0014.4820.5115.0026.2312.938.0014.93
1.822.6917.527.2214.4411.516.4017.8912.0020.8911.698.0011.6216.158.0021.0310.156.0011.98
1.918.461422.6411.61913.3212.916.0017.239.456.009.1911.649.0017.338.034.009.68
215.3211.519.019.587.511.1411.125.0014.658.105.007.5710.008.0014.526.644.007.90
M S S ( 9 ) 1370.46313.5388.17369.15306.5375.18370.32250.00392.33369.71250.00376.49369.39253.00377.28370.61252.00374.99
1.1219.78180.5232.80193.73160.625198.72198.35130.00221.03168.32113.00178.71197.35130.00210.57168.15120.00177.68
1.2138.21113149.28112.0592.5116.81118.2375.00135.9389.9860.0096.09116.8977.00127.2088.0464.0095.91
1.392.8975.5101.9870.8258.574.9174.3744.0087.7555.0336.0061.2273.1343.0082.2954.1639.0057.19
1.465.6652.573.5547.1938.550.8649.8429.0060.1634.9922.0040.0348.5228.0060.2433.9625.0038.15
1.548.3438.555.4333.4327.536.6136.2721.0044.3823.4314.0026.8235.5820.0044.0222.6317.0025.64
1.636.652943.0424.5119.527.4126.4614.0033.0917.4510.0020.3525.7913.0033.6516.2813.0019.19
1.728.732234.6318.721521.2720.2910.0026.1212.877.0015.1819.2710.0026.7211.0010.0014.87
1.822.8817.528.1414.6211.516.8014.429.0021.6310.136.0012.1113.978.0020.859.028.0011.85
1.918.6514.523.1211.839.513.7112.146.0017.927.974.009.7412.016.0017.536.606.009.53
215.6611.519.729.687.511.369.635.0015.226.533.007.798.865.0014.975.135.007.94
M S S ( 10 ) 1370.20306.5386.05369.79309376.85369.95246.00400.58370.58251.00384.72369.99248.00381.01370.15251.00364.44
1.1215.30177.5229.16193.10161197.78195.87126.00220.30167.96114.00179.01194.50137.00211.75164.12123.00177.52
1.2136.48111.5148.08113.4394118.53116.6072.00137.3990.6361.0096.49114.8983.00127.0988.8165.0094.08
1.391.6774102.0971.115975.4873.6543.0088.0256.0237.0062.1972.6352.0087.2355.0239.0058.06
1.465.055273.8147.943951.9948.6028.0060.0034.7722.0040.1546.1935.0061.2233.6126.0038.40
1.547.893855.9533.5627.536.8535.9219.0045.0423.8514.0027.6334.5925.0044.0222.5718.0026.62
1.636.2228.543.4924.992027.9825.9413.0033.2617.4510.0020.2024.8119.0034.0217.3613.0019.51
1.728.3221.534.5818.881521.4920.4214.0025.5513.9810.0015.1720.1710.0026.7513.088.0015.47
1.822.611728.1214.8511.517.1816.3112.0021.5310.088.0012.0415.817.0021.4310.216.0012.36
1.918.3913.523.3711.85913.9412.579.0017.747.807.009.6812.485.0017.417.974.009.67
215.151119.549.737.511.5111.158.0014.886.265.008.1510.674.0014.706.653.008.05
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hyder, M.; Raza, S.M.M.; Mahmood, T.; Abbas, N. Enhanced Dispersion Monitoring Structures Based on Modified Successive Sampling: Application to Fertilizer Production Process. Symmetry 2023, 15, 1108. https://doi.org/10.3390/sym15051108

AMA Style

Hyder M, Raza SMM, Mahmood T, Abbas N. Enhanced Dispersion Monitoring Structures Based on Modified Successive Sampling: Application to Fertilizer Production Process. Symmetry. 2023; 15(5):1108. https://doi.org/10.3390/sym15051108

Chicago/Turabian Style

Hyder, Mehvish, Syed Muhammad Muslim Raza, Tahir Mahmood, and Nasir Abbas. 2023. "Enhanced Dispersion Monitoring Structures Based on Modified Successive Sampling: Application to Fertilizer Production Process" Symmetry 15, no. 5: 1108. https://doi.org/10.3390/sym15051108

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop