Bayesian Test for Lifetime Performance Index of Ailamujia Distribution Under Squared Error Loss Function
Lanping Li
Department of Basic Subjects, Hunan University of Finance and Economics, Changsha, China
Email address:
To cite this article:
Lanping Li. Bayesian Test for Lifetime Performance Index of Ailamujia Distribution Under Squared Error Loss Function. Pure and Applied Mathematics Journal. Vol. 5, No. 6, 2016, pp. 181-185. doi: 10.11648/j.pamj.20160506.12
Received: October 5, 2016; Accepted: October 14, 2016; Published: November 7, 2016
Abstract: The aim of this paper is to estimate the lifetime performance index of Ailamujia distribution. A Bayesian test procedure is developed under squared error loss function. Firstly, Bayes estimation of life performance index is derived, then a Bayesian test procedure for lifetime performance index using posteriori probability ratio test method. Finally, an example is used to illustrate the effectiveness and feasibility of the method.
Keywords: Bayes Estimation, Lifetime Performance Index, Bayes Test, Ailamujia Distribution, Squared Error Loss Function
1. Introduction
For manufacture industry, process capability index is an effective and convenient tool for quality assessment. Many process capability indices have been put forward. For example, Juran [1] introduced the first process capability index , Kane [2] proposed the process capability index , which are the two most commonly used indices, Boyles [3] and Pearn et al. [4] introduce two more-advanced indices Cpm and Cpmk respectively. Those four process capability indices have been defined explicitly as:
(i) ,
where USL and LSL are the upper and the lower specification limits and is the process standard deviation.
(ii) ,
where is the process mean.
(iii) ,
where is the target value.
(iv)
The statistical inferences of process capability indices have drawn great attention by many authors. For example, Shiau et al. [5] applied Bayesian method to the estimation of Cpm and Cpk under the restriction that the process mean equals to the midpoint of the two specification limits. Pearn and Wu [6] studied the test of Cpk for general situation without restriction on the process mean based on Bayesian approach. Chen and Hsu [7] proposed a test for Cpk that is asymptotically equivalent to the likelihood ratio test. Baral and Anis [8] developed method of generalized confidence interval to measure the process capability index Cpm in presence of measurement errors. Macintyre [9] using Bayes approach to estimate the process capability indices and system availability for the inverse Rayleigh lifetime distribution.
To assess the products with the larger-the-better type of the quality characteristics, Montgomery [10] proposed a special unilateral specification process capability index, named as lifetime performance index ,
(1)
where is the lower bound of the specifications.
The statistical analysis of lifetime performance index for products whose lifetime distributed various distributions have been widely studied. For example, Wu et al. [11] discussed the maximum likelihood estimation, minimum variance unbiased estimation lifetime performance index of Rayleigh distribution product under progressively type II censored test. Lee et al. [12] derived maximum likelihood estimator and constructed a hypothesis testing procedure for lifetime performance index based on censored samples, which products’ lifetime comes from the normal distribution but sample data modeled by fuzzy numbers. Lee et al. [13] studied the Bayesian estimation and Bayesian testing procedures for lifetime performance index under squared error loss function in which products’ lifetime distributed with Rayleigh distribution. Liu and Ren [14] studied the Bayesian estimation and Bayesian test of lifetime performance index for exponential product under progressively type II censored samples.
The statistical inferences about various lifetime distributions, such as exponential distribution, Weibull distribution, normal distribution and Pareto distribution, etc. have been studied a lot [15-20]. Recent years, many new distributions are proposed for various engineer application. Ailamujia (Эрланга) distribution is one of these distribution proposed by Lv et al. [21]. Pan et al. [22] studied the interval estimation and hypothesis test of Ailamujia distribution based on small sample. Long [23] studied the Bayes estimation of Эрлангa distribution under type-II censored samples based on three different prior distributions. Li [24] discussed the minimax estimation of the parameter ofЭрлангa distribution under a non-informative prior distribution with three different loss functions.
Assume that is the products’s lifetime, and it distributed the Ailamujia distribution whose probability density function (pdf) and cumulative distribution function (cdf) respectively:
(2)
(3)
Here, is the unknown parameter.
In this paper, we study the estimation of lifetime performance index using Bayesian approach, and the Bayes test procedure of will also be constructed. Section 2 introduces some properties of the lifetime of product with the Ailamujia distribution. Moreover, the relationship between and conforming rate is also discussed. Furthermore, the Bayesian estimator of based on the conjugate Gamma prior distribution is also obtained under squared error loss function. A new Bayesian hypothesis testing procedure is developed in Section 3, and a practical example is given in Section 4. Finally, a conclusion is given in Section 5.
2. Lifetime Performance Index
Let be the lifetime of such a product whose lifetime distribution is Ailamujia distribution with pdf (2).
It is easily to verify that the process mean and the process standard deviation .
Then the parameter is often called the mean time. Then the lifetime performance index of Ailamujia distribution can be rewritten as follows
(4)
The failure rate function is defined by
(5)
The derivative of the failure rate function with respect to is
(6)
From Eq. (6), we see that is always strictly bigger than zero. Then the failure rate function is strictly increasing function with respect to the parameter . Thus we can see that the smaller , i.e. the larger the mean , the smaller the failure rate and larger the lifetime performance index . Therefore, the lifetime performance index is a reasonably and accurately index for assess the performance of products.
Moreover, conforming rate of the product can be defined as a probability of the lifetime of a product exceeding the lower specification limit , that is
(7)
It is easily to verify that the conforming rateand lifetime performance index have a strictly increasing relationship. The relationships of and can be calculated using Matlab software for a list of various values. For the values which are not listed in Table 1, the conforming rate can be obtained through interpolation. Tong et al. [25] pointed out that the conforming rate can be calculated by dividing the number of conforming products by total number of products sampled, while Montgomery in 1985 suggested that the sample size must be large to accurately estimation. However, a large sample size is usually not practical considering with the perspective of cost. Since there exist a one-to-one mathematical relationship between the conforming rate and the lifetime performance index . Then lifetime performance index can be a flexible and effective tool for estimating the conforming rate .
3. Estimation
This section will discuss the maximul likelihood estimation and Bayes estimation of lifetime performance index of Ailamujia distribution
3.1. Maximum Likelihood Estimation
Let represent the lifetime of sample from the Ailamujia distribution with pdf (2), is the observation of and is the observation of . Then the likelihood function corresponding to pdf (2) is given by
(8)
Then the log- likelihood function can be obtained as follows:
The maximum likelihood estimator of can be easily derived from the log-likelihood equation
.
The maximum likelihood estimator of is
(9)
Where .
Then by invariance of maximum likelihood estimation, we can get the maximum likelihood estimator of lifetime performance indexas follows:
(10)
We can also easily show that is a random variable distributed with the Gamma distribution, which has the following probability density function:
(11)
3.2. Bayes Estimation
In this section, we shall discuss the Bayes estimation of the lifetime performance index of Ailamujia distribution with pdf (2) using Bayesian approach. Squared error loss function is one of the most important loss functions in Bayesian statistical analysis, and the formula of it is
(12)
Assume that the conjugate prior distribution of is the Gamma prior distribution, with pdf
(13)
Where are two prior hyper parameters.
Combining the likelihood function (8) with the prior probability density function in Eq. (13), the posterior pdf of can be derived using Bayes Theorem as follows
(14)
That is
(15)
Then under the squared error loss function (12), Bayes estimator of is the posterior mean, i.e.
(16)
Further the Bayes estimator of is
(17)
4. Bayes Test of Lifetime Performance Index
This section will construct a Bayesian testing procedure to assess whether the lifetime performance index adheres to the required level. Assume that the required value of lifetime performance is larger than the target value. . First, we establish the following hypothesis:
(18)
Then the new proposed Bayesian testing procedure of is as follows:
Step 1. Determine the lower lifetime limit and sample size .
Step 2. Calculate the Bayesian estimator
(19)
Where .
Step 3. Calculate the posterior probability odds ratio
(20)
where ,
and is the observation of .
Step 4. The decision rules are provided as follows:
If and , we reject to the null hypothesis , then it is concluded that the lifetime performance index or conforming rate of the products meets the required level;
If and , we accept the null hypothesis , then it is concluded that the lifetime performance index or conforming rate of the products does not meet the required level.
5. Numerical Example
To illustrate the practicability and feasibility of the proposed testing method, a Monte Carlo simulation is used to generate a sample of Ailamujia distribution with and n=20. The data set is:
5.4833, 18.6857, 3.2724, 7.2842, 3.9008, 19.4737, 7.0814, 4.2670, 6.0806, 22.5402, 7.6014, 14.8026, 14.3133, 1.9649, 7.7563, 4.6101, 4.1231, 9.7974, 8.4994, 3.6257
Now we give the steps of the proposed Bayesian testing procedure about as follows:
Step1. Calculate and here we assume the lower lifetime limit =4.078. To deal with the lifetime performances, the conforming rate is required to exceed 0.80. According to Eq. (7), the value of is required to exceed 0.8375. Thus, the target value of performance index is set at =0.8375, then we establish the following testing hypothesis
.
Step 2. Under squared loss function, we get the Bayesian estimate;
Step 3. Suppose that the prior parameter values and, then
.
Therefore, the posterior odds ratio
.
Step 4. Obviously, and , then we can accept the null hypothesis . That is, we conclude that the lifetime performance index does not meet the required level.
6. Conclusions
Process capability indices are well effective tools to assess the performance and potentiality of their process and widely employed by manufactures. Lifetime performance index is a the-larger-the-better index which is especially useful for non-normal distributions. This paper studied the Bayesian estimation and Bayesian test of life performance index under squared error loss function. The new proposed Bayesian testing method is easier than other classical approaches and the test process is easy to operate by using ordinary programming software as Matlab, Excel, etc. This testing method can be similar used to other life distributions. The testing procedure can provide reference for the enterprise engineers to assess whether the true lifetime performance of products meets the requirements.
Acknowledgement
This study is partially supported by Natural Science Foundation of Hunan Province (No. 2015JJ3030 and No. 2016JJ4012) and Foundation of Hunan Educational Committee (No.15C0228). The author also gratefully acknowledges the helpful comments and suggestions of the reviewers, which have improved the presentation.
References