Next Article in Journal
Type-3 Fuzzy Control of Robotic Manipulators
Next Article in Special Issue
Estimation of the Modified Weibull Additive Hazards Regression Model under Competing Risks
Previous Article in Journal
Halo Orbits under Some Perturbations in cr3bp
Previous Article in Special Issue
Computational Analysis for Fréchet Parameters of Life from Generalized Type-II Progressive Hybrid Censored Data with Applications in Physics and Engineering
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Improved Charting Scheme to Monitor the Process Mean Using Two Supplementary Variables

1
Department of Mathematics and Statistics, Institute of Southern Punjab, Multan 59300, Pakistan
2
Department of Mathematics, College of Arts and Science, Prince Sattam bin Abdul Aziz University, AL-Kharj 16242, Saudi Arabia
3
Department of Basic Science, College of Science and Theoretical Studies, Saudi Electronic University, Riyadh 11673, Saudi Arabia
4
Department of Business Administration, College of Business Administration, Princes Nourah bint Abdulrahman University, P.O. Box 84428, Riyadh 11671, Saudi Arabia
5
Department of Mathematics, Women University of Azad Jammu and Kashmir, Bagh 12500, Pakistan
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(2), 482; https://doi.org/10.3390/sym15020482
Submission received: 23 December 2022 / Revised: 6 February 2023 / Accepted: 6 February 2023 / Published: 11 February 2023

Abstract

:
A control chart is the most well-known statistical monitoring tecnique to address unfavourable process parameter (s) changes. Quality practitioners always desire a charting device that promptly identifies the undesired changes in the process. This study intends to design a sensitive homogeneously weighted moving average chart using two supplementary variables (hereafter, TAHWMA). The two supplementary variables are correlated with the study variable in the form of a regression estimator, which is an efficient and unbiased estimator for the process mean. The suggested TAHWMA charting structure is checked out and compared in terms of appearance and non-appearance of multicollinearity amidst the two additional variables. Average run length-related measures are taken as performance measures. It is observed that the proposed TAHWMA scheme performs effectively when the two supplementary variables have no collinearity. A comprehensive comparison between the proposed TAHWMA and existing charts is also carried out, showing the proposed’s supremacy over existing counterparts. For execution purposes, two illustrative examples, one belonging to carbon fibre manufacturing-related data and the other using a simulated dataset and where our simulated dataset belongs to symmetrical distribution, are also presented for the application of the recommended TAHWMA chart.

1. Introduction

Quality control is an essential aspect of production management. Many management and engineering techniques are widely used to maintain the quality of goods and services that fulfil increasing customer demand. Companies can monitor their processes with the application of control charts for producing high-quality products. In an ongoing process, change/variation is an inevitable output factor. Statistical process control (SPC) is handy for controlling the variation of methods. Variations are divided into two significant categories: natural and unnatural variations, respectively. Suppose natural changes occur in any running process. In that case, the process is assumed to be statistically in-control (IC), while unusual changes lead the running process to an out-of-control (OOC) state (Montgomery [1]).
Control charts are customarily utilized when unnatural variations exist in the process. The charting mechanism based on the Shewhart [2] structure is a type of memoryless chart because it accepts only recent sample information. The cumulative sum (CUSUM) charting scheme is an example of a memory-type charting mechanism that was originated by Page [3]. Another memory-type control chart, namely the exponentially weighted moving average (EWMA), was initiated by Roberts [4], and the homogenously weighted moving average proposed by Abbas [5] is also a memory-type chart. The above memory-type charts utilize the previous information along with recent sample information. In any manufacturing industry, SPC quality inspectors have adopted several supplementary techniques to improve the output of any continuing process.
To spot a shift in the process parameter (s) efficiently and enhance the control chart’s sensitivity in SPC literature, utilising the supplementary information and different sampling schemes is considered the best option. Recently, the concept of supplementary information has been tested for enhancing the performance of the existing charting schemes by many researchers. The charting situation when the study characteristic is observed correlated with another supplementary feature; such a structure is called an AIB-based charting design (cf. Haq and Khoo [6]). The extraneous information at the estimation phase or information other than the sample is called supplementary or supplementary information.
The regression-based control charts were designed by Mandel [7] and Zhang [8] to screen the process. The utilization of supplementary information in the control chart was initiated by Riaz [9,10] for supervising the process dispersion and location parameters, respectively. Riaz et al. [11] suggested a new AIB-Shewhart control chart under normal and non-normal scenarios for monitoring the process. Furthermore, the (AIB-GWMCV) chart was proposed by Nuriman et al. [12]. Haq and Khoo [13] developed an AIB-synthetic charting scheme to improve the efficiency of the synthetic mean model. Recently, the monitoring of the coefficient of variation using supplementary information has been initiated by Abbasi [14]. In addition, Nuriman et al. [12] proposed a supplementary information-based control chart for efficiently monitoring the process CV. Muhammad Arslan et al. [15] studied a mixed EWMA Dual Crosier CUSUM chart with and without supplementary information. The interested readers are referred to Abbas et al. [16], Chen, J. H., & Lu, S. L. [17], Rasheed et al. [18], Aslam et al. [19], Anwar et al. [20], Rasheed et al. [21], and Zhang et al. [22] for more recent work on HWMA structure and supplementary information-based charting schemes.
Recently, Zichuan et al. [23] designed two AIB EWMA (name hereafter; TAEWMA) control charts to monitor small shifts in the process mean promptly. Taking inspiration from Zichuan et al. [23], this study proposes a new HWMA charting scheme based on two supplementary variables (TAHWMA) to monitor small changes in the process mean quickly. The proposed TAHWMA chart is also investigated under the appearance and non-appearance of the multicollinearity behaviour among the two supplementary variables. The performance of the proposed design has been evaluated using run-length (RL) characteristics, where run length is defined as the number of samples before a chart signal. The expected value of the run lengths is recognized as average run length (ARL). A R L 0 and A R L 1 known as IC ARL and OOC ARL values, respectively.
The rest of the article is organized as follows: The design structures of the existing control charts are briefly described in Section 2. The designed structure of the proposed TAHWMA chart and RL evaluation is presented in Section 3. The comparison of the proposed TAHWMA chart against the existing charts is presented in Section 4. Illustrative examples are conferred in Section 5, and the conclusions of this article are provided in Section 6.

2. Design Structures of Some Existing Control Charts

In this section of the article, existing counterparts are briefly described. The classical HWMA, AHWMA, classical EWMA, and AEWMA control charts are the competitors of the proposed TAHWMA control chart.

2.1. Design Structure of the HWMA Control Chart

Assume the quality characteristic Y i j ~ N ( μ , σ 2 ) to be monitored; when the process is IC, it is assumed that both process parameters μ and σ are known, i.e., μ = μ 0   and   σ = σ y , respectively. To monitor the process mean, the classical HWMA chart was recently designed by Abbas [5]. The plotting statistic of the HWMA chart is given as:
H i = λ Y ¯ i + ( 1 λ ) Y ¯ ¯ i 1 ,
where λ is known as smoothing parameter that is parameter of the HWMA charting scheme; its value lies between 0 and 1. The HWMA chart is more efficient at small choices of the smoothing parameter λ . Abbas [5] showed that the Shewhart control chart becomes the particular case of the HWMA chart at λ = 1. Where Y ¯ ¯ i 1 = k = 1 i 1 Y ¯ k i 1 is the mean of the remaining ( i 1 ) sample means, and The Y ¯ i is the sample average of i t h the sample. Control limits of the HWMA are given as:
L C L = μ 0 L λ 2 σ y 2 n , if   i = 1 μ 0 L λ 2 σ y 2 n + ( 1 λ ) 2 σ y 2 n ( i 1 ) , if   i   >   1 C L = μ 0 , U C L = μ 0 + L λ 2 σ y 2 n , if   i = 1 μ 0 + L λ 2 σ y 2 n + ( 1 λ ) 2 σ y 2 n ( i 1 ) , if   i   >   1
where the width of the control chart is known as L. The HWMA chart produces an OOC signal if the H i statistic presented in Equation (1) goes beyond the control limits described in Equation (2).

2.2. Design Structure of the AHWMA Control Chart

For monitoring the process location, the supplementary-based HWAM (AHWMA) chart was proposed by Adegoke et al. [20]. According to Adegoke et al. [24], the study variable Y i j is correlated with the supplementary variable X i j (Cochran [25]). The regression can be expressed as:
R i = Y ¯ i + b Y X ( μ X X ¯ i )
where b Y X = ρ Y X σ Y σ X is expressed as the regression coefficient. The mean and the variance of the regression estimator are μ R = μ Y and σ R 2 = σ Y 2 n ( 1 ρ Y X 2 ) , respectively. The AHWMA charting scheme is given below,
D i = λ R i + ( 1 λ ) R ¯ i 1 .
In Equation (4) R i is the regression estimate of the process variable while R ¯ i 1 is the average of all previous samples. The control limits of the AHWMA chart are expressed as:
L C L = μ 0 L σ Y λ 2 n ( 1 ρ y x 2 ) , if   i = 1 μ 0 L σ Y ( λ 2 n + ( 1 λ ) 2 n ( i 1 ) ) ( 1 ρ y x 2 ) , if   i   >   1 C L = μ 0 , U C L = μ 0 + L σ Y λ 2 n ( 1 ρ y x 2 ) , if   i = 1 μ 0 + L σ Y ( λ 2 n + ( 1 λ ) 2 n ( i 1 ) ) ( 1 ρ y x 2 ) , if   i   >   1
where L is the width of control limits of the AHWMA control chart and λ is chosen to achieve a desired IC ARL for the chart. The AHWMA chart produces an OOC signal if the D i statistic presented in Equation (4) goes outside the control limits given in Equation (5).

2.3. Design Structure of the Classical EWMA Control Chart

The classical EWMA charting scheme to examine the process location was suggested by Roberts [4]. The classical EWMA charting statistic is given below:
S i = λ Y ¯ i + ( 1 λ ) S i 1 ,
where S i 1 is the information at (i − 1)th time E ( S i ) = μ 0 , and V a r ( S i ) = σ Y ¯ 2 ( λ 2 λ ( 1 ( 1 λ ) 2 i ) ) are the mean and variance for the IC process. The control limits for the classical EWMA charts are:
L C L = μ 0 L σ Y ¯ λ 2 λ ( 1 ( 1 λ ) 2 i ) , C L = μ 0 , U C L = μ 0 + L σ Y ¯ λ 2 λ ( 1 ( 1 λ ) 2 i ) .
It is interesting to note that the Shewhart chart is a special case of EWMA at λ = 1. The classical EWMA chart detects an OOC state if any plotting statistic S i falls beyond the control limits described in Equation (7).

2.4. Design Structure of the AEWMA Control Chart

The supplementary variable-based EWMA chart for monitoring the process location was suggested by Abbas et al. [12]. Let us assume X i known as a supplementary variable and is associated with the variable of interest Y i . The term ρ Y X is known as the correlation between the two variables. The bivariate symmetrical distribution can be expressed as ( Y , X ) ~ N 2 ( μ Y , μ X , σ Y 2 , σ X 2 , ρ Y X ) . For monitoring the population mean μ 0 , the regression estimator based on supplementary information is given as (Cochran [25]):
C Y = Y ¯ + b Y X ( μ X X ¯ ) ,
where b Y X = ρ Y X σ Y σ X is presenting change the measures Y due to a one-unit change in X, and in Equation (9) mean and variance of the regression estimator is given below
E ( C Y ) = μ 0 , V a r ( C Y ) = σ Y 2 = σ Y 2 n ( 1 ρ Y X 2 ) = σ Y 2 b Y X 2 σ X 2 n ,
C Y is an unbiased estimator of μ 0 and V a r ( C Y ) < V a r ( Y ) ¯ for ρ Y X 2 > 0 . The AEWMA statistic based on the regression estimator is defined as
E i = λ C Y + ( 1 λ ) E i 1
The value of EWMA statistics E i 1 presents past information and is taken from the initial to (i − 1) sample group. The value E 0 is commonly accepted as equal to the target mean μ 0 . The variance V a r ( C Y ) = σ Y 2 = σ Y 2 n ( 1 ρ Y X 2 ) = σ Y 2 b Y X 2 σ X 2 n and its means are the target value of the process. The control limits of the AHWMA control chart are given as;
L C L = μ 0 L σ Y ( 1 ρ Y X 2 ) ( λ 2 λ ( 1 ( 1 λ ) 2 i ) ) , C L = μ 0 , U L C = μ 0 + L σ Y ( 1 ρ Y X 2 ) ( λ 2 λ ( 1 ( 1 λ ) 2 i ) ) .
The AEWMA charting scheme detects OOC if any plotting statistic falls beyond the control limits.

3. Proposed TAHWMA Control Chart

The design structure of the proposed TAHWMA control chart and performance metrics are discussed in this section. The proposed TAHWMA structure is designed under the presence and absence of multicollinearity.

3.1. Design Structure of the Proposed TAHWMA Control Chart

This section of the article provides the detailed design structure of the proposed TAHWMA charting scheme. Here three variables are selected from a trivariate symmetrical distribution such as Y ,   X , and Z . X and Z are the supplementary variables, and Y presents the study variable. The matrix form of the variables is organized below:
Y X Z ~ N 3 μ Y μ X μ Z , σ y y σ y x σ y z σ x y σ x x σ x z σ z y σ z x σ z z
The regression-based estimator initiated by Kadilar and Cingi [26] is given as,
G i = y ¯ + b y x ( μ x x ¯ ) + b y z ( μ z z ¯ ) ,
where b y x = s y x s x x and b y z = s y z s z z are the regression coefficient, s y x and s y z are sample covariances of Y ,   X and Y ,   Z . While s x x and s z z are the sample variances of X and Z , respectively. The mean and variance of Equation (12) are E ( G ) = μ G = μ 0   V a r ( G ) = σ G 2 = ( 1 ρ y x 2 ρ y z 2 + 2 ρ y x ρ y z ρ x z ) σ Y 2 n , respectively (cf. Kadilar and Cingi [27]). The proposed TAHWMA charting scheme based on the regression estimator (in Equation (12)) is given as,
H i = λ G i + ( 1 λ ) G ¯ ¯ i 1
L C L = μ 0 L λ 2 σ G 2 n , if   i = 1 μ 0 L λ 2 σ G 2 n + ( 1 λ ) 2 σ G 2 n ( i 1 ) , if   i > 1   CL = μ 0 , U C L = μ 0 + L λ 2 σ G 2 n . if   i = 1 μ 0 + L λ 2 σ G 2 n + ( 1 λ ) 2 σ G 2 n ( i 1 ) . if   i > 1
Control limits of the proposed TAHWMA control chart are presented in Equation (14), where L represents the width of the control limits. Simulation codes are established in R software for the performance evaluation of the proposed chart. The amount of shift in the process mean can be mathematically expressed as δ = μ 1 μ 0 σ y / n , where μ 1 denotes the shifted mean of the study variable and n = 1, has taken without loss generality. We use 50,000 iterations in simulation to find the desired average run length. The proposed TAHWMA control chart has two designed parameters λ , and L . λ is identified as the smoothing parameter, and the various values of the λ is λ { 0.03 , 0.05 , 0.1 , 0.25 } are considered in this study.

3.2. Performance Metrics

In this section, a detailed discussion of the run length properties of the proposed TAHWMA control chart is carried out using run-length (RL) characteristics. The run length (RL) is defined as the number of sample observations before a chart alarm describes as run-length (RL). Extra quadratic loss (EQL) is also taken as the performance measure of run-length properties, describing the charting schemes’ overall effectiveness. EQL is mathematically defined as E Q L = 1 δ max δ min δ min δ max σ 2 A R L ( δ ) d δ . Another performance measure is the relative mean index (RMI). Mathematically it is defined as R M I = 1 N i = 1 N A R L ( δ ) A R L * ( δ ) A R L * ( δ ) . Where N represents the number of shifts to be considered. For the specific shift δ , A R L ( δ ) is the A R L 1 value of a control chart, while A R L * ( δ ) represents the smallest value of the A R L 1 . The percentage decrease in ARL is A R L d , where A R L d = A R L 0 A R L 1 A R L 0 × 100 % . A R L 0 shows the ARL when the process is working in stable conditions and A R L 1 presents ARL values when the process is in an OOC situation. For each value of λ , the A R L 0 is fixed at 500 using the Monte Carlo simulation method. A comprehensive discussion on the RL distribution of the newly suggested TAHWMA charting scheme is examined under the appearance and non-appearance of the multicollinearity among the two supplementary characteristics.

3.2.1. Performance of the Proposed TAHWMA Control Chart under the Non-Appearance of Multicollinearity

In this section, we examined the performance of the suggested TAHWMA charting scheme under the appearance and non-appearance of multicollinearity among the two supplementary variables. If there is no relationship between the supplementary variable such as ρ x z = 0, both supplementary variables have a partial effect on the study variable. Under the appearance of multicollinearity, the ARL values of the suggested TAHWMA charting scheme at several choices of the parameters are given in Table 1. For tracing the small shift in the process mean, it is obvious that small choices of smoothing parameter, high values of the ρ YX   and ρ YZ   , and the suggested TAHWMA charting scheme becomes more sensitive. The suggested TAHWMA charting scheme with designed parameters λ = 0.03, L = 2.272, ρ YX     = 0.25 and ρ YZ = 0.50 , provides ARL 1 = 330.14 and A R L d = 33.97% at δ = 0.05. At λ = 0.03, L = 2.272, ρ YX     = 0.75 and ρ YZ = 0.50 , at 3% increase in the process mean the suggested TAHWMA charting scheme yields A R L d = 65.05% (see Table 1). The effect of the ρ YX   and ρ YZ on the performance of the proposed TAHWMA model can be seen in Figure 1 and Figure 2. From Figure 1 and Figure 2, it can be observed that the values of ρ YX   and ρ YZ increase, and the performance of the suggested TAHWMA charting scheme also becomes highly sensitive.

3.2.2. Performance of the Proposed TAHWMA Control Chart under the Appearance of Multicollinearity

If some linear relationship occurs between the two supplementary variables, this term is expressed as multicollinearity (cf. Hocking and Pendleton [26]). Unfortunately, multicollinearity happens due to mistakes or a lack of understanding of the model. It is vital to know about the effect of multicollinearity between two supplementary variables on the performance of the suggested TAHWMA charting scheme. The performance of the suggested TAHWMA charting scheme is reported at various choices of correlation between two supplementary characteristics, i.e., ρ x z   0.05, 0.15, and 0.25 (cf. Table 2). In Table 2, it can be noticed that as the value ρ x z increases, the performance of the suggested TAEWMA charting structure decreases. For example, at λ = 0.05, L = 2.633 and ρ x z = 0.05 at δ = 0.05, the proposed scheme provides A R L 1 = 349.49 and at ρ x z = 0.25, λ = 0.05 it gives A R L 1 = 354.30 respectively (cf. Table 2).

4. Comparative Study

In this section, a comparison of the proposed TAHWMA chart is provided against some existing control charts. The ARL is a comparative index for classical EWMA, AEWMA, HWMA, AHWMA, and proposed TAHWMA charts.

4.1. Proposed Versus EWMA and AEWMA Control Charts

The classical EWMA charting scheme is famous for examining the small shifts in the process; the classical EWMA chart was initiated by Roberts [4]. Abbas et al. [12] suggested a new charting scheme, AEWMA, for the efficient monitoring of the process mean. In Table 3, ARL values of both existing charting schemes such as classical EWMA and AEWMA are described. At λ = 0.05 and δ = 0.075 classical EWMA yields A R L 1 = 338.65 and for AEWMA at λ = 0.05 and δ = 0.075 ρ y x = 0.25, A R L 1 = 328.86 respectively. The proposed TAHWMA chart at λ = 0.05, δ = 0.075, ρ y x = 0.25 and ρ y z = 0.50 produces A R L 1 = 254.06, respectively (cf. Table 1). From Table 1 and Table 3, we noticed that the suggested TAHWMA charting scheme performs more efficiently than AEWMA and classical EWMA charting schemes to detect small and moderate shifts in the process mean level.

4.2. Proposed Versus HWMA and AHWMA Control Charts

The HWMA scheme was recently developed by Abbas [5] to address small changes in the process mean. To increase the sensitivity of the HWMA model, Adegoke et al. [24] designed the AHWMA chart using a regression estimator that is based on the single supplementary variable. The values of ARL for both classical HWMA and AHWMA control charts are given in Table 4 at various combinations of designed parameters. For example, At λ = 0.03 the classical HWMA scheme yields A R L 1 = 205.43 at δ = 0.1. For the AHWMA chart, at ρ y x = 0.25, λ = 0.3 and δ = 0.1, the value of A R L 1 = 199.38 (cf. Table 4). The suggested TAHWMA charting scheme at λ = 0.03, δ = 0.1 and ρ y x = 0.25 and ρ y z = 0.50 yields A R L 1 = 165.07 (cf. Table 1). The supremacy of the suggested TAHWMA charting scheme is obvious over the classical HWMA and AHWMA schemes (cf. Table 1 and Table 4).
In Table 5, the relative mean index (RMI) and extra quadratic loss (EQL) are also considered performance measures. For the proposed TAHWMA control chart λ ( 0.03 , 0.05 , 0.1 , 0.25 ) is used when ρ y x = 0.75, ρ y z = 0.50, and ρ x z = 0. Other existing charting schemes such as HWMA, EWMA, AEWMA, and AHWMA are also considered when λ = 0.03. For the proposed TAHWMA chart, when λ = 0.03, RMI has a minimum value of zero and the proposed TAHWMA chart also has the minimum value of the EQL, which is 2.12. The minimum values of the RMI and EQL show the supremacy of the suggested TAHWMA charting scheme against its existing charting structure.
In Figure 3, the supremacy of the proposed TAHWMA charting structure (in terms of smallest ARL1) is evident compared to all the competitors at small, moderate and large shifts in the process mean. The performance of the AHWMA scheme is the second best.

5. Illustrative Example

Along with exploring the several properties of the suggested TAHWMA charting scheme hypothetically, it is necessary to offer an illustrative example of the developed design. Two examples are presented in this section to show the supremacy of the proposed TAHWMA chart; one consists of the real dataset and the second is based on simulated data. The first real dataset belongs to the carbon fibre tube (manufacturing data), and the second dataset is taken from the trivariate normal distribution.

5.1. Real-Life Application

In this section, a real dataset example related to manufacturing carbon fibre (graphite fibre) tubes is used to show the application of the proposed TAWHMA chart and its competitors. This carbon dataset is used from the R statistical package “MSQC” and is also used by Abbasi and Adegoke [28]. Carbon fibre tubes are frequently comprised of carbon atoms. It has numerous assets for illustration, high tensile strength, low weight, etc. (cf. Zhang et al. [29]). It is a polymer that is solid but lightweight. It is used to manufacture aircraft, cars, and many other machinery parts. Carbon fibre is produced from polyacrylonitrile (PAN) and a small quantity of petroleum pitch. Carbon fibre is much stronger than steel; it is a fibre having many characteristics, and the size of the fibre is nearly equal to 5 to 10 micrometres. It is low-density material but has high-temperature tolerance, high chemical resistance, high stiffness, and high tensile strength.
In this analysis, Y represents the inner diameter, X represents thickness, and Z represents length. These are the first and second supplementary variables, respectively. These variables belong to the manufacturing process of the fiber tubes. A dataset has 30 samples, where 10 samples are taken with 0.05 amount of shift in the process mean, and 20 samples are taken as the IC process. For the execution of the proposed TAHWMA plotting scheme and prevailing schemes, planned parameters are set at the desired value of ARL = 500 for all the charts included in this example. Some descriptive statistics of the carbon fibre data of the first 20 samples are; X ¯ is 1.022, and Z ¯ is 50.1005 and Y ¯ is 0.993. where s y 2 = 0.001706316, s x 2 = 0.01139579, s z 2 = 0.05572079, ρ y x = 0.2778609, ρ y z = 0.617334, ρ x z = 0.1219354. Figure 4 and Figure 5 show that the HWMA and AHWMA charts are not detecting a 0.05 amount of shift at any sample, whereas the proposed TAHWMA charting structure drops this change at the 25th to 30th sample observations (cf. Figure 6).
In detecting the small changes in the process mean compared to its counterparts, the real-life application indicates that the suggested TAHWMA charting structure is the most powerful tool.

5.2. Simulation Study

This section describes the application using the hypothetical dataset to implement the planned design. In this simulated part of the study, the dataset consists of 50 samples. These 50 samples belong to trivariate symmetrical distribution; for the IC process, 25 samples are considered, and 25 samples are incorporated with a shift of 0.5 in the process mean for the OOC process. The execution of the recommended charting structure is executed considering the appearance and non-appearance of multicollinearity among the two supplementary variables. Further details on the selection of designed parameters for the proposed and existing counterparts are given in Figure 7, Figure 8, Figure 9 and Figure 10. From Figure 7, Figure 8, Figure 9 and Figure 10, it can be observed that all the charts remain in the IC situation for the first 25 samples. In Figure 7, the application of the proposed TAHWMA chart is provided when there is no multicollinearity between the supplementary variables, and it traces a shift at the 39th samples. In Figure 8, the presentation of the proposed TAHWMA control chart is delivered when multicollinearity exists between the supplementary variables. The proposed TAHWMA chart traces the shift at the 41st samples in this case. In Figure 9, the existing chart, namely the HWMA chart, outlines the shift at the 49th samples. In Figure 10, the existing chart, namely the AHWMA chart, detects the shift at 49th samples. Based on the discussion, it is observed that the proposed TAHWMA charting scheme has performed efficiently against existing charting schemes, particularly in the non-appearance of multicollinearity among the two supplementary variables.

6. Summary, Conclusions and Recommendations

The single supplementary information-based (AIB) homogeneously weighted moving (HWMA) chart is an innovative version of the HWMA charting scheme to monitor the process mean shift. This study aims to enhance the HWMA and AIB HWMA charts and propose an HWMA chart based on two supplementary variables, indicated as the TAHWMA charting scheme, to improve the process mean shift monitoring. The run-length characteristic of the proposed TAHWMA charting scheme has been discussed in the absence and presence of multicollinearity between the two supplementary variables. It has been observed that the proposed structure performs better for various correlation coefficient choices, particularly when both supplementary variables consume a fractional effect on the study variable.
To evaluate the performance of the proposed TAHWMA chart against other control existing charts, an algorithm is developed in R software using the Monte Carlo simulation technique to obtain numerical results. Based on numerical results, performance evaluation measures such as ARL, EQL, and RMI are calculated. The analysis based on performance evaluation measures and visual presentation reveals that the proposed charting scheme outperformed existing counterparts. Furthermore, two illustrative examples, one related to the hypothetical dataset and another with regard to a Carbon fibre tube (manufacturing process), are also provided to show the practical implementation of the proposed charting scheme. The proposed charting scheme can be extended for multivariate structures to monitor the process location, dispersion, or both parameters.

Author Contributions

Conceptualization, M.A., S.A.L. and S.A.; methodology, S.A.L. and S.S.; software, N.M.G. and M.A.; validation, S.A., N.M.G. and M.A.; formal analysis, S.A.L. and M.A.; investigation, S.A. and Z.R.; resources, M.A. and S.A.L.; data curation, S.S. and N.M.G.; writing—original draft preparation, M.A. and S.S.; writing—review and editing, Z.R. and S.A.L.; visualization, S.A.L. and Z.R.; supervision, S.A.L. and Z.R.; project administration, S.A.L. and Z.R.; funding acquisition, M.A. and S.A.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are openly available, and a reference is provided.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Montgomery, D.C. Introduction to Statistical Quality Control, 7th ed.; John Wiley & Sons: New York, NY, USA, 2012. [Google Scholar]
  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, 1, 239–250. [Google Scholar] [CrossRef]
  5. Abbas, N. Homogeneously weighted moving average control chart with an application in substrate manufacturing process. Comput. Ind. Eng. 2018, 120, 460–470. [Google Scholar] [CrossRef]
  6. Haq, A.; Khoo, M.B.C. A new double sampling control chart for monitoring process mean using auxiliary information. J. Stat. Comput. Simul. 2017, 88, 869–899. [Google Scholar] [CrossRef]
  7. Mandel, B. The regression control chart. J. Qual. Technol. 1969, 1, 1–9. [Google Scholar] [CrossRef]
  8. Zhang, G. Cause-selecting control charts—A new type of quality control charts. QR J. 1985, 12, 221–225. [Google Scholar]
  9. Riaz, M. Monitoring process mean level using auxiliary information. Stat. Neerl. 2008, 62, 458–481. [Google Scholar] [CrossRef]
  10. Riaz, M. Monitoring process variability using auxiliary information. Comput. Stat. 2008, 23, 253–276. [Google Scholar] [CrossRef]
  11. Riaz, M.; Mehmood, R.; Ahmad, S.; Abbasi, S.A. On the Performance of Auxiliary-based Control Charting under Normality and Nonnormality with Estimation Effects. Qual. Reliab. Eng. Int. 2013, 29, 1165–1179. [Google Scholar] [CrossRef]
  12. Nuriman, M.A.; Mashuri, M.; Ahsan, M. Auxiliary information based generally weighted moving coefficient of variation (AIB-GWMCV) control chart. IOP Conf. Ser. Mater. Sci. Eng. 2021, 1115, 012033. [Google Scholar] [CrossRef]
  13. Haq, A.; Khoo, M.B.C. A new synthetic control chart for monitoring process mean using auxiliary information. J. Stat. Comput. Simul. 2016, 86, 3068–3092. [Google Scholar] [CrossRef]
  14. Abbasi, S.A. Efficient Control Charts for Monitoring Process CV Using Auxiliary Information. IEEE Access 2020, 8, 46176–46192. [Google Scholar] [CrossRef]
  15. Arslan, M.; Ashraf, M.A.; Anwar, S.M.; Rasheed, Z.; Hu, X.; Abbasi, S.A. Novel Mixed EWMA Dual-Crosier CUSUM Mean Charts without and with Auxiliary Information. Math. Probl. Eng. 2022, 1362193. [Google Scholar] [CrossRef]
  16. Abbas, Z.; Nazir, H.Z.; Akhtar, N.; Riaz, M.; Abid, M. On designing a progressive mean chart for efficient monitoring of process location. Qual. Reliab. Eng. Int. 2020, 36, 1716–1730. [Google Scholar] [CrossRef]
  17. Chen, J.H.; Lu, S.L. A New Sum of Squares Exponentially Weighted Moving Average Control Chart Using Auxiliary In-formation. Symmetry 2020, 12, 1888. [Google Scholar] [CrossRef]
  18. Rasheed, Z.; Zhang, H.; Anwar, S.M.; Zaman, B. Homogeneously Mixed Memory Charts with Application in the Substrate Production Process. Math. Probl. Eng. 2021, 2021, 2582210. [Google Scholar] [CrossRef]
  19. Aslam, M.; Anwar, S.M.; Khan, M.; Abiodun, N.L.; Rasheed, Z. Efficient Auxiliary Information–Based Control Charting Schemes for the Process Dispersion with Application of Glass Manufacturing Industry. Math. Probl. Eng. 2022, 1265204. [Google Scholar] [CrossRef]
  20. 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]
  21. Rasheed, Z.; Khan, M.; Abiodun, N.L.; Anwar, S.M.; Khalaf, G.; Abbasi, S.A. Improved Nonparametric Control Chart Based on Ranked Set Sampling with Application of Chemical Data Modelling. Math. Probl. Eng. 2022, 7350204. [Google Scholar] [CrossRef]
  22. Zhang, H.; Rasheed, Z.; Khan, M.; Namangale, J.J.; Anwar, S.M.; Hamid, A. A Distribution-Free THWMA Control Chart under Ranked Set Sampling. Math. Probl. Eng. 2022, 2022, 3823013. [Google Scholar] [CrossRef]
  23. Zichuan, M.; Arslana, M.; Abbasb, Z.; Abbasic, S.A.; Zafar, H. Improving the Performance of EWMA mean chart using Two Auxiliary Variables. Rev. Argent. Clin. Psicol. 2020, 29, 2016–2024. [Google Scholar]
  24. Adegoke, N.A.; Smith, A.N.H.; Anderson, M.J.; Sanusi, R.A.; Pawley, M.D.M. Efficient Homogeneously Weighted Moving Average Chart for Monitoring Process Mean Using an Auxiliary Variable. IEEE Access 2019, 7, 94021–94032. [Google Scholar] [CrossRef]
  25. Cochran, W.G. Sampling Techniques, 3rd ed.; Wiley: New York, NY, USA, 1977. [Google Scholar]
  26. Abbas, N.; Ahmad, S.; Riaz, M. Reintegration of auxiliary information based control charts. Comput. Ind. Eng. 2022, 171, 108479. [Google Scholar] [CrossRef]
  27. Kadilar, C.; Cingi, H. A new estimator using two auxiliary variables. Appl. Math. Comput. 2005, 162, 901–908. [Google Scholar] [CrossRef]
  28. Abbasi, S.A.; Adegoke, N.A. Auxiliary—information—based efficient variability control charts for Phase I of SPC. Qual. Reliabil. Eng. Int. 2020, 36, 2322–2337. [Google Scholar] [CrossRef]
  29. Zhang, G. A new type of control charts and a theory of diagnosis with control charts. In World Quality Congress Transactions; American Society for Quality Control: Milwaukee, WI, USA, 1984; pp. 175–185. [Google Scholar]
Figure 1. The performance of the proposed TAHWMA control chart at fixed ρ y x and various choices of ρ y z .
Figure 1. The performance of the proposed TAHWMA control chart at fixed ρ y x and various choices of ρ y z .
Symmetry 15 00482 g001
Figure 2. The performance of the proposed TAHWMA control chart at several choices of ρ y x and fixed ρ y z .
Figure 2. The performance of the proposed TAHWMA control chart at several choices of ρ y x and fixed ρ y z .
Symmetry 15 00482 g002
Figure 3. The OOC performance of the proposed TAHWMA and existing control charts.
Figure 3. The OOC performance of the proposed TAHWMA and existing control charts.
Symmetry 15 00482 g003
Figure 4. The real-life application of the AHWMA control chart when ( ρ Y X = 0.25 ,   L = 2.608 ).
Figure 4. The real-life application of the AHWMA control chart when ( ρ Y X = 0.25 ,   L = 2.608 ).
Symmetry 15 00482 g004
Figure 5. The real-life application of the HWMA control chart when (L = 2.608).
Figure 5. The real-life application of the HWMA control chart when (L = 2.608).
Symmetry 15 00482 g005
Figure 6. The real-life application of the proposed TAHWMA control chart when ( λ = 0.05; L = 2.608).
Figure 6. The real-life application of the proposed TAHWMA control chart when ( λ = 0.05; L = 2.608).
Symmetry 15 00482 g006
Figure 7. The real-life application of the proposed TAHWMA control chart without multicollinearity with design parameters L = 2.608 ,   λ = 0.05 ,     ρ Y X = 0.25 ,   ρ Y Z = 0.50   and   ρ X Z = 0 .
Figure 7. The real-life application of the proposed TAHWMA control chart without multicollinearity with design parameters L = 2.608 ,   λ = 0.05 ,     ρ Y X = 0.25 ,   ρ Y Z = 0.50   and   ρ X Z = 0 .
Symmetry 15 00482 g007
Figure 8. The real-life application of the proposed TAHWMA control chart with multicollinearity with designed parameters L = 2.633 ,   λ = 0.05 ,     ρ Y X = 0.25 ,   ρ Y Z = 0.50   and   ρ X Z = 0.05 .
Figure 8. The real-life application of the proposed TAHWMA control chart with multicollinearity with designed parameters L = 2.633 ,   λ = 0.05 ,     ρ Y X = 0.25 ,   ρ Y Z = 0.50   and   ρ X Z = 0.05 .
Symmetry 15 00482 g008
Figure 9. The real-life of the HWMA control chart design parameters ( L = 2.272 ,   λ = 0.05 ).
Figure 9. The real-life of the HWMA control chart design parameters ( L = 2.272 ,   λ = 0.05 ).
Symmetry 15 00482 g009
Figure 10. The real-life application of the AHWMA control chart with design parameters ( L = 2.272 ,   λ = 0.05 ,     ρ Y X = 0.25 ).
Figure 10. The real-life application of the AHWMA control chart with design parameters ( L = 2.272 ,   λ = 0.05 ,     ρ Y X = 0.25 ).
Symmetry 15 00482 g010
Table 1. ARL values of the proposed TAHWMA chart in absence of multicollinearity at various choices of design parameters.
Table 1. ARL values of the proposed TAHWMA chart in absence of multicollinearity at various choices of design parameters.
δ Small ShiftsModerate ShiftsLarge Shifts
λ = 0.03, L = 2.272 λ = 0.05, L = 2.608 λ = 0.10, L = 2.938 λ = 0.25, L = 3.075
ρ y x = 0.25 ρ y x = 0.5 ρ y x = 0.75 ρ y x = 0.25 ρ y x = 0.5 ρ y x = 0.75 ρ y x = 0.25 ρ y x = 0.5 ρ y x = 0.75 ρ y x = 0.25 ρ y x = 0.5 ρ y x = 0.75
ρ y z = 0.50, ρ x z = 0 ρ y z = 0.50, ρ x z = 0 ρ y z = 0.50, ρ x z = 0 ρ y z = 0.50, ρ x z = 0
0499.92500.59500.50500.87500.48500.00499.70499.68500.22500.96499.81500.92
0.03422.12393.17291.19432.70408.11315.37439.43421.68335.05466.80455.76398.38
0.05330.14288.40174.72345.68310.76194.21362.74330.82218.70419.96392.24290.20
0.075226.50193.27104.73254.06217.28120.66272.73237.35136.16346.73310.04189.49
0.1165.07136.6969.23186.99157.0681.87208.11173.6591.86274.62239.11127.87
0.125125.40100.8348.69144.15118.1059.37160.82132.5366.29223.38184.6589.11
0.17579.3662.0428.5192.9374.4135.17103.9283.4539.74145.17114.7250.07
0.264.2250.3222.9276.3260.7028.6985.9668.5832.40117.8792.2839.13
0.2545.8235.6915.8055.4043.7619.8162.2948.8422.8182.9863.0226.01
0.514.7411.355.3418.5314.136.4821.1916.437.4424.2018.197.52
0.757.696.093.089.407.333.6111.028.574.0811.538.763.95
15.014.022.116.004.752.466.885.442.717.095.442.61
1.52.892.361.213.402.761.353.813.081.503.682.931.45
21.971.581.022.311.841.032.582.071.062.441.951.06
Table 2. ARL values of proposed TAHWMA control chart under the presence of multicollinearity when ρ y x = 0.25 and ρ y z = 0.50.
Table 2. ARL values of proposed TAHWMA control chart under the presence of multicollinearity when ρ y x = 0.25 and ρ y z = 0.50.
δ ρ x z = 0.05 ρ x z = 0.15 ρ x z = 0.25
λ = 0.05, L = 2.633 λ = 0.05, L = 2.679 λ = 0.05, L = 2.719
ARLSDRLMDRLARLSDRLMDRLARLSDRLMDRL
0501.18372.02438499.75372.26437495.53368.77436
0.03434.94329.12376434.31331.05373429.45329.72367
0.05349.49274.23291353.02278.69295354.30279.99294
0.075253.89200.40210261.87204.59216262.13205.74218
0.1190.50147.65160192.57150.72161196.71154.48164
0.125145.68110.16123150.41114.68127151.95116.03127
0.17594.2469.618196.5271.678298.1073.0983
0.278.0656.576779.2557.856881.4759.9569
0.2556.4240.244957.6241.415058.6841.9551
0.518.5912.441619.2412.931719.5813.2717
0.759.535.9589.876.20910.086.339
16.173.4966.233.5166.413.676
1.53.441.7433.521.7733.601.853
22.341.2632.401.2832.441.313
Table 3. The ARL values of the EWMA and AEWMA control charts.
Table 3. The ARL values of the EWMA and AEWMA control charts.
δ EWMAAEWMA ( ρ y x = 0.25)
λ = 0.03
L = 2.483
λ = 0.05
L = 2.639
λ = 0.1
L = 2.824
λ = 0.25
L = 3.001
λ = 0.03
L = 2.483
λ = 0.05
L = 2.639
λ = 0.1
L = 2.824
λ = 0.25
L = 3.001
0500.64499.68500.94499.79500.33500.36500.39500.78
0.03456.41462.91480.81488.19453.72466.58475.58484.32
0.05388.18412.31438.70466.30382.72408.43438.29462.55
0.075304.78338.65376.08435.67301.97328.86372.44427.63
0.1232.24266.78318.91390.56229.75262.54314.47384.82
0.125181.85211.65264.66346.91177.55202.55251.93336.03
0.175113.49135.20177.02262.29110.97130.36168.29254.75
0.293.56111.53148.70224.7790.34106.50139.50219.07
0.2567.2976.75102.97169.3363.2673.5897.39161.22
0.521.1923.7429.1247.8520.2522.2626.8144.31
0.7510.7711.9313.5819.2710.1411.1812.7818.02
16.607.408.2510.386.346.927.719.76
1.53.423.774.174.793.293.593.984.49
22.252.412.652.932.152.312.522.80
Table 4. The ARL values of the HWMA and AHWMA control charts.
Table 4. The ARL values of the HWMA and AHWMA control charts.
δ HWMAAHWMA ( ρ y x = 0.25)
λ = 0.03
L = 2.272
λ = 0.05
L = 2.608
λ = 0.1
L = 2.938
λ = 0.25
L = 3.075
λ = 0.03
L = 2.272
λ = 0.05
L = 2.608
λ = 0.1
L = 2.938
λ = 0.25
L = 3.075
0500.70499.35499.48499.69499.98500.18500.96500.34
0.03440.12449.24456.35473.57442.47448.12453.61473.88
0.05359.03382.47397.11441.31359.39380.54393.94445.27
0.075274.16298.00313.35382.72266.64286.42310.21377.80
0.1205.43229.66250.83329.77199.38222.00242.77317.69
0.125158.89179.95201.22270.99151.54174.35191.82265.10
0.175101.43119.55133.63187.6297.86115.19127.79179.99
0.284.5299.66111.81158.7680.3296.17106.49148.73
0.2561.1973.0681.59113.0158.5969.1578.16106.86
0.520.0825.2628.5733.9619.0923.7627.2731.84
0.7510.3612.8014.8916.129.8212.2514.0715.24
16.647.999.379.746.287.638.819.14
1.53.724.414.974.923.564.224.754.66
22.553.003.333.182.442.853.183.05
Table 5. ARL comparisons of the proposed TAHWMA and the existing control chart.
Table 5. ARL comparisons of the proposed TAHWMA and the existing control chart.
TAHWMAEWMAAEWMAHWMAAHWMA
Shift λ = 0.03 λ = 0.05 λ = 0.1 λ = 0.25 λ = 0.03 λ = 0.03 λ = 0.03 λ = 0.03
0500.5500500.22500.92500.64500.33500.7499.98
0.03291.19315.37335.05398.38456.41453.72440.12442.47
0.05174.72194.21218.7290.2388.18382.72359.03359.39
0.075104.73120.66136.16189.49304.78301.97274.16266.64
0.169.2381.8791.86127.87232.24229.75205.43199.38
0.12548.6959.3766.2989.11181.85177.55158.89151.54
0.17528.5135.1739.7450.07113.49110.97101.4397.86
0.222.9228.6932.439.1393.5690.3484.5280.32
0.2515.819.8122.8126.0167.2963.2661.1958.59
0.55.346.487.447.5221.1920.2520.0819.09
0.753.083.614.083.9510.7710.1410.369.82
12.112.462.712.616.66.346.646.28
1.51.211.351.51.453.423.293.723.56
21.021.031.061.062.252.152.552.44
EQL2.122.372.612.606.286.016.486.18
RMI0.000.170.300.522.212.102.031.91
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

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. https://doi.org/10.3390/sym15020482

AMA Style

Arslan M, Anwar S, Gunaime NM, Shahab S, Lone SA, Rasheed Z. An Improved Charting Scheme to Monitor the Process Mean Using Two Supplementary Variables. Symmetry. 2023; 15(2):482. https://doi.org/10.3390/sym15020482

Chicago/Turabian Style

Arslan, Muhammad, Sadia Anwar, Nevine M. Gunaime, Sana Shahab, Showkat Ahmad Lone, and Zahid Rasheed. 2023. "An Improved Charting Scheme to Monitor the Process Mean Using Two Supplementary Variables" Symmetry 15, no. 2: 482. https://doi.org/10.3390/sym15020482

APA Style

Arslan, M., Anwar, S., Gunaime, N. M., Shahab, S., Lone, S. A., & Rasheed, Z. (2023). An Improved Charting Scheme to Monitor the Process Mean Using Two Supplementary Variables. Symmetry, 15(2), 482. https://doi.org/10.3390/sym15020482

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