In this paper, acceptance sampling plans are developed for the odds exponential log logistic distribution (OELLD) introduced by Rosaiah et al. [1K. Rosaiah, G.S. Rao, K. Kalyani, and D.C.U. Sivakumar, "Odds exponential log logistic distribution: Properties and estimation", J. Math. Stat., vol. 12, no. 4, 2016. (in press).] based on lifetime percentiles when the life test is truncated at a predetermined time. The minimum sample size necessary to ensure the specified lifetime percentile is obtained under a given customer’s risk. The operating characteristic values of the sampling plans as well as the producer’s risk are presented. One example with real data set is also given as an illustration.
Open Peer Review Details | |||
---|---|---|---|
Manuscript submitted on 24-05-2016 |
Original Manuscript | A New Acceptance Sampling Plans Based on Percentiles for Odds Exponential Log Logistic Distribution |
In most of the statistical quality control experiment, it is not possible to perform 100% inspection due to various reasons. The acceptance sampling plans are concerned with accepting or rejecting a submitted lot of a large size of products on the basis of the quality of the products inspected in a sample taken from the lot. An acceptance sampling plan is a specified plan that establishes the minimum sample size to be used for testing. In most acceptance sampling plans for a truncated life test, the major issue is to determine the sample size from a lot under consideration. Traditionally, when the life test indicates that the mean life of products exceeds the specified one, the lot of products is accepted, otherwise it is rejected. For the purpose of reducing the test time and cost, a truncated life test may be conducted to determine the smallest sample size to ensure a certain mean life time/percentile life time of products when the life test is terminated at a time, t0 and the number of failures observed does not exceed a given acceptance number c.
A common practice in life testing is to terminate the life test by a predetermined time t0 and note the number of failures. One of the objectives of these experiments is to set a lower confidence limit on the mean life/percentile lifetime. It is then to establish a specified mean life with a given probability, which provides protection to consumers. The test may be terminated before the time is reached or when the number of failures exceeds the acceptance number c in which case the decision is to reject the lot.
In the literature of acceptance sampling, there are many methods of testing. Epstein [2B. Epstein, "Truncated life tests in the exponential case", Ann. Math. Stat., vol. 25, pp. 555-564, 1954.
[http://dx.doi.org/10.1214/aoms/1177728723] ] was the first who considered truncated life tests in the exponential distribution. Truncated life tests are considered by many authors for various distributions: for example Sobel and Tischendrof [3M. Sobel, and J.A. Tischendrof, "Acceptance sampling with new life test Objectives", In: Proceedings of 5th National Symposium on Reliability and Quality Control, Philadelphia, Pennsylvania, 1959, pp. 108-118.], Gupta and Groll [4S.S. Gupta, and P.A. Groll, "Gamma distribution in acceptance sampling based on life tests", J. Am. Stat. Assoc., vol. 56, pp. 942-970, 1961.
[http://dx.doi.org/10.1080/01621459.1961.10482137] ], Kantam and Rosaiah [5R.R.L. Kantam, and K. Rosaiah, "Half logistic distribution in acceptance sampling based on life tests", IAPQR Trans., vol. 23, no. 2, pp. 117-125, 1998.], Baklizi and El Masri [6A. Baklizi, and A.E.K. El Masri, "Acceptance sampling based on truncated life tests in the Birnbaum Saunders model", Risk Anal., vol. 24, no. 6, pp. 1453-1457, 2004.
[http://dx.doi.org/10.1111/j.0272-4332.2004.00541.x] [PMID: 15660603] ], Tsai and Wu [7T.R. Tsai, and S.J. Wu, "Acceptance sampling based on truncated life tests for generalized Rayleigh distribution", J. Appl. Stat., vol. 33, pp. 595-600, 2006.
[http://dx.doi.org/10.1080/02664760600679700] ], Balakrishnan et al. [8N. Balakrishnan, V. Leiva, and J. Lopez, "Acceptance sampling plans from truncated life tests based on the generalized Birnbaum-Saunders distribution", Commun. Stat. Simul. Comput., vol. 36, pp. 643-656, 2007.
[http://dx.doi.org/10.1080/03610910701207819] ], Kantam et al. [9R.R.L. Kantam, K. Rosaiah, and G.S. Rao, "Acceptance sampling based on life tests: Log-logistic models", J. Appl. Stat., vol. 28, pp. 121-128, 2001.
[http://dx.doi.org/10.1080/02664760120011644] ], Rao et al. [10G.S. Rao, M.E. Ghitany, and R.R.L. Kantam, "Acceptance sampling plans for Marshall-Olkin extended Lomax distribution", Int. J. Appl. Math., vol. 21, pp. 315-325, 2008.] and the references there in.
The above authors considered the design of acceptance sampling plans based on the population mean under a truncated life test. Gupta [11S.S. Gupta, "Life test sampling plans for normal and log-normal distribution", Technometrics, vol. 4, pp. 151-175, 1962.
[http://dx.doi.org/10.1080/00401706.1962.10490002] ] suggested that for a skewed distribution the median represents a better quality parameter than the mean. On the other hand, for a symmetric distribution, mean is preferable to use as a quality parameter. Lio et al. [12Y.L. Lio, T.-R. Tsai, and S.-J. Wu, "Acceptance sampling plan based on the truncated life test in the Birnbaum Saunders distribution for percentiles", Commun. Stat. Simul. Comput., vol. 39, pp. 119-136, 2009.
[http://dx.doi.org/10.1080/03610910903350508] , 13Y.L. Lio, T.-R. Tsai, and S.-J. Wu, "Acceptance sampling plans from truncated life tests based on the Burr type XII percentiles", J. Chinese Inst. Ind. Eng., vol. 27, no. 4, pp. 270-280, 2010.
[http://dx.doi.org/10.1080/10170661003791029] ] considered acceptance sampling plans based on the truncated life tests to Birnbaum-Saunders distribution and Burr type XII for percentiles and they proposed that the acceptance sampling plans based on mean may not satisfy the requirement of engineering on the specific percentile of strength or breaking stress. When the quality of a specified low percentile is concerned, the acceptance sampling plans based on the population mean could pass a lot which has the low percentile below the required standard of consumers. Furthermore, a small decrease in the mean with a simultaneous small increase in the variance can result in a significant downward shift in small percentiles of interest. This means that a lot of products could be accepted due to a small decrease in the mean life after inspection. But the material strengths of products are deteriorated significantly and may not meet the consumer’s expectation. Therefore, engineers pay more attention to the percentiles of lifetime than the mean life in life testing applications. Moreover, most of the employed life distributions are not symmetric. Actually, percentiles provide more information regarding a life distribution than the mean life does. When the life distribution is symmetric, the 50th percentile or the median is equivalent to the mean life. Hence, developing acceptance sampling plans based on percentiles of a life distribution can be treated as a generalization of developing acceptance sampling plans based on the mean life of items. Balakrishnan et al. [8N. Balakrishnan, V. Leiva, and J. Lopez, "Acceptance sampling plans from truncated life tests based on the generalized Birnbaum-Saunders distribution", Commun. Stat. Simul. Comput., vol. 36, pp. 643-656, 2007.
[http://dx.doi.org/10.1080/03610910701207819] ] proposed the acceptance sampling plans could be used for the quantiles and derived the formulae whereas Lio et al. [12Y.L. Lio, T.-R. Tsai, and S.-J. Wu, "Acceptance sampling plan based on the truncated life test in the Birnbaum Saunders distribution for percentiles", Commun. Stat. Simul. Comput., vol. 39, pp. 119-136, 2009.
[http://dx.doi.org/10.1080/03610910903350508] , 13Y.L. Lio, T.-R. Tsai, and S.-J. Wu, "Acceptance sampling plans from truncated life tests based on the Burr type XII percentiles", J. Chinese Inst. Ind. Eng., vol. 27, no. 4, pp. 270-280, 2010.
[http://dx.doi.org/10.1080/10170661003791029] ] developed for the acceptance sampling plans for any other percentiles of the Birnbaum-Saunders (BS) and Burr Type XII models. They have developed the acceptance sampling plans for percentile by replacing the scale parameter by the 100qth percentile. Rao and Kantam [14G.S. Rao, and R.R.L. Kantam, "Acceptance sampling plans from truncated life tests based on the log-logistic distributions for percentiles", Econ. Qual. Control, vol. 25, no. 2, pp. 153-167, 2010.
[http://dx.doi.org/10.1515/eqc.2010.008] ] and Rao et al. [15G.S. Rao, R.R.L. Kantam, K. Rosaiah, and J. Pratapa Reddy, "Acceptance sampling plans for percentiles based on the inverse Rayleigh distribution", Electron. J. App. Stat. Anal., vol. 5, no. 2, pp. 164-177, 2012.] developed acceptance sampling plans from truncated life tests based on the log-logistic and inverse Rayleigh distributions for Percentiles. Rao et al. [16B. Rao, M.Ch. Priya, and R.R.L. Kantam, "Acceptance sampling plans for percentiles assuming the linear failure rate distribution", Econ. Qual. Control, vol. 29, no. 1, pp. 1-9, 2013.] developed acceptance sampling plans for percentiles assuming linear failure rate distribution. Balamurali et al. [17S. Balamurali, M. Aslam, and M.S. Fallah Nezhad, "An acceptance sampling plan under frechet distribution assuring median life", Res. J. Appl. Sci., Eng. Technol., vol. 6, no. 24, pp. 4519-4523, 2013.] developed acceptance sampling plans based on median life for Fréchet distribution. These reasons motivate to develop acceptance sampling plans based on the percentiles, since odds exponential log logistic distribution (OELLD) is a skewed distribution, we prefer to use the percentile point as the quality parameter, and it will be denoted by. The rest of the paper is organized as follows. In Section 2, we describe concisely the OELLD. The design of proposed acceptance sampling plan for lifetime percentiles under a truncated life test is discussed in Section 3. In Section 4, we present the description of the proposed plan and obtain the necessary results. An example with real data set is also given as an illustration. Finally, conclusions are made in Section 5.
In this section, we provide a brief summary about the odds exponential log logistic distribution. The OELLD was introduced and studied quite extensively by Rosaiah et al. [1K. Rosaiah, G.S. Rao, K. Kalyani, and D.C.U. Sivakumar, "Odds exponential log logistic distribution: Properties and estimation", J. Math. Stat., vol. 12, no. 4, 2016. (in press).]. The probability density function (pdf) and cumulative distribution function (cdf) of OELLD respectively are given as follows:
(1) |
(2) |
where, λ, σ are the scale parameters and θ respectively shape parameters.
The 100q-th quantile of the OELLD is given as:
(3) |
Hence, for the fixed values of λ = λ 0 and θ = θ 0 , the quantile tq given in Equation (3) is the function of scale parameter σ = σ 0, that is , where
(4) |
Note that σ 0 also depends on λ 0 = θ 0, to build up acceptance sampling plans for the OELLD ascertain tq ≥ tq 0, equivalently that σ exceeds σ 0.
In this paper, we obtain the minimum sample size necessary to ensure a percentile life when the life test is terminated at a pre-assigned time t 0q and when the observed number of failures does not exceed a given acceptance number c. The decision procedure is to accept a lot only if the specified percentile of lifetime is established with a pre-assigned high probability α, which provides protection to the consumer. The life test experiment gets terminated at the time at which (c + 1)th failure is observed or at quantile time tq, whichever is earlier. Acceptance sampling plan based on truncated life tests for the OELLD using average life was discussed by Rosaiah et al. [1K. Rosaiah, G.S. Rao, K. Kalyani, and D.C.U. Sivakumar, "Odds exponential log logistic distribution: Properties and estimation", J. Math. Stat., vol. 12, no. 4, 2016. (in press).].
The probability of accepting lot based on the number of failures from all groups under a truncated life test at the test time schedule is given by
(5) |
Where n is the sample size, c is the acceptance number, and p is the probability of getting a failure within the life test schedule t 0,. If the product lifetime follows an OELLD, then p = F(t 0). Usually, it would be convenient to determine the experiment termination time t 0 as t 0 = δq 0 tq 0 for a constant δq 0 and the targeted 100q-th lifetime percentile, t q 0. Suppose tq is the true 100q-th lifetime percentile. Then, p can be rewritten as
(6) |
In order to find the design parameters of the proposed acceptance sampling plan, we prefer the approach based on two points on the OC curve by considering the producer’s and consumer’s risks. In our approach, the quality level is measured through the ratio of its percentile lifetime to the lifetime, tq/tq0. These percentile ratios are very helpful for the producer to enhance the quality of products. From the producer’s perspective, the probability of lot acceptance should be at least 1 - α at the acceptable reliability level (ARL), p1. So the producer demands that a lot should be accepted at various levels, say tq/tq0 = 2, 4, 6, 8, 10 in Equation (5). On the other hand, from the consumer’s viewpoint, the lot rejection probability should be at most β at the lot tolerance reliability level (LTRL), p2. In this way, the consumer considers that a lot should be rejected when tq/tq0 = 1. From Equation (5), we have
(7) |
(8) |
where p1 and p2 are given by
(9) |
The plan parametric quantities for distinct values of parameters λ are θ constructed. Given the producer’s risk α = 0.5 and termination time schedule t 0 = δq 0tq 0 with δq 0 = 1.0,1.5,2.5,2.0,3.0, the two parameters of the proposed sampling plan under the truncated life test at the pre-specified time, t 0, with λ = 2 and θ = 0.1,1.0,2.0 are obtained according to the consumer’s confidence levels β = 0.25,0.10,0.05,0.01 for 50th percentile and the operating characteristic values are also calculated and the results are presented in Tables 1-3. The plan parameters are presented in Tables 1-3 for λ = 2 and θ = 1.0,1.5,2.0 with 50th percentiles, whereas Table 4 shows the plan parameters for and are the MLE’s from the data given in Subsection 4.2 at 50th percentile. We noticed from Tables 1-4 that the percentile ratio increases, the sample size n decreases for all the cases. The sampling plans for OELLD were developed by Rosaiah et al. [1K. Rosaiah, G.S. Rao, K. Kalyani, and D.C.U. Sivakumar, "Odds exponential log logistic distribution: Properties and estimation", J. Math. Stat., vol. 12, no. 4, 2016. (in press).].
Assume that the producer wants to implement a single sampling plan for assuring that the 50th percentile life of the products under inspection is at least 1000 hours when = β 0.10 at the percentile ratio tq/tq0. He wants to run this experiment for 1000 hours. From the past data, if it is observed that the lifetime of the product follows OELLD with λ = 2 and θ = 2. The optimal plan from Table 3 for specified requirements such as β = 0.10, λ = 2, θ = 2, tq/tq0 = 2 and δq = 1.0 is obtained as n = 17 and c = 5 with the acceptance probability is 0.9587. To the best our knowledge, most of the life testing with single sampling plans for various life time distributions available in the literature is based on one point on the OC curve approach for assuring mean or percentile life time (except Balamurali et al. [17S. Balamurali, M. Aslam, and M.S. Fallah Nezhad, "An acceptance sampling plan under frechet distribution assuring median life", Res. J. Appl. Sci., Eng. Technol., vol. 6, no. 24, pp. 4519-4523, 2013.]). But in this study, we have designed sampling plans based on two-points on the OC curve approach for assuring percentile life time of the products under OELLD.
In this section, we present the application of the proposed OELLD distribution considered by Lemonte [18A.J. Lemonte, "The beta log-logistic distribution", Braz. J. Probab. Stat., vol. 28, no. 3, pp. 313-332, 2014.
[http://dx.doi.org/10.1214/12-BJPS209] ] for a real data set to illustrate its potentiality. The following real data set corresponds to an uncensored data set from Nichols and Padgett [19M.D. Nichols, and W.J. Padgett, "A Bootstrap control chart for Weibull percentiles", Qual. Reliab. Eng. Int., vol. 22, pp. 141-151, 2006.
[http://dx.doi.org/10.1002/qre.691] ] on breaking stress of carbon fibres (in Gba):
The goodness of fit for our model by plotting the superimposed for the data shows that the OELLD is a good fit in the Figure and also goodness of fit is emphasized with Q-Q plot, displayed in the Fig. (1). The maximum likelihood estimates of the two-parameters of OELLD for the breaking stress of carbon fibres are and and the Kolmogorov-Smirnov test we found that the maximum distance between the data and the fitted of the OELLD is 0.0604 with p-value is 0.8582. Therefore, the three-parameter OELLD provides reasonable fit for the breaking stress of carbon fibres.
Fig. (1) Estimated density and Q-Q plot for OELLD. |
Suppose that it is desired to develop the single acceptance sampling plan to satisfy that the 50th percentile lifetime is greater than breaking stress of carbon fibres 0.35 through the experiment to be completed by breaking stress of carbon fibres/measurement 0.35. Let us fix that the consumer's risk is at 25% when the true 50th percentile is breaking stress of carbon fibres 0.35 and the producer's risk is 5% when the true 50th percentile is breaking stress of carbon fibres 0.70. Since and , the consumer's risk is 25%, δq 0 = 1.0 and tq/tq 0 = 2, the minimum sample size and acceptance number given by n = 7 and c = 2 from Table 4. Thus the design can be implemented as follows. Select a sample of 7 runoff amounts, we will accept the lot when no more than two failures occurs before breaking stress of carbon fibres 0.70. According to this plan, the breaking stress of carbon fibres could have been accepted because there is only one failure before the termination time, breaking stress of carbon fibres 0.70.
In this article, we developed the single acceptance sampling plans based on the OELLD percentiles, when the life test is truncated at a pre-fixed time. To ensure that the life quality of products exceeds a specified one in terms of the percentile lifetime, the acceptance sampling plans based on percentiles can be used. We have designed sampling plans based on two-points on the OC curve approach for assuring percentile lifetime of the products. Some tables are provided for practical use in industry and also proposed plan illustrated with real data set. We fitted the proposed OELLD curve for the above data which is shown in Fig. (1).
The authors confirm that this article content has no conflict of interest.
Declared none.
[1] | K. Rosaiah, G.S. Rao, K. Kalyani, and D.C.U. Sivakumar, "Odds exponential log logistic distribution: Properties and estimation", J. Math. Stat., vol. 12, no. 4, 2016. (in press). |
[2] | B. Epstein, "Truncated life tests in the exponential case", Ann. Math. Stat., vol. 25, pp. 555-564, 1954. [http://dx.doi.org/10.1214/aoms/1177728723] |
[3] | M. Sobel, and J.A. Tischendrof, "Acceptance sampling with new life test Objectives", In: Proceedings of 5th National Symposium on Reliability and Quality Control, Philadelphia, Pennsylvania, 1959, pp. 108-118. |
[4] | S.S. Gupta, and P.A. Groll, "Gamma distribution in acceptance sampling based on life tests", J. Am. Stat. Assoc., vol. 56, pp. 942-970, 1961. [http://dx.doi.org/10.1080/01621459.1961.10482137] |
[5] | R.R.L. Kantam, and K. Rosaiah, "Half logistic distribution in acceptance sampling based on life tests", IAPQR Trans., vol. 23, no. 2, pp. 117-125, 1998. |
[6] | A. Baklizi, and A.E.K. El Masri, "Acceptance sampling based on truncated life tests in the Birnbaum Saunders model", Risk Anal., vol. 24, no. 6, pp. 1453-1457, 2004. [http://dx.doi.org/10.1111/j.0272-4332.2004.00541.x] [PMID: 15660603] |
[7] | T.R. Tsai, and S.J. Wu, "Acceptance sampling based on truncated life tests for generalized Rayleigh distribution", J. Appl. Stat., vol. 33, pp. 595-600, 2006. [http://dx.doi.org/10.1080/02664760600679700] |
[8] | N. Balakrishnan, V. Leiva, and J. Lopez, "Acceptance sampling plans from truncated life tests based on the generalized Birnbaum-Saunders distribution", Commun. Stat. Simul. Comput., vol. 36, pp. 643-656, 2007. [http://dx.doi.org/10.1080/03610910701207819] |
[9] | R.R.L. Kantam, K. Rosaiah, and G.S. Rao, "Acceptance sampling based on life tests: Log-logistic models", J. Appl. Stat., vol. 28, pp. 121-128, 2001. [http://dx.doi.org/10.1080/02664760120011644] |
[10] | G.S. Rao, M.E. Ghitany, and R.R.L. Kantam, "Acceptance sampling plans for Marshall-Olkin extended Lomax distribution", Int. J. Appl. Math., vol. 21, pp. 315-325, 2008. |
[11] | S.S. Gupta, "Life test sampling plans for normal and log-normal distribution", Technometrics, vol. 4, pp. 151-175, 1962. [http://dx.doi.org/10.1080/00401706.1962.10490002] |
[12] | Y.L. Lio, T.-R. Tsai, and S.-J. Wu, "Acceptance sampling plan based on the truncated life test in the Birnbaum Saunders distribution for percentiles", Commun. Stat. Simul. Comput., vol. 39, pp. 119-136, 2009. [http://dx.doi.org/10.1080/03610910903350508] |
[13] | Y.L. Lio, T.-R. Tsai, and S.-J. Wu, "Acceptance sampling plans from truncated life tests based on the Burr type XII percentiles", J. Chinese Inst. Ind. Eng., vol. 27, no. 4, pp. 270-280, 2010. [http://dx.doi.org/10.1080/10170661003791029] |
[14] | G.S. Rao, and R.R.L. Kantam, "Acceptance sampling plans from truncated life tests based on the log-logistic distributions for percentiles", Econ. Qual. Control, vol. 25, no. 2, pp. 153-167, 2010. [http://dx.doi.org/10.1515/eqc.2010.008] |
[15] | G.S. Rao, R.R.L. Kantam, K. Rosaiah, and J. Pratapa Reddy, "Acceptance sampling plans for percentiles based on the inverse Rayleigh distribution", Electron. J. App. Stat. Anal., vol. 5, no. 2, pp. 164-177, 2012. |
[16] | B. Rao, M.Ch. Priya, and R.R.L. Kantam, "Acceptance sampling plans for percentiles assuming the linear failure rate distribution", Econ. Qual. Control, vol. 29, no. 1, pp. 1-9, 2013. |
[17] | S. Balamurali, M. Aslam, and M.S. Fallah Nezhad, "An acceptance sampling plan under frechet distribution assuring median life", Res. J. Appl. Sci., Eng. Technol., vol. 6, no. 24, pp. 4519-4523, 2013. |
[18] | A.J. Lemonte, "The beta log-logistic distribution", Braz. J. Probab. Stat., vol. 28, no. 3, pp. 313-332, 2014. [http://dx.doi.org/10.1214/12-BJPS209] |
[19] | M.D. Nichols, and W.J. Padgett, "A Bootstrap control chart for Weibull percentiles", Qual. Reliab. Eng. Int., vol. 22, pp. 141-151, 2006. [http://dx.doi.org/10.1002/qre.691] |