# A New Generalized Range Control Chart for the Weibull Distribution

Authors: Osama H. Arif; Muhammad Aslam
Journal: Complexity
 (2018)

---

## Abstract

In this study, a generalized range control chart is designed for the Weibull distribution using generally weighted moving average statistics. The proposed chart is based on minimum generally weighted moving average statistic and maximum generally weighted moving average statistics. We utilize the inverse erf function to transform the Weibull data to normal data. The necessary measures are given to assess the performance of the proposed control chart. The comparison study shows that the proposed control chart outperforms the existing control charts based on exponentially weighted moving average statistic in terms of the average run length. A real example is given for applying the proposed chart in the industry.


---

## Body

## 1. Introduction

The control charts are an integral part of statistical process control which is used to monitor and improve the manufacturing process in the industry. Among the control charts, the Shewhart control charts are preferred due to ease in assumptions and use in the industry. These charts are usually designed under the assumption that the quality of interest follows the normal distribution. But the data obtained from the chemical, semiconductor, cutting tool wear, and failure time processes are often skewed [1, 2]. So, monitoring the skewed data using the Shewhart control chart may mislead to observe the manufacturing process.The Weibull distribution has been widely used to characterize material breaking strength, describe the lifetime of electronic items, and monitor the reliability data obtained from life test [3, 4]. For the skewed data, it is not a proper approach to monitor the process using Shewhart X-bar, standard deviation, and range charts. According to [5], these charts are not effective to detect a shift in the process in low population percentiles of these distributions and increase the false signal out-of-control. Several people worked on the designing of control chart for the Weibull distribution. Nelson [6] designed various control charts for the Weibull process. Bai and Choi [1] designed X-bar and range charts to monitor the skewed data. Hawkins and Olwell [7] proposed cumulative sum (CUSUM) control chart to monitor the scale parameter of the Weibull distribution. Ramalhoto and Morais [8] focused on the development of Shewhart for the Weibull distribution with fixed and variable sampling intervals. Chang and Bai [9] designed various control chart including exponentially weighted moving average (EWMA) chart using weighted standard deviations for the skewed distribution. More details about these charts can be seen in [10–12].Transforming the Weibull distributed data to normal data may make ease for industrial engineering. Therefore, several transformation methods including power transformation, inverse erf function, and Box–Cox transformation are available in the literature. Batson et al. [13] designed moving average control chart using the power transformation. Pascual [14] studied moving average chart using natural logarithm of the Weibull distribution. More details about control chart using transformation can be read in [3, 15–17] and in [18].Recently, Wang [19] proposed the MaxEWMA control chart for the Weibull distribution using three different transformations. Wang [19] recommended that inverse erf function performs better than power transformation and Box–Cox transformation when shape parameter of the Weibull distribution is greater than one. Woodall and Montgomery [20] pointed that development of control chart for the Weibull distribution has a potential for future research.In this study, a new generalized range control chart is designed for the Weibull distribution using inverse erf function which has been commonly used to transform the Weibull distributed data to normal data. Some tables for practical use are also given. The comparison study shows that the proposed control chart outperforms the existing control chart in terms of average run length. A real example is also given to apply the proposed chart in the industry.


## 2. Design of Proposed Control Chart

The Weibull distribution with shape parameterβ and scale parameter θ has the following probability density function:
(1)fx;θ,β=βxβ−1θβe−x/θβ,x>0,θ>0,β>0.Usually, industrial engineers have the information about the parametersβ and θ. These parameters can be estimated from the monitoring data if unknown in the practice. This density function is skewed, so symmetric type of control limits may not work for a Weibull distribution. But there are several methods of transforming Weibull distributed data to a normal distributed data. Wang [19] discussed and analyzed three transformation methods and recommended the method in [17] when β≥1.Faraz et al. [17] suggested the use of erf function to convert the Weibull distributed random variable to standard normal variable. When x is a Weibull distributed data, the following transformed data were shown to follow a standard normal data:
(2)u=2erf−11−2e−x/θβ,where erfz=1/π∫−zze−u2du.As suggested by Wang [19], we utilized another statistic from a sequence of above transformed data ui’s:
(3)vi=ϕ−1χ12ui−ui−12,i=1,2,3,…n.Note here thatΦ−1. presents the inverse cumulative standard normal distribution and χ12. cumulative chi-square distribution with one degree of freedom.The generally weighted moving average (GWMA) statistic is the generalization of traditional EWMA statistic. The GWMA statistic empowers the decision about the state of process by utilizing the weight assigned to current and past information; for more details, see [21–24]. Suppose that M denotes the number of samples in the sequence of independent samples until the 1st occurrence of event A since the previous occurrence of event A[25]. Suppose that PM=1,PM=2,…,PM=i are assigned weight for current and previous sample such that PM>i shows the weight to target of the process.Two GWMA statistic to monitor the mean and variability of the process are defined.The GWMA statistics applied to two transformed data are given by(4)Gi=PM=1Ui+PM=2Ui−1+⋯+PM=2U1+PM>iG0,Hi=PM=1Vi+PM=2Vi−1+⋯+PM=iV1+PM>1H0,where i=1,2,3,… and the starting values of G0 and H0 are usually set at 0.The MaxEWMA control chart has the ability to detect a small shift in the mean or variability than the traditional EWMA control charts [25]. The range chart (Ri) which is computed from the difference of maximum quantity, say maxGWMA, and minimum quantity, say minGWMA, has been commonly used to monitor the system variability.We define the maxGWMA statistic as(5)MGi,max=maxGi,Hi,i=1,2,….The minGWMA statistic is defined as follows(6)MGi,min=minGi,Hi,i=1,2,….Following [26], Sheu et al. [25] suggested that when Ui and Vi are independently normally distributed, the quantity PM>i can be evaluated as follows for the design parameter q and adjustment parameter α.
(7)PM=i=qi−1α−qiα,where i=1,2,3,…, 0≤q≤1, and α>0.Now, the plotting statisticR can be defined as follows
(8)Ri=MGi,max−MGi,min.The proposed control chart is stated as follows.Step 1.
Obtain the Weibull dataXi’s at each subgroup of size n and transform it to the standard normal data using Ui=2erf−11−2e−Xi/θβ and compute Vi=Φ−1χ12Ui−Ui−12.Step 2.
Obtain GWMA statisticsGi and Hi. Determine MGi,max=maxGi,Hi and MGi,min=minGi,Hi and compute statistic Ri=MGi,max−MGi,min.Step 3.
Declare the process out-of-control ifRi>UCL or Ri<LCL.
The proposed control chart is a generalization of several control charts. The proposed control chart reduces the chart in [27] when α=1 and q=1−λ, where λ is smoothing constant in EWMA statistic. The proposed control chart becomes the chart in [19] when α=1 and R=0.The operational process of the proposed control chart is based on two control limits which are defined as(9)UCL=R¯+kSDR,LCL=R¯−kSDR,where R¯ is the average and SD is the standard deviation of statistic Ri, respectively. Also, k is the control chart coefficient which can be determined by closed-form solution, the Markov chain approach, and Monte Carlo simulation (MCS). These are alternative methods for obtaining the properties of the run-length distribution. In this paper, the control chart coefficient will be determined through the following MCS.
(1)
Specify the shape parameterβ and scale parameter θ of the Weibull distribution. Generate a random sample of size n from the Weibull distribution by assuming in-control state
(2)
Compute statisticRi and plot it on the control limits. Note that the first of out-of-control is called run length (RL). Repeat this process 10,000 times
(3)
Compute average run length (ARL) and standard deviation (SD) using these RL. Determinek for which ARL0≥r0, where r0 is specified ARL
(4)
Repeat Steps1-3 to find the ARL for the shifted process, say ARL1 when scale parameter is shifted to θ1=dθ0, where d is a constantUsing the above MCS, the values ofARL1 and SD are reported for various control chart parameters in Tables 1–5 as well as the mean and standard deviation when q=0.95, n=5, and various values of α, β, and θ. The ARL1, k, and SD for various values of shape and scale parameters of the Weibull distribution and shift constant d are reported in Tables 2–4.Table 1
Mean and standard deviation at different level.
q
n
α
β
θ
Mean
SD
0.95
5
1.1
1
50
6.6191
3.3631
0.95
5
1.1
1.5
50
6.6187
3.3612
0.95
5
1.1
2
50
6.6254
3.3679
0.95
5
1
1
50
6.6876
3.4146
0.95
5
1
1.5
50
6.6847
3.4188
0.95
5
1
2
50
6.6786
3.4004
0.95
5
0.8
1
50
6.7956
3.4937
0.95
5
0.8
1.5
50
6.7950
3.4979
0.95
5
0.8
2
50
6.7872
3.4946
0.95
5
1.1
1
100
6.6293
3.3689
0.95
5
1.1
1.5
100
6.6285
3.3737
0.95
5
1.1
2
100
6.6204
3.3653
0.95
5
1
1
100
6.6840
3.4112
0.95
5
1
1.5
100
6.6935
3.4197
0.95
5
1
2
100
6.6841
3.4160
0.95
5
0.8
1
100
6.7972
3.4950
0.95
5
0.8
1.5
100
6.8013
3.5018
0.95
5
0.8
2
100
6.7910
3.4972
0.95
5
0.8
3.2
2.4
6.7891
3.4936
0.95
5
0.8
0.92
53.41
6.7850
3.4944Table 2
The values of ARLs whenq=0.95, α=1.10, θ=50, and β=1.0,1.5,2.0.
n=5, ARL0=370
β
1
1.5
2
k
4.25
4.24
4.226
d
ARL
SD
ARL
SD
ARL
SD
1
371.78
366.58
370.39
366.12
370.95
365.76
0.99
347.09
340.70
337.02
331.35
338.41
332.47
0.98
312.59
310.27
303.72
301.29
305.73
303.73
0.95
228.88
227.29
223.44
222.09
224.57
222.97
0.92
169.19
170.02
165.97
166.05
167.07
168.33
0.90
136.61
136.92
132.96
132.62
133.32
132.19
0.85
83.46
82.18
82.53
82.55
82.91
83.09
0.80
52.48
52.28
51.06
49.90
51.11
50.30
0.75
31.03
29.71
32.18
31.50
26.90
26.56
0.7
17.31
16.40
21.65
19.03
21.16
21.81
0.65
12.57
11.75
12.38
11.30
12.85
12.89
0.60
8.44
7.31
8.30
7.32
7.68
6.71Table 3
The values of ARLs whenq=0.95, α=1.0, θ=50, and β=1.0,1.5,2.0.
n=5, ARL0=370
β
1
1.5
2
k
4.245
4.24
4.245
d
ARL
SD
ARL
SD
ARL
SD
1
370.89
362.40
371.97
366.53
369.51
359.22
0.99
346.87
346.27
344.74
338.86
336.47
328.97
0.98
313.37
304.43
311.20
308.88
301.17
300.46
0.95
225.28
229.94
227.23
228.12
223.06
224.30
0.92
168.89
164.14
169.99
167.66
165.44
162.40
0.9
139.04
140.38
136.79
136.48
134.74
134.46
0.85
85.76
85.73
83.86
83.05
81.63
80.44
0.8
52.48
51.16
52.38
51.41
51.54
50.79
0.75
32.75
32.81
32.61
32.48
32.30
31.65
0.7
21.20
18.72
19.37
17.65
19.57
19.11
0.65
13.58
13.05
13.86
13.09
12.28
12.72
0.6
8.22
7.40
8.20
8.20
9.65
7.91Table 4
The values of ARLs whenq=0.95, α=0.80, θ=50, and β=1.0,1.5,2.0.
n=5, ARL0=370
β
1
1.5
2
k
4.255
4.245
4.245
d
ARL
SD
ARL
SD
ARL
SD
1
373.74
367.11
371.71
369.38
373.00
366.31
0.99
345.51
336.44
337.18
333.44
345.25
336.30
0.98
310.97
309.32
308.68
307.83
310.29
308.20
0.95
227.81
228.32
229.18
227.59
227.54
227.96
0.92
170.82
168.79
169.74
164.95
170.77
168.88
0.9
137.30
136.83
135.39
135.90
136.80
136.52
0.85
83.51
82.64
85.79
85.58
83.72
83.06
0.8
52.97
52.01
52.46
52.45
52.57
51.64
0.75
32.66
32.38
32.93
32.33
32.78
32.35
0.7
19.55
16.39
21.76
21.24
22.85
22.10
0.65
15.26
13.62
14.51
13.51
14.75
13.55
0.6
8.68
9.23
9.40
8.39
9.50
9.80Table 5
The values of ARLs whenq=0.95and α=0.80.
n=5, ARL0=370
β
β=0.92, θ=53.41
β=3.2, θ=2.4
k
4.25
4.245
d
ARL
SD
ARL
SD
1
372.81
366.44
368.06
365.01
0.99
332.58
329.81
334.46
331.46
0.98
303.79
306.22
307.32
305.61
0.95
230.89
229.60
225.06
223.30
0.92
166.86
168.52
169.12
164.21
0.9
137.81
138.77
134.44
136.91
0.85
83.46
81.42
83.10
81.47
0.8
52.10
50.96
52.71
52.80
0.75
32.85
32.67
32.67
32.51
0.7
23.11
24.92
21.06
19.98
0.65
13.65
12.36
14.18
16.64
0.6
8.89
9.63
8.92
7.88From Tables1–5, we note the following trend in the control chart parameters.
(1)
For same values of all other parameters, the values ofARL1 and SD decrease as α values decrease. For example, when α=1.10 and d=0.99, the value of ARL1 from Table 2 is 347, and it is 346 when α=1.0
(2)
For the same values of all other parameters, the values ofARL1 and SD decrease for β≤1.5 and decrease when β≥2
(3)
For the same values of all other parameters andq=0.95, α=0.80, the values of ARL1 and SD increase when β increases and θ decreases


## 3. Comparison and Simulation Study

In this section, we compare the performance of the proposed control chart with chart based on EWMA statistic and the chart in [19].


### 3.1. Proposed Chart vs EWMA Chart

As mentioned earlier, the proposed control chart is the extensions of several charts. The proposed control chart reduces to control chart based on EWMA statistic whenα=1 and q=1−λ. The performance of the proposed control chart will be compared with chart based on EWMA statistic proposed by [27] using the simulation data. In this study, it is assumed that ARL0=370, n=5, and α=1.10. Suppose that the process is declared to be θ0=50 and β=2.0. The first 20 observations of subgroup size n=5 are generated from the Weibull distribution with parameters β=2.0 and θ0=50. The next 30 observations of subgroup size n=5 are generated from the Weibull distribution with parameters β=2.0 and shift in scale parameter that is d=0.75. The simulation data, Gi and Hi, of the proposed control chart are shown in Table 6.Table 6
The simulation data whenα=1.10 and β=2.0.
Sr#
xi
Gi
Hi
Ri
1
32.24
17.55
48.11
33.38
43.82
−1.11
−3.38
−2.68
−3.22
−2.85
−3
−2.51
−0.29
−0.42
−1.15
3.08
2
30.48
35.51
78.97
43.77
15.71
−1.08
−1.53
1.26
1.38
−1.27
−0.74
−2.47
0.18
1.87
3.71
3.53
3
41.35
59.57
19.31
90.87
43.95
0.26
1.61
−0.6
2.92
2.95
−0.45
−0.31
2.53
7.56
9.76
9.49
4
58.55
35.03
37.58
33.85
26
1.13
0.52
0.15
−0.52
−1.9
0.02
0.81
−1.85
−3.94
−4.77
4.75
5
48.06
80.37
38.48
27.83
42.64
0.7
3.46
3.01
1.62
1.59
−3.19
−1.64
0.79
0.15
0.11
3.34
6
78.87
33.23
47.25
16.96
39.57
2.69
1.83
2.16
−0.37
−0.54
2.69
5.28
4.73
6.52
7.44
7.07
7
42.42
28.86
64.51
55.99
43.07
0.25
−0.88
0.88
1.94
1.95
−2.44
−2.53
−0.33
−1.63
−2.09
2.28
8
33
22.6
24.1
35.3
10
−1.5
−3.2
−4.6
−4.9
−8
−1
−1.5
−4.3
−4.5
−2.1
6.98
9
12.41
34.61
57.29
54.83
28.37
−3.29
−3.7
−2.29
−1.12
−2.2
1.9
3.37
3.89
0.83
2
3.05
10
37.5
33.98
51.95
32.21
32.07
−0.52
−1.14
−0.27
−1.05
−1.81
−5.82
−7.76
−7.1
−6.22
−10.81
10.54
11
53.78
86.47
24.49
26.85
59.38
1.2
4.34
2.55
1.08
2.35
−1.58
−0.16
4.1
1.36
3.12
4.18
12
98.47
49.62
61.55
67.63
87.55
5.08
5.44
6.63
8.17
10.96
0.67
3.29
2.35
0.37
0.37
10.58
13
9.196
70.773
48.55
49.672
23.734
−3.761
−1.407
−0.767
−0.078
−1.69
3.083
8.214
8.201
4.188
5.243
8.13
14
48.605
11.519
26.189
43.754
57.548
0.297
−2.894
−4.118
−3.714
−2.27
−0.079
3.024
3.58
3.733
3.061
4.03
15
35.87
31.84
46.78
32.71
37.42
−0.68
−1.48
−0.99
−1.69
−1.94
−4.38
−6.21
−5.96
−5.84
−7.36
6.68
16
63.75
8.42
39.62
88.5
57.94
1.89
−1.94
−2.01
1.44
2.62
0.93
5.81
8.41
10.77
11.25
10.31
17
22.75
60.78
28.48
20.93
33.67
−2.13
−0.56
−1.67
−3.51
−3.98
−0.43
2.07
3.73
2.56
2.31
3.55
18
52.27
78.73
16.52
46.44
48.71
1.02
3.66
1
1.31
1.78
−1.55
−0.66
4.06
5.91
2.73
5.25
19
58.07
38.71
71.19
37.11
79.1
1.8
1.46
3.55
2.97
5.51
−2.73
−2.3
−0.62
1.14
3.46
4.88
20
52.24
69.83
26.96
34.98
34.17
0.87
2.91
1.45
0.8
0.11
−1.35
−1.36
1.43
0.32
−3.4
3.29
21
88.6
24.8
136.5
115.1
90.4
5
3.2
9.3
13.8
16.5
−4.2
0.4
7.9
7.5
7.5
16.06
22
27.76
30.53
86.21
46.54
31.06
−1.18
−2.07
1.23
1.57
0.57
0.31
−2.13
1.56
3.49
3.28
3.17
23
34.647
76.426
146.008
41.612
17.992
0.052
2.583
9.353
8.839
6.039
2.325
4.598
8.163
14.244
14.787
14.73
24
13
11.1
22.8
32.7
33.7
−3.9
−6.9
−8.2
−8.5
−8.7
−1
−3.2
−2.7
−3.1
−6.5
7.64
25
59.43
78.15
4.81
27.22
5.87
0.75
3.37
−1.41
−2.63
−6.79
3.07
2.87
9.6
11.7
13.3
12.54
26
51.14
48.06
43.06
65.66
25.54
0.91
1.37
1.41
3.12
1.5
−2.24
−4.67
−6.35
−5.43
−2.56
5.433
27
32.7
16.11
17.84
22.11
3.98
−1.95
−4.36
−6.41
−7.85
−12.25
0.49
1.11
−1.57
−3.14
−0.64
11.76
28
11.86
23.35
43.18
50.06
73.61
−3.28
−4.78
−4.38
−3.45
−0.89
2.19
2.29
2.87
1.17
1.64
3.886
29
4.9
71.5
24.3
56.9
9.8
−5.2
−2.7
−4.1
−2.7
−6
4.4
10.6
13.1
14.3
17.6
14.92
30
89.34
34.78
51.53
13.8
36.46
3.42
2.66
3.28
0.25
−0.22
3.66
6.84
6.51
9.08
10.06
9.83
31
37.77
16.82
59.81
40.03
12.46
−0.76
−3.14
−1.59
−1.62
−4.54
−0.89
0.29
3.47
3.59
5.51
5.22
32
10.8
49.9
8
57.8
75
−3.6
−2.7
−6.4
−4.8
−2
3.3
6.4
10
14
13
11.95
33
88.95
53.4
56.02
68.26
41.26
4.08
4.77
5.6
7.25
6.79
0.85
2.38
−0.5
−1.27
−0.2
7.05
34
45.28
6.59
65.19
30.65
58.42
0.12
−4.01
−2.03
−2.85
−1.4
−0.17
3.7
8.97
10.38
11.08
10.95
35
57.48
33
6.17
27.99
25.49
0.59
−0.17
−4.39
−5.36
−6.5
2.14
2.95
5.61
7.58
4.65
7.40
36
13.44
23.78
20.2
26.33
96.45
−3.28
−4.72
−6.47
−7.46
−3.17
1.34
1.18
−0.84
−2.04
2.81
6.61
37
10.8
4.5
21.7
65.3
31.7
−4.2
−8.6
−10
−7.6
−8
1
1.1
3.1
5.9
7.4
8.97
38
33.52
119.48
110.18
25.36
2.49
−0.46
4.86
9.32
7.34
1.41
1.03
6.54
4.76
10.25
13.09
12.62
39
111.87
20.28
12.98
60.19
30.5
4.85
2.59
−0.52
0.89
−0.13
5.75
11.9
10.59
13.76
14.43
14.29
40
18.74
50.75
25.15
39.65
25.39
−2.63
−1.78
−3.15
−3.13
−4.4
0.25
2.37
3.5
3.31
3.08
4.14
41
28.76
41.14
23.85
104.3
20.06
−1.17
−1.15
−2.7
1.8
−0.31
0.6
0.19
0.57
6.09
11.74
11.55
42
44.23
30.62
8.55
58.3
4.8
−0.55
−1.47
−5.09
−3.53
−7.91
1.87
1.51
3.37
7.67
12.68
12.13
43
14.4
14.36
54.04
27.03
41.23
−3.32
−5.89
−4.59
−5.61
−5.3
0.31
−5.32
−1.96
−0.5
−0.54
5.57
44
34.48
183.41
69.23
10.86
52.78
0.19
9.3
10.85
6.94
7.37
3.83
13.08
19.09
22.86
25.09
24.90
45
19.69
24.96
10.02
97.08
32.18
−2.46
−3.83
−7.04
−2.75
−3.37
0.99
−0.5
0.42
7.36
11.15
10.73
46
71.24
53.7
59.17
25.16
10.36
2.15
2.97
4.14
2.42
−1.11
1.55
1.36
−0.63
1.43
2.19
3.50
47
50.3
7.8
69.5
63.8
80
1
−2.9
−0.7
1
3.8
−4.6
−0.3
5.3
3.1
2.5
4.95
48
56.22
29.3
117.62
50.95
27.29
1.58
0.42
5.57
5.99
4.37
−1.9
−0.55
5.35
8.98
9.47
9.05
49
24.66
74.53
64.5
14.03
77.6
−1.3
1.18
2.84
−0.11
2.49
2.2
5.42
3.96
7.7
12.21
12.09
50
74.3
99.2
21.5
39.2
127.2
3.4
7.2
4.9
4.4
9.9
−1
−0.5
5
5.3
10.5
10.02The calculations of statisticR is also shown in Table 6. The values of statistic R are plotted when α=1.10; β=2.0, α=1.10; β=1.0 (exponential distribution) and α=1.0; β=2.0 ([27] statistic) are plotted on control charts given in Figure 1. The UCL and LCL of the three charts are also shown in Figure 1. From Figure 1, it can be noted that the proposed chart using GWMA statistic detects shift in the process at around 45th values while the proposed control charts using the statistic in [27] and exponential distributions are unable to detect shift in the process. From Figure 1, it can be concluded that the proposed control chart using GWMA statistic performs better than the control chart using the statistic in [27] in detecting earlier shift in the process.Figure 1
The proposed control chart and existing charts using simulation data.
(a)
(b)
(c)


### 3.2. Proposed Chart vs the Chart in [19]

The performance of the proposed chart in terms of ARL and SD will be compared with the MaxEWMA control chart proposed by [19] for the same values of specified control chart parameters. To save the space, the comparison will be given when q=0.95, α=1.10, β=2, and θ=50. The ARLs of both control charts are depicted in Figure 2.Figure 2
The ARL curve between the proposed chart and the chart in [19].From Figure2, it can be noted that the proposed control chart has lower curve of ARLs as compared to the chart proposed by [19]. From Figure 2, it can be observed that the proposed chart has smaller values of ARLs than the chart in [19] for small-shift level that is d≤0.95. For example, when d=0.99, the ARL and SD from the proposed control charts are 338 and 332, respectively, while the ARL and SD from the chart in [19] are 344 and 336, respectively. The two control charts have the same performance when shift is large that is when d≥0.95. So, by comparison, it is concluded that the proposed control chart is more sensitive to detect small shift in the process.


## 3.1. Proposed Chart vs EWMA Chart

As mentioned earlier, the proposed control chart is the extensions of several charts. The proposed control chart reduces to control chart based on EWMA statistic whenα=1 and q=1−λ. The performance of the proposed control chart will be compared with chart based on EWMA statistic proposed by [27] using the simulation data. In this study, it is assumed that ARL0=370, n=5, and α=1.10. Suppose that the process is declared to be θ0=50 and β=2.0. The first 20 observations of subgroup size n=5 are generated from the Weibull distribution with parameters β=2.0 and θ0=50. The next 30 observations of subgroup size n=5 are generated from the Weibull distribution with parameters β=2.0 and shift in scale parameter that is d=0.75. The simulation data, Gi and Hi, of the proposed control chart are shown in Table 6.Table 6
The simulation data whenα=1.10 and β=2.0.
Sr#
xi
Gi
Hi
Ri
1
32.24
17.55
48.11
33.38
43.82
−1.11
−3.38
−2.68
−3.22
−2.85
−3
−2.51
−0.29
−0.42
−1.15
3.08
2
30.48
35.51
78.97
43.77
15.71
−1.08
−1.53
1.26
1.38
−1.27
−0.74
−2.47
0.18
1.87
3.71
3.53
3
41.35
59.57
19.31
90.87
43.95
0.26
1.61
−0.6
2.92
2.95
−0.45
−0.31
2.53
7.56
9.76
9.49
4
58.55
35.03
37.58
33.85
26
1.13
0.52
0.15
−0.52
−1.9
0.02
0.81
−1.85
−3.94
−4.77
4.75
5
48.06
80.37
38.48
27.83
42.64
0.7
3.46
3.01
1.62
1.59
−3.19
−1.64
0.79
0.15
0.11
3.34
6
78.87
33.23
47.25
16.96
39.57
2.69
1.83
2.16
−0.37
−0.54
2.69
5.28
4.73
6.52
7.44
7.07
7
42.42
28.86
64.51
55.99
43.07
0.25
−0.88
0.88
1.94
1.95
−2.44
−2.53
−0.33
−1.63
−2.09
2.28
8
33
22.6
24.1
35.3
10
−1.5
−3.2
−4.6
−4.9
−8
−1
−1.5
−4.3
−4.5
−2.1
6.98
9
12.41
34.61
57.29
54.83
28.37
−3.29
−3.7
−2.29
−1.12
−2.2
1.9
3.37
3.89
0.83
2
3.05
10
37.5
33.98
51.95
32.21
32.07
−0.52
−1.14
−0.27
−1.05
−1.81
−5.82
−7.76
−7.1
−6.22
−10.81
10.54
11
53.78
86.47
24.49
26.85
59.38
1.2
4.34
2.55
1.08
2.35
−1.58
−0.16
4.1
1.36
3.12
4.18
12
98.47
49.62
61.55
67.63
87.55
5.08
5.44
6.63
8.17
10.96
0.67
3.29
2.35
0.37
0.37
10.58
13
9.196
70.773
48.55
49.672
23.734
−3.761
−1.407
−0.767
−0.078
−1.69
3.083
8.214
8.201
4.188
5.243
8.13
14
48.605
11.519
26.189
43.754
57.548
0.297
−2.894
−4.118
−3.714
−2.27
−0.079
3.024
3.58
3.733
3.061
4.03
15
35.87
31.84
46.78
32.71
37.42
−0.68
−1.48
−0.99
−1.69
−1.94
−4.38
−6.21
−5.96
−5.84
−7.36
6.68
16
63.75
8.42
39.62
88.5
57.94
1.89
−1.94
−2.01
1.44
2.62
0.93
5.81
8.41
10.77
11.25
10.31
17
22.75
60.78
28.48
20.93
33.67
−2.13
−0.56
−1.67
−3.51
−3.98
−0.43
2.07
3.73
2.56
2.31
3.55
18
52.27
78.73
16.52
46.44
48.71
1.02
3.66
1
1.31
1.78
−1.55
−0.66
4.06
5.91
2.73
5.25
19
58.07
38.71
71.19
37.11
79.1
1.8
1.46
3.55
2.97
5.51
−2.73
−2.3
−0.62
1.14
3.46
4.88
20
52.24
69.83
26.96
34.98
34.17
0.87
2.91
1.45
0.8
0.11
−1.35
−1.36
1.43
0.32
−3.4
3.29
21
88.6
24.8
136.5
115.1
90.4
5
3.2
9.3
13.8
16.5
−4.2
0.4
7.9
7.5
7.5
16.06
22
27.76
30.53
86.21
46.54
31.06
−1.18
−2.07
1.23
1.57
0.57
0.31
−2.13
1.56
3.49
3.28
3.17
23
34.647
76.426
146.008
41.612
17.992
0.052
2.583
9.353
8.839
6.039
2.325
4.598
8.163
14.244
14.787
14.73
24
13
11.1
22.8
32.7
33.7
−3.9
−6.9
−8.2
−8.5
−8.7
−1
−3.2
−2.7
−3.1
−6.5
7.64
25
59.43
78.15
4.81
27.22
5.87
0.75
3.37
−1.41
−2.63
−6.79
3.07
2.87
9.6
11.7
13.3
12.54
26
51.14
48.06
43.06
65.66
25.54
0.91
1.37
1.41
3.12
1.5
−2.24
−4.67
−6.35
−5.43
−2.56
5.433
27
32.7
16.11
17.84
22.11
3.98
−1.95
−4.36
−6.41
−7.85
−12.25
0.49
1.11
−1.57
−3.14
−0.64
11.76
28
11.86
23.35
43.18
50.06
73.61
−3.28
−4.78
−4.38
−3.45
−0.89
2.19
2.29
2.87
1.17
1.64
3.886
29
4.9
71.5
24.3
56.9
9.8
−5.2
−2.7
−4.1
−2.7
−6
4.4
10.6
13.1
14.3
17.6
14.92
30
89.34
34.78
51.53
13.8
36.46
3.42
2.66
3.28
0.25
−0.22
3.66
6.84
6.51
9.08
10.06
9.83
31
37.77
16.82
59.81
40.03
12.46
−0.76
−3.14
−1.59
−1.62
−4.54
−0.89
0.29
3.47
3.59
5.51
5.22
32
10.8
49.9
8
57.8
75
−3.6
−2.7
−6.4
−4.8
−2
3.3
6.4
10
14
13
11.95
33
88.95
53.4
56.02
68.26
41.26
4.08
4.77
5.6
7.25
6.79
0.85
2.38
−0.5
−1.27
−0.2
7.05
34
45.28
6.59
65.19
30.65
58.42
0.12
−4.01
−2.03
−2.85
−1.4
−0.17
3.7
8.97
10.38
11.08
10.95
35
57.48
33
6.17
27.99
25.49
0.59
−0.17
−4.39
−5.36
−6.5
2.14
2.95
5.61
7.58
4.65
7.40
36
13.44
23.78
20.2
26.33
96.45
−3.28
−4.72
−6.47
−7.46
−3.17
1.34
1.18
−0.84
−2.04
2.81
6.61
37
10.8
4.5
21.7
65.3
31.7
−4.2
−8.6
−10
−7.6
−8
1
1.1
3.1
5.9
7.4
8.97
38
33.52
119.48
110.18
25.36
2.49
−0.46
4.86
9.32
7.34
1.41
1.03
6.54
4.76
10.25
13.09
12.62
39
111.87
20.28
12.98
60.19
30.5
4.85
2.59
−0.52
0.89
−0.13
5.75
11.9
10.59
13.76
14.43
14.29
40
18.74
50.75
25.15
39.65
25.39
−2.63
−1.78
−3.15
−3.13
−4.4
0.25
2.37
3.5
3.31
3.08
4.14
41
28.76
41.14
23.85
104.3
20.06
−1.17
−1.15
−2.7
1.8
−0.31
0.6
0.19
0.57
6.09
11.74
11.55
42
44.23
30.62
8.55
58.3
4.8
−0.55
−1.47
−5.09
−3.53
−7.91
1.87
1.51
3.37
7.67
12.68
12.13
43
14.4
14.36
54.04
27.03
41.23
−3.32
−5.89
−4.59
−5.61
−5.3
0.31
−5.32
−1.96
−0.5
−0.54
5.57
44
34.48
183.41
69.23
10.86
52.78
0.19
9.3
10.85
6.94
7.37
3.83
13.08
19.09
22.86
25.09
24.90
45
19.69
24.96
10.02
97.08
32.18
−2.46
−3.83
−7.04
−2.75
−3.37
0.99
−0.5
0.42
7.36
11.15
10.73
46
71.24
53.7
59.17
25.16
10.36
2.15
2.97
4.14
2.42
−1.11
1.55
1.36
−0.63
1.43
2.19
3.50
47
50.3
7.8
69.5
63.8
80
1
−2.9
−0.7
1
3.8
−4.6
−0.3
5.3
3.1
2.5
4.95
48
56.22
29.3
117.62
50.95
27.29
1.58
0.42
5.57
5.99
4.37
−1.9
−0.55
5.35
8.98
9.47
9.05
49
24.66
74.53
64.5
14.03
77.6
−1.3
1.18
2.84
−0.11
2.49
2.2
5.42
3.96
7.7
12.21
12.09
50
74.3
99.2
21.5
39.2
127.2
3.4
7.2
4.9
4.4
9.9
−1
−0.5
5
5.3
10.5
10.02The calculations of statisticR is also shown in Table 6. The values of statistic R are plotted when α=1.10; β=2.0, α=1.10; β=1.0 (exponential distribution) and α=1.0; β=2.0 ([27] statistic) are plotted on control charts given in Figure 1. The UCL and LCL of the three charts are also shown in Figure 1. From Figure 1, it can be noted that the proposed chart using GWMA statistic detects shift in the process at around 45th values while the proposed control charts using the statistic in [27] and exponential distributions are unable to detect shift in the process. From Figure 1, it can be concluded that the proposed control chart using GWMA statistic performs better than the control chart using the statistic in [27] in detecting earlier shift in the process.Figure 1
The proposed control chart and existing charts using simulation data.
(a)
(b)
(c)


## 3.2. Proposed Chart vs the Chart in [19]

The performance of the proposed chart in terms of ARL and SD will be compared with the MaxEWMA control chart proposed by [19] for the same values of specified control chart parameters. To save the space, the comparison will be given when q=0.95, α=1.10, β=2, and θ=50. The ARLs of both control charts are depicted in Figure 2.Figure 2
The ARL curve between the proposed chart and the chart in [19].From Figure2, it can be noted that the proposed control chart has lower curve of ARLs as compared to the chart proposed by [19]. From Figure 2, it can be observed that the proposed chart has smaller values of ARLs than the chart in [19] for small-shift level that is d≤0.95. For example, when d=0.99, the ARL and SD from the proposed control charts are 338 and 332, respectively, while the ARL and SD from the chart in [19] are 344 and 336, respectively. The two control charts have the same performance when shift is large that is when d≥0.95. So, by comparison, it is concluded that the proposed control chart is more sensitive to detect small shift in the process.


## 4. Case Study

Monitoring the breaking strength fibrous composite is very important in the industries which guarantee the safety of material used in the aero industry and construction of bridges. The application of the proposed control chart is given on the breaking strength of carbon fiber data which is modeled by the Weibull distribution withβ=4.8 and θ=3.2. The data were originally reported by [5] and used by [17] in the area of the control chart. Let ARL0=370 and the control chart coefficient for this real data is 4.27. The data for the breaking strength of carbon fibers for n=5 are reported in Table 7 whereas the first 10 samples are recorded when the process is in-control. The values of statistic Ri are also reported in Table 7.Table 7
The real example data.
Sr#
xi
Gi
Hi
Ri
1
3.7
2.74
2.73
2.5
3.6
2.31
1.58
0.85
−0.44
1.41
0.72
2.65
−1.97
−3.24
−0.72
2.80
2
3.11
3.27
2.87
1.47
3.11
0.16
1.01
0.7
−3.21
−2.64
−0.35
−2.08
−2.3
0.8
4.49
4.32
3
4.42
2.41
3.19
3.22
1.69
4.73
3.01
3.46
3.97
0.43
4.85
10.26
10.81
6.64
9.62
10.37
4
3.28
3.09
1.87
3.15
4.9
1.41
1.68
−1.25
−0.68
5.85
−2.13
−3.53
−0.87
1.86
7.33
6.65
5
3.75
2.43
2.95
2.97
3.39
2.52
0.97
0.87
0.82
1.98
0.85
3.91
3.79
−0.35
−0.55
3.56
6
2.96
2.53
2.67
2.93
3.22
−0.14
−1.29
−2.01
−1.99
−1.16
−3.12
−3.23
−5.08
−5.93
−6.47
6.32
7
3.4
2.8
4.2
3.3
2.5
1.7
1.1
4.9
5.7
4.2
−1.8
−1.2
2.6
4.4
5.4
4.51
8
3.31
2.85
3.31
2.56
3.56
1.2
0.82
1.75
0.57
2.25
−1.44
−1.4
−1.37
−0.22
1.82
2.02
9
3.15
2.35
2.55
2.59
2.38
0.053
−1.567
−2.587
−3.442
−4.784
−0.258
0.931
−0.703
−4.026
−5.323
5.26
10
2.8
2.8
2.2
2.8
1.9
−1
−1.5
−3.5
−3.6
−6.1
−1.1
−4.4
−3.8
−3
−1.3
5.12
11
1.41
3.68
2.97
1.36
0.98
−5.14
−2.77
−2.58
−6.58
−11.48
1.71
7.33
8.03
11.17
10.18
9.77
12
2.76
4.91
3.68
1.84
1.59
−0.4
6.15
7.92
4.56
0.71
0.5
7.2
10.7
14.61
12.5
14.21
13
3.19
1.57
0.81
5.56
1.73
0.38
−3.25
−8.92
0.98
−2.21
2.86
6.37
7.32
21.97
33.14
32.76
14
1.6
2
1.2
1.1
1.7
−5.4
−7.6
−11.7
−15.9
−18.2
−3.2
−3.5
−2.2
−4.4
−3.7
16.00
15
2.17
1.17
5.08
2.48
1.18
−2.61
−7.17
0.47
−0.8
−5.41
1.75
3.45
14.98
22.39
23.82
23.34
16
3.51
2.17
1.69
1.25
4.38
1.2
−0.94
−4.19
−8.43
−3.5
2.67
5.53
5.06
4.53
12.9
11.96
17
1.84
0.39
3.68
2.48
0.85
−4.42
−11.86
−9.11
−9.84
−14.94
−1.08
3.24
12.58
14.59
17.62
16.54
18
1.6
2.8
4.7
2
1.8
−3.9
−4.1
1.8
−0.7
−3.7
3.6
5.6
11.1
18.2
15.8
17.53
19
1.57
1.08
2.03
1.61
2.12
−5.23
−9.91
−11.79
−14.6
−15.93
−2.52
−2.31
−0.49
−0.85
−0.82
15.44
20
1.89
2.88
2.82
2.05
3.65
−3.152
−3.2133
−3.4189
−5.5926
−3.2656
1.4784
3.1056
−0.0038
0.9536
4.7244
5.588The values ofRi are also plotted in Figure 3. From Figure 3, it can be noted that 13th values are out-of-control which lead the experimenter to take action to bring back the process to the in-control state.Figure 3
The proposed chart using the real data.


## 5. Concluding Remarks

In this study, a new generalized range control chart is designed for the Weibull distribution using inverse erf function which has been commonly used to transform the Weibull distributed data to normal data. Some tables for practical use are also given. From comparison study, it is concluded that the proposed control chart is more efficient than the existing charts in detecting quick shift at small-shift levels. The proposed control chart can be applied in the industry for monitoring of the process that follows the Weibull distribution. The proposed control chart using some other sampling scheme and cost models can be considered as future research.



---

*Source: 9453589-2018-11-05.xml*