- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Journal of Quality and Reliability Engineering
Volume 2013 (2013), Article ID 190437, 13 pages
On the Mean Residual Life Function and Stress and Strength Analysis under Different Loss Function for Lindley Distribution
Department of Decision Sciences, Bocconi University, via Roenthen 1, 20136 Milan, Italy
Received 4 November 2012; Accepted 4 February 2013
Academic Editor: Shey-Huei Sheu
Copyright © 2013 Sajid Ali. 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.
Purpose. Mathematical properties of Lindley distribution are derived under different loss functions. These properties include mean residual life function, Lorenz curve, stress and strength characteristic, and their respective posterior risk via simulation scheme. Methodology. Bayesian approach is used for the reliability characteristics. Results are compared on the basis of posterior risk. Findings. Using prior information on the parameter of Lindley distribution, Bayes estimates for reliability characteristics are compared under different loss functions. Practical Implications. Since Lindley distribution is a mixture of gamma and exponential distribution, so Bayesian estimation of reliability characteristics will have a great implication in reliability theory. Originality. A real life application to waiting time data at the bank is also described for the developed procedures. This study is useful for researcher and practitioner in reliability theory.
Exponential distribution is frequently used as a lifetime distribution in statistics and applied areas; the Lindley distribution has been ignored in the literature since 1958. Lindley distribution originally developed by Lindley  and some classical statistic properties are investigated by Ghitany et al. . Sankaran  introduced a discrete version of Lindley distribution known as discrete Poisson-Lindley distribution, and Ghitany and Al-Mutairi  described some estimation methods. The distribution of zero-truncated Poisson-Lindley was introduced by Ghitany et al.  who used the distribution for modeling count data in the case where the distribution has to be adjusted for the count of missing zeros. Zamani and Ismail  introduced negative binomial distribution as an alternative to zero-truncated Poisson-Lindley distribution. Recently, Ghitany et al.  introduced a two-parameter weighted Lindley distribution and pointed that Lindley distribution is particularly useful in modelling biological data from mortality studies.
The rest of the study is organized as follows. Section 2 deals with the derivation of posterior distribution using different noninformative and informative priors. Using different loss functions, the Bayes estimators and their respective posterior risks are discussed in Section 3. Elicitation of hyperparameter is also discussed in Section 3. Simulation study of Bayes estimates of mean residual life and their posterior risks is performed in Section 4. Lorenz curve discussion for Lindley distribution is given in Section 5 while stress and strength reliability characteristics and simulation study under different loss functions is discussed/performed in Section 6. Real life application is illustrated in Section 7. Finally, Section 8 deals with a conclusion and some future remarks.
2. Likelihood Function and Posterior Distributions
The posterior distribution summarizes available probabilistic information on the parameters in the form of prior distribution and the sample information contained in the likelihood function. The likelihood principle suggests that the information on the parameter should depend only on its posterior distribution. Bayesian scientist’s job is to assist the investigator to extract features of interest from the posterior distribution. In this section, we will use the Lindley model as sampling distribution mingles with noninformative priors for the derivation of posterior distribution. A random variable is said to possess a Lindley distribution if it has the following form:
It is obvious from Figure 1 that behavior of Lindley distribution is close to exponential or gamma distribution.
The likelihood function for a random sample which is taken from Lindley distribution is
2.1. Posterior Distribution Using Noninformative Priors
2.1.1. Uniform Prior
An argument in the favor of uniform prior is that when the data are sufficiently informative, so that likelihood function is sharply peaked, then it really does not matter what prior is used, since all reasonably smooth prior densities will lead to approximately the same posterior density. The uniform density in most cases is convenient to simplify calculations of the posterior. This argument supports the uniform prior only in those cases where it produces approximately the same conclusions as the highly imprecise prior constructed from a sufficiently large class of prior densities. If the data are highly informative, the uniform prior may produce reasonable inferences. The uniform prior for is defined as
2.1.2. Jeffreys Prior
Jeffreys was motivated by invariance requirements and suggested a solution to provide a noninformative prior. He used differential geometry method. The requirements are invariance under 1-1 transformations and invariance under sufficient statistics. One-dimensional version of Jeffreys prior has been justified from many different viewpoints. Jeffreys  proposed a formal rule for obtaining a noninformative prior as follows. If is a k-vector valued parameter, then, JP of is where is a kxk Fisher’s (information) matrix whose element is . Fisher’s information matrix is not directly related to the notation of lack of information. The connection comes from the role of Fisher’s matrix in asymptotic theory. Jeffreys noninformative priors based on Fisher’s information matrix often lead to a family of improper priors. The Jeffreys prior of the parameter is
2.2. Posterior Distribution Using Informative Prior
In case of an informative prior, the use of prior information is equivalent to adding a number of observations to a given sample size and, therefore, leads to a reduction of the variance/posterior risk of the Bayes estimates. Bansal  discussed a method to evaluate the relevance of a prior information in terms of the number of additional observation supposed to be added to a given sample size. We used the gamma and conjugate informative prior for analysis.
2.2.1. Gamma Prior
Posterior distribution using informative gamma prior (GP) is
2.2.2. Likelihood Matching Prior (Conjugate Prior)
The likelihood matching prior (LMP) for Lindley distribution is
and the posterior distribution using LMP is
3. Bayes Estimators and Posterior Risk under Different Loss Functions
This section spotlight is on the derivation of the Bayes estimator (BE) under different loss functions and their respective posterior risk (PR). The results are compared for noninformative as well as informative priors. If the decision is a choice of an estimator, then, the Bayes decision is a Bayes estimator. The Bayes estimators are evaluated under squared error loss function (SELF), weighted squared error loss function (WSELF), precautionary loss function (PLF), modified (quadratic) squared error loss function (M/Q SELF), logarithmic loss function (SLLF), entropy loss function (ELF), and K-Loss function. K-loss function proposed by Wasan  is well fitted for a measure of inaccuracy for an estimator of a scale parameter of a distribution defined on = ; this loss function is called K-loss function (KLF). Kanefuji and Iwase  used KLF for the estimation of a scale parameter with a known coefficient of variations. Table 1 (by ) will show the Bayes estimators and their posterior risks for the above-mentioned loss function.
3.1. Elicitation of Hyperparameter(s)
Even though many authors have pointed a need for a formal and comprehensive process for elicitation of hyperparameters, there is no standard method. For elicitation, mainly two points are considered; the functional form of the prior distribution and hyperparameter(s), that is why a natural conjugate prior distribution has been generally recommended because its functional form is identical to likelihood function and posterior distribution can be determined by the way of conjugancy. To determine hyperparameter, we adopted the method discussed by Ali et al. .
4. Mean Residual Life Function
For a continuous distribution with the density and cumulative distribution function , the mean residual life function is defined as
Bayramoglu and Gurler  study the mean residual life function of out of system with nonidentical components while Govil and Aggarwal  and Abdous and Berred  compared for different distributions like gamma, exponential, Pareto, uniform, truncated normal, Maxwell, and lognormal distributions. The mean residual life function for Lindley distribution is
From , one can easily observe that mean residuals life is a diminishing function of time because the distribution belongs to exponential family and for larger parameter value it is close to zero. To evaluate the Bayes estimates and their risk, since the integral appears in both numerator and denumerator, we required a suitable approximate method to obtain Bayes estimates and respective posterior risks. The simplest method is Lindley’s  approximation method, which approaches the ratio of the integrals as a whole and produces a single numerical results. Thus, Lindley approximation (LA) given by Lindley  for obtaining the Bayes estimator and posterior risk of (Mathematica can be used for the solution of integral but takes large time as compared to LA). Many researchers have used this approximation for solving the ratio of integrals for different numbers of parameters for lifetime distributions; see among others Howlader and Hossain , Singh et al. , and Preda et al. .
If is sufficiently large, the ratio of the integral of the form according to Lindley  can be computed as
whereof likelihood;log of prior of, can be evaluated as
Using SELF for , Bayes estimates are overestimated while for and results are underestimated and by increasing sample size these approaches to true parameter values. This behaviour can be observed in all loss functions except PLF. Making comparison between symmetric and asymmetric loss functions, one can easily observes that SELF has smaller posterior risk. Since there is defect in symmetric loss function, that is, SELF assign equal weight to over and under estimation, so we have to look for an alternative choice. In case of asymmetric loss functions; WSELF, SLLF, and PLF can be alternative choices.
5. Lorenz Curve (See Figure 2)
For a positive random variable , the Lorenz curve is defined by the graph of the ratio against with the properties , and for . If represents annual income, then, is the proportion of total income that accrues to individuals having the lowest incomes; see Gail and Gastwirth  for details of Lorenz curves. For the exponential distribution, it is well known that Lorenz curve is given by
For the Lindley distribution, the Lorenz curve is
The comparison of Lorenz curve for exponential and Lindley distribution is given in Figures 2 and 3. Lindley distribution performance for Lorenz curve is better as compared to exponential distribution for different values of .
6. Stress and Strength Analysis for Lindley Distribution
Let and be two random variables such that represents “strength” and “stress”; then, reliability of the stress and strength model is presented as
where is a relationship which represents the probability that the strength exceeds the stress and is joint p.d.f. of and for example, receptor in human eye operates only if it is simulated by a source where magnitude is greater than a random lower threshold for the eye, so here is the probability that the receptor operates. In mechanical reliability of a system gives the probability of a system failure, if the system fails whenever the strength is less than the applied stress.
This quandary has a long narration starting with the revolutionary work of Birnbaum  and Birnbaum and McCarty . The stress strength was first introduced by Church and Harris . A comprehensive handling of the different stress strength models can be found in outstanding monograph by Kotz et al. . Some most fresh work on the stress strength model can be found in Kundu and Gupta [25, 26], Raqab and Kundu , Kundu and Raqab , Krishnamoorthy et al. , Eryilmaz , and the references cited therein. Recently, Al-Mutairi et al.  considered the stress and strength analysis of Lindley distribution using SELF but we are generalizing their work by considering different types of loss functions.
Suppose that and are two independent Lindley rv’s with respective parameters and , then,
Remark. If , then, .
where function of and only; log of likelihood; log of prior of and, can be evaluated as
The Bayes estimates of stress and strength under SELF, MSELF, and ELF are underestimated. By comparing symmetric and asymmetric loss functions, it is noted that posterior risk of SELF is smaller than asymmetric loss functions. In case of asymmetric loss functions, WSELF and PLF have smaller posterior risk than other available loss functions.
Evaluating the performance of informative and noninformative priors, one can easily observe that informative priors have smaller posterior risk due to the availability of compact information. LM and gamma priors both have approximately the same behaviour depending upon the choice of hyperparameters value. More compact information will lead to correct hyperparameters which will lead to definitely better results and smaller posterior risk than noninformative priors. Although there are some depicts where informative priors have posterior risks greater than noninformative priors which is just due to random generation. Increasing sample size in case of SLLF has an inverse effect.
7. Real Life Application
Ghitany et al.  provide waiting times (in minutes) before service of 100 bank customers data set for Lindley distribution. They fitted both Lindley and exponential distributions (both have the same number of parameters) by method of maximized likelihood method and found Lindley distribution provides better fit. The data is given in Tables 16 and 17.
We fit on both data sets Kolmogorov-Smirnov test and found that Lindley distribution is good fitted. The values of K-S test along value are given in Table 18.
Since the Lindley distribution belongs to the exponential family so the natural conjugate prior is Gamma distribution.
The posterior risks of LMP and GP are approximately the same as compared to noninformative priors. There are some posterior risk values which are greater than noninformative priors. These are just due to hyperparameters value effect that is, more accurate values will lead to the smaller posterior risk. PLF and WSELF loss functions have smaller posterior risk as compared to other loss functions.
8. Conclusion and Suggestions
We consider the Bayesian analysis of the Lindley model via informative and informative priors under different loss functions. Based on posterior distribution, different properties, we conclude that informative priors (LMP, GP) performance approximately equal and have smaller posterior risk’s as compared to the noninformative priors; also Jeffreys prior results are more precised than uniform prior. In other words, we can summarize result as
The choice of loss function as concerned, one can easily observe based on evidence (different properties as discussed above) that PLF, SLLF, and WSELF are suitable than other asymmetrical loss functions. One thing is common as we increase sample size posterior risk comes down. In future, this work can be extended using censored data.
- D. V. Lindley, “Fiducial distributions and Bayes' theorem,” Journal of the Royal Statistical Society. Series B, vol. 20, pp. 102–107, 1958.
- M. E. Ghitany, B. Atieh, and S. Nadarajah, “Lindley distribution and its application,” Mathematics and Computers in Simulation, vol. 78, no. 4, pp. 493–506, 2008.
- M. Sankaran, “The discrete Poisson-Lindley distribution,” Biometrics, vol. 26, pp. 145–149, 1970.
- M. E. Ghitany and D. K. Al-Mutairi, “Estimation methods for the discrete Poisson-Lindley distribution,” Journal of Statistical Computation and Simulation, vol. 79, no. 1, pp. 1–9, 2009.
- M. E. Ghitany, D. K. Al-Mutairi, and S. Nadarajah, “Zero-truncated Poisson-Lindley distribution and its application,” Mathematics and Computers in Simulation, vol. 79, no. 3, pp. 279–287, 2008.
- H. Zamani and N. Ismail, “Negative binomial-Lindley distribution and its application,” Journal of Mathematics and Statistics, vol. 6, no. 1, pp. 4–9, 2010.
- M. E. Ghitany, F. Alqallaf, D. K. Al-Mutairi, and H. A. Husain, “A two-parameter weighted Lindley distribution and its applications to survival data,” Mathematics and Computers in Simulation, vol. 81, no. 6, pp. 1190–1201, 2011.
- H. Jeffreys, Theory of Probability, Oxford University Press, 3rd edition, 1964.
- A. K. Bansal, Bayesian Parametric Inference, Narosa Publishing House, New Delhi, India, 2007.
- M. Wasan, Parametric Estimation, McGraw-Hill, New York, NY, USA, 1970.
- K. Kanefuji and K. Iwase, “Estimation for a scale parameter with known coefficient of variation,” Statistical Papers, vol. 39, no. 4, pp. 377–388, 1998.
- S. Ali, M. Aslam, and S. M. A. Kazmi, “A study of the effect of the loss function on Bayes Estimate, posterior risk and hazard function for Lindley distribution,” Applied Mathematical Modelling, vol. 37, no. 8, pp. 6078–6078, 2013.
- I. Bayramoglu and S. Gurler, “On the mean residual life function of the k-out-of-n system with nonidentical components,” in International Conference on Mathematical and Statistical Modeling in Honor of Enrique Castillo, pp. 28–30, Ciudad Real, Spain, June 2006.
- K. K. Govil and K. K. Aggarwal, “Mean residual life function for normal, gamma and lognormal densities,” Reliability Engineering, vol. 5, no. 1, pp. 47–51, 1983.
- B. Abdous and A. Berred, “Mean residual life estimation,” Journal of Statistical Planning and Inference, vol. 132, no. 1-2, pp. 3–19, 2005.
- D. V. Lindley, “Approximate Bayesian methods,” Trabajos de Estadistica Y de Investigacion Operativa, vol. 31, no. 1, pp. 223–245, 1980.
- H. A. Howlader and A. M. Hossain, “Bayesian survival estimation of Pareto distribution of the second kind based on failure-censored data,” Computational Statistics and Data Analysis, vol. 38, no. 3, pp. 301–314, 2002.
- R. Singh, S. K. Singh, U. Singh, and G. P. Singh, “Bayes estimator of generalized-exponential parameters under Linex loss function using Lindley's approximation,” Data Science Journal, vol. 7, pp. 65–75, 2008.
- V. Preda, E. Panaitescu, and A. Constantinescu, “Bayes estimators of modified-Weibull distribution parameters using Lindley's approximation,” WSEAS Transactions on Mathematics, vol. 9, no. 7, pp. 539–549, 2010.
- M. H. Gail and J. L. Gastwirth, “A scale-free goodness-of-fit test for the exponential distribution based on the Lorenz curve,” Journal of the American Statistical Association, vol. 73, no. 364, pp. 787–793, 1978.
- Z. W. Birnbaum, “On a use of the Mann-Whitney statistic,” in Proceedings of the 3rd Berkeley Symposium on Mathematical Statistics and Probability, vol. 1, pp. 13–17, University of California Press, Berkeley, Calif, USA, 1956.
- Z. W. Birnbaum and R. C. McCarty, “A distribution-free upper confidence bound for , based on independent samples of and ,” Annals of Mathematical Statistics, vol. 29, pp. 558–562, 1958.
- J. D. Church and B. Harris, “The estimation of reliability from stress strength relationships,” Technometrics, vol. 12, pp. 49–54, 1970.
- S. Kotz, Y. Lumelskii, and M. Pensky, The Stress-Strength Model and Its Generalizations: Theory and Applications, World Scientific, River Edge, NJ, USA, 2003.
- D. Kundu and R. D. Gupta, “Estimation of for Weibull distributions,” IEEE Transactions on Reliability, vol. 55, no. 2, pp. 270–280, 2006.
- D. Kundu and R. D. Gupta, “Estimation of for generalized exponential distribution,” Metrika, vol. 61, no. 3, pp. 291–308, 2005.
- M. Z. Raqab and D. Kundu, “Comparison of different estimators of for a scaled Burr type distribution,” Communications in Statistics: Simulation and Computation, vol. 34, no. 2, pp. 465–483, 2005.
- D. Kundu and M. Z. Raqab, “Estimation of for three-parameter Weibull distribution,” Statistics and Probability Letters, vol. 79, no. 17, pp. 1839–1846, 2009.
- K. Krishnamoorthy, S. Mukherjee, and H. Guo, “Inference on reliability in two-parameter exponential stress-strength model,” Metrika, vol. 65, no. 3, pp. 261–273, 2007.
- S. Eryilmaz, “On system reliability in stress-strength setup,” Statistics and Probability Letters, vol. 80, no. 9-10, pp. 834–839, 2010.
- D. K. Al-Mutairi, M. E. Ghitany, and D. Kundu, “Inferences on stress-strength reliability from Lindley distribution,” to appear in Communications in Statistics—Theory and Methods, available from http://home.iitk.ac.in/~kundu/lindley-ss.pdf.