Submit manuscript...
eISSN: 2378-315X

Biometrics & Biostatistics International Journal

Research Article Volume 8 Issue 1

A two-parameter power Rama distribution with properties and applications

Berhane Abebe, Mussie Tesfay, Tesfalem Eyob, Rama Shanker

Department of Statistics, College of Science, Eritrea Institute of Technology, Eritrea

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

Received: December 29, 2018 | Published: January 21, 2019

Citation: Abebe B, Tesfay M, Eyob T, et al. A two-parameter power Rama distribution with properties and applications. Biom Biostat Int J. 2019;8(1):6-11. DOI: 10.15406/bbij.2019.08.00262

Download PDF

Abstract

In this paper a two-parameter Power Rama distribution, which includes one parameter Rama distribution as a particular case, has been proposed. Its statistical and reliability properties including shapes of density for varying values of parameters, the moments about origin, the mean and variance, survival function, hazard rate function, mean residual life function have been discussed. The maximum likelihood estimation for estimating the parameters has been discussed. Finally, the goodness of fit of the proposed distribution has been discussed with two real lifetime dataset.

Keywords: Rama distribution, Hazard rate function, moments, maximum likelihood estimation, goodness of fit

Introduction

Shanker1 proposed a continuous lifetime distribution named Rama distribution defined by its probability density function (pdf) and cumulative distribution function (cdf)

f(y;θ)=θ4θ3+6(1+y3)eθy;​​​​​​​​   y0,θ0   (1.1)

F(y,θ)=1[1+θ3y3+3θ2y2+6θyθ3+6]eθy   ;x>0,θ>0   (1.2)

The pdf (1.1) is a convex combination of exponential (θ) and gamma (4,θ) distribution with mixing proportion θ3θ3+6 . We have

f(y,θ)=pg1(y)+(1p)g2(y)

where

p=θ3θ3+6,  g1(y)=θeθy,  and g2(y)=θ4y3eθy6.

The statistical properties, estimation of parameter using maximum likelihood estimation and applications for modeling lifetime data of Rama distribution are available in Shanker.1

The pdf and the cdf of Power Lindley distribution (PLD) introduced by Ghitany et al.,2 are given by

f(x;θ,α)=αθ2(θ+1)xα1(1+xα)eθxα  ;x>0,θ>0,α>0   (1.3)

F(x;θ,α)=1[1+θxαθ2+1]eθxα;x>0,θ>0,α>0   (1.4)

Obviously at α=1, PLD reduces to Lindley distribution, introduced by Lindley3 having pdf and cdf

f(x;θ)=θ2θ+1(1+x)eθx;x>0,θ>0   (1.5)

F(x;θ)=1[1+θxθ+1]eθx;x>0,θ>0   (1.6)

The statistical properties, estimation of parameter and application of Lindley distribution are discussed in Ghitany et al.,4 Shanker et al.,5 have detailed study on applications of Lindley distribution and exponential distribution to model real lifetime datasets from engineering and biomedical sciences. Since Rama distribution has only one parameter, it has less flexibility to model data of varying natures. In the present paper an attempt has been made to derive two-parameter power Rama distribution which includes one parameter Rama distribution as particular cases as power transformation of Rama distribution. The shapes of the density, moments, hazard rate function, and mean residual life function of the distribution have been discussed. The maximum likelihood estimation has been explained. The goodness of fit of the proposed distribution has been discussed with two real lifetime dataset and fit shows quite satisfactory fit over other one parameter and two-parameter lifetime distributions.

A two-parameter power Rama distribution

Taking the power transformation X=Y1/α in (1.1), the pdf of the random variable X can be obtained as

f1(x;  θ,α)=αθ4θ3+6(1+x3α)xα1eθxα;x0,  θ0,  α0   (2.1)

=pg1(x;θ,α)+(1p)g2(x;θ,α)   (2.2)

where

p=θ3θ3+6,g1(x;θ,α)=αθeθxαxα1;       x>0,θ>0,α>0

and

g2(x;θ,α)=αθ4eθxαx4α16;     x>0,θ>0,α>0

We would call the pdf in (2.1) power Rama distribution (PRD). At α=1, (2.1) reduces to one Rama distribution. Like Rama distribution, the PRD is also a convex combination of Weibull distribution (with shape parameter α and scale parameter θ) and a generalized gamma distribution (with shape parameter 4 and α and scale α parameter θ) with mixing proportion p=θ3θ3+6.

The corresponding cdf of PRD can be obtained as

F1(x;θ,α)=1θ3eθxα+Γ(4,θxα)θ3+6  x>0,θ>0  ,α>0   (2.3)

Various graphs of the pdf of PRD for varying values of parameters have been drawn and presented in Figure 1. As the value of the shape parameter alpha increases the graph of PRD approaches normal distribution. Various graphs of the cdf of PRD for varying values of its parameters have been drawn and shown in Figure 2.

Figure 1 Behavior of the pdf of PRD for varying values of parameters θ and α.

Figure 2 Behavior of the cdf of PRD for varying values of parameters θ and α.

Statistical constants

Using (2.2), the rth moment about the origin, of PRD can be obtained as

μr=E(Xr)=p0xrg1(x;θ,α)dx+(1p)0xrg2(x;θ,α)dx

=θ3Γ(rα+1)+Γ(rα+4)θrα(θ3+6)  ;r=1,2,3,...   (3.1)

At α=1, the above expression reduces to the rth moment about origin of Rama distribution given by

μr=θ3Γ(r+1)+Γ(r+4)θr(θ3+6)  ;r=1,2,3,...

Thus the first four moments about origin of the PRD are obtained as

μ1=1αΓ(1α)[θ3+(1α+3)(1α+2)(1α+1)]θ1α(θ3+6)

μ2=2αΓ(2α)[θ3+(2α+3)(2α+2)(2α+1)]θ2α(θ3+6)

μ3=3αΓ(3α)[θ3+(3α+3)](3α+2)(3α+1)θ3α(θ3+6)

μ4=4αΓ(4α)[θ3+(4α+3)](4α+2)(4α+1)θ4α(θ3+6)

The variance of the PRD thus can be expressed as

σ2=1θ2α(θ3+6)[2αΓ(2α)[θ3+(2α+3)(2α+2)(2α+1)]{1αΓ(1α)[θ3+(1α+3)(1α+2)(1α+1)]}2(θ3+3)]

Reliability properties

The survival function, S(x;θ,α) of PRD can be obtained as

S(x;θ,α)=P(X>x)=xf(t;θ,α)dt

=θ3eθxα+Γ(4,θxα)θ3+6  ;    x>0,θ>0,α>0

The hazard rate function, h(x;θ,α) and the mean residual function, m(x;θ,α) of PRD are obtained as

h(x;θ,α)=f(x;θ,α)S(x;θ,α)=αθ4(1+x3α)xα1eθxαθ3eθxα+Γ(4,θxα)

and

m(x;θ,α)=1S(x;θ,α)xtf(t;θ,α)dtx

=1αΓ(1α,θxα)[θ3+(1α+1)(1α+2)](1α+3)+eθxα(θxα)[θ3+(θxα)3(1α+3)(θxα)2+(1α+3)(1α+2)(θxα)]θ1α[Γ(4,θxα)+θ3eθxα]x

Obviously m(x;θ,α)=1αΓ(1α)[θ3+(1α+3)(1α+2)(1α+1)]θ1α(θ3+6)=μ1 for x=0, the mean of PRD .

The behaviors of hazard rate function of PRD for varying values of parameters have been shown graphically in Figure 3. The graphs of hazard rate function of PRD are monotonically increasing for varying values of parameters. The behaviors of mean residual life function of PRD for varying values of parameters are shown in Figure 4. The graphs of mean residual life function of PRD are monotonically deceasing for varying values of parameters.

Figure 3 Behavior of the hazard rate function of PRD for varying values of parameters θ and α.

Figure 4 Behavior of the hazard rate function of PRD for varying values of parameters θ and α.

Estimation of parameters

The natural log likelihood function based on random samples (x1,x2,...,xn) from PRD (θ,α) can be expressed as

lnL=i=1nlnf(xi;θ,α)

=n[4lnθ+lnαln(θ3+6)]+(α1)i=1nlnxi+i=1n(1+xi)θi=1nxiα

Thus, the maximum likelihood estimates, (θ^,α^) of parameters (θ,α) of PRD is the solution of the following log-likelihood equations

lnLθ=4nθ3nθ2θ3+6i=1nxi=0

lnLα=nαi=1nlnxiθi=1nxiα=0

Since these two natural log likelihood equations cannot be expressed in closed form, these two log-likelihood equations are not directly solvable. The of parameters can be obtained directly from the natural log likelihood equation using Newton-Raphson iteration using R-software till sufficiently close values of and are obtained.

Applications

The application and the goodness of fit of PRD based on maximum likelihood estimates of parameters have been explained with two real datasets from engineering and biomedical Sciences. The fit given by one parameter exponential, Lindley and Rama distributions and two-parameter power Lindley distribution (PLD) and Weibull distribution, introduced by Weibull,6 have also given for ready comparison. The following two real lifetime datasets have been considered for showing the superiority of PRD over other lifetime distributions.

Data set 1: The data set is from Smith & Naylor7 relating to the strength of 1.5cm glass fibers measured at the National Physical Laboratory, England.

0.55 0.93 1.25 1.36 1.49 1.52 1.58 1.61 1.64 1.68 1.73 1.81 2.00 0.74 1.04 1.27 1.39 1.49 1.53 1.59 1.61 1.66 1.68 1.76 1.82 2.01 0.77 1.11 1.28 1.42 1.50 1.54 1.60 1.62 1.66 1.69 1.76 1.84 2.24 0.81 1.13 1.29 1.48 1.50 1.55 1.61 1.62 1.66 1.70 1.77 1.84 0.84 1.24 1.30 1.48 1.51 1.55 1.61 1.63 1.67 1.70 1.78 1.89

In order to compare considered distributions, values of MLE(θ^,α^) along with their standard errors, -2ln L, AIC (Akaike Information Criterion), K-S statistic and p-value for the dataset have been computed and presented in Table 1. The best distribution is the distribution corresponding to lower values of -2ln L, AIC and K-S statistic.

Model

ML estimate

se(θ^,α^)

-2ln L

AIC

K-S

p-value

PRD

θ^=0.8436
α^=3.2697

0.0942
0.2241

26.19

30.19

0.1297

0.2398

PLD

θ^=0.2224
α^=4.4583

0.0466
0.3871

29.38

33.38

0.1442

0.1457

Weibull

θ^=0.0598
α^=5.7800

0.0205
0.5752

30.41

34.31

0.1523

0.1075

Rama

θ^=1.7313

0.0991

169.72

171.72

0.3580

< 0.0001

Lindley

θ^=0.9961

0.0948

162.56

164.56

0.3864

< 0.0001

Exponential

θ^=0.6636

0.0836

177.66

179.66

0.4180

< 0.0001

Table 1 MLE"s, -2ln L, se(θ^,α^), AIC, K-S statistic, and p-value of the fitted distribution of dataset 1

The Variance-Covariance matrix and 95% confidence interval (CI"s) for the parameters and of PRD for the dataset 1 has been presented in Table 2.

Parameters

Variance-covariance matrix

95% CI

θ^

α^

Lower

Upper

θ^

0.0089

-0.0179

0.0668

1.0367

α^

-0.0179

0.0504

2.8450

3.7278

Table 2 Variance-Covariance matrix and 95% confidence interval (CI"s) for the parameters and of PRD

In order to see the closeness of the fit given by one parameter exponential, Lindley and Rama distributions and two-parameter PLD and Weibull distribution, the fitted plot of pdfs of these distributions for the dataset 1 have been shown in Figure 5. Obviously the fitted plots of the distributions and the histogram of the original dataset shows that PRD gives much closer fit over the considered distributions.

Figure 5 Fitted pdf plots of the considered distributions for the given dataset 1.

Data set 2: The following dataset represents waiting time (in minutes) of 154 patients, waiting before OPD (Out Patient Diagnosis) from the 25th - 30th December, 2017(in the 4th week of December) at Halibet Hospital., available in the master thesis of Berhane Abebe,8 Department of Statistics, College of Science, Eritrea Institute of Technology, Eritrea.

2(13)

3(29)

4(32)

5(29)

6(18)

7(15)

8(6)

9(6)

10(2)

11(2)

12

17

The maximum likelihood estimates of parameters, standard error of estimates of parameters, AIC, K-S and p-values of the fitted distributions for dataset 2 are shown in Table 3. The Variance-Covariance matrix and 95% confidence interval (CI"s) for the parameters θ^ and α^ of PRD for the given dataset 2 has been presented in Table 4.

Model

ML Estimate

se(θ^,α^)

-2ln L

AIC

K-S

P-value

PRD

θ^=0.5378
α^=1.2173

0.0555
0.0576

659.57

663.57

0.1329

0.0087

PLD

θ^=0.1054
α^=1.7223

0.0173
0.0873

665.69

669.69

0.1368

0.0062

Weibull

θ^=0.0189
α^=2.2892

0.0049
0.1271

673.14

677.14

0.1402

0.0047

Rama

θ^=0.7575

0.0292

675.87

677.87

0.1581

0.0009

Lindley

θ^=0.3475

0.0201

746.74

748.74

0.2903

< 0.0001

Exponential

θ^=0.1995

0.0161

804.51

806.51

0.3659

< 0.0001

Table 3 MLE"s, -2ln L, se(θ^,α^), AIC, K-S statistic, and P-value of the fitted distribution of dataset 2

Parameters

Variance-covariance matrix

95% CI

θ^

α^

Lower

Upper

θ^

0.0031

-0.0030

0.4355

0.6526

α^

-0.0030

0.0033

1.1080

1.3339

Table 4 Variance-Covariance matrix and 95% confidence interval (CI"s) for the parameters θ^ and α^ of PRD for dataset 2

In order to see the closeness of the fit given by one parameter exponential, Lindley and Rama distributions and two-parameter PLD and Weibull distribution the fitted plot of pdfs of these distributions for the dataset 2 have been shown in Figure 6. Obviously the fitted plots of the distributions and the histogram of the dataset 2 shows that PRD gives much closer fit over the considered distributions.

Figure 6 Fitted pdf plots of the considered distributions for the given dataset 2.

Concluding remarks

A two-parameter Power Rama distribution (PRD) has been proposed which includes Rama distribution, introduced by Shanker (2017), as a particular case. The statistical and reliability properties including moments, survival function, hazard rate function, mean residual function of PRD have been discussed. Maximum likelihood estimation has been explained for estimating the parameters. The goodness of fit of PRD has been discussed with two real lifetime datasets and the fit has been found quite satisfactory over one parameter Rama, Lindley and Exponential distributions and two-parameter Power Lindley distribution (PLD) and Weibull distribution. Hence, PRD can be considered an important two-parameter lifetime distribution.

Acknowledgements

None.

Conflict of interest

The author declares there is no conflicts of interest.

References

Creative Commons Attribution License

©2019 Abebe, 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.