Journal of Industrial Mathematics

Volume 2013 (2013), Article ID 489285, 6 pages

http://dx.doi.org/10.1155/2013/489285

Research Article

Truncated Life Test Plans for Economic Reliability Based on Four-Parametric Burr Distribution

Ramkumar Balan and Sajana Kunjunni

Department of Statistics, St. Thomas' College, Thrissur, Kerala, India

Received 1 June 2013; Revised 1 October 2013; Accepted 8 October 2013

Academic Editor: Tamer Eren

Copyright © 2013 Ramkumar Balan and Sajana Kunjunni. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Burr distribution is considered as a probability model for the lifetime of products. Reliability test plans are those sampling plans in which items from a lot are put to test to make conclusions on the estimate of life, and hence acceptance or rejection of the submitted lot is done. A test plan designs the termination time of the experiment and the termination number for a given sample size and producer’s risk. Tables and graphs were provided for certain specific values of designs, and it is useful to verify the optimum reliability test plan realized by Burr distributions.

1. Introduction

Reliability study plays a vital role in the quality control, and it can save time and money by realizing early conclusions. If a genuine product (reaching defined life) is rejected on the basis of sample information it is called Type-1 error. On the other hand, if an ingenuine product (not reaching defined life) is accepted by the consumer, then it is called Type-2 error. The decision to accept or reject a lot is subjected to the risks associated with these errors. The procedure is termed as “reliability test plan” or “acceptance sampling based on life test” [1].

Tailed probability distributions are the basis of reliability test plans. These distributions are used to find the reliability sampling plans which will be more economical for the experimenter. Kantam developed a detailed study on life tests based on log-logistic distribution [2]. Rosaiah and Kantam used Inverse Rayleigh distribution to present acceptance sampling [3]. Kantam et al. introduced economic reliability test plan for log-logistic distribution [4], and Aslam and Shahbaz considered generalized exponential distribution to explain economic reliability test plan [5]. Aslam also presented economic reliability test for a generalized Rayleigh Distribution [6]. Rao et al. explained the economic reliability test plan on the basis of Marshall-Olkin extended exponential distribution [7]. Rao et al. considered reliability test plans for type-II exponentiated log-logistic distribution [8]. Mugahal et al. introduced economic reliability group acceptance sampling plans for lifetimes, a Marshall-Olkin extended distribution [9]. Aslam et al. considered generalized exponential distribution to explain time truncated acceptance sampling [10].

It is found that a null hypothesis about scale parameter such as “the scale parameter is greater than or equal to a specified value” is equivalent to establish that the “average life of a product realized by a given scaled density exceeds the specified average life.” Acceptance of this hypothesis by a test procedure is equivalent to getting the sample lifetime echoing the lot quality, and the rejection of hypothesis implies that the lot is a bad.

Ramkumar and Sajana [11] had designed a reliability test plan on four-parametric Burr distribution showing the minimum termination time at a prefixed risk, and this study tries to develop an economic sampling plan under the same assumptions. To verify the efficacy of economic sampling plans, a comparative study of the plans was made by finding the minimum termination time and the acceptance numbers.

General definition of four-parameter Burr distribution and the method of finding the minimum sample size are described in Section 2. The rest of the paper is organized as follows. Section 3 describes the method of finding operating characteristic. An example for illustrating the difference between these two approaches is given in Section 4. Comparative studies between the two sampling plans were discussed in Section 5, and conclusions were drawn and presented in Section 6.

2. The Sampling Plan

Assume that the lifetime of a product follows four-parametric Burr distribution. The cumulative distribution function and probability distribution function of the Burr distribution are given by

Let , then the CDF becomes

Also the probability density function of four-parameter Burr distribution is

The graphical representation of probability density function of four-parameter Burr distribution for specified parameters , , , and is shown in Figure 1.

Figure 1: The graph of the probability density function of four-parameter Burr distribution of parameters , , , and .

Consider a null hypothesis

Let the lifetime of a product follows the above Burr distribution with eventual failures; the above hypothesis becomes the average life of the product. If is accepted on the basis of sample lifetime data, one may conclude that the submitted lot has a better average lifetime than the specified quality. Hence, the lot can be termed as good and accepted.

In reliability test plans based on Burr distribution, the minimum sample size to make a decision about the lot for a given waiting time in terms of (i.e., ), acceptance number , and some risk probability (say ) for a specified of which is the probability of detecting is required.

If the number of failures is less than or equal to the acceptance number in a sample of size , the lot is accepted and hence the probability of accepting the lot:

where.

For , the above probability of acceptance will increase. Therefore, if is a prefixed risk, probability of minimum risk is

For a given and , this is a single inequality in two unknowns and for the assumed parameters and . The above inequality can be solved for with successive values of from zero. The minimum values of which satisfies the inequality (6) are for = 0.01, 0.05, 0.10, and 0.25 and = 0.306, 0.413, 0.521, 0.841, 1.055, 1.483, and 1.911 at , , and , given in Table 1. The minimum value of satisfies the binomial approximations which are obtained for the combination producers risk () and operating ratio .

Table 1: Minimum sample size for the specified ratio , confidence level , acceptance number , , and using binomial approximation.

0.306
0.413
0.521
0.627
0.841
1.055
1.483
1.911

0.75
0
7
4
3
2
1
1
1
1

0.75
1
15
9
6
4
3
2
—
—

0.75
2
22
13
9
7
4
3
—
—

0.75
3
30
17
12
9
6
5
—
—

0.75
4
37
22
15
11
8
6
—
—

0.75
5
44
26
18
13
9
8
—
—

0.75
6
50
30
21
16
11
9
—
—

0.75
7
57
34
23
18
13
10
—
—

0.75
8
64
38
26
20
14
12
—
—

0.75
9
71
42
29
22
16
13
—
—

0.75
10
77
46
32
25
18
14
—
—

0.9
0
12
7
4
3
2
1
1
1

0.9
1
22
12
8
6
4
3
2
2

0.9
2
30
17
12
9
6
4
3
3

0.9
3
38
22
15
11
7
6
4
4

0.9
4
46
27
18
14
9
7
6
5

0.9
5
54
32
22
16
11
9
7
6

0.9
6
61
36
25
19
13
10
8
7

0.9
7
69
41
28
21
15
12
9
8

0.9
8
76
45
31
24
16
13
10
9

0.9
9
84
49
34
26
18
14
12
10

0.9
10
91
54
37
28
20
16
13
12

0.95
0
16
9
6
4
2
2
1
1

0.95
1
27
15
10
7
5
3
2
2

0.95
2
36
21
14
10
7
5
4
3

0.95
3
44
26
17
13
8
6
5
4

0.95
4
53
31
21
16
10
8
6
5

0.95
5
61
36
24
18
12
10
7
6

0.95
6
69
40
27
21
14
11
9
8

0.95
7
77
45
31
23
16
13
10
9

0.95
8
84
50
34
26
18
14
11
10

0.95
9
92
54
37
28
19
15
12
11

0.95
10
99
59
40
31
21
17
13
12

0.99
0
24
14
9
6
4
3
1
1

0.99
1
37
21
14
10
6
5
3
2

0.99
2
47
27
18
13
8
6
4
4

0.99
3
57
33
22
16
11
8
6
5

0.99
4
66
38
26
19
13
10
7
6

0.99
5
75
44
29
22
15
12
8
7

0.99
6
84
49
33
25
17
13
10
8

0.99
7
92
54
37
27
18
14
11
9

0.99
8
100
59
40
30
20
16
12
11

0.99
9
109
64
43
33
22
17
13
12

0.99
10
117
69
47
35
24
19
15
13

Fixing “” and “” be a natural number less than sample size (). The acceptance criteria can be restated as follows.

(1)Put “” items on test.(2)Stop the process if failure occurs before and the lot is rejected.

One may be interested that the probability of acceptance should be as large as possible for a given , and specify “” as a multiple of “” (). Then,

Given the values of (where ), , and (integer), the above inequality can be solved for “” using cumulative probabilities of binomial distribution. Then the values of “” can be used in the cumulative density function for to find the values of . With the choices of () and , the termination ratio could be found satisfying the above inequality. The termination ratios are given in Table 2 for various values of and .

Table 2: Life test termination time as a ratio of scale parameter (/β ₀) for Burr distribution with and .

= 0.05

1
0.10989
0.09147
0.07928
0.07091
0.06463
0.05989
0.05632
0.05196
0.04978

2
0.22796
0.17984
0.15388
0.13652
0.12391
0.11442
0.10656
0.10062
0.09496

3
0.29303
0.22842
0.19398
0.17164
0.15483
0.14292
0.13294
0.12564
0.11873

4
0.33525
0.25927
0.21873
0.19372
0.17483
0.16096
0.14978
0.13998
0.13286

5
0.36551
0.28128
0.23661
0.20887
0.18892
0.17342
0.16182
0.15172
0.14368

6
0.38747
0.29698
0.24976
0.21963
0.19851
0.18285
0.17047
0.16013
0.14999

7
0.40514
0.30944
0.25994
0.22871
0.20673
0.19054
0.17634
0.16642
0.15713

8
0.42012
0.31963
0.26846
0.23654
0.21382
0.19647
0.18266
0.17157
0.16216

9
0.43194
0.32842
0.27569
0.24252
0.21906
0.20119
0.18698
0.17564
0.16574

10
0.44252
0.33537
0.28147
0.24763
0.22352
0.20516
0.19096
0.17914
0.16962

1
0.02010
0.01360
0.01200
0.01098
0.00916
0.00848
0.00791
0.00747
0.00710

2
0.09164
0.06510
0.05642
0.04896
0.04479
0.04107
0.03884
0.03604
0.03422

3
0.14044
0.11112
0.09430
0.08341
0.07556
0.06919
0.06479
0.06069
0.05834

4
0.18949
0.14985
0.12321
0.10834
0.09862
0.09019
0.08455
0.07952
0.07486

5
0.22011
0.17143
0.14521
0.12812
0.11541
0.10727
0.09892
0.09345
0.08813

6
0.24937
0.19142
0.16248
0.14337
0.13001
0.11901
0.11075
0.10499
0.09841

7
0.27425
0.21021
0.17672
0.15591
0.140412
0.12999
0.12033
0.11296
0.10732

8
0.29307
0.22386
0.18802
0.16573
0.14988
0.13776
0.12825
0.12049
0.11429

9
0.31030
0.23574
0.19901
0.17457
0.15773
0.14524
0.13537
0.12726
0.12006

10
0.32198
0.24562
0.20684
0.18268
0.16501
0.15201
0.14098
0.13232
0.13011

3. Operating Characteristic Function

If unknown lifetime of the product deviates from the specified lifetime, it should result in a considerable change in the probability of acceptance of the lot based on the sampling plan. Hence, the probability of acceptance can be regarded as a function of the deviation of specified average lifetime from the true average lifetime. This function is called operating characteristic function of the sampling plan. Specifically if is the cumulative distribution function of the lifetime random variable of the product, “” corresponds to a specified life.

Thus,

where “” corresponds to true but unknown average life. The ratio in the right hand side of the above equation can be taken as a measure of changes between true and specified life. For instance, implies that the true mean life is less than the declared life leading to less acceptance probability or more failure risk. Similarly, implies more acceptance probability or less failure risk. Hence, giving a set of hypothetical values, say , the corresponding acceptance probability for the given sampling plan is determined. For selected sampling plans, the operating characteristic (OC) values were given in Table 3. The OC curves for different values of as given in Table 3 are shown in Figure 2.

Table 3: Operating characteristic values for the proposed plan based on Burr distribution with and for specified values of terminating number and producer’s risk , 0.01, 0.10, and 0.25.

(OC curve 1)
(OC curve 2)
(OC curve 3)
(OC curve 4)

,
, α* = 0.01
,
,

,
,
,
,

0.1
1
1
1
1

0.2
1
0.958
1
1

0.3
0.99
0.612
0.999
1

0.4
0.9
0.258
0.978
0.998

0.5
0.688
0.095
0.862
0.97

0.6
0.451
0.035
0.652
0.872

0.7
0.269
0.014
0.436
0.706

0.8
0.154
0.006
0.272
0.521

0.9
0.087
0.003
0.163
0.362

1.0
0.049
0.001
0.097
0.242

1.1
0.028
0.0006
0.058
0.159

1.2
0.017
0.0003
0.035
0.103

1.3
0.009
0.0002
0.022
0.068

1.4
0.006
0.0001
0.013
0.044

1.5
0.004
0.00005
0.009
0.029

1.6
0.002
0.00004
0.006
0.02

1.7
0.002
0.00002
0.004
0.013

1.8
0.001
0.00001
0.002
0.009

1.9
0.001
0.000009
0.002
0.006

Figure 2: The operating characteristic curve of Table 3.

4. Example

Consider a simulated problem associated with failure times of bulbs whose lifetime follows four-parameter Burr distribution with , , and . This data can be regarded as an ordered sample of size with observations 794, 1284, 1087, 1369, 2206, 1346, 939, 1578, 953, and 1092. Goodness of fit details for this sample is shown in Table 5.

4.1. Illustration of Example Using Reliability Test Plan

Let the required average life time be 1000 hours and the testing time be , leading to ratio = 0.521 with a corresponding sample size and acceptance number obtained from Table 1, for . Therefore, the sampling plan for the above sample data is (, , = 0.521). Based on the data, it is to decide whether the lot is accepted or rejected. Accept the product only if the number of failures before 521 hours is less than or equal to 1.

In the above sample of size 10, there are no failures before the termination time hours and the acceptance number of the plan . Therefore, accept the product.

4.2. Illustration of Example Using Economic Sampling Plan

From the life test termination in Table 2, corresponding to (), for column 5, the value is 0.13652. As the acceptable average lifetime is 1000 hours (given) for four-parameter Burr distribution, then the termination time is obtained as (approximately).

According to the above sampling plan, if the first failure is realized before 137th hour of the test, reject the lot; otherwise, accept the lot. In either case terminating the experiment as soon as the first failure is reached or 137th hour of the test time is realized whichever is earlier. In the case of acceptance, the assurance is that the average life of the submitted products is at least 1000 hours.

There are no failures before 137th hour, and accept the lot (in the example) by this approach. Thus, for both approaches, the sample size, acceptance number (termination number), the risk probability, and the decision about the lot are the same. The decision on the first approach can be reached only at the 521st hour while in the second approach decision is realized at the 137th hour, showing that the economic sampling plan requires less waiting time.

5. Comparative Study

In order to compare the economic sampling plan with reliability test plan, the values are presented for , in Table 4. The entries given in the first row represent economic sampling plan, and the values of the second row show reliability sampling plan. Termination values in the example are taken for , , and from Table 4. For any , , and , the present test plan gives minimum termination time than reliability test plan saving money and time of experiment

Table 4: Proportions of life test termination time for the present sampling plans and reliability test plan with producer’s risk , 0.01.

1
0.109
0.091
0.079
0.071
0.065
0.059
0.056
0.052
0.049
1
0.020
0.013
0.012
0.011
0.009
0.008
0.008
0.007
0.007

1.055
0.792
0.627
0.591
0.521
0.483
0.447
0.413
0.392
1.357
1.055
0.841
0.755
0.627
0.613
0.548
0.521
0.504

2
0.228
0.179
0.154
0.137
0.124
0.114
0.107
0.101
0.094
2
0.092
0.065
0.056
0.049
0.044
0.041
0.038
0.036
0.034

0.985
0.734
0.598
0.521
0.478
0.428
0.406
0.379
0.357
1.214
0.841
0.714
0.627
0.564
0.521
0.484
0.447
0.427

3
0.293
0.228
0.194
0.172
0.155
0.143
0.133
0.126
0.119
3
0.140
0.111
0.094
0.083
0.076
0.069
0.065
0.061
0.058

0.957
0.684
0.566
0.497
0.448
0.413
0.378
0.356
0.334
1.055
0.824
0.676
0.568
0.521
0.478
0.426
0.413
0.379

4
0.335
0.259
0.219
0.194
0.175
0.161
0.149
0.139
0.133
4
0.189
0.149
0.123
0.108
0.099
0.090
0.085
0.079
0.075

0.841
0.665
0.548
0.482
0.431
0.396
0.367
0.339
0.325
1.055
0.749
0.627
0.559
0.493
0.447
0.412
0.394
0.369

5
0.366
0.281
0.237
0.209
0.189
0.173
0.162
0.152
0.144
5
0.220
0.171
0.145
0.128
0.115
0.107
0.098
0.093
0.088

0.841
0.645
0.536
0.468
0.415
0.382
0.356
0.334
0.315
1.055
0.736
0.618
0.523
0.478
0.426
0.397
0.365
0.353

6
0.387
0.297
0.249
0.219
0.199
0.183
0.170
0.160
0.149
6
0.249
0.191
0.162
0.143
0.130
0.119
0.111
0.105
0.098

0.841
0.627
0.521
0.457
0.413
0.376
0.349
0.326
0.309
1.055
0.726
0.594
0.511
0.456
0.417
0.385
0.357
0.334

7
0.405
0.309
0.259
0.229
0.207
0.191
0.176
0.166
0.157
7
0.274
0.210
0.177
0.156
0.140
0.129
0.120
0.113
0.107

0.841
0.627
0.518
0.451
0.401
0.366
0.341
0.322
0.302
0.967
0.714
0.562
0.493
0.446
0.413
0.376
0.354
0.338

8
0.420
0.319
0.268
0.237
0.214
0.196
0.183
0.172
0.162
8
0.293
0.224
0.188
0.166
0.149
0.138
0.128
0.120
0.114

0.841
0.617
0.509
0.445
0.399
0.365
0.338
0.317
0.298
0.932
0.696
0.563
0.494
0.429
0.398
0.364
0.348
0.327

9
0.432
0.328
0.276
0.243
0.219
0.201
0.187
0.176
0.166
9
0.310
0.236
0.199
0.175
0.158
0.145
0.135
0.127
0.120

0.841
0.612
0.502
0.439
0.391
0.359
0.331
0.313
0.294
0.928
0.662
0.543
0.477
0.421
0.394
0.367
0.339
0.316

10
0.443
0.335
0.281
0.248
0.224
0.205
0.191
0.179
0.169
10
0.322
0.246
0.207
0.183
0.165
0.152
0.141
0.132
0.130

0.834
0.602
0.498
0.434
0.386
0.355
0.331
0.308
0.291
0.919
0.667
0.534
0.473
0.429
0.384
0.358
0.335
0.321

Table 5: Goodness of fit details.

Kolmogorov-Smirnov

Sample size
10

Statistic
0.19705

P value
0.76382

Rank
33

0.2
0.1
0.05
0.02
0.01

Critical value
0.3226
0.36866
0.40925
0.45662
0.48893

Reject?
No
No
No
No
No

Anderson-Darling

Sample size
10

Statistic
0.43355

Rank
30

0.2
0.1
0.05
0.02
0.01

Critical value
1.3749
1.9286
2.5018
3.2892
3.9074

Reject?
No
No
No
No
No

Chi-squared

Deg. of freedom
1

Statistic
0.07696

P value
0.78146

Rank
28

0.2
0.1
0.05
0.02
0.01

Critical value
1.6424
2.7055
3.8415
5.4119
6.6349

Reject?
No
No
No
No
No

6. Conclusion

This paper presents the economic test plan for lifetime of products following four-parameter Burr distribution. The proposed plan suggests a way for achieving optimum reliability test plan for tailed Burr distribution. Comparative study affirms the conclusion that the present test plan works with minimum termination time (waiting time) ensuring minimum experimental cost at more or less same sampling design.

The continuous improvement and review of acceptance sampling plan are important to improve the quality of the products and to ensure customer satisfaction.

Acknowledgment

The authors hereby extent sincere gratitude to the reviewers for their valuable suggestions for rectification of errors as well as for modification of the format.

References

A. J. Ducan, Quality Control and Industrial Statistics, Irwin Home Wood III, Toronto, Canada, 5th edition, 1986.

R. R. L. Kantam, K. Rosaiah, and G. S. Rao, “Acceptance sampling based on life tests: log-logistic model,” Journal of Applied Statistics, vol. 28, no. 1, pp. 121–128, 2001.

K. Rosaiah and R. R. L. Kantam, “Acceptance sampling based on the inverse Rayleigh distribution,” Economic Quality Control, vol. 20, no. 2, pp. 277–286, 2005.

R. R. L. Kantam, G. Srinivasa Rao, and B. Sriram, “An economic reliability test plan: log-logistic distribution,” Journal of Applied Statistics, vol. 33, no. 3, pp. 291–296, 2006.

M. Aslam and M. Q. Shahbaz, “Economic reliability test plan, generalized exponential distribution,” Journal of Statistics, vol. 14, pp. 52–59, 2007.

M. Aslam, “Economic reliability test generalized rayleigh distribution,” Journal of Statistics, pp. 26–35, 2008.

G. S. Rao, M. E. Ghitany, and R. R. L. Kantam, “An economic reliability test plan for Marshall-Olkin extended exponential distribution,” Applied Mathematical Sciences, vol. 5, no. 3, pp. 103–112, 2011.

G. S. Rao, R. R. L. Kantam, K. Rosaiah, and S. V. Prasad, “Reliability test plans for type-II exponentiated log-logistic distribution,” Journal of Reliability and Statistical Studies, vol. 5, no. 1, pp. 55–64, 2012.

A. Mugahal, M. Aslam, J. Hussain, and A. Rehman, “Economic relibility group acceptance sampling plans for lifetimes a Marshall-Olkin extended distribution,” Middle Eastern Finance and Economics, no. 7, 2010.

M. Aslam, D. Kundu, and M. Ahmad, “Time truncated acceptance sampling plans for generalized exponential distribution,” Journal of Applied Statistics, vol. 37, no. 4, pp. 555–566, 2010.

T. B. Ramkumr and O. K. Sajana, “Reliability test plan based on Burr distribution,” International Journal of Mathematics and Computer Applications Research, vol. 1, no. 2, pp. 28–40, 2011.