Submit manuscript...
eISSN: 2378-315X

Biometrics & Biostatistics International Journal

Research Article Volume 7 Issue 2

A quasi exponential distribution

Rama Shanker,1 Kamlesh Kumar Shukla,1 Shambhu Sharma,2 Ravi Shanker3

1lDepartment of Statistics, Eritrea Institute of Technology, Eritrea
2Department of Mathematics, Dayalbagh Educational Institute, India
3Department of Mathematics, N.P University, India

Correspondence: Rama Shanker, Department of Statistics, College of Science, Eritrea Institute of Technology, Asmara, Eritrea

Received: February 23, 2018 | Published: March 9, 2018

Citation: Shanker R, Shukla KK, Sharma S. A quasi exponential distribution. Biom Biostat Int J. 2018;7(2):90–94. DOI: 10.15406/bbij.2018.07.00194

Download PDF

Abstract

In this paper, a quasi exponential distribution has been proposed and studied. Its statistical properties including moments, survival function, hazard rate function, mean residual life function and stochastic ordering have been discussed. Maximum likelihood estimation has been discussed for estimating the parameter of the distribution. Application of the distribution has been explained with a real lifetime data from biomedical science and its goodness of fit shows quite satisfactory over exponential distribution.

Keywords: exponential distribution, moments, hazard rate function, mean residual life function, stochastic ordering, parameter estimation, application.

Abbreviations

QED, quasi exponential distribution; PDF, probability density function; CDF, cumulative distribution function; MLE, maximum likelihood estimate; AIC, Akaike information criterion

Introduction

The exponential distribution having scale parameter  (also known as rate parameter for failure, death, arrival etc) is defined by its probability density function (pdf) and cumulative distribution function (cdf)

f1(x;θ)=θe-θx;x>0,θ>0f1(x;θ)=θeθx;x>0,θ>0           (1.1)

F1(x;θ)=1-e-θx;x>0,θ>0F1(x;θ)=1eθx;x>0,θ>0           (1.2)

The exponential distribution provides a simple, elegant and closed form solution to many problems in lifetime testing and reliability studies. One reason for the wide popularity of the exponential distribution in reliability is that it is the limiting lifetime distribution of a series system of substantially similar components. Also the exponential distribution is important due to its memory loss property. However, the exponential distribution does not provide a significant fit for some real lifetime applications, where the failure rates are not constant. Epstein1 has detailed study on exponential distribution and its role in life testing. Sato et al.2 obtained a discrete exponential distribution and discussed its properties and applied the distribution to model defect count distribution in semi-conductor deposition equipment and defect count distribution per chips. Gupta & Kundu3 have introduced generalized exponential distribution and discussed its properties and applications. Most of the works done on exponential distribution, its extension, mixture and applications in different fields of knowledge are available in Ahsanullah &Hamedani.4 Although exponential distribution was the first one parameter lifetime distribution for modeling lifetime data from engineering and biomedical sciences, but due to one parameter it is not always a suitable lifetime distribution. Also the hazard rate function and the mean residual life function of exponential distribution are always constant.

The statistical modeling and analysis of lifetime data are crucial in almost all branches of knowledge especially engineering and biomedical sciences. The one parameter exponential distribution is popular in statistics literature for modeling lifetime data. It has been observed that exponential distribution is not always a suitable model either due to theoretical or applied point of view. An attempt has been made to find a one parameter lifetime distribution which competes well with exponential distribution.

In this paper a quasi exponential distribution has been proposed. Its statistical properties including shapes of pdf for varying values of parameter, moments, hazard rate function, mean residual life function, stochastic ordering have been discussed. Estimation of parameter has been discussed using maximum likelihood estimation. Finally, application to a real lifetime data from biomedical science has been presented to test its goodness of fit over exponential distribution.

A quasi exponential distribution

A quasi exponential distribution (QED) having scale parameter θθ is defined by its pdf

&f2(x;θ)=2θ1/2Γ(12)e-θx2;x>0,θ>0f2(x;θ)=2θ1/2Γ(12)eθx2;x>0,θ>0           (2.1)

The survival (reliability) function of QED (2.1) can be obtained as

S(x;θ)=P(X>x)=2θ1/2Γ(12)xe-θt2dt=Γ(12,θx2)Γ(12);x>0,θ>0S(x;θ)=P(X>x)=2θ1/2Γ(12)xeθt2dt=Γ(12,θx2)Γ(12);x>0,θ>0         (2.2)

Where the function  is the upper incomplete gamma function defined as

Γ(α,z)=zeyyα1dy;α>0,z0Γ(α,z)=zeyyα1dy;α>0,z0,          (2.3)

Thus the cdf of QED (2.1) can be defined as

F2(x;θ)=1-S(x;θ)=1-Γ(12,θx2)Γ(12);x>0,θ>0F2(x;θ)=1S(x;θ)=1Γ(12,θx2)Γ(12);x>0,θ>0     (2.4)

The behavior of the pdf and cdf of QED and exponential distribution (ED) for varying values of parameter have been shown in figures 1 and 2 respectively.

  • Figure 1 Behavior of the pdf of the QED and ED for varying varying values of parameter θθ .

  • Figure 2Behavior of the cdf of the QED and ED for varying values of parameter . θθ .

Moments

The rth moment about origin of QED can be obtained as

μr=E(Xr)=2θ1/2Γ(12)0e-θx2xrdx=Γ(r+12)θr/2Γ(12);r=1,2,3,....μr=E(Xr)=2θ1/2Γ(12)0eθx2xrdx=Γ(r+12)θr/2Γ(12);r=1,2,3,....      (3.1)

&μr=E(Xr)=2θ1/2Γ(12)0e-θx2xrdx=Γ(r+12)θr/2Γ(12);r=1,2,3,....μr=E(Xr)=2θ1/2Γ(12)0eθx2xrdx=Γ(r+12)θr/2Γ(12);r=1,2,3,....

Taking r=1,2,3,and4r=1,2,3,and4;in (3.1), the first four moments about origin of QED are obtained as

μ1=1θ1/2Γ(12)μ1=1θ1/2Γ(12), μ2=12θμ2=12θ, μ3=1θ3/2Γ(12)μ3=1θ3/2Γ(12);and μ4=34θ2μ4=34θ2

The variance of QED can be obtained as

μ2=(Γ(12))2-22θ(Γ(12))2μ2=(Γ(12))222θ(Γ(12))2.

The third and fourth central moments are not being given here because their expressions are lengthy. However, they can be calculated using relationshipμr=E(Y-μ1)r=rk=0(rk)μk(-μ1)r-k ;between central moments and moments about mean.

Hazard rate function and mean residual life function

Let X be a continuous random variable with pdf f(x) and cdf F(x) . The hazard rate     function (also known as the failure rate function) and the mean residual life                 function of are respectively defined as

h(x)=limΔx0P(X<x+Δx|X>x)Δx=f(x)1-F(x)          (4.1)

and  m(x)=E[Xx|X>x]=11F(x)x[1F(t)]dt =1S(x)xtf(t)dtx        (4.2)

The corresponding hazard rate function,h(x) and the mean residual life function,m(x) of QED (2.1) are obtained as

h(x)=2θ1/2e-θx2Γ(12,θx2)        (4.3)

and

m(x)=1S(x;θ)xtf2(t;θ)dtx=Γ(1,θx2)θ1/2Γ(12,θx2)x =eθx2θ1/2Γ(12,θx2)x       (4.4)

The hazard rate function and the mean residual life function of the QED show flexibility over exponential distribution because in case of exponential distribution these measures are always constant. The behavior of  and  of QED for varying values of parameter are shown in Figure 3 &Figure 4.

Figure 3 Behavior of the hazard rate function of QED for varying values of parameter θ .

Figure 4Behavior of the mean residual life function of QED for varying values of parameter θ .

Stochastic ordering

Stochastic ordering of positive continuous random variables is an important tool for judging their comparative behavior.6 A random variable X is said to be smaller than a random variable Y in the

  1. stochastic order (XstY) if FX(x)FY(x) for all x
  2. hazard rate order (XhrY) if hX(x)hY(x) for all x
  3. mean residual life order (XmrlY) ifmX(x)mY(x) for all x
  4. likelihood ratio(XlrY) order if fX(x)fY(x) decreases inx .

The following results due to Shaked & Shanthikumar5 are well known for establishing stochastic ordering of distributions

XlrYXhrYXmrlYXstY

The QED is ordered with respect to the strongest ‘likelihood ratio ordering’ as established in the following theorem

Theorem: Let XQED(θ1,α1) and YQED(θ2,α2) If θ1>θ2 , then XlrY and henceXhrY XmrlY , andXSTY .

Proof: We have

fX(x;θ1,α1)fY(x;θ2,α2)=(θ1θ2)1/2e-(θ1-θ2)x2;x>0

Now, lnfX(x;θ1,α1)fY(x;θ2,α2)=12ln(θ1θ2)-(θ1-θ2)x2

This gives

ddx{lnfX(x;θ1,α1)fY(x;θ2,α2)}=-2(θ1-θ2)x

Thus ifθ1>θ2 ,ddx{lnfX(x;θ1,α1)fY(x;θ2,α2)}<0 . This means that XlrY and hence XhrY , XmrlY and XstY . This shows flexibility of QED over exponential distribution.

Estimation of parameter

In this section, the estimation of parameters of QED has been discussed using method of maximum likelihood. Let (x1,x2,...,xn) be a random sample of size from the QED (2.1). The natural log likelihood function of QED (2.1) can be given by

lnL=n[ln2+12lnθ-lnΓ(12)]-θni=1xi2

The maximum likelihood estimate (MLE) of the parameter  of QED (2.1) is the solutions of log likelihood equation

dlnLdθ=n2θ-ni=1xi2=0

which gives

θ=n2ni=1xi2=12nni=1xi2

Note that the method of moment estimate θ of the parameter θ is also the same as the MLE of the parameter.

Goodness of fit

The goodness of fit of the QGD has been explained with a real dataset from medical science. The following data represents the lifetime’s data relating to relief times (in minutes) of 20 patients receiving an analgesic and reported by Gross & Clark.6

1.1          1.4          1.3          1.7          1.9          1.8          1.6          2.2          1.7          2.7          4.1          1.8

1.5          1.2          1.4          3.0          1.7          2.3          1.6          2.0

For this data set, QED has been fitted along with one parameter exponential distribution.

The ML estimates,-2lnL values of , Akaike Information criteria (AIC), and K-S statistics of the fitted distributions are presented in Table 1. The AIC and K-S Statistics are computed using the following formulae:

AIC=-2lnL+2k ;and K-S=Supx|Fn(x)-F0(x)| ,

Where k= the number of parameters, n= the sample size,Fn(x) is the empirical (sample) cumulative distribution function, and F0(x) is the theoretical cumulative distribution function. The best distribution is the distribution corresponding to lower values of , AIC, and K-S statistics. It is clear from the goodness of fit in Table 1 that QED gives better fit than exponential distribution and hence it can be considered an important lifetime distribution to model lifetime data from medical science.

The profiles of the likelihood estimate θ of QED for the given dataset is presented in Figure 5. The fitted plots of the considered distributions for the given dataset are presented in Figure 6.

 

Model

ML estimate(θ)

S.E(θ)

 

- 2 log L

 

AIC

 

K-S

QED

0.122518

0.03874

57.16

59.16

0.461

Exponential

0.526314

0.11768

65.67

67.67

0.471

Table 1 Summary of ML estimates, -2lnL , AIC, and K-S of the fitted distributions

Figure 5 Profiles of the likelihood estimates of of QED for the considered dataset.

Figure 6 Fitted plots of the quasi exponential and exponential distributions for the given dataset.

Concluding remarks

This paper proposes a quasi exponential distribution for modeling real lifetime data from medical science. Its statistical properties including moments, survival function, hazard rate function, mean residual life function and stochastic ordering have been studied. The estimation of parameter has been discussed using maximum likelihood. Profile of the likelihood estimate of the parameter for QED has been shown graphically. Application of the distribution has been explained with a real lifetime data from biomedical science and its goodness of fit shows quite satisfactory over exponential distribution.

Acknowledgement

None.

Conflict of interest

Authors declares there is no conflict of interest.

References

Creative Commons Attribution License

©2018 Shanker, et al. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.