arxiv: v1 [stat.ap] 12 Jun 2014

Transcription

1 On Lindley-Exponential Distribution: Properties and Application Deepesh Bhati 1 and Mohd. Aamir Malik 2 arxiv: v1 [stat.ap] 12 Jun Assistant Professor, Department of Statistics, Central University of Rajasthan, deepesh.bhati@curaj.ac.in 2 Department of Statistics, Central University of Rajasthan aamirmalik.stats@gmail.com November 25, 216 Abstract In this paper, we introduce a new distribution generated by Lindley random variable which offers a more flexible model for modelling lifetime data. Various statistical properties like distribution function, survival function, moments, entropy, and iting distribution of extreme order statistics are established. Inference for a random sample from the proposed distribution is investigated and maximum likelihood estimation method is used for estimating parameters of this distribution. The applicability of the proposed distribution is shown through real data sets. Keyword: Lindley Distribution, Entropy, Stress-Strength Reliability Model, Maximum Likelihood Estimator. AMS 21 Subject Classification: 6E5 1 Introduction Lifetime distribution represents an attempt to describe, mathematically, the length of the life of a system or a device. Lifetime distributions are most frequently used in the fields like medicine, engineering etc. Many parametric models such as exponential, gamma, Weibull have been frequently used in statistical literature to analyze lifetime data. But there is no clear motivation for the gamma and Weibull distributions. They only have more general mathematical closed form than the exponential distribution with one additional parameter. Recently, one parameter Lindley distribution has attracted the researchers for its use in modelling lifetime data, and it has been observed in several papers that this distribution has performed excellently. The Lindley distribution was originally proposed by Lindley [11] in the context of Bayesian statistics, as a counter example of fudicial statistics which can be seen that as a mixture of expθ) and gamma2, θ). More details on the Lindley distribution can be found in Ghitany et al. [7]. A random variable X is said to have Lindley distribution with parameter θ if its probability density function is defined as: f X x; θ) = θ2 θ + 1) 1 + x)e θx ; x >, θ > 1) 1

2 with cumulative distribution function F x) = 1 e θx 1 + θ + θx) 1 + θ Some of the advances in the literature of Lindley distribution are given by Ghitany et al. [5] who has introduced a two-parameter weighted Lindley distribution and has pointed that Lindley distribution is particularly useful in modelling biological data from mortality studies. Mahmoudi et. al. [12] have proposed generalized Poisson Lindley distribution. Bakouch et al. [3] have come up with extended Lindley EL) distribution, Adamidis and Loukas [1] have introduced exponential geometric EG) distribution. Shanker et. al. [14] have introduced a two-parameter Lindley distribution. Zakerzadeh et al.[15] have proposed a new two parameter lifetime distribution: model and properties. M.K. Hassan [7] has introduced convolution of Lindley distribution. Ghitany et al.[6] worked on the estimation of the reliability of a stressstrength system from power Lindley distribution. Elbatal et al.[4] has proposed a new generalized Lindley distribution. Ristić [13] has introduced a new family of distributions with survival function given by F X x) = 1 Γα) loggx)) t α 1 e tdt, x R, θ > In this paper we introduce a new family of distribution generated by a random variable T which follows one parameter Lindley distribution. The survival function of this new family is given as: F X x) = θ2 1 + θ loggx)) 1 + t)e θt dt, x R, θ > 2) where θ > and Gx) is a cumulative distribution functioncdf) which we use to generate a new distribution. The cdf Gx) is referred to as a transformer and the corresponding probability density function pdf) is given by fx) = θ2 1 + θ 1 log Gx)) Gx)θ 1 gx), 3) We consider the transformer to follow exponential distribution with cdf Gx) = 1 e λx ). Hence the survival function of the new distribution is given by F x) = 1 with corresponding density given by 1 e λx ) θ 1 + θ θ log 1 e λx )) 1 + θ fx) = θ2 λe λx 1 e λx) θ 1 1 log 1 e λx )) 1 + θ) We refer the random variable with survival function 4) as Lindley-ExponentialL-E) distribution with parameters θ and λ which we denote by L-Eθ, λ). The aim of this paper is to study the mathematical properties of the L-E distribution and to illustrate its applicability. The contents are organized as follows. The analytical shapes of the pdfin equations 5) are established in section 2. The quantile function presented in section 3. The expressions for the moment generating function and moments corresponding to Equation 5) are given in Section 4. Limiting distribution of sample statistics like maximum and minimum has 4) 5) 2

3 been shown in section 5. In section 6, entropy of L-E distribution is presented. The maximum likelihood estimation procedure is considered in Section 7. The performance of the maximum likelihood estimators for small samples is assessed by simulation in Section 8. Section 9 gives estimation of stress-strength parameter R by using maximum likelihood estimation method. Finally we conclude the paper by showing applicability of the model to the real data sets. 2 Shape of the density Here, the shape of pdf 5) follows from theorem 1. Theorem 1: The probability density function of the L-E distribution is decreasing for < θ < 1 and unimodel for θ > 1. In the latter case, mode is a root of the following equation: Proof: The first order derivative of logfx)) is 1 θe λx ) 1 log1 e λx ) ) + e λx = d logfx)) dx = λ rx) 1 e λx ) 1 log 1 e λx )) 6) where, rx) = θe λx + 1 θe λx) log 1 e λx) e λx 1 ). For < θ < 1, the function rx) is negative. So f x) < for all x >. This implies that f is decreasing for < θ < 1. Also note that, log f) ) = and log f) ) <. This implies that for θ > 1, gx) has a unique mode at x such that rx) > for x < x and rx) < for x > x. So, g is unimodal function with mode at x = x. The pdf for various values of λ and θ are shown in Figure 1. We, now, consider the hazard rate function hrf) of the L-E distribution, which is given by hx) = θ2 λe λx 1 e λx) θ 1 1 log 1 e λx )) 1 + θ 1 e λx ) θ 1 + θ θ log 1 e λx )) 7) Proposition 1: For θ > the hazard rate function follows relation hx) = λ. x Proof: The proof is straight forward and is omitted. In Figure 2, hazard function for different values of parameters θ and λ. 3 The Quantile Function of L-E distribution The cdf, F X x) = 1 F x), can be obtained by using eq.4). Further, it can be noted that F X is continuous and strictly increasing so the quantile function of X is Q X γ) = F 1 X γ), < γ < 1. In the following theorem, we give an explicit expression for Q X in terms of the Lambert W function. For more details on Lambert W function we refer the reader to Jodrá [9]. Theorem 2: For any θ, λ >, the quantile function of the L-E distribution X is x γ = F 1 γ) = 1 λ log 1 W 1 e θ 1 θ + 1)γ ) ) 1/θ 8) θ + 1)γ 3

4 g HxL x g HxL x g HxL x g HxL x g HxL.3 g HxL x x Figure 1: PDF plot for various values of λ and θ. 4

5 h x) h x) x x Figure 2: Hazard function for various values of λ and θ. where W 1 denotes the negative branch of the Lambert W function. Proof: By assuming p = 1 e λx, the cdf can be written as F X x) = pθ 1 + θ θ logp)) 1 + θ for fixed θ, λ > and γ, 1), the γ th quantile function is obtained by solving F X x) = γ. By re-arranging the above, we obtain logp θ ) γp θ 1 + θ) = 1 + θ) taking exponential and multiplying 1 + θ)γ on both sides, we get p θ 1 + θ)γe p θ 1+θ)γ = 1 + θ)γe 1+θ) 9) By using definition of Lambert-W function W z)e W z) = z, where z is a complex number), we see that p θ 1 + θ)γ is the Lambert W function of the real argument 1 + θ)γe 1+θ). Thus, we have W 1 + θ)γe 1+θ)) = p θ 1 + θ)γ 1) Moreover, for any θ, λ > it is immediate that p θ 1 + θ)γ > 1, and it can also be checked that 1 + θ)γe 1+θ) 1/e, ) since γ, 1). Therefore, by taking into account the properties of the negative branch of the Lambert W function, we have W θ)γe 1+θ)) = p θ 1 + θ)γ Also by substituting p = 1 e λx in cdf and solving it for x, we get x γ = F 1 γ) = 1 λ log 1 W 1 e θ 1 θ + 1)γ ) ) 1/θ 11) θ + 1)γ Further the first three quantiles we obtained by substituting γ = 1 4, 1 2, 3 4 in equation 11). 5

6 ) Q 1 = F 1 1/4) = 1 λ log 4W e θ 1 1+θ) θ) ) MedianM d ) = Q 2 = F 1 1/2) = 1 λ log 2W e θ 1 1+θ) θ) ) Q 3 = F 1 3/4) = 1 λ log 4W 3e θ 1 1+θ) θ) 1 θ 1 θ 1 θ 12) 4 Moments The moment generating function of the random variable X follow L-E distribution is given as M X t) = Ee tx θ 2 λe xλ t) 1 e λx) θ 1 )) 1 log 1 e λx ) = dx θ + 1 ) ψ ) θ) + ψ ) t λ + θ + 1) + 1 ) = θγθ + 1)Γ 1 t λ θ + 1)Γ t 13) λ + θ + 1) where, ψ n) z) = dn ψz) dz n and ψz) = Γ z) Γz) known as digamma function. ) Hence the first and second raw moments can be obtained by respectively. dmx t) dx t= and d 2 M X t) dx 2 )t= EX) = γθ + 1 θ + θ ψ1) θ + 1)) + θ + 1)ψ ) θ) + γ + 1 θ + 1)λ 14) EX 2 ) = 6γ 2 θ+π 2 θ+6θψ ) θ+1) 3 +62γθ+θ)ψ ) θ+1) 2 6θψ ) θ)ψ ) θ+1) 2 12γθψ ) θ) γθ 1)ψ ) θ) 62γθ+θ+3)ψ 1) θ+1) 12θψ ) θ)ψ 1) θ+1)+6θψ 2) θ+1)+6γ 2 +π 2 +12γ 6θ + 1)λ 2 15) where γ is Eulergamma constant = Table 1 displays the mode, mean and median for L-E distribution for different choices of parameter λ and θ. It can be observed from the table that all the three measures of central tendency decrease with increase in λ and increase with an increase in θ. Also for any choice of λ and θ it is observed that Mean > Median > Mode, which is an indication of positive skewness. 5 Limiting Distribution of Sample Minima and Maxima We can derive the asymptotic distribution of the sample minimum X 1:n by using theorem of Arnold t. al.[2], it follows that the asymptotic distribution of X 1:n is Weibull type with shape parameter θ > if F tx) t + F t) = xθ 6

7 Table 1: Mode, Mean and Median for various value of parameter θ λ Mode Mean Median Mode Mean Median Mode Mean Median Mode Mean Median Mode Mean Median Mode Mean Median Mode Mean Median for all x >. Then, by using L Hópitals rule, it follows that Since and F tx) t + F t) = x ftx) t + ft) = x 1 e λtx ) θ 1 1 log1 e λtx )) t + 1 e λt ) θ 1 1 log1 e λt = x θ )) 1 e λtx ) θ 1 t + 1 e λt ) θ 1 = xθ 1 1 log1 e λtx )) t + 1 log1 e λt )) = 1 Hence, we obtain that the asymptotic distribution of the sample minima X 1:n is of the Weibull type with shape parameter θ. Further, it can be seen that 1 F t + x) 1 + θ) 1 e λx+t) ) θ 1 + θ θ log1 e λx+t) )) = t 1 F t) t 1 + θ) 1 e λt ) θ 1 + θ θ log1 e λt )) by using L-Hópitals rule, t e λx 1 e λt+x)) θ 1 1 log 1 e λt+x) )) 1 e λt ) θ 1 1 log = e λθx 1 e λt))) 7

8 Since and t 1 e λt+x)) ) θ 1 1 e λt ) θ 1 = e λxθ 1) 1 log 1 e λ t+x))) t 1 log 1 e λ t)) = 1 Hence, it follows from Theorem in Leadbetter et al. 1983) that there must be norming constants a n,b n, c n > and d n such that and Pr{a n M n b n ) x} e e θλx 16) Pr{c n m n d n ) x} 1 e xθ 17) as n. By following Corollary in Leadbetter et al. 1983), we can determine the form of the norming constants. As an illustration, one can see that a n = θ and b n = F 1 1 1/n), where F 1.) denotes the inverse function of F.). 6 Entropy In many field of science such as communication, physics and probability, entropy is an important concept to measure the amount of uncertainty associated with a random variable X. Several entropy measures and information indices are available but among them the most popular entropy measure called Rényi entropy is defined as Jζ) = 1 1 ζ log f ζ x)dx, for ζ > 1 and ζ 1 18) In our case f ζ dx = R + R + R + θ 2 ) ζ λ e xλζ 1 e xλ) ζθ 1) 1 log 1 e xλ )) ζ dx θ + 1 substituting x = 1 λ log1 e u ) and using power series expansion 1 z) α = 1) j α j) z j, the above expression reduces to f ζ θ 2ζ λ ζ 1 dx = θ + 1) ζ u + 1) ζ 1 e u) ζ 1 e u ζθ 1)+1)) du R + = θ2ζ λ ζ 1 ) ζ 1 θ + 1) ζ 1) j e u ζθ 1)+j+1)) u + 1) ζ du j j= = θ2ζ λ ζ 1 ) ζ 1 θ + 1) ζ 1) j u + 1) ζ e u ζθ 1)+j+1)) du j = θ2ζ λ ζ 1 θ + 1) ζ j= j= 1) j e ζθ 1)+j+1 E ζ [j + θ 1)ζ + 1)] 19) j= 8

9 where E n z) = e zt t n dt known as exponential integral function. For more details 1 see Thus according to 18) the Rényi entropy of L-Eθ, λ) distribution is given by Jζ) = 1 1 ζ log θ2ζ λ ζ 1 θ + 1) ζ 1) j e ζθ 1)+j+1 E ζ [j + θ 1)ζ + 1)] 2) j= Moreover, the Shannon entropy is defined by E[ logfx))]. This is a special case derived from ζ 1 Jζ) 7 Maximum Likelihood function In this section we shall discuss the point and interval estimation on the parameters that index the L-Eθ, λ). Let the log-likelihood function L = lθ) of single observationsay x i ) for the vector of parameter Θ = θ, λ) can be written as l n x, Θ) = 2 logθ)+logλ) λx i logθ+1)+θ 1) log1 e λxi )+log1 log1 e λxi )), x i > The associated score function is given by U n Θ) = l n θ, ), ln λ where l n θ = 2 θ 1 θ log 1 e λxi) 21) l n λ = 1 λ x i + θ 1) x ie λxi 1 e x i e λx λxi 1 e λxi ) 1 log 1 e λxi )) 22) As we know th xpctd value of score function equals zero, i.e. EUΘ) = ), which implies E log 1 e λx)) = 1 θ+1 2 θ The total log-likelihood of the random sample x = x 1,, x n ) of size n from X is given by l n = n l i) and th total score function is given by U n = n U i), where l i) is the log-likelihood of 1 i th observation. The maximum likelihood estimator ˆΘ of Θ is obtained by solving equation21) and 22) numerically or this can also be obtained easily by using nlm) function in R. Moreover the Fisher information matrix is given by [ ] κθ,θ κ K = K n Θ) = n θ,λ 23) κ λ,θ κ λ,λ i=1 where 2 κ θ,θ = θ + 1) 2 2 θ 2 ) Xe λx κ λ,θ = E 1 e λx κ λ,λ = 1 λ 2 + θ 1)E X 2 e λx 1 e λx ) 2 ) e λx X 2 log 1 e λx) + e λx) ) E 1 e λx ) 2 log 1 e λx )) 2 24) The above expressions depend on some expectations which easily computed using numerical integration. Under the usual regularity conditions, the asymptotic distribution of n ˆΘ Θ ) is N 2, KΘ) 1 ) 25) 9

10 where n = K nθ) 1 = KΘ) 1. The asymptotic multivariate normal N 2, KΘ) 1 ) distribution of ˆΘ can be usd to construct approximate confidence intervals. An asymptotic confidence interval with significance level α for each parameter θ and λ is ACI θ, 11 α)%) =ˆθ z α/2 κθ,θ), ˆθ + z α/2 κθ, θ)) ACI λ, 11 α)%) =ˆλ z α/2 κλ,λ), ˆλ + z α/2 κλ, λ)) 26) where z 1 α/2 denotes 1 α/2 quantile of standard normal random variable. 8 Simulation In this section, we investigate the behavior of the ML estimators for a finite sample size n). Simulation study based on different L-Eθ, λ) distribution is carried out. The random variable are generated by using cdf technique presented in section 4 from L-Eθ, λ) are generated. A simulation study consisting of following steps is being carried out for each triplet θ, λ, n), where θ =.5, 1, 2, λ =.5, 1, 2, 3 and n = 2, 5, 75, Choose the initial values of θ, λ for the corresponding elements of the parameter vector Θ = θ, λ), to specify L-Eθ, λ) distribution; 2. choose sample size n; 3. generate N independent samples of size n from L-Eθ, λ); 4. compute the ML estimate ˆ Θ n of Θ for each of the N samples; 5. compute the mean of the obtained estimators over all N samples, the average bias = N N 1 N Θ i Θ ) and the average mean square error MSE Θ) = 1 N Θ i Θ ) 2, of i=1 simulated estimates. i=1 θ =.5 θ = 1 θ = 2 Table 2: Average bias of the simulated estimates. n λ =.5 λ = 1 λ = 2 λ = 3 bais θ) bais λ) bais θ) bais λ) bais θ) bais λ) bais θ) bais λ)

11 θ =.5 θ = 1 θ = 2 Table 3: Average MSE of the simulated estimates. n λ =.5 λ = 1 λ = 2 λ = 3 MSEθ) MSEλ) MSEθ) MSEλ) MSEθ) MSEλ) MSEθ) MSEλ) Application to Real Datasets In this section, we illustrate, the applicability of L-E Distribution by considering two different datasets used by different researchers. We also fit L-E distribution, Power-Lindley distribution [6], New Generalized Lindley Distribution [4], Lindley Distribution, Weibull distribution and Exponential distribution. Namely i) Power-Lindley distribution PLα, β)): f 1 x) = αβ2 1 + β 1 + xα )x α 1 e βx α, x, α, β >. ii) New Generalized Lindley distribution NGLDα, β, θ)): iii) Lindley Distribution Lθ)) f 2 x) = e θx θ α+1 x α 1 + θβ x β 1 ), x, α, θ, β > 1 + θ Γα) Γβ) f 3 x) = θ2 1 + θ 1 + x)e θx, x, α, β >. In each of these distributions, the parameters are estimated by using the maximum likelihood method, and for comparison we use negative log-likelihood values LL), the Akaike information criterion AIC) and Bayesian information criterion BIC) which are defined by 2LL + 2q and 2LL+q logn), respectively, where q is the number of parameters estimated and n is the sample size. Further K-SKolmogorov-Smirnov) test statistic defined as K S = sup x F n x) F x), n where F n x) = 1 n I xi x is empirical distribution function and F x) is cumulative distribution i=1 function is calculated and shown for all the datasets. 9.1 Illustration 1 We consider an uncensored data set corresponding to remission times in months) of a random sample of 128 bladder cancer patientslee and Wang[1]) as presented in Appendix A.1.The data sets are presented in appendix A.1 in Table 6). The results for these data are presented in Table 4. We observe that the L-E distribution is a competitive distribution as compared with 11

12 LE NGLD PL L Fn x x Figure 3: PDF plot for various values of λ and θ. other distributions. In fact, based on the values of the AIC, BIC and as well as the value of the K-S test statistic, we observe that the L-E distribution provides the best fit for these data among all the models considered. In Figure 2, we have plotted probability density function and empirical distribution function for all considered distributions for these data. Table 4: The estimates of parameters and goodness-of-fit statistics for Illustration 1. Model Parameter -LL AIC BIC K-S statistic L-E ˆλ=.962, ˆθ= PL ˆθ=.385, ˆβ= L ˆθ= NGLD ˆθ=.18, ˆα=4.679, ˆβ= Illustration 2 As second example, we consider 1 observations on waiting time in minutes) before the customer service in a bank see Ghitany et al.[7]). The data sets are presented in appendix A.2 in Table 7). The results for these data are presented in Table 5. From these results we can observe that L-E distribution provide smallest AIC and BIC values as compare to Power lindley, new generalized Lindley distribution, Lindley and exponential and hence best fits the data among all the models considered. The results are presented in Table 5 and probability density function and empirical distribution function are shown in Figure 3. Table 5: The estimates of parameters and goodness-of-fit statistics for Illustration 2. Model Parameter -LL AIC BIC K-S L-E ˆθ=2.65; ˆλ= PL ˆθ=.153; ˆβ= L ˆθ= E ˆθ= NGLD ˆθ=.233; ˆβ=2.8; ˆα=

13 LE NGLD PL L E Fn x x Figure 4: PDF plot for various values of λ and θ. 1 Estimation of the Stress-Strength Parameter The stress-strength parameter R) plays an important role in the reliability analysis as it measures the system performance. Moreover, R provides the probability of a system failure, the system fails whenever the applied stress is greater than its strength i.e. R = P X > Y ). Here X L-Eθ 1, λ) denotes the strength of a system subject to stress Y, and Y L-Eθ 2, λ), X and Y are independent of each other. In our case, the stress-strength parameter R is given by R = P X > Y ) = = S X y)f Y y)dy P X > Y Y = y) f Y y)dy θ = 1 2λ 2 1 e λy) θ 1+θ 2 1 e λy 1 log 1 e λy)) 1 + θ 1 1 log 1 e λy ) )) dy 1 + θ 1 ) 1 + θ 2 ) ) =1 θ2 2 θ θ 2 + 3) θ1 2 + θ 2 + 1) θ 2 + 3) θ 1 + θ2 2 + θ 2 θ 1 + 1) θ 2 + 1) θ 1 + θ 2 ) 3 27) Remarks: i) R is independent of λ ii) When θ 1 = θ 2, R=.5. This is intuitive that X and Y are i.i.d. and there is an equal chance that X is bigger than Y. Since R in equation 27) is a function of stress-strength parameters θ 1 and θ 2 we need to obtain the maximum likelihood estimators MLEs) of θ 1 and θ 2 to compute the MLE of R under invariance property of the MLE. Suppose that X 1, X 2,, X n and Y 1, Y 2,, Y m are independent random samples from L-E θ 1, λ) and L-E θ 2, λ) respectively. Thus, the likelihood function based on the observed sample is given by 13

14 Lθ 1, θ 2 ) = θ2m 2 θ1 2n λ m+n n e λxi 1 e λxi) θ log 1 e λx i )) θ 1 + 1) θ 2 + 1) i=1 m e λyj 1 e ) λyj θ log 1 e λy j )) The log - likelihood function is given by j=1 log Lθ 1, θ 2 ) =2n log θ 1 ) + 2m log θ 2 ) + m + n) logλ) n log θ 1 + 1) m log θ 2 + 1) n n m θ 1 1) s 1 λ x i + log1 log1 e λxi )) θ 2 1) s 2 λ + i=1 i=1 m log1 log1 e λyj )) j=1 where s 1 = n log ) 1 e λxi and s2 = m log ) 1 e λyj. i=1 j=1 The MLE of θ 1 and theta 2, say ˆθ 1 and ˆθ 2 respectively, can be obtained as the solutions of the following equations logl = 2n n θ 1 θ 1 θ s 1 logl = 2m m 28) θ 2 θ 2 θ s 2 from above equations ˆθ 1 = s 1 n) + s 1 n) 2 + 8ns 1 2s 1 ˆθ 2 = s 2 m) + s 2 m) 2 + 8ms 2 2s 2 Hence, using the invariance property of the MLE, the maximum likelihood estimator ˆR mle of R can be obtained by substituting ˆθ k for k=1,2 in equation 27). ) ˆR mle = 1 θ2 2 θ θ 2 + 3) θ1 2 + θ 2 + 1) θ 2 + 3) θ 1 + θ2 2 + θ 2 θ 1 + 1) θ 2 + 1) θ 1 + θ 2 ) 3 3) θ1=ˆθ 1,θ 2=ˆθ Asymptotic Confidence For an estimator ˆθ k to be asymptotically efficient for estimating θ k for large samples, we should have n k ˆθk θ k ) D N, Iθk ) 1) 2 ) logl Iθ k ) = E θk 2 = E n I 1 = Iθ 1 ) = E θ 1 + 1) 2 2n ) θ1 2 = 2n θ1 2 m I 2 = Iθ 2 ) = E θ 2 + 1) 2 2m ) θ2 2 = 2m θ2 2 n k θ k + 1) 2 2n k θ 2 k ) n θ 1 + 1) 2 m θ 2 + 1) 2 = 2n k θ 2 1 n k θ 1 + 1) 2 j=1 y j 29) 14

15 where n 1 = n and n 2 = m therefore, as n and m ˆR R d 2 2 mi 2 + d2 1 ni 1 D N, 1) where d 1 = R = θ 1θ2 2 θ θ 2 + 3) θ1 2 + θ θ ) θ θ θ )) θ 1 θ 1 + 1) 2 θ 2 + 1) θ 1 + θ 2 ) 4 d 2 = R = θ2 1θ 2 θ θ θ 2 + θ1 2 θ 2 + 2) + 2θ 1 θ θ ) + 6 ) θ 2 θ 1 + 1) θ 2 + 1) 2 θ 1 + θ 2 ) 4 Interval estimators and 11 α)% confidence interval for R can be obtained by using the asymptotic distribution of ˆR, and obtained as ˆd2 2 ˆd ˆR 2 1 ± Z α + 2 mî2 nî1 Conclusion: We have proposed the new distribution Lindley-Exponential L-E) distribution generated by Lindley distribution. We have derived important properties of the L-E distribution like moments, entropy, asymptotic distribution of sample maximum and sample Minimum. We have illustrated the application of L-E distribution to two real data sets used by researchers earlier. By comparing L-E distribution with other popular models we conclude that L-E distribution performs satisfactorily or better. References [1] Adamidis K., and Loukas S.,1998) A lifetime distribution with decreasing failure rate, Statistics and Probability Letters, Vol39), [2] Arnold B.C., Balakrishnan N. and Nagaraja H.N.213): A First Course in Order Statistics, Wiley, New York, [3] Bakouch H. S., Al-Zahrani B. M., Al-Shomrani A. A., Marchi V. A., and Louzada F.,212) An extended Lindley distribution, Journal of the Korean Statistical Society, Vol41), [4] Elbatal I., Merovci F., and Elgarhy M.213): A new generalized Lindley distribution, Mathematical Theory and Modeling, Vol3) no. 13. [5] Ghitany M. E., Alqallaf F., Al-Mutairi D. K., and Husain H. A., 211) A two-parameter weighted Lindley distribution and its applications to survival data, Mathematics and Computers in Simulation, Vol. 81), no. 6, [6] Ghitany M. E., Al-Mutairi D. K., and Aboukhamseen S. M., 213) Estimation of the reliability of a stress-strength system from power Lindley distributions, Communications in Statistics - Simulation and Computation, Vol 78), [7] Ghitany M. E., Atieh B., and Nadarajah, S., 28) Lindley distribution and its application, Mathematics and Computers in Simulation, Vol 78),

16 [8] Hassan M.K.,214), On the Convolution of Lindley Distribution, Columbia International Publishing Contemporary Mathematics and Statistics, Vol. 2) No. 1, [9] Joŕda P.,21) Computer generation of random variables with Lindley or Poisson-Lindley distribution via the Lambert W function, Mathematics and Computers in Simulation, Vol81), [1] Lee E.T., and Wang J.W.,23) Statistical methods for survival data analysis, John Wiley & Sons, inc., Hoboken, New Jersey, 3 edition. [11] Lindley D. V.,1958) Fiducial distributions and Bayes theorem, Journal of the Royal Statistical Society, Series B Methodological), [12] Mahmoudi E., and Zakerzadeh H., 21) Generalized Poisson Lindley distribution, Communications in Statistics: Theory and Methods, Vol 39), [13] Miroslav M. Risti and Narayanaswamy Balakrishnan 212): The gammaexponentiated exponential distribution, Journal of Statistical Computation and Simulation, 82:8, , DOI: 1.18/ [14] Shanker R., Sharma S., and Shanker R., 213) A Two-Parameter Lindley Distribution for Modeling Waiting and Survival Times Data, Applied Mathematics, Vol 4), [15] Zakerzadeh, H. and Mahmoudi, E.,212) A new two parameter lifetime distribution: model and properties. arxiv: v1 [stat.co]. 16

17 11 Appendix 11.1 A.1- Dataset used in Illustration 1: Table 6: The remission times in months) of bladder cancer patients A.2- Dataset used in Illustration 2: Table 7: Waiting times min.) of 1 bank customers