# Monitoring Mortality Caused by COVID-19 Using Gamma-Distributed Variables Based on Generalized Multiple Dependent State Sampling

Authors: Muhammad Aslam; G. Srinivasa Rao; Muhammad Saleem; Rehan Ahmad Khan Sherwani; Chi-Hyuck Jun
Journal: Computational and Mathematical Methods in Medicine
 (2021)

---

## 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.


---

## Body

## 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 [2–5]). 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 t-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, 20–29, 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. LetX 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(1)pX≤x=1−∑k=0α−1e−x/βx/βkk!.Wilson and Hilferty [37] recommended that if X follows a gamma distribution with specific parameters, then the transformed variable X∗=X1/3 can be distributed approximately as normal with mean μX∗ and variance σX∗2, where(2)μX∗=β1/3Γα+1/3Γα,σX∗2=β2/3Γα+2/3Γα−μX∗2.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 by(3)UCL1=μX∗+k1σX∗=β1/3Γα+1/3Γα+k1β2/3Γα+2/3Γα−μX∗2,(4)LCL1=μX∗−k1σX∗=β1/3Γα+1/3Γα−k1β2/3Γα+2/3Γα−μX∗2,(5)UCL2=μX∗+k2σX∗=β1/3Γα+1/3Γα+k2β2/3Γα+2/3Γα−μX∗2,(6)LCL2=μX∗−k2σX∗=β1/3Γα+1/3Γα−k2β2/3Γα+2/3Γα−μX∗2,where k1 and k2 are the chart constants to be found when the in-control ARL is approximately equal to preassigned value r0. The convenient form of the above control limits is given as follows: UCL1=β1/3UL1, UCL2=β1/3UL2, LCL1=β1/3LL1, and LCL2=β1/3LL2, where(7)LL1=Γα+1/3Γα−k1Γα+2/3Γα−Γα+1/3Γα2,(8)LL2=Γα+1/3Γα−k2Γα+2/3Γα−Γα+1/3Γα2,(9)UL1=Γα+1/3Γα+k1Γα+2/3Γα−Γα+1/3Γα2,(10)UL2=Γα+1/3Γα+k2Γα+2/3Γα−Γα+1/3Γα2.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 byX. Compute the transformed variable X∗ as X∗=X1/3(2)
The process can be considered under control ifLCL2≤X∗≤UCL2, and the process can be considered out-of-control if X∗≥UCL1orX∗≤LCL1. Or else, go to Step 3(3)
The process can be considered under control wheneverk out of m proceeding subgroups have been declared as under control, that is, LCL2≤X∗≤UCL2; otherwise, the output of the product can be considered out-of-control and go back to Step 1The probability of declaring as in-control for the proposed control chart when the process is actually in-control is given as follows:(11)Pin.0=Pa.0+Ps.0∑j=kmmjPa.0j1−Pa.0m−j,where(12)Pa.0=pLCL2≤X∗≤UCL2∣β=β0=pX∗≤UCL2∣β=β0−pX∗≤LCL2∣β=β0=∑k=0α−1e−LL23LL23kk!−∑k=0α−1e−UL23UL23kk!,(13)Ps.0=pLCL1≤X∗≤LCL2∣β=β0+pUCL2≤X∗≤UCL1∣β=β0=∑k=0α−1e−LL13LL13kk!−∑k=0α−1e−LL23LL23kk!+∑k=0α−1e−UL23UL23kk!−∑k=0α−1e−UL13UL13kk!.Therefore, the in-control average run length (ARL) when the process is under control is given by(14)ARL0=11−Pin.0.Assume the gamma scale parameter has been changed fromβ=β0 to β=β1=sβ0, where s 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:(15)Pin.1=Pa.1+Ps.1∑j=kmmjPa.1j1−Pa.1m−j,where(16)Pa.1=pLCL2≤X∗≤UCL2∣β=β1=pX∗≤UCL2∣β=β1−pX∗≤LCL2∣β=β1=∑k=0α−1e−LL23/sLL23/skk!−∑k=0α−1e−UL23/sUL23/skk!,(17)Ps.1=pLCL1≤X∗≤LCL2∣β=β1+pUCL2≤X∗≤UCL1∣β=β1=∑k=0α−1e−LL13/sLL13/skk!−∑k=0α−1e−LL23/sLL23/skk!+∑k=0α−1e−UL23/sUL23/skk!−∑k=0α−1e−UL13/sUL13/skk!.The out-of-control average run length (ARL) when the process is out-of-control is given as(18)ARL1=11−Pin.1.The proposed control chart parametersk1 and k2 along with ARL1 are obtained using the following algorithm:(1)
Decide the predetermined in-control ARL asr0(2)
Fix the known values form and k(3)
Obtain the ARL0 using Equation (14), which consists of chart parameters k1 and k2(4)
Determine the most possible values of chart parametersk1 and k2 such that, ARL0≥r0(5)
In the above step, we get more values ofk1 and k2 to satisfy the condition. Choose the best values of k1 and k2 for which the value of ARL0 is almost equal to r0(6)
Using the best parametric values ofk1 and k2 determined in the previous step, work out the ARL1 using Equation (18) and hence obtain standard deviation (SD) of run-length (SDRL) for various shift (s) valuesThe 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 k1 and k2 are obtained. The out-of-control ARLs and SDRL are computed for a choice of shift values, s 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 m considered are 4, 5, and 6 and α0=5, 10, and 20. Table 1 is for r0=370 and α0=5, Table 2 is for r0=370 and α0=10, Table 3 is for r0=370 and α0=20, Table 4 is for r0=500 and α0=5, Table 5 is for r0=500 and α0=10, and Table 6 is for r0=500 and α0=20.Table 1
ARLs and SDRLs of proposed control chart for gamma distribution whenα0=5 and ARL0=370.
sk1=3.1125k1=3.0045k1=3.0025k1=3.0240k1=3.0045k1=3.0060k1=3.0615k1=3.0215k1=3.0225k2=1.5025k2=2.0305k2=2.5235k2=1.6985k2=2.0955k2=2.5455k2=1.7260k2=2.1075k2=2.5325m=4, k=2m=4, k=3m=4, k=4m=5, k=3m=5, k=4m=5, k=5m=6, k=4m=6, k=5m=6, k=6ARLSDRLARLSDRLARLSDRLARLSDRLARLSDRLARLSDRLARLSDRLARLSDRLARLSDRL1.0370.05369.55370.01369.51370.02369.52370.00369.50370.01369.51370.04369.54370.04369.54370.03369.53370.02369.521.1206.61206.11207.63207.13208.25207.75205.73205.23206.33205.83207.32206.82200.51200.01203.00202.50204.91204.411.2100.3699.86102.57102.07104.12103.6299.9399.43100.88100.38103.04102.5493.7793.2796.9796.47100.3799.871.350.2349.7352.5152.0154.5654.0649.9649.4550.9950.4953.6853.1845.0744.5647.8747.3651.6251.121.427.1726.6729.0028.5031.1030.6026.8626.3627.8227.3130.4629.9623.4822.9725.6325.1329.0328.531.516.0015.5017.3716.8619.2718.7715.6815.1716.4915.9918.8318.3313.4412.9315.0414.5317.8617.351.610.219.7011.2210.7112.8712.369.919.4010.6010.0812.5612.058.437.929.639.1211.9111.391.77.006.487.767.249.158.646.746.227.326.808.958.435.755.236.676.158.507.981.85.114.585.695.176.876.354.904.375.384.856.736.214.223.694.944.416.425.901.93.933.394.403.865.394.873.763.224.173.645.314.783.292.743.873.345.094.562.03.162.613.543.004.403.873.032.483.392.844.343.812.692.133.182.634.183.652.51.661.041.861.262.301.731.631.011.841.242.301.731.560.931.801.212.271.703.01.280.601.410.761.671.061.290.611.420.771.681.071.270.591.420.771.671.053.51.140.401.230.541.400.741.160.431.250.551.400.751.160.431.250.561.400.754.01.080.301.150.411.250.561.100.331.160.421.260.571.100.331.160.431.260.57Table 2
ARLs and SDRLs of proposed control chart for gamma distribution whenα0=10 and ARL0=370.
sk1=3.0575k1=3.0145k1=3.0135k1=3.1515k1=3.0785k1=3.0765k1=3.1775k1=3.1225k1=3.1145k2=1.5790k2=2.0750k2=2.5765k2=1.6005k2=2.0085k2=2.4645k2=1.6655k2=2.0190k2=2.4505m=4,k=2m=4,k=3m=4,k=4m=5,k=3m=5,k=4m=5,k=5m=6,k=4m=6,k=5m=6,k=6ARLSDRLARLSDRLARLSDRLARLSDRLARLSDRLARLSDRLARLSDRLARLSDRLARLSDRL1.0370.01369.51370.02369.52370.02369.52370.01369.51370.00369.50370.05369.55370.01369.51370.04369.54370.01369.511.1176.21175.71178.51178.01179.79179.29167.40166.90168.84168.33171.65171.15161.99161.49163.28162.77167.95167.451.266.3265.8268.9968.4971.5071.0057.8457.3460.0259.5264.2263.7153.3952.8955.4654.9661.1960.691.327.0126.5128.9528.4531.6131.1121.8921.3823.5623.0527.1126.6119.5119.0021.1320.6325.4224.921.412.6512.1413.9113.4016.1215.619.899.3710.9810.4713.5313.028.708.199.799.2812.6412.131.56.866.347.677.159.378.855.354.826.085.567.857.334.754.225.484.967.376.851.64.253.724.804.276.085.563.392.843.903.365.164.633.062.523.593.054.904.371.72.942.393.342.804.323.792.431.862.812.263.733.192.251.682.652.093.593.051.82.231.662.541.983.302.751.911.322.211.642.912.361.821.222.131.552.832.271.91.821.222.071.482.662.101.621.001.861.262.401.831.570.941.811.212.351.782.01.560.941.771.162.251.671.430.791.631.022.061.481.410.761.610.992.031.442.51.130.381.220.511.400.751.110.351.190.481.350.681.120.371.200.481.340.683.01.040.201.080.301.160.431.040.201.070.281.140.401.040.211.080.281.140.393.51.020.121.030.191.070.281.020.131.030.181.060.251.020.131.030.181.060.254.01.010.081.020.121.030.181.010.081.010.121.030.171.010.091.010.121.030.17Table 3
ARLs and SDRLs of proposed control chart for gamma distribution whenα0=20 and ARL0=370.
sk1=3.2715k1=4.1315k1=4.4505k1=3.3265k1=4.5695k1=4.7925k1=3.6385k1=4.6190k1=4.5345k2=1.4485k2=1.7555k2=2.2170k2=1.5435k2=1.8270k2=2.2580k2=1.5880k2=1.8835k2=2.2915m=4,k=2m=4,k=3m=4,k=4m=5,k=3m=5,k=4m=5,k=5m=6,k=4m=6,k=5m=6,k=6ARLSDRLARLSDRLARLSDRLARLSDRLARLSDRLARLSDRLARLSDRLARLSDRLARLSDRL1.0370.02369.52370.04369.54370.05369.55370.01369.51370.02369.52370.05369.55370.04369.54370.05369.55370.05369.551.1119.55119.05120.53120.02121.50120.99113.23112.73114.29113.79118.35117.84108.88108.38109.14108.64115.61115.111.225.9225.4126.7426.2429.4428.9322.7622.2624.2623.7628.1827.6721.5721.0722.4421.9327.1626.651.37.937.418.317.799.949.437.416.907.517.009.549.036.375.856.986.469.248.731.42.902.363.763.224.694.163.312.763.492.944.574.042.942.393.322.784.493.961.52.071.502.271.702.842.292.011.432.171.602.822.261.871.282.131.552.822.261.61.521.001.671.052.051.461.510.881.641.022.051.471.460.821.641.032.071.491.71.330.661.380.731.641.031.280.601.380.731.661.051.270.591.400.751.691.071.81.180.461.230.541.420.771.170.441.240.551.440.791.170.451.260.571.460.811.91.110.341.150.411.280.601.110.341.160.431.300.621.110.351.170.451.310.642.01.060.261.100.331.190.481.070.271.110.341.200.491.070.281.120.361.210.512.51.010.091.010.121.030.181.010.091.020.121.030.191.010.101.020.131.040.193.01.000.031.000.041.010.071.000.031.000.051.010.081.000.041.000.051.010.083.51.000.011.000.021.000.031.000.011.000.021.000.031.000.011.000.021.000.034.01.000.001.000.011.000.011.000.011.000.011.000.011.000.011.000.011.000.01Table 4
ARLs and SDRLs of proposed control chart for gamma distribution whenα0=5 and ARL0=500.
sk1=3.2215k1=3.7115k1=4.5980k1=3.3615k1=3.9440k1=4.4575k1=3.5845k1=4.3465k1=4.4760k2=1.5255k2=1.8115k2=2.2670k2=1.5835k2=1.8755k2=2.3070k2=1.6325k2=1.9275k2=2.3390m=4,k=2m=4,k=3m=4,k=4m=5,k=3m=5,k=4m=5,k=5m=6,k=4m=6,k=5m=6,k=6ARLSDRLARLSDRLARLSDRLARLSDRLARLSDRLARLSDRLARLSDRLARLSDRLARLSDRL1.0500.03499.53500.04499.54500.05499.55500.01499.51500.02499.52500.05499.55500.04499.54500.01499.51500.05499.551.1270.28269.78271.57271.07273.31272.81263.92263.42266.57266.07268.09267.59259.99259.49262.25261.75264.22263.711.2126.46125.96123.19122.69126.16125.66117.99117.49117.61117.11121.35120.84112.09111.58113.21112.71117.87117.371.361.1360.6357.0556.5560.1559.6554.2253.7253.0952.5957.0856.5849.4948.9950.1849.6854.9454.431.432.0831.5828.9928.4831.5931.0927.2926.7926.5126.0029.7729.2724.1823.6824.7824.2828.5428.031.518.4117.9116.3715.8618.4017.8915.2414.7314.8414.3317.3116.8013.2912.7813.8213.3116.6016.091.611.5010.9910.209.6811.7511.249.398.889.238.7111.0910.588.157.638.618.0910.6710.151.77.747.236.916.398.127.606.315.796.285.757.707.185.494.975.895.377.456.931.85.565.045.024.505.995.464.564.034.604.075.715.194.013.474.353.825.565.041.94.233.693.873.334.654.123.502.963.583.044.483.943.122.573.422.884.393.862.03.362.823.132.583.783.242.832.272.922.373.663.122.562.002.822.273.613.072.51.701.091.691.082.011.421.560.931.661.052.011.421.510.871.661.052.021.443.01.290.621.320.651.510.881.250.561.330.661.520.891.240.551.340.671.540.913.51.150.411.180.461.300.621.140.391.190.481.310.641.140.401.200.491.320.654.01.090.301.110.351.190.481.080.301.120.371.200.491.090.311.130.381.210.50Table 5
ARLs and SDRLs of proposed control chart for gamma distribution whenα0=10 and ARL0=500.
sk1=3.2635k1=3.8515k1=4.5495k1=3.3245k1=5.0695k1=4.5285k1=3.6205k1=4.4735k1=4.5660k2=1.5185k2=1.8085k2=2.2730k2=1.6015k2=1.8730k2=2.3135k2=1.6325k2=1.9290k2=2.3460m=4,k=2m=4,k=3m=4,k=4m=5,k=3m=5,k=4m=5,k=5m=6,k=4m=6,k=5m=6,k=6ARLSDRLARLSDRLARLSDRLARLSDRLARLSDRLARLSDRLARLSDRLARLSDRLARLSDRL1.0500.01499.51500.02499.51500.05499.55500.04499.54500.03499.53500.05499.55500.05499.55500.05499.55500.05499.551.1221.13220.63222.17221.66223.48222.98215.54215.04217.15216.65218.55218.05208.30207.80209.08208.58214.77214.271.271.1670.1672.2171.7175.7575.2569.5269.0268.1567.6572.5372.0362.5262.0163.2062.7070.1669.661.325.3124.8126.1225.6229.1128.6124.8224.3224.0823.5827.6127.1021.1920.6821.9521.4526.5526.041.411.0610.5511.5711.0613.5913.0810.7610.2510.6010.0912.8912.389.058.539.679.1612.4311.921.55.895.376.215.697.557.035.675.145.735.207.226.704.824.295.294.767.016.491.63.212.883.903.364.824.303.522.983.643.104.664.133.072.533.432.894.584.051.72.502.052.762.213.442.892.491.932.632.073.362.822.241.672.521.963.342.791.82.101.322.151.572.662.101.951.362.081.502.632.071.811.212.031.452.632.071.91.591.151.791.182.191.611.631.021.751.152.181.611.550.931.741.132.201.622.01.440.911.560.941.881.291.440.801.550.921.891.301.400.751.550.931.911.322.51.120.361.150.421.270.591.110.361.160.431.290.611.120.361.170.451.300.623.01.040.201.060.241.110.341.040.201.060.251.110.361.040.211.070.271.120.363.51.010.121.020.151.050.221.020.131.020.161.050.231.020.131.030.171.050.234.01.010.081.010.101.020.151.010.081.010.101.020.151.010.081.010.111.020.15Table 6
ARLs and SDRLs of proposed control chart for gamma distribution whenα0=20 and ARL0=500.
sk1=4.6725k1=4.2825k1=4.5495k1=3.4165k1=4.4505k1=4.5120k1=3.5885k1=4.6790k1=4.7185k2=2.2755k2=1.8025k2=2.2730k2=1.5815k2=1.8740k2=2.3165k2=1.6365k2=1.9295k2=2.3490m=4,k=2m=4,k=3m=4,k=4m=5,k=3m=5,k=4m=5,k=5m=6,k=4m=6,k=5m=6,k=6ARLSDRLARLSDRLARLSDRLARLSDRLARLSDRLARLSDRLARLSDRLARLSDRLARLSDRL1.0500.04499.54500.04499.54500.05499.55500.05499.55500.05499.55500.05499.55500.04499.54500.05499.55500.05499.551.1153.53153.03154.54154.04155.79155.29143.46142.96145.72145.22151.15150.65138.53138.03139.18138.68147.91147.411.230.9330.4232.1131.6035.5535.0527.6027.1028.8728.3633.8033.3025.7425.2326.6726.1732.6432.131.38.237.729.468.9511.4410.938.087.068.477.9510.8910.387.196.677.857.3310.5510.041.43.292.754.123.585.194.673.583.043.773.235.024.493.172.633.583.044.934.401.52.121.552.411.843.052.512.111.542.291.713.012.461.961.382.241.663.002.451.61.691.081.731.132.151.571.550.931.691.082.151.571.500.871.691.082.171.591.71.360.701.420.771.701.101.310.631.410.761.721.111.290.621.430.781.741.141.81.200.491.250.561.460.811.180.461.260.571.470.841.180.471.280.601.490.861.91.120.361.160.431.310.631.110.361.170.451.320.651.120.371.190.471.330.672.01.070.281.110.341.210.501.070.281.120.361.220.521.080.291.130.381.230.532.51.010.091.010.121.030.191.010.101.020.131.040.201.010.101.020.141.040.203.01.000.031.000.051.010.081.000.041.000.051.010.081.000.041.000.051.010.083.51.000.011.000.021.000.031.000.011.000.021.000.031.000.021.000.021.000.044.01.000.011.000.011.000.011.000.011.000.011.000.021.000.011.000.011.000.02We pointed out the following several noteworthy comments from Tables1–6 for the developed control charts:(1)
The out-of-control ARL and SDRL values decline speedily when the shift (s) of the manufacturing process increases(2)
It is detected that the chart coefficientk2 shows an increasing tendency for increased value of k for a fixed value of m when other parametric combinations are fixed(3)
From the tables, it is noticed that ARL1 and SDRL values decrease when m values increase. In addition, ARL1 and SDRL values increased with the increase of k value (i.e., m-2 to m-0). It also observed the same inclination over the other parametric combinations and ARL0=370 and 500(4)
It is interesting to observe from the results that the values of ARL1 and SDRL are small for k=m-2 and these values are increasing from k=m-2 to k=m for fixed values of m. In addition, noticed that ARL1 and SDRL values are large at k=m as compared to the values at k=m-1 and k=m-2 (we know that if k=m, 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 whenARL0=370 and ARL0=500; the shape parameter of gamma distribution is given as α0=5, 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 ARL0=370 and m=4 and in Table 8 for ARL0=500 and m=5 at various shape parameters of the gamma distribution.Table 7
Developed gamma control chart ARL comparison with existing control charts whenARL0=370 and m=4.
sα0=5α0=10α0=20Proposedk1=3.1125k2=1.5025MDSk1=3.0025k2=2.5235ShewhartL=2.9605Proposedk1=3.0575k2=1.5790MDSk1=3.0135k2=2.5765ShewhartL=2.9821Proposedk1=3.2715k2=1.4485MDSk1=4.4505k2=2.2170ShewhartL=2.99171.0370.05370.02370.96370.01370.02370.96370.02370.05370.961.1206.61208.25217.16176.21179.79188.41119.55121.50144.071.2100.36104.12114.7366.3271.5080.5425.9229.4446.781.350.2354.5663.7827.0131.6138.477.939.9419.111.427.1731.1038.4412.6516.1220.972.904.699.591.516.0019.2724.986.869.3712.772.072.845.641.610.2112.8717.304.256.088.501.522.053.751.77.009.1512.642.944.326.071.331.642.731.85.116.879.642.233.304.601.181.422.151.93.935.397.631.822.663.641.111.281.782.03.164.406.221.562.253.001.061.191.552.51.662.303.071.131.401.651.011.031.113.01.281.672.071.041.161.271.001.011.023.51.141.401.631.021.071.121.001.001.004.01.081.251.401.011.031.061.001.001.00Table 8
Developed gamma control chart ARLs comparison with existing control charts when ARL0 =500 and m =5.
sα0=5α0=10α0=20Proposedk1=3.3615k2=1.5835MDSk1=4.4575k2=2.3070ShewhartL=3.0458Proposedk1=3.3245k2=1.6015MDSk1=4.5285k2=2.3135ShewhartL=3.0701Proposedk1=3.4165k2=1.5815MDSk1=4.5120k2=2.3165ShewhartL=3.08041.0500.01500.05500.94500.04500.05500.93500.05500.05500.041.1263.92268.09283.07215.54218.55244.74143.46151.15185.051.2117.99121.35144.6069.5272.53100.6327.6033.8057.411.354.2257.0878.2624.8227.6146.638.0810.8922.641.427.2929.7746.1510.7612.8924.793.585.0211.041.515.2417.3129.455.677.2214.782.113.016.331.69.3911.0920.093.524.669.661.552.154.121.76.317.7014.472.493.366.801.311.722.961.84.565.7110.911.952.635.081.181.472.291.93.504.488.541.632.183.981.111.321.882.02.833.666.901.441.893.241.071.221.612.51.562.013.291.111.291.721.011.041.123.01.251.522.171.041.111.291.001.011.033.51.141.311.691.021.051.141.001.001.014.01.081.201.441.011.021.071.001.001.00It is noticed that from the results on the basis of Tables7 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 (s) values and various parametric values studied in this article. At a glance, when ARL0=370, α0=5 and s=1.4 from Table 7, for the developed control chart ARL1=27.17 whereas ARL1=31.10 for MDS scheme and ARL1=38.44 from the Shewhart type control chart. Similarly, for ARL0=500, α0=10, and s=1.5 from Table 8, we sense that the developed control chart gives ARL1=5.67 while ARL1=7.22 for the MDS control chart and ARL1=14.28 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.Figure 1
ARL curves of gamma control chart for three charts form=4 and ARL0=370.Figure 2
ARL curves of gamma control chart for three charts form=5 and ARL0=500.


### 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α0=5 and in-control scale parameter 1 and last 30 random samples are generated from a gamma distribution with shape parameter α0=5 and out-of-control scale parameter 1.4 (i.e., the shift of s=1.4). The data is reported in Table 9 alongside computed statistic Xi∗=Xi1/3. The control chart coefficients at m=5, α0=5, and ARL0=500 are available in Table 8. The Shewhart type gamma control chart is given in Figure 3, and MDS gamma control chart when m=5 and k=5 is provided in Figure 4. The gamma control chart using GMDS sampling when m=5 and k=3 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 m=5) Xi∗ 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 k=3 and m=5) Xi∗ values are within the interior control limits.Table 9
The simulated data whenm=5, α0=5, and ARL0=500.
S. no.XX∗S. no.XX∗S. no.XX∗S. no.XX∗17.60631.9666162.36081.3315318.04372.0036464.62611.666222.87431.4218172.58001.3715324.66801.6713475.49351.764533.23011.4782184.96371.7058338.64532.0524481.99991.259944.66711.6712194.10351.6010349.13952.0908496.44791.861357.03981.9165202.76471.4035354.02001.59005011.10012.230764.36211.6339213.87281.5704366.27341.8443517.49441.956973.51451.5204228.00952.0008372.75841.4024528.13112.010987.88311.9902238.31952.0263386.49971.8662533.43481.508897.23281.9339243.28211.4861397.64331.9698543.28601.4867107.06401.9187255.59561.77534010.10042.1616554.86311.6942113.71531.5488263.46081.51264111.09292.2302567.67221.9723125.00591.7106274.24621.6193425.99051.8162576.28981.8459132.62491.3794286.65231.8807433.94661.5803586.14691.8318143.69221.5456296.57041.8730449.11402.0888599.31272.1039154.54751.6567305.25601.73874511.00672.2244603.22131.4769Figure 3
Gamma control chart for Shewhart type for simulated data whenm=5, α0=5, and ARL0=500.Figure 4
Gamma control chart using MDS sampling for simulated data whenm=5, α0=5, and ARL0=500.Figure 5
Gamma control chart using GMDS sampling for simulated data whenm=5, α0=5, and ARL0=500.From Figures3 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.0026≈2.0, and the scale parameter is 3.9185. The Kolmogorov-Smirnov test is 0.1197, and p 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.Table 10
Death cause by COVID-19 in Wuhan city of Hubei province in China.
Subgroup no.Days from ICU intake to death (X)X∗Subgroup no.Days from ICU intake to death (X)X∗141.587418172.5713261.81711951.7100351.71002082.0000471.91292182.0000551.71002282.0000641.58742311.0000721.259924122.2894861.817125102.15449102.154426122.28941011.00002741.58741171.91292841.58741292.08012921.259913222.80203041.587414112.22403111.00001561.817132112.22401682.000033252.924017142.4101Figure 6
The empirical and theoretical cdfs and Q-Q plots for the GD for the days from ICU intake to death.The chart constants at the estimated shape parameter of 2 are obtained using a simulation procedure given in Section2. The control chart constants at α0=2 and m=4 are k1=3.1035 and k2=1.4645 for a developed control chart, k1=3.7525 and k2=2.1935 for MDS control chart, and L=2.8828 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 7–9. In Figures 7–9, the proposed control charts for GMDS, MDS, and Shewhart type control charts are displayed by plotting the control limits and chat statisticsX∗. 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 X∗ 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 X∗ fall in the inner control limits.Figure 7
Shewhart type gamma control chart for days from ICU intake to death data.Figure 8
Gamma control chart using MDS sampling for days from ICU intake to death data.Figure 9
Gamma control chart using GMDS sampling for days from ICU intake to death data.It is clearly noticeable from Figure7 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.


## 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α0=5 and in-control scale parameter 1 and last 30 random samples are generated from a gamma distribution with shape parameter α0=5 and out-of-control scale parameter 1.4 (i.e., the shift of s=1.4). The data is reported in Table 9 alongside computed statistic Xi∗=Xi1/3. The control chart coefficients at m=5, α0=5, and ARL0=500 are available in Table 8. The Shewhart type gamma control chart is given in Figure 3, and MDS gamma control chart when m=5 and k=5 is provided in Figure 4. The gamma control chart using GMDS sampling when m=5 and k=3 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 m=5) Xi∗ 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 k=3 and m=5) Xi∗ values are within the interior control limits.Table 9
The simulated data whenm=5, α0=5, and ARL0=500.
S. no.XX∗S. no.XX∗S. no.XX∗S. no.XX∗17.60631.9666162.36081.3315318.04372.0036464.62611.666222.87431.4218172.58001.3715324.66801.6713475.49351.764533.23011.4782184.96371.7058338.64532.0524481.99991.259944.66711.6712194.10351.6010349.13952.0908496.44791.861357.03981.9165202.76471.4035354.02001.59005011.10012.230764.36211.6339213.87281.5704366.27341.8443517.49441.956973.51451.5204228.00952.0008372.75841.4024528.13112.010987.88311.9902238.31952.0263386.49971.8662533.43481.508897.23281.9339243.28211.4861397.64331.9698543.28601.4867107.06401.9187255.59561.77534010.10042.1616554.86311.6942113.71531.5488263.46081.51264111.09292.2302567.67221.9723125.00591.7106274.24621.6193425.99051.8162576.28981.8459132.62491.3794286.65231.8807433.94661.5803586.14691.8318143.69221.5456296.57041.8730449.11402.0888599.31272.1039154.54751.6567305.25601.73874511.00672.2244603.22131.4769Figure 3
Gamma control chart for Shewhart type for simulated data whenm=5, α0=5, and ARL0=500.Figure 4
Gamma control chart using MDS sampling for simulated data whenm=5, α0=5, and ARL0=500.Figure 5
Gamma control chart using GMDS sampling for simulated data whenm=5, α0=5, and ARL0=500.From Figures3 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.0026≈2.0, and the scale parameter is 3.9185. The Kolmogorov-Smirnov test is 0.1197, and p 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.Table 10
Death cause by COVID-19 in Wuhan city of Hubei province in China.
Subgroup no.Days from ICU intake to death (X)X∗Subgroup no.Days from ICU intake to death (X)X∗141.587418172.5713261.81711951.7100351.71002082.0000471.91292182.0000551.71002282.0000641.58742311.0000721.259924122.2894861.817125102.15449102.154426122.28941011.00002741.58741171.91292841.58741292.08012921.259913222.80203041.587414112.22403111.00001561.817132112.22401682.000033252.924017142.4101Figure 6
The empirical and theoretical cdfs and Q-Q plots for the GD for the days from ICU intake to death.The chart constants at the estimated shape parameter of 2 are obtained using a simulation procedure given in Section2. The control chart constants at α0=2 and m=4 are k1=3.1035 and k2=1.4645 for a developed control chart, k1=3.7525 and k2=2.1935 for MDS control chart, and L=2.8828 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 7–9. In Figures 7–9, the proposed control charts for GMDS, MDS, and Shewhart type control charts are displayed by plotting the control limits and chat statisticsX∗. 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 X∗ 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 X∗ fall in the inner control limits.Figure 7
Shewhart type gamma control chart for days from ICU intake to death data.Figure 8
Gamma control chart using MDS sampling for days from ICU intake to death data.Figure 9
Gamma control chart using GMDS sampling for days from ICU intake to death data.It is clearly noticeable from Figure7 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.



---

*Source: 6634887-2021-04-29.xml*