Abstract

More recently in statistical quality control studies, researchers are paying more attention to quality characteristics having nonnormal distributions. In the present article, a generalized multiple dependent state (GMDS) sampling control chart is proposed based on the transformation of gamma quality characteristics into a normal distribution. The parameters for the proposed control charts are obtained using in-control average run length (ARL) at specified shape parametric values for different specified average run lengths. The out-of-control ARL of the proposed gamma control chart using GMDS sampling is explored using simulation for various shift size changes in scale parameters to study the performance of the control chart. The proposed gamma control chart performs better than the existing multiple dependent state sampling (MDS) based on gamma distribution and traditional Shewhart control charts in terms of average run lengths. A case study with real-life data from ICU intake to death caused by COVID-19 has been incorporated for the realistic handling of the proposed control chart design.

1. Introduction

One of the important techniques for improving manufactured product quality and for reducing the manufacturing costs is statistical quality control (SQC). Since the pioneer work by Shewhart A. Walter during 1920s in Bell Telephone Laboratories, wide varieties of control chart techniques have been constructed and extensively implemented in SQC. The main feature of control charting is to identify the amount of assignable cause(s) and hence rectify it by taking necessary action on the production process before sending the outcome of the products into the market. This control charting helps to avoid nonconforming products from being manufactured by the company. More details about Shewhart control charts can be seen in Montgomery [1].

Usually, control charts are being designed and operating under the assumption of the normality for the variable of interest. Nevertheless, these assumptions may not be true for various realistic situations and other distributions away from normality had been considered and discussed by many authors in the literature (e.g., see [25]). The waiting time of an event, for example, can be represented by a gamma distribution as in [6]. Numerous researchers concentrate on quality characteristic understudy which follows a nonnormal distribution or transformed into normality to apply Shewhart type control charts. For skewed data, the gamma distribution is widely used. The works on the control charts for the gamma distribution are presented by Al-Oraini and Rahim [7], Jearkpaporn et al. [8], Sheu and Lin [9], Aslam et al. [10], and Zhang et al. [11]. Santiago and Smith [5] used transformation given by Johnson and Kotz [12] and Nelson [13]. Mohammed [14], Mohammed and Laney [15], and Aslam et al. [16] discussed the application of the -chart.

Several researchers have developed diversified sampling designs to obtain more efficient control charts. Recently, researchers focused on multiple dependent state (MDS) sampling in the creation of a control chart. Wortham and Baker [17] proposed the MDS sampling in quality control charts. MDS design is more competent than the existing single sampling plans because it considers the previous lot information along with the current lot to make a decision whether the process is under control or not (see [18]). Aslam et al. [19] developed a control chart for gamma distribution using MDS sampling. The control chart scheme using MDS sampling was studied by different authors for various schemes (see [16, 19, 2029, 30, 31]).

More recently, Raza and Aslam [32], Rao et al. [33], Rao et al. [34], and Aslam et al. [35, 36] formulated generalized MDS (GMDS) sampling for various schemes. GMDS is more flexible and efficient than MDS to design the control chart using the gamma distribution. The aim of this article is to construct a gamma control chart for monitoring the process mean based on GMDS sampling design. The application of the proposed chart will be given using the COVID-19 mortality data. It is expected that the proposed chart will perform better than the existing Shewhart control chart and control charts using MDS in terms of average run length and standard deviation of run length.

2. Design of Control Chart for Gamma Distribution Based on GMDS Sampling

The proposed control chart for a gamma distribution using gamma to normal transformation is discussed. Let be a random variable from a gamma distribution with shape parameter and scale parameter . The cumulative distribution function (cdf) of the gamma distribution is given by

Wilson and Hilferty [37] recommended that if follows a gamma distribution with specific parameters, then the transformed variable can be distributed approximately as normal with mean and variance , where

The proposed gamma control chart using GMDS sampling comprises the two pairs of control chart limits. The inner lower control limit (LCL) and upper control limit (UCL) are denoted by subscript 1, and the outer lower control limit (LCL) and upper control limit (UCL) are denoted by subscript 2. The four control limits are given bywhere and are the chart constants to be found when the in-control ARL is approximately equal to preassigned value . The convenient form of the above control limits is given as follows: , , , and , where

The operation of the proposed control chart using GMDS scheme is described as follows:(1)Obtain quality measurement from the manufacturing process, and denote the quality characteristic by . Compute the transformed variable as (2)The process can be considered under control if , and the process can be considered out-of-control if . Or else, go to Step 3(3)The process can be considered under control whenever out of proceeding subgroups have been declared as under control, that is, ; otherwise, the output of the product can be considered out-of-control and go back to Step 1

The probability of declaring as in-control for the proposed control chart when the process is actually in-control is given as follows:where

Therefore, the in-control average run length (ARL) when the process is under control is given by

Assume the gamma scale parameter has been changed from to , where is the shift value.

The probability of process is declared as in-control while the scale parameter which has been changed can be obtained as follows:where

The out-of-control average run length (ARL) when the process is out-of-control is given as

The proposed control chart parameters and along with ARL1 are obtained using the following algorithm:(1)Decide the predetermined in-control ARL as (2)Fix the known values for and (3)Obtain the ARL0 using Equation (14), which consists of chart parameters and (4)Determine the most possible values of chart parameters and such that, (5)In the above step, we get more values of and to satisfy the condition. Choose the best values of and for which the value of is almost equal to (6)Using the best parametric values of and determined in the previous step, work out the using Equation (18) and hence obtain standard deviation (SD) of run-length (SDRL) for various shift () values

The R codes to find the design parameters of the control chart are given in the appendix.

3. Numerical Results and Discussion

The performance of the proposed gamma control chart using GMDS sampling is considered based on ARL, such as ARL0 and ARL1. These ARL values are used to know the effectiveness of the developed control chart. The developed chart is said to be efficient if it shows larger in-control ARL and smaller out-of-control ARL. Using the aforementioned algorithm in Section 2, the chart coefficients and are obtained. The out-of-control ARLs and SDRL are computed for a choice of shift values, from 1.0 to 2.0 with an interval of 0.1 and 2.0 to 4.0 with an interval of 0.5. The values of considered are 4, 5, and 6 and , 10, and 20. Table 1 is for and , Table 2 is for and , Table 3 is for and , Table 4 is for and , Table 5 is for and , and Table 6 is for and .

We pointed out the following several noteworthy comments from Tables 16 for the developed control charts:(1)The out-of-control ARL and SDRL values decline speedily when the shift () of the manufacturing process increases(2)It is detected that the chart coefficient shows an increasing tendency for increased value of for a fixed value of when other parametric combinations are fixed(3)From the tables, it is noticed that ARL1 and SDRL values decrease when values increase. In addition, ARL1 and SDRL values increased with the increase of value (i.e., -2 to -0). It also observed the same inclination over the other parametric combinations and and 500(4)It is interesting to observe from the results that the values of ARL1 and SDRL are small for -2 and these values are increasing from -2 to for fixed values of . In addition, noticed that ARL1 and SDRL values are large at as compared to the values at -1 and -2 (we know that if , the developed plan becomes MDS design). Hence, it is concluded from the results that gamma control chart using GMDS sampling is an enormous amount of accurate than gamma control chart using MDS sampling

4. Comparison with Existing Charts

In this part, a comparison is made between the developed control chart and the existing Shewhart type control chat and MDS control chart for gamma distribution. Also, the application of developed control chart and its dominance over available control chart schemes studied using real data set is presented. In addition, through a simulation study, the supremacy of the developed control chart when compared with the existing control charts is examined. The performance of the developed control chart is studied through ARL values and we know that a control chart with smaller ARL values is more desirable. In this investigation, we studied when and ; the shape parameter of gamma distribution is given as , 10, and 20 to compare the developed gamma control chart under GMDS with the existing MDS and Shewhart type control chart at various shift values. These comparisons are presented in Table 7 for and and in Table 8 for and at various shape parameters of the gamma distribution.

It is noticed that from the results on the basis of Tables 7 and 8, the developed gamma control charts show smaller quantity ARL1 values as compared with the MDS and Shewhart type control charts at various shifts () values and various parametric values studied in this article. At a glance, when , and from Table 7, for the developed control chart whereas for MDS scheme and from the Shewhart type control chart. Similarly, for , , and from Table 8, we sense that the developed control chart gives while for the MDS control chart and from the Shewhart type control chart. The graphical presentation is given to show the performance of developed control chat over the existing MDS and Shewhart type control charts along with various shift values (see Figures 1 and 2). From these two figures, it is articulated that the developed gamma control chart based on GMDS is certified extra sensitive as compared to the MDS and the Shewhart-type control charts. To draw attention to this conclusion, a real data illustration and a simulation study are also carried out in the following subsections.

4.1. Simulation Analysis

In order to investigate the implementation of the planned control chart over the available control charts, a simulation study is conducted. In this investigation, 30 samples are generated from the gamma distribution with shape parameter and in-control scale parameter 1 and last 30 random samples are generated from a gamma distribution with shape parameter and out-of-control scale parameter 1.4 (i.e., the shift of ). The data is reported in Table 9 alongside computed statistic . The control chart coefficients at , , and are available in Table 8. The Shewhart type gamma control chart is given in Figure 3, and MDS gamma control chart when and is provided in Figure 4. The gamma control chart using GMDS sampling when and is depicted in Figure 5. According to the above scheme for gamma control charts, we implement the MDS chart as follows: if previous 5 (since ) values are displayed between the inner control limits, then the process is considered to be under control while for proposed gamma control chart under GMDS, the process is said to be declared as under control if no less than 3 out of 5 previous (since and ) values are within the interior control limits.

From Figures 3 and 4, it is apparent that the gamma control charts based on Shewhart type and MDS scheme are unable to notice the shift. On the other hand, in Figure 5, it can be found that using gamma control charts under the GMDS scheme detects out-of-control at sample numbers 34, 40, 41, 44, 45, 50, and 59. Present simulation examines that the developed gamma control chart based on GMDS sampling is more efficient than the gamma control charts as compared to MDS and Shewhart type design.

4.2. Application of the Proposed Chart for COVID-19 Data

In the present section, the developed gamma control chart using GMDS sampling is applied to the monitoring of Coronavirus (COVID-19) outbreak in China. The data set is borrowed from Li et al. [38], and they discussed mortality caused by COVID-19 based on collected data of 33 death cases in Wuhan city of Hubei province during an early outbreak as well as confirmed cases and death toll. They have studied the COVID-19 outbreak in China by considering some specific regions as representative samples. The days from ICU intake to death caused by COVID-19 data are presented in Table 10 along with the transformation of the variable. To fit the gamma distribution, the parameters are estimated using the maximum likelihood approach, and the shape is 2.00262.0, and the scale parameter is 3.9185. The Kolmogorov-Smirnov test is 0.1197, and value is 0.7322. Figure 6 displays the histogram and Q-Q plot to highlight the goodness of fit of the gamma distribution. Hence, gamma distribution furnishes a good fit for the days from ICU intake to death caused by COVID-19 data.

The chart constants at the estimated shape parameter of 2 are obtained using a simulation procedure given in Section 2. The control chart constants at and are and for a developed control chart, and for MDS control chart, and for Shewhart type control chart. The chart limits of the developed control chart, MDS, and Shewhart type control charts for days from ICU intake to death data are given in Figures 79. In Figures 79, the proposed control charts for GMDS, MDS, and Shewhart type control charts are displayed by plotting the control limits and chat statistics. Using days from ICU intake to death data, the developed control chart can be exemplified as follows: declare the process as in-control when 4 earlier values of fall in the inner control limits of the developed gamma control chart using the MDS plan. While in the case of the developed gamma control chart using GMDS scheme, the process can be expected as in-control when at least 2 out of 4 earlier values of fall in the inner control limits.

It is clearly noticeable from Figure 7 that all sample statistics are inside the upper and lower control limits. Hence, it shows that Shewhart-type gamma control chart fails to detect a change in the process. On the other hand, if we consider gamma control chart, using MDS design also failed to identify changes in the process due to all sample statistics inside the inner upper and inner lower control limits (see Figure 8), whereas the gamma control chart based on GMDS design experiences the out-of-control signals at sample numbers 10, 13, 23, 31, and 33, respectively (see Figure 9). Therefore, from the methodology explained in Section 2, it reveals that the developed gamma control chart based on GMDS design is faster in detecting process variation as compared with the existing control chart based MDS and Shewhart type. Hence, from real application in the field of medical sciences, data for days from ICU intake to death revels that the gamma control chart based on GMDS remains superior methodology as compared to the existing control charts considered in this study.

5. Conclusions

We developed the gamma control chart based on GMDS sampling. The computational methodology is also discussed for ARL and SDRL when the process is in-control and out-of-control. The control chart parameters for the proposed control chart are obtained using in-control average run length (ARL0) and out-of-control ARL (ARL1) of the proposed gamma control chart using GMDS sampling. The performance of the proposed control chart is investigated using simulation for the various shifts in the scale parameter. Tables of chart parameters alongside out-of-control ARL (ARL1) for various shift values for specified shape parameters are displayed. Furthermore, a comparative study is also carried out with the developed gamma control chart based on GMDS sampling over the existing MDS and Shewhart type gamma control charts using the ARLs. The results display that the developed gamma control chart based on GMDS sampling shows reduced ARL1 values as compared with the existing two gamma control chart discussed in this article.

The implementation of the developed control chart is demonstrated using simulation study as well as tangible data from ICU intake to death cause by COVID-19 and it is shown that designed gamma control chart using GMDS sampling detected out-of-control samples whereas the MDS and Shewhart type gamma control charts failed to detect the out-of-control signal. Hence, we conclude that the gamma control chart based on GMDS is a superior methodology as compared to the existing control charts considered in this study to detect a shift in the parameter. The developed control chart method in this paper can be used in different industrial and medical situations specifically when the researcher would like to discover a small and moderate shift in quality characteristics. Future research maybe considered as control charts for some nonnormal distributions and cost consideration using GMDS sampling design.

Appendix

R Coding to Run the Algorithm

## compute mean and standard deviation

mts<-b^(1/3)gamma(a+1/3)/gamma(a)

vts<-b^(2/3)gamma(a+2/3)/gamma(a)-mts^2

sdts<-sqrt(vts)

## fixing the range of k1 and k2

q1<-seq(2.61,6.96, by =0.0005)

 q2<-seq(1.21,6.86, by =0.0005)

vv<-gamma(a+2/3)/gamma(a)-(gamma(a+1/3)/gamma(a))^2

   ll1<-gamma(a+1/3)/gamma(a)-k1sqrt(vv)

   ll2<-gamma(a+1/3)/gamma(a)-k2sqrt(vv)

   u1<-gamma(a+1/3)/gamma(a)+k1sqrt(vv)

   u2<-gamma(a+1/3)/gamma(a)+k2sqrt(vv)

   lcl1<-max(0,b^(1/3)ll1)

   ucl1<-b^(1/3)ul1

   lcl2<-max(0,b^(1/3)ll2)

   ucl2<-b^(1/3)ul2

pa<-sum(dpois(0:(a-1),ll2^3))-sum(dpois(0:(a-1),ul2^3))

ps<-sum(dpois(0:(a-1),ll1^3))-sum(dpois(0:(a-1),ll2^3))+sum(dpois(0:(a-1), ul2^3))

      -sum(dpois(0:(a-1),ul1^3))

s1<-pbinom(m, m, pa, lower.tail = TRUE, log.p = FALSE)

      -pbinom(k-1, m, pa, lower.tail = TRUE, log.p = FALSE)

   Pin0<-pa+pss1

   ARL0=1.0/(1.0-Pin0)

   if ((ARL0>=370) && (ARL0<=370.5)) then choose those k1 and k2 values

Data Availability

The data is given in the paper.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This article was funded by the Deanship of Scientific Research (DSR) at King Abdulaziz University, Jeddah. The authors therefore acknowledge DSR with thanks for the technical and financial support. The work by Chi-Hyuck Jun was supported by Korea Institute for Advancement of Technology (KIAT) grant funded by the Korea Government (MSIT) (N0008691, the Competency Development Program for Industry Specialist).