Submit manuscript...
eISSN: 2378-315X

Biometrics & Biostatistics International Journal

Research Article Volume 11 Issue 5

Uma distribution with properties and applications

Rama Shanker

Department of Statistics, Assam University, Silchar, Assam, India

Correspondence: Rama Shanker, Department of Statistics, Assam University, Silchar, Assam, India

Received: December 09, 2022 | Published: December 22, 2022

Citation: Shanker R. Uma distribution with properties and applications. Biom Biostat Int J. 2022;11(5):165-169. DOI: 10.15406/bbij.2022.11.00372

Download PDF

Abstract

The stochastic natures of lifetime data are really a challenge for statistician to search a suitable distribution for modeling and analysis of lifetime data. Keeping in mind the stochastic natures of lifetime data, a new lifetime distribution named Uma distribution has been suggested. Its several statistical properties, estimation of parameter and applications have been discussed. Applications of the distribution have been presented with three datasets and the goodness of fit of Uma distribution has been compared with exponential, Lindley, Shanker, Akash and Sujatha distributions.

Keywords: lifetime distributions, statistical properties, estimation of parameter, applications

Introduction

Due to stochastic nature of lifetime data, the search for lifetime distribution in the field of lifetime data analysis is expanding exponentially and getting popularity among policy makers to model data. In recent decades several lifetime distributions have been suggested in statistics literature. For example, Lindley distribution by Lindley,1 Shanker distribution by Shanker,2 Akash distribution by Shanker,3 and Sujatha distribution by Shanker,4 is some among others. Shanker et al.5 discussed the modeling of lifetime data using exponential and Lindley distributions and observed that there are some datasets in which these two distributions do not give good fit. Further, Shanker et al.6 put an effort to have comparative study on modeling of lifetime data using exponential, Lindley and Akash distribution and found that Akash distributions gives much better fit than both exponential and Lindley distribution but still there are some data sets in which these three distributions do not give good fit. Then, Shanker and Hagos,7 tried to model the real lifetime datasets using exponential, Lindley, Shanker and Akash and observed that still there are some datasets in which these distributions do not give good fit. Flexibility and tractability are the two important characteristics of a lifetime distributions and if the existing distributions are not flexible or tractable for the given dataset, then the search for a new distribution starts. Sometimes, data are being transformed to satisfy some assumptions of the distribution so that distribution fits well. But this is not useful practice because the original nature of the dataset is lost. Therefore, the most preferable is to search a distribution which fits the given data well than to modify the existing distributions.

While testing the goodness of fit of some well-known one parameter lifetime distributions available in literature, it has been observed that the existing distributions do not fit the data well. In this paper, in the search for a new distribution, we propose a new distribution named Uma distribution which fits the data well over the existing distributions. The statistical properties, estimation of

parameter and applications of the distribution has been presented systematically. It is hoped and expected that the distribution will draw attention of researchers to model lifetime data and preferred over the existing one parameter lifetime distributions.

Uma distribution

Taking the convex combination of exponential (θ)(θ) , gamma (2,θ)(2,θ) and gamma (4,θ)(4,θ) with respective mixing proportions θ3θ3+θ2+6,θ2θ3+θ2+6θ3θ3+θ2+6,θ2θ3+θ2+6 and , 6θ3+θ2+66θ3+θ2+6 a probability density function (pdf) can be expressed as

f(x;θ)=θ4θ3+θ2+6(1+x+x3)eθx;x>0,θ>0f(x;θ)=θ4θ3+θ2+6(1+x+x3)eθx;x>0,θ>0

We would call this distribution as ‘Uma distribution’. Since it is a convex combination of exponential and gamma distributions, it is expected to give better fit over exponential and gamma distribution and other distributions developed using convex combinations of exponential and gamma distribution. The cumulative distribution function (cdf) of Uma distribution can be obtained as

F(x;θ)=1[1+θx(θ2x2+3θx+θ2+6)θ3+θ2+6]eθx;x>0,θ>0F(x;θ)=1[1+θx(θ2x2+3θx+θ2+6)θ3+θ2+6]eθx;x>0,θ>0

The behaviour of the pdf and the cdf of Uma distribution for varying values of parameter θθ have been presented in Figures 1,2 respectively.

Figure 1 Graphs of the pdf of Uma distribution for selected values of the parameter.

Figure 2 Graphs of the cdf of Uma distribution for selected values of the parameter.

Reliability properties

  • Hazard rate function

The hazard rate function of a random variable XX having pdf f(x;θ)f(x;θ) and cdf F(x;θ)F(x;θ) is defined as

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

Thus, the hazard rate function of Uma distribution can be obtained as

h(x)=θ4(1+x+x3)θ3(x3+x+1)+θ2(3x2+1)+6θx+6h(x)=θ4(1+x+x3)θ3(x3+x+1)+θ2(3x2+1)+6θx+6

This gives. h(0)=θ4θ3+θ2+6=f(0)h(0)=θ4θ3+θ2+6=f(0) The behaviour of the hazard rate function of Uma distribution for various values of parameter is shown in the following Figure 3.

Figure 3 Graphs of the hazarad rate function of Uma distribution for selected values of the parameter.

  • Mean residual life function

Let XX be a random variable over the support (0,)(0,) representing the lifetime of a system. Mean Residual life (MRL) function measures the expected value of the remaining lifetime of the system, provided it has survived up to time. Let us consider the conditional random variable Xx=(Xx|X>x);x>0Xx=(Xx|X>x);x>0 . Then, the MRL function, denoted by, m(x)m(x) is defined as

m(x)=E(Xx)=1S(x)x[1F(t)]dtm(x)=E(Xx)=1S(x)x[1F(t)]dt

The MRL function of Uma distribution can thus be obtained as

m(x)=1{θ3(x3+x+1)+θ2(3x2+1)6θx+6}eθxx[θ3(t3+t+1)+θ2(3t2+1)6θt+6]eθtdtm(x)=1{θ3(x3+x+1)+θ2(3x2+1)6θx+6}eθxx[θ3(t3+t+1)+θ2(3t2+1)6θt+6]eθtdt =θ3(x3+x+1)+θ2(6x2+2)+18θx+24θ{θ3(x3+x+1)+θ2(3x2+1)+6θx+6}=θ3(x3+x+1)+θ2(6x2+2)+18θx+24θ{θ3(x3+x+1)+θ2(3x2+1)+6θx+6}

This givesm(0)=θ3+2θ2+24θ(θ3+θ2+6)=μ1 . The behaviour of the mean residual life function of Uma distribution for various values of parameter θ is shown in the following Figure 4.

Figure 4 Graphs of the mean residual life function of Uma distribution for selected values of the parameter.

Reverse hazard rate and Mill’s ratio

The reverse hazard rate of a random variable X having pdf f(x;θ) and cdf F(x;θ) is defined as

hr(x)=f(x;θ)F(x;θ)

Thus, the reverse hazard rate function of Uma distribution can be obtained ashr(x)=θ4(1+x+x3)eθx(θ3+θ2+6)[(θ3+θ2+6)+θx(θ2x2+3θx+θ2+6)]eθx

 Mill’s ratio of a random X xmlns='http://www.w3.org/1998/Math/MathML'> X MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=Mj0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFf ea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadIfaaaa@37EB@ variable having pdf f(x;θ) and cdf F(x;θ) is defined as

 Mill’s ratio=1h(x)=1F(x;θ)f(x;θ)

Thus, the Mill’s ratio of Uma distribution can be obtained as1h(x)=θ3(x3+x+1)+θ2(3x2+1)+6θx+6θ4(1+x+x3)

Stochastic ordering

In Probability theory and statistics, a stochastic order quantifies the concept of one random variable being bigger than another. A random variable  is said to be smaller than a random variable  in the

  1. Stochastic order (XstY) for allx
  2. Hazard rate order (XhrY) if hX(x)hY(y) for allx
  3. Mean residual life order (XmrlY) if  mX(x)mY(y) for allx
  4. Likelihood ratio order (XlrY) if fX(x)fY(y) decrease inx

The following results due to Shaked and Shantikumar,8 are well known for establishing stochastic ordering of distributionsX<lrYX<hrYX<mrlYX<stY

Theorem: Let X~ Uma distribution (θ1) and Y~ Uma. (θ2) If θ1>θ2 , then X<lrY hence X<hrY ,X<mrlY and X<stY .

Proof: We have

fX(x)fY(x)==θ14(θ23+θ22+6)θ24(θ13+θ12+6)e(θ1θ2)x

We have ,log[fX(x)fY(x)]=log[θ14(θ23+θ22+6)θ24(θ13+θ12+6)](θ1θ2)x

Therefore,ddxlog[fX(x)fY(x)]=(θ1θ2)

Thus, for θ1>θ2 ,ddxlog[fX(x)fY(x)]<0 . this means that X<lrY hence X<hrY ,and X<stY .

Moments based descriptive measures

The r th moment about origin μr of Uma distribution can be obtained as μr=E(Xr)=θ4θ3+θ2+60xr(1+x+x2)eθxdx

=r!{θ3+(r+1)θ2+(r+1)(r+2)(r+3)}θr(θ3+θ2+6);r=1,2,3,

Substituting r=1,2,3,4 in the above equation, the first four moments about origin of Uma distribution can be obtained as

μ1=θ3+2θ2+24θ(θ3+θ2+6) ,μ2=2(θ3+3θ2+60)θ2(θ3+θ2+6)

μ3=6(θ3+4θ2+120)θ3(θ3+θ2+6) ,μ4=24(θ3+5θ2+210)θ4(θ3+θ2+6) .

The moments about the mean, using relationship between moments about the mean and the moments about the origin, can thus be obtained as

μ2=θ6+4θ5+2θ4+84θ3+60θ2+144θ2(θ3+θ2+6)2 μ3=2(θ9+6θ8+6θ7+200θ6+270θ5+108θ4+324θ3+432θ2+864)θ3(θ3+θ2+6)3 μ4=3(3θ12+24θ11+44θ10+968θ9+2336θ8+2016θ7+7488θ6+13248θ5+5760θ4+31104θ3+24192θ2+31104)θ4(θ3+θ2+6)4

The descriptive constants including coefficient of variation (CV), coefficient of skewness (CS), coefficient of kurtosis (CK) and the index of dispersion (ID) of Uma distribution are thus obtained as

CV=μ2μ1=θ6+4θ5+2θ4+84θ3+60θ2+144θ3+2θ2+24 CS=μ32μ23=4(θ9+6θ8+6θ7+200θ6+270θ5+108θ4+324θ3+432θ2+864)2(θ6+4θ5+2θ4+84θ3+60θ2+144)3 CK=μ4μ22=3(3θ12+24θ11+44θ10+968θ9+2336θ8+2016θ7+7488θ6+13248θ5+5760θ4+31104θ3+24192θ2+31104)(θ6+4θ5+2θ4+84θ3+60θ2+144)2 ID=μ2μ1=θ6+4θ5+2θ4+84θ3+60θ2+144θ(θ3+θ2+6)(θ3+2θ2+24)

Behaviour of coefficient of variation, coefficient of skewness, coefficient of kurtosis and index of dispersion for changing values of parameter are shown in the Figure 5.

Figure 5 Graph of CV, CS, CK and ID of Uma distribution for different values of the parameter.

Deviations from mean and median

Mean deviation about the mean and the mean deviation about median of a random variable X having pdf f(x) and cdf F(x) are defined byδ1(x)=0|xμ|f(x)dx=2μF(μ)2μ0xf(x)dx

 and δ2(x)=0|xM|f(x)dx=μ+2Mxf(x)dx respectively, where μ=E(X) andM=Median(X) .

Using pdf and expressions for the mean of Uma distribution, we get

μ0xf(x;θ)dx=μ[θ4(μ4+μ2+μ)+θ3(4μ3+2μ+1)+2θ2(6μ2+1)+24(θμ+1)]eθμθ(θ3+θ2+6) M0xf(x;θ)dx=μ[θ4(M4+M2+M)+θ3(4M3+2M+1)+2θ2(6M2+1)+24(θM+1)]eθMθ(θ3+θ2+6)

Using above expressions some algebraic simplifications, the meanδ1(x) deviation about the mean, and the mean deviation about the median δ2(x) of Uma distribution are obtained asδ1(x)=2[θ3μ3+6θ2μ2+θ3μ+18θμ+(θ3+2θ2+24)]eθ(θ3+θ2+6)θμ δ2(x)=2[θ4(M4+M2+M)+θ3(4M3+2M+1)+2θ2(6M2+1)+24(θM+1)]eθ(θ3+θ2+6)θMμ

Parameter estimation of Uma distribution

Suppose (x1,x2,x3,...,xn) be a random sample of size n from Uma distribution. The log likelihood function, L of Uma distribution is given by

logL=ni=1logf(xi;θ)=n{4logθlog(θ3+θ2+6)}+ni=1log(1+xi+xi3)nθˉx

The maximum likelihood estimate (MLE) (ˆθ) of the parameters (θ) of Uma distribution is the solution of the following log likelihood equation

dlogLdθ=4nθ(3θ2+2θ)nθ3+θ2+6nˉx=0

This gives

ˉxθ4+(ˉx1)θ32θ2+6ˉxθ24=0.

This is a fourth degree polynomial equation in θ . It should be noted that the method of moment estimate is also the same as that of the MLE. The above equation can easily be solved using Newton-Raphson method, taking the initial value of the parameterθ=0.5 .

Applications and goodness of fit

The applications and the goodness of fit of Uma distribution has been discussed with three datasets. Keeping in mind the flexibility and tractability of the distribution with the dataset following three datasets have been considered.

Data set 1: This data set represents the lifetime data relating to relief times (in minutes) of 20 patients receiving an analgesic and reported by Gross and Clark.9

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, 1.7, 2.3, 1.6, 2.0

Data Set 2: This data set is the strength data of glass of the aircraft window reported by Fuller et al.10:

18.83, 20.80, 21.657, 23.03, 23.23, 24.05, 24.321, 25.5, 25.52, 25.80, 26.69, 26.77, 26.78, 27.05, 27.67, 29.90, 31.11, 33.2, 33.73, 33.76, 33.89, 34.76, 35.75, 35.91, 36.98, 37.08, 37.09, 39.58, 44.045, 45.29, 45.381

Data Set 3: The following data represent the tensile strength, measured in GPa, of 69 carbon fibers tested under tension at gauge lengths of 20 mm(Bader and Priest, 1982)11:  

1.312, 1.314, 1.479, 1.552, 1.700, 1.803, 1.861, 1.865, 1.944, 1.958, 1.966, 1.997, 2.006, 2.021, 2.027, 2.055, 2.063, 2.098, 2.140, 2.179, 2.224, 2.240, 2.253, 2.270, 2.272, 2.274, 2.301, 2.301, 2.359, 2.382, 2.382, 2.426, 2.434, 2.435, 2.478, 2.490, 2.511, 2.514, 2.535, 2.554, 2.566, 2.570, 2.586, 2.629, 2.633, 2.642, 2.648, 2.684, 2.697, 2.726, 2.770, 2.773, 2.800, 2.809, 2.818, 2.821, 2.848, 2.880, 2.954, 3.012, 3.067, 3.084, 3.090, 3.096, 3.128, 3.233, 3.433, 3.585, 3.585 .

The values ML estimates of parameter, , AIC (Akaike Information Criterion), AICC (Akaike Information Criterion corrected), BIC (Bayesian Information criterion), K-S (Kolmogorov-Smirnov) for the considered distributions for the given datasets have been computed and presented in Tables 1–3 respectively.

It is clear from the goodness of fit in the Tables 1 to 3 that Uma distribution gives much better fit over exponential, Lindley, Shanker, Akash and Sujatha distributions.

Sl. No

Distributions                                                ˆθ

  2logL


 

AIC

AICC

BIC

K-S

1

Uma

1.6024

     38.61

40.61

40.83

41.60

0.238

2

Sujatha

1.1367

     57.50

59.50

59.72

60.49

0.309

3

Akash

1.1569

     59.52

61.52

61.74

62.51

0.320

4

Shanker

0.8038

     59.78

61.78

61.22

62.51

0.315

5

Lindley

0.8161

     60.50

62.50

62.72

63.49

0.341

6

Exponential

0.5263

     65.67

67.67

67.90

68.67

0.389

Table 1 ML estimates, , AIC, AICC, BIC, K-S of the distribution for the dataset-1

Sl. No

Distributions

   ˆθ

2logL

 

AIC

AICC

BIC

K-S

1

Uma

0.1299

232.54

234.54

234.67

235.97

0.233

2

Sujatha

0.0956

241.50

243.50

243.64

244.94

0.27

3

Akash

0.0971

240.68

242.68

242.82

244.11

0.266

4

Shanker

0.0647

252.35

254.35

254.49

255.78

0.326

5

Lindley

0.0629

253.99

255.99

256.13

257.42

0.333

6

Exponential

0.0325

274.53

276.53

276.67

277.96

0.426

Table 2 ML estimates, , AIC, AICC, BIC, K-S of the distributions for the dataset-2

Sl. No

Distributions

   ˆθ

2logL

AIC

AICC

BIC

K-S

1

Uma

1.3828

156.41

158.41

158.47

160.64

0.312

2

Sujatha

0.9361

221.61

223.61

223.67

225.84

0.348

3

Akash

0.9647

224.28

226.28

226.34

228.51

0.348

4

Shanker

0.6580

233.01

235.01

235.06

237.24

0.355

5

Lindley

0.6590

238.38

240.38

240.44

242.61

0.390

6

Exponential

0.4079

261.74

263.74

263.80

265.97

0.434

Table 3 ML estimates, , AIC, AICC, BIC, K-S of the distributions for the dataset-3

Conclusion and future works

A new lifetime distribution named Uma distribution has been suggested. Statistical properties, estimation of parameter and applications of the distribution has been presented. As the distribution is new one, it is expected and hoped that it will be of great use to statisticians working in the field of modeling lifetime data from different fields of knowledge.

Being a new lifetime distribution with flexibility, tractability and practicability, a lot of future works can be done on Uma distribution.

Acknowledgments

Author is really grateful to the Editor-In-Chief of the Journal and the anonymous reviewer for quick and valuable comments on the paper.

Conflicts of interest

There aren't any conflict of interests.

Funding

None.

References

  1. Lindley DV. Fiducial distributions and Bayes’ theorem. Journal of the Royal Statistical Society, Series B. 1958;20:102–107.
  2. Shanker R. Shanker distribution and its applications. International Journal of Statistics and Applications. 2015a;5(6):338–348.
  3. Shanker R. Akash Distribution and its applications. International Journal of Probability and Statistics. 2015b;4(3):65–75.
  4. Shanker R. Sujatha Distribution and its applications. Statistics in Transition-New series. 2016a;17(3):391–410.
  5. Shanker R, Hagos F, Sujatha S. On Modeling of Lifetimes data using Exponential and Lindley distributions. Biom Biostat Int J. 2015;2(5):140–147.
  6. Shanker R, Hagos F, Sujatha S. On modeling of lifetime data using one parameter Akash, Lindley and exponential distributions. Biom Biostat Int J. 2016;3(2):54–62.
  7. Shanker R, Hagos F. On modeling of lifetime data using Akash, Shanker, Lindley and exponential distributions. Biom & Biostat Int J. 2016;3(6):214–224.
  8. Shaked M, Shanthikumar JG. Stochastic Orders and Their Applications. Academic Press, New York: 1994.
  9. Gross AJ, Clark VA. Survival Distributions: Reliability Applications in the Biometrical Sciences. New York: John Wiley. 1975.
  10. Fuller EJ, Frieman S, Quinn J, et al. Fracture mechanics approach to the design of glass aircraft windows: A case study, SPIE Proc 1994;2286:419–430.
  11. Bader MG, Priest AM. Statistical aspects of fiber and bundle strength in hybrid composites. In: Hayashi T, Kawata K, Umekawa S, editors. Progress in Science in Engineering Composites, ICCM-IV. Tokyo; 1982:1129–1136.
Creative Commons Attribution License

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