Power Hamza Distribution with Application to Lifetime Data

Samuel U. Enogwe^1*, Chike H. Nwankwo², Eric U. Oti³

¹Department of Statistics, Michael Okpara University of Agriculture, Umudike, Nigeria.
²Department of Statistics, Nnamdi Azikiwe University, Awka, Nigeria.
³Department of Statistics, Federal Polytechnic, Ekowe, Nigeria.

DOI: 10.4236/jamp.2022.101004                

Abstract

In this paper, a three-parameter lifetime distribution named power Hamza distribution (PH) is proposed. The PH distribution is a useful generalization of the Hamza distribution which accommodates heavy-tailed, upside-down bathtub and J-shaped hazard rates making it more flexible than the Hamza distribution for modelling various kinds of lifetime data. A comprehensive account of the properties of this distribution is presented. The maximum likelihood estimators of the unknown model parameters are discussed. Finally, a real-life data is analyzed for illustrative purpose proving that the PH outperforms the Hamza distribution and several other lifetime distributions.

Keywords

Hamza Distribution, Hazard Rate Function, Conditional Moments, Stochastic Ordering, Maximum Likelihood Method

1. Introduction

The Hamza distribution was introduced by [1] with the cumulative distribution function (cdf) and the corresponding probability density function (pdf) respectively given by

$F_H\left(y\right)=1-\left[1+\frac{\theta x\left({\theta }^5{y}^5+6{\theta }^4{y}^4+30{\theta }^3{y}^3+120{\theta }^2{y}^2+360\theta y+720\right)}{\left(\alpha {\theta }^5+120\right)}\right]{\text{e}}^{-\theta y}$ (1)

and

$f_H\left(x\right)=\frac{{\theta }^6}{\alpha {\theta }^5+120}\left(\alpha +\frac{\theta }{6}{y}^6\right){\text{e}}^{-\theta y}$ (2)

for $y>0$, $\theta >0$ and $\alpha >0$.

It may be observed that the Hamza distribution was obtained by compounding the exponential distribution having scale parameter $\theta$ and the gamma distribution having shape parameter 7 and scale parameter $\theta$, with mixing proportions $p=\alpha {\theta }^5/\left(\alpha {\theta }^5+120\right)$ and $q=\alpha {\theta }^5/\left(\alpha {\theta }^5+120\right)$ such that $p+q=1$.

[1] studied the properties and applications of this distribution in the context of lifetime analysis, showing that the distribution is superior to Lindley distribution due to [2], Ishita distribution by [3] and Pranav distribution by [4], respectively.

The aim of this paper is to introduce a new distribution, called the power Hamza distribution, which is a direct generalization of the Hamza distribution. Some of the distributions proposed using the power transformation include the power Lindley due to [5], power Akash and Shanker proposed by [6] [7], power Ishita and power Aradhana due to [8] [9], power Rama and power Garima due to Abebe et al. [10] [11], power Pranav due to [12], power Sujatha by [13], power Prakaamy by [14]. From the literature reviewed in this paper, all the power transformed distributions were shown to be more flexible than their corresponding baseline distributions and more useful for analyzing complex data structures in various fields of life.

The rest of the paper is organized as follows. The pdf, cdf and hazard rate function of the new distribution is given in Section 2. Section 3 provides a comprehensive account of the properties of the distribution including the moment generating function, moments, skewness, kurtosis, mean residual lifetime, mean deviations, Bonferroni and Lorenz curves, stochastic ordering, entropy measure, stress-strength reliability, distributions and moments of order statistics. In Section 4, the maximum likelihood estimates of the parameters of the distribution are given. Also, Section 5 gives the asymptotic confidence intervals. Section 6 illustrates the proposed model in two real-datasets. The paper is concluded in Section 7.

2. The Power Hamza (PH) Distribution

The probability density function, cumulative distribution function and hazard function of the power Hamza distribution having parameters $\theta$, $\alpha$ and $\beta$ are provided in as Propositions 1, 2 and 3.

Proposition 1. A random variable X is said to have a PH distribution if its pdf is of the form

$f_{PH}\left(x\right)=\frac{\beta {\theta }^6}{\alpha {\theta }^5+120}\left(\alpha +\frac{\theta }{6}{x}^{6\beta }\right)x^{\beta -1}{\text{e}}^{-\theta x^{\beta }}$ ; $x>0$, $\theta >0$, $\alpha >0$, $\beta >0$ (3)

Proof. Given the distribution of the Hamza random variable Y defined in (2). Assume that another random variable X is related to Y by the power function $X=g\left(Y\right)=Y^{1/\beta }$. Then the distribution of X is the power Hamza distribution. To derive the distribution of X, we notice that X is a one-to-one function of Y and so, $x=0$ when $y=0$ and $x=\infty$ when $y=\infty$, which implies that the support of the distribution of X is $\left(0,\infty \right)$.

Letting $x=g\left(y\right)=y^{1/\beta }$, and $y=g^{-1}\left(x\right)=x^{\beta }$ in (2), gives

$f_H\left(g^{-1}\left(x\right)\right)=\frac{{\theta }^6}{\alpha {\theta }^5+120}\left(\alpha +\frac{\theta }{6}{x}^{6\beta }\right){\text{e}}^{-\theta x^{\beta }}$, $x>0$ (4)

According to [15], the probability density function of a continuous random variable $X=g\left(Y\right)$ is gotten by

$f\left(x\right)=f\left(g^{-1}\left(x\right)\right)\left|\frac{\text{d}y}{\text{d}x}\right|$ (5)

Substituting (4) and $\text{d}y/\text{d}x=\beta x^{\beta -1}$ into (5), one obtains

$f_{PH}\left(x\right)=\frac{{\theta }^6}{\alpha {\theta }^5+120}\left(\alpha +\frac{\theta }{6}{x}^{6\beta }\right){\text{e}}^{-\theta x^{\beta }}\times \left|\beta x^{\beta -1}\right|$ ; $x>0$ (6)

Further simplification of (6) yields the pdf of the power Hamza random variable X defined in (3), and the proof of Proposition 1 is complete.

Corollary 1. Let X be a power Hamza random variable, then the function defined by Equation (3) is a pdf.

Proof. We show that $f_{PH}\left(x\right)\ge 0\forall x$ and ${\displaystyle \underset{0}{\overset{\infty }{\int }}{f}_{PH}\left(x\right)\text{d}x}=1$ are satisfied.

1) $f_{PH}\left(x\right)\ge 0\forall x$

2) $\begin{array}{c}{\displaystyle \underset{0}{\overset{\infty }{\int }}{f}_{PH}\left(x\right)\text{d}x}={\displaystyle \underset{0}{\overset{\infty }{\int }}\frac{\beta {\theta }^6}{\alpha {\theta }^5+120}\left(\alpha +\frac{\theta }{6}{x}^{6\beta }\right)x^{\beta -1}{\text{e}}^{-\theta x^{\beta }}\text{d}x}\\ =\frac{\beta {\theta }^6}{\alpha {\theta }^5+120}\left[\alpha {\displaystyle \underset{0}{\overset{\infty }{\int }}{x}^{\beta -1}{\text{e}}^{-\theta x^{\beta }}\text{d}x}+\frac{\theta }{6}{\displaystyle \underset{0}{\overset{\infty }{\int }}{x}^{7\beta -1}{\text{e}}^{-\theta x^{\beta }}\text{d}x}\right]\end{array}$

By setting $x=y^{1/\beta }$ in the above integral and noting that ${\displaystyle \underset{0}{\overset{\infty }{\int }}{y}^{\alpha -1}{\text{e}}^{-\theta y}\text{d}y}={\Gamma \left(\alpha \right)/\theta }^{\alpha }$ one obtains

$\begin{array}{c}{\displaystyle \underset{0}{\overset{\infty }{\int }}{f}_{PH}\left(x\right)\text{d}x}=\frac{\beta {\theta }^6}{\alpha {\theta }^5+120}\left[\frac{\alpha }{\beta }{\displaystyle \underset{0}{\overset{\infty }{\int }}{\text{e}}^{-\theta y}\text{d}y}+\frac{\theta }{6\beta }{\displaystyle \underset{0}{\overset{\infty }{\int }}{y}^6{\text{e}}^{-\theta y}\text{d}y}\right]\\ =\frac{\beta {\theta }^6}{\alpha {\theta }^5+120}\left[\frac{\alpha }{\beta }\frac{\Gamma \left(1\right)}{\theta }+\frac{\theta }{6\beta }\frac{\Gamma \left(7\right)}{{\theta }^7}\right]\\ =\frac{\beta {\theta }^6}{\alpha {\theta }^5+120}\left[\frac{\alpha }{\beta \theta }+\frac{120}{\beta {\theta }^6}\right]\\ =\frac{\beta {\theta }^6}{\alpha {\theta }^5+120}\left[\frac{\alpha {\theta }^5+120}{\beta {\theta }^6}\right]=1\end{array}$

Henceforth, a random variable X that follows the distribution in (3) is symbolized by $X~PH\left(\theta ,\alpha ,\beta \right)$. The power Hamza distribution reduces to the Hamza distribution when $\beta =1$.

Proposition 2. For $x>0$, $\theta >0$, $\alpha >0$ and $\beta >0$, the cdf of $X~PH\left(\theta ,\alpha ,\beta \right)$ is given by

$F_{PH}\left(x\right)=1-\left[1+\frac{\theta x^{\beta }\left({\theta }^5{x}^{5\beta }+6{\theta }^4{x}^{4\beta }+30{\theta }^3{x}^{3\beta }+120{\theta }^2{x}^{2\beta }+360\theta x^{\beta }+720\right)}{\left(\alpha {\theta }^5+120\right)}\right]{\text{e}}^{-\theta x^{\beta }}$ (7)

Proof. [15] defines the cdf of a continuous random variable X as

$F\left(x\right)=P\left(X\le x\right)={\displaystyle \underset{0}{\overset{x}{\int }}f\left(u\right)\text{d}u}$ (8)

Substituting (3) into (8) leads to

$F_{PH}\left(x\right)=\frac{\beta {\theta }^6}{\alpha {\theta }^5+120}{\displaystyle \underset{0}{\overset{x}{\int }}\left(\alpha +\frac{\theta }{6}{u}^{6\beta }\right)u^{\beta -1}{\text{e}}^{-\theta u^{\beta }}\text{d}u}$ (9)

Letting $x=u^{1/\beta }$, $u=x^{\beta }$, $\text{d}u/\text{d}x=\beta x^{\beta -1}$, the integral (9) becomes

$F_{PH}\left(x\right)=\frac{\beta {\theta }^6}{\alpha {\theta }^5+120}\left[\frac{\alpha }{\theta }{\displaystyle \underset{0}{\overset{\theta x^{\beta }}{\int }}{\text{e}}^{-x}\text{d}x}+\frac{1}{6{\theta }^6}{\displaystyle \underset{0}{\overset{\theta x^{\beta }}{\int }}{x}^6{\text{e}}^{-x}\text{d}x}\right]$ (10)

Applying direct integration to the first part of the square bracket in (10) and integration by parts to the second part the square bracket in (10) gives

${\displaystyle \underset{0}{\overset{\theta x^{\beta }}{\int }}{\text{e}}^{-x}\text{d}x}=1-{\text{e}}^{-\theta x^{\beta }}$ (11)

and

$\begin{array}{l}{\displaystyle \underset{0}{\overset{\theta x^{\beta }}{\int }}{x}^6{\text{e}}^{-x}\text{d}x}=-{\left(\theta x^{\beta }\right)}^6{\text{e}}^{-\theta x^{\beta }}-6{\left(\theta x^{\beta }\right)}^5{\text{e}}^{-\theta x^{\beta }}-30{\left(\theta x^{\beta }\right)}^4{\text{e}}^{-\theta x^{\beta }}\\ -120{\left(\theta x^{\beta }\right)}^3{\text{e}}^{-\theta x^{\beta }}-360{\left(\theta x^{\beta }\right)}^2{\text{e}}^{-\theta x^{\beta }}-720\theta x^{\beta }{\text{e}}^{-\theta x^{\beta }}+720\end{array}$ (12)

respectively. Substituting (11) and (12) into (10), we obtain

$\begin{array}{l}{F}_{PH}\left(x\right)=\frac{\beta {\theta }^6}{\alpha {\theta }^5+120}[\frac{\alpha }{\theta }\left(1-{\text{e}}^{-\theta x^{\beta }}\right)+\frac{1}{6{\theta }^6}(-{\left(\theta x^{\beta }\right)}^6{\text{e}}^{-\theta x^{\beta }}-6{\left(\theta x^{\beta }\right)}^5{\text{e}}^{-\theta x^{\beta }}\\ -30{\left(\theta x^{\beta }\right)}^4{\text{e}}^{-\theta x^{\beta }}-120{\left(\theta x^{\beta }\right)}^3{\text{e}}^{-\theta x^{\beta }}-360{\left(\theta x^{\beta }\right)}^2{\text{e}}^{-\theta x^{\beta }}\\ \stackrel{}{}-720\theta x^{\beta }{\text{e}}^{-\theta x^{\beta }}+720)]\end{array}$ (13)

Further simplification of (13) gives (7), which completes the proof of Proposition 2.

Proposition 3. Let $X~PH\left(\theta ,\alpha ,\beta \right)$, then the hazard rate function of X is given by

$h_{PH}\left(x\right)=\frac{\beta {\theta }^6\left(\alpha +\frac{\theta }{6}{x}^{6\beta }\right)x^{\beta -1}}{\left[\theta x^{\beta }\left({\theta }^5{x}^{5\beta }+6{\theta }^4{x}^{4\beta }+30{\theta }^3{x}^{3\beta }+120{\theta }^2{x}^{2\beta }+360\theta x^{\beta }+720\right)+6\left(\alpha {\theta }^5+120\right)\right]}$ (14)

Proof. The proof of Proposition 3 follows from using (4) and (7) in the relation

$h_{PH}\left(x\right)=\frac{f_{PH}\left(x\right)}{1-F_{PH}\left(x\right)}$ (15)

Figure 1 and Figure 2 demonstrate the graphs of the pdf and hazard function of the PH distribution for different values of $\theta$, $\alpha$ and $\beta$.

(a) (b)

Figure 1. (a) pdf plot of power hamza distribution; (b) pdf plot of power hamza distribution.

(a) (b)

Figure 2. (a) hazard plot of power hamza distribution; (b) hazard plot of power hamza distribution.

3. Properties of the Power Hamza (PH) Distribution

3.1. Moment Generating Function

Proposition 4. Let $X~PH\left(\theta ,\alpha ,\beta \right)$, then the moment generating function of X is given by

$M_X\left(t\right)={\displaystyle \underset{r=0}{\overset{\infty }{\sum }}\frac{6\alpha {\theta }^5\Gamma \left(r/\beta +1\right)+\Gamma \left(r/\beta +7\right)}{6\left(\alpha {\theta }^5+120\right)r!{\theta }^{r/\beta }}{\left(\frac{t}{{\theta }^{1/\beta }}\right)}^r}$ (16)

Proof. The moment generating function of $X~PH\left(\theta ,\alpha ,\beta \right)$ is obtained as follows

$\begin{array}{c}{M}_X\left(t\right)=E\left({\text{e}}^{tX}\right)={\displaystyle \underset{0}{\overset{\infty }{\int }}{\text{e}}^{tx}{f}_{PH}\left(x\right)\text{d}x}\\ =\frac{\beta {\theta }^6}{\alpha {\theta }^5+120}{\displaystyle \underset{0}{\overset{\infty }{\int }}{\text{e}}^{tx}\left(\alpha +\frac{\theta }{6}{x}^{6\beta }\right)x^{\beta -1}{\text{e}}^{-\theta x^{\beta }}\text{d}x}\\ =\frac{\beta {\theta }^6}{\alpha {\theta }^5+120}{\displaystyle \underset{r=0}{\overset{\infty }{\sum }}\frac{t^r}{r!}}\left[\alpha {\displaystyle \underset{0}{\overset{\infty }{\int }}{x}^{\beta +r-1}{\text{e}}^{-\theta x^{\beta }}\text{d}x}+\frac{\theta }{6}{\displaystyle \underset{0}{\overset{\infty }{\int }}{x}^{7\beta +r-1}{\text{e}}^{-\theta x^{\beta }}\text{d}x}\right]\end{array}$ (17)

Letting $x=y^{1/\beta }$, $y=x^{\beta }$ and $\text{d}x=\left(1/\beta \right)y^{\beta -1}\text{d}y$, (17) reduces to

$\begin{array}{c}{M}_X\left(t\right)=\frac{\beta {\theta }^6}{\alpha {\theta }^5+120}{\displaystyle \underset{r=0}{\overset{\infty }{\sum }}\frac{t^r}{r!}}\left[\frac{\alpha }{\beta }{\displaystyle \underset{0}{\overset{\infty }{\int }}{y}^{r/\beta }{\text{e}}^{-\theta y}\text{d}y}+\frac{\theta }{6\beta }{\displaystyle \underset{0}{\overset{\infty }{\int }}{x}^{r/\beta +6}{\text{e}}^{-\theta y}\text{d}y}\right]\\ =\frac{{\theta }^6}{\alpha {\theta }^5+120}{\displaystyle \underset{r=0}{\overset{\infty }{\sum }}\frac{t^r}{r!}}\left[\frac{\alpha \Gamma \left(r/\beta +1\right)}{{\theta }^{r/\beta +1}}+\frac{\theta \Gamma \left(r/\beta +7\right)}{6{\theta }^{r/\beta +7}}\right]\end{array}$ (18)

Further simplification of (18) gives (16), which completes the proof of Proposition 4.

3.2. Non-Central Moment

Proposition 5. Let $X~PH\left(\theta ,\alpha ,\beta \right)$, then the rth non-moment of X is given by

${{\mu }^{\prime }}_r=\frac{6\alpha {\theta }^5\Gamma \left(r/\beta +1\right)+\Gamma \left(r/\beta +7\right)}{6{\theta }^{r/\beta }\left(\alpha {\theta }^5+120\right)};r=1,2,3,\cdots$ (19)

Proof. The rth moment of $X~PH\left(\theta ,\alpha ,\beta \right)$ is obtained as follows

$\begin{array}{c}{{\mu }^{\prime }}_r=E\left(X^r\right)={\displaystyle \underset{0}{\overset{\infty }{\int }}{x}^r{f}_{PH}\left(x\right)\text{d}x}=\frac{\beta {\theta }^6}{\alpha {\theta }^5+120}{\displaystyle \underset{0}{\overset{\infty }{\int }}{x}^r\left(\alpha +\frac{\theta }{6}{x}^{6\beta }\right)x^{\beta -1}{\text{e}}^{-\theta x^{\beta }}\text{d}x}\\ =\frac{\beta {\theta }^6}{\alpha {\theta }^5+120}\left[\alpha {\displaystyle \underset{0}{\overset{\infty }{\int }}{x}^{\beta +r-1}{\text{e}}^{-\theta x^{\beta }}\text{d}x}+\frac{\theta }{6}{\displaystyle \underset{0}{\overset{\infty }{\int }}{x}^{7\beta +r-1}{\text{e}}^{-\theta x^{\beta }}\text{d}x}\right]\end{array}$ (20)

Putting $x=y^{1/\beta }$, $y=x^{\beta }$ and $\text{d}x=\left(1/\beta \right)y^{\beta -1}\text{d}y$ into (20), yields

$\begin{array}{c}{{\mu }^{\prime }}_r=\frac{{\theta }^6}{\alpha {\theta }^5+120}\left[\alpha {\displaystyle \underset{0}{\overset{\infty }{\int }}{y}^{\frac{r}{\beta }-1}{\text{e}}^{-\theta y}\text{d}y}+\frac{\theta }{6}{\displaystyle \underset{0}{\overset{\infty }{\int }}{y}^{\left(\frac{r}{\beta }+6\right)-1}{\text{e}}^{-\theta y}\text{d}y}\right]\\ =\frac{{\theta }^6}{\alpha {\theta }^5+120}\left[\frac{\alpha \Gamma \left(r/\beta +1\right)}{{\theta }^{r/\beta +1}}+\frac{\theta \Gamma \left(r/\beta +7\right)}{6{\theta }^{r/\beta +7}}\right]\end{array}$ (21)

Simplifying (21) a little bit, we get (19), hence the proof of Proposition 5.

Corollary 2. The first four non-central moments of $X~PH\left(\theta ,\alpha ,\beta \right)$ are

${{\mu }^{\prime }}_1=\frac{6\alpha {\theta }^5\Gamma \left(1/\beta +1\right)+\Gamma \left(1/\beta +7\right)}{6{\theta }^{1/\beta }\left(\alpha {\theta }^5+120\right)}=\mu$ (22)

${{\mu }^{\prime }}_2=\frac{6\alpha {\theta }^5\Gamma \left(2/\beta +1\right)+\Gamma \left(2/\beta +7\right)}{6{\theta }^{2/\beta }\left(\alpha {\theta }^5+120\right)}$ (23)

${{\mu }^{\prime }}_3=\frac{6\alpha {\theta }^5\Gamma \left(3/\beta +1\right)+\Gamma \left(3/\beta +7\right)}{6{\theta }^{3/\beta }\left(\alpha {\theta }^5+120\right)}$ (24)

${{\mu }^{\prime }}_4=\frac{6\alpha {\theta }^5\Gamma \left(4/\beta +1\right)+\Gamma \left(4/\beta +7\right)}{6{\theta }^{4/\beta }\left(\alpha {\theta }^5+120\right)}$ (25)

Proof. The proof of (22)-(25) follows directly from Proposition 5 by substituting $r=1,2,3$ and 4 into (19).

3.3. Variance

Apart from the non-central moments, variance of a distribution is always relevant for measuring the spread from the mean. So, we provide the variance of the PH distribution in Proposition 6.

Proposition 6. Let $X~PH\left(\theta ,\alpha ,\beta \right)$, then the variance of X is given by

${\sigma }^2=\frac{6\left(\alpha {\theta }^5+120\right)\left(6\alpha {\theta }^5\Gamma \left(2/\beta +1\right)+\Gamma \left(2/\beta +7\right)\right)-{\left(6\alpha {\theta }^5\Gamma \left(1/\beta +1\right)+\Gamma \left(1/\beta +7\right)\right)}^2}{36{\theta }^{2/\beta }{\left(\alpha {\theta }^5+120\right)}^2}$ (26)

Proof. The variance of a random variable X can be computed using the relation

${\sigma }^2={{\mu }^{\prime }}_2-{\left({{\mu }^{\prime }}_1\right)}^2$ (27)

The proof of (26) follows directly from substituting (22) and (23) into (27)

${\sigma }^2=\frac{6\alpha {\theta }^5\Gamma \left(2/\beta +1\right)+\Gamma \left(2/\beta +7\right)}{6{\theta }^{2/\beta }\left(\alpha {\theta }^5+120\right)}-{\left(\frac{6\alpha {\theta }^5\Gamma \left(1/\beta +1\right)+\Gamma \left(1/\beta +7\right)}{6{\theta }^{1/\beta }\left(\alpha {\theta }^5+120\right)}\right)}^2$ (28)

3.4. Coefficient of Variation and Index of Dispersion

The coefficient of variation is given by

$\begin{array}{c}cv=\frac{\sigma }{\mu }=\frac{{\left[\frac{6\left(\alpha {\theta }^5+120\right)\left(6\alpha {\theta }^5\Gamma \left(2/\beta +1\right)+\Gamma \left(2/\beta +7\right)\right)-{\left(6\alpha {\theta }^5\Gamma \left(1/\beta +1\right)+\Gamma \left(1/\beta +7\right)\right)}^2}{36{\theta }^{2/\beta }{\left(\alpha {\theta }^5+120\right)}^2}\right]}^{1/2}}{\frac{6\alpha {\theta }^5\Gamma \left(1/\beta +1\right)+\Gamma \left(1/\beta +7\right)}{6{\theta }^{1/\beta }\left(\alpha {\theta }^5+120\right)}}\\ =\frac{{\left[6\left(\alpha {\theta }^5+120\right)\left(6\alpha {\theta }^5\Gamma \left(2/\beta +1\right)+\Gamma \left(2/\beta +7\right)\right)-{\left(6\alpha {\theta }^5\Gamma \left(1/\beta +1\right)+\Gamma \left(1/\beta +7\right)\right)}^2\right]}^{1/2}}{6\alpha {\theta }^5\Gamma \left(1/\beta +1\right)+\Gamma \left(1/\beta +7\right)}\end{array}$ (29)

The index of dispersion is given by

$\begin{array}{c}\gamma =\frac{{\sigma }^2}{\mu }=\frac{\left[\frac{6\left(\alpha {\theta }^5+120\right)\left(6\alpha {\theta }^5\Gamma \left(2/\beta +1\right)+\Gamma \left(2/\beta +7\right)\right)-{\left(6\alpha {\theta }^5\Gamma \left(1/\beta +1\right)+\Gamma \left(1/\beta +7\right)\right)}^2}{36{\theta }^{2/\beta }{\left(\alpha {\theta }^5+120\right)}^2}\right]}{\frac{6\alpha {\theta }^5\Gamma \left(1/\beta +1\right)+\Gamma \left(1/\beta +7\right)}{6{\theta }^{1/\beta }\left(\alpha {\theta }^5+120\right)}}\\ =\frac{6\left(\alpha {\theta }^5+120\right)\left(6\alpha {\theta }^5\Gamma \left(2/\beta +1\right)+\Gamma \left(2/\beta +7\right)\right)-{\left(6\alpha {\theta }^5\Gamma \left(1/\beta +1\right)+\Gamma \left(1/\beta +7\right)\right)}^2}{6{\theta }^{1/\beta }\left(\alpha {\theta }^5+120\right)\left(6\alpha {\theta }^5\Gamma \left(1/\beta +1\right)+\Gamma \left(1/\beta +7\right)\right)}\end{array}$ (30)

3.5. Central Moment

Proposition 7. Let $X~PH\left(\theta ,\alpha ,\beta \right)$, then the central moment of X is given by

${\mu }_r={\displaystyle \underset{i=0}{\overset{r}{\sum }}\left(\begin{array}{c}r\\ i\end{array}\right){\left(-1\right)}^{r-i}{\left(\frac{6\alpha {\theta }^5\Gamma \left(1/\beta +1\right)+\Gamma \left(1/\beta +7\right)}{6{\theta }^{1/\beta }\left(\alpha {\theta }^5+120\right)}\right)}^{r-i}\times \left(\frac{6\alpha {\theta }^5\Gamma \left(i/\beta +1\right)+\Gamma \left(i/\beta +7\right)}{6{\theta }^{i/\beta }\left(\alpha {\theta }^5+120\right)}\right)}$ (31)

Proof. The rth moment of the $X~PH\left(\theta ,\alpha ,\beta \right)$ can be obtained from the relation

${\mu }_r=E{\left(X-\mu \right)}^r={\displaystyle \underset{i=0}{\overset{r}{\sum }}\left(\begin{array}{c}r\\ i\end{array}\right){\left(-\mu \right)}^{r-i}E\left(X^i\right)}$ (32)

Using (22) and (19) in (32) gives (31) and the proof of Proposition 7 is complete.

3.6. Conditional Moment

A function that is useful in deriving the mean residual life function of a component as well as the mean deviations is the conditional moment. Given that X follows a power Hamza distribution with parameters $\theta$, $\alpha$ and $\beta$, then

$E\left(X^n|X>x\right)=\frac{1}{\left(1-F\left(x\right)\right)}{J}_n\left(x\right)$ (33)

where

$\begin{array}{c}{J}_n\left(x\right)={\displaystyle \underset{x}{\overset{\infty }{\int }}{t}^n{f}_{BEI}\left(t\right)\text{d}t}=\frac{\beta {\theta }^6}{\alpha {\theta }^5+120}{\displaystyle \underset{x}{\overset{\infty }{\int }}{t}^n\left(\alpha +\frac{\theta }{6}{t}^{6\beta }\right)t^{\beta -1}{\text{e}}^{-\theta t^{\beta }}\text{d}t}\\ =\frac{\beta {\theta }^6}{\alpha {\theta }^5+120}\left[\alpha {\displaystyle \underset{x}{\overset{\infty }{\int }}{t}^{\beta +n-1}{\text{e}}^{-\theta t^{\beta }}\text{d}t}+\frac{\theta }{6}{\displaystyle \underset{x}{\overset{\infty }{\int }}{t}^{7\beta +n-1}{\text{e}}^{-\theta t^{\beta }}\text{d}t}\right]\end{array}$ (34)

Letting $x=\theta t^{\beta }$ in (34) leads to

$\begin{array}{c}{J}_n\left(x\right)=\frac{\beta {\theta }^6}{\alpha {\theta }^5+120}\left[\frac{\alpha }{\beta {\theta }^{\beta /n+1}}{\displaystyle \underset{\theta x^{\beta }}{\overset{\infty }{\int }}{x}^{\left(\beta /n+1\right)-1}{\text{e}}^{-x}\text{d}x}+\frac{\theta }{6\beta {\theta }^{\beta /n+7}}{\displaystyle \underset{\theta x^{\beta }}{\overset{\infty }{\int }}{x}^{\left(\beta /n+7\right)-1}{\text{e}}^{-x}\text{d}x}\right]\\ =\frac{\beta {\theta }^6}{\alpha {\theta }^5+120}\left[\frac{\alpha }{\beta {\theta }^{\beta /n+1}}\Gamma \left(\frac{\beta }{n}+1,\theta x^{\beta }\right)+\frac{\theta }{6\beta {\theta }^{\beta /n+7}}\Gamma \left(\frac{\beta }{n}+7,\theta x^{\beta }\right)\right]\end{array}$

$J_n\left(x\right)=\frac{\alpha {\theta }^5\Gamma \left(\frac{\beta }{n}+1,\theta x^{\beta }\right)+\Gamma \left(\frac{\beta }{n}+7,\theta x^{\beta }\right)}{6\left(\alpha {\theta }^5+120\right){\theta }^{\beta /n}}$ (35)

where $\Gamma \left(a,z\right)={\displaystyle \underset{x}{\overset{\infty }{\int }}{t}^{a-1}{\text{e}}^{-t}\text{d}t}$ is the complementary incomplete gamma function.

Also from (7), we get

$\left(1-F\left(x\right)\right)=\left[1+\frac{\theta x^{\beta }\left({\theta }^5{x}^{5\beta }+6{\theta }^4{x}^{4\beta }+30{\theta }^3{x}^{3\beta }+120{\theta }^2{x}^{2\beta }+360\theta x^{\beta }+720\right)}{\left(\alpha {\theta }^5+120\right)}\right]{\text{e}}^{-\theta x^{\beta }}$ (36)

Plugging (35) and (36) into (33), we have the conditional moment as

$E\left(X^n|X>x\right)=\frac{\left[\alpha {\theta }^5\Gamma \left(\frac{\beta }{n}+1,\theta x^{\beta }\right)+\Gamma \left(\frac{\beta }{n}+7,\theta x^{\beta }\right)\right]{\text{e}}^{\theta x^{\beta }}}{\left[\theta x^{\beta }\left({\theta }^5{x}^{5\beta }+6{\theta }^4{x}^{4\beta }+30{\theta }^3{x}^{3\beta }+120{\theta }^2{x}^{2\beta }+360\theta x^{\beta }+720\right)+6\left(\alpha {\theta }^5+120\right)\right]{\theta }^{\beta /n}}$ (37)

3.7. Mean Residual Life Function

In many life testing experiments, it is always of interest to know the additional lifetime given that a component has survived until a certain amount of time. To achieve this purpose, the mean residual life function (MRL), which refers to the expected remaining life, $X-x$, given that the item has survived up to time x, is required. It may be observed from (37), that the MRL function is derived from the conditional as follows

$m_X\left(x\right)=E\left(X-x|X>x\right)-x$ (38)

Putting $n=1$ in (37) and substituting the result into (38), we obtain MRL function as

$m_X\left(x\right)=\frac{\left[\alpha {\theta }^5\Gamma \left(\beta +1,\theta x^{\beta }\right)+\Gamma \left(\beta +7,\theta x^{\beta }\right)\right]{\text{e}}^{\theta x^{\beta }}}{\left[\theta x^{\beta }\left({\theta }^5{x}^{5\beta }+6{\theta }^4{x}^{4\beta }+30{\theta }^3{x}^{3\beta }+120{\theta }^2{x}^{2\beta }+360\theta x^{\beta }+720\right)+6\left(\alpha {\theta }^5+120\right)\right]{\theta }^{\beta }}-x$ (39)

3.8. Mean Deviations

In statistical modelling, it is often an interest to measure the spread in a population from either the mean or the median. To achieve this, two indices, namely mean deviation about the mean $\left({\mu }_d\right)$ and mean deviation about the median $\left(M_d\right)$ are used. Let $\mu$ denote the mean and M, the median of a power Hamza distributed random variable X. The values of ${\mu }_d$ and $M_d$ can be calculated using the relationships

${\mu }_d={\displaystyle \underset{0}{\overset{\infty }{\int }}\left|x-\mu \right|f_{PH}\left(x\right)\text{d}x}=2\mu F\left(\mu \right)-2\mu +2J_1\left(\mu \right)$ (40)

and

$M_d={\displaystyle \underset{0}{\overset{\infty }{\int }}\left|x-M\right|f_{PH}\left(x\right)\text{d}x}=2J_1\left(M\right)-\mu$ (41)

respectively. By replacing x with $\mu$ and M in (7) and (35), yields the following

$\begin{array}{l}{\mu }_d=2\mu \left[1-\left(1+\frac{\theta {\mu }^{\beta }\left({\theta }^5{\mu }^{5\beta }+6{\theta }^4{\mu }^{4\beta }+30{\theta }^3{\mu }^{3\beta }+120{\theta }^2{\mu }^{2\beta }+360\theta {\mu }^{\beta }+720\right)}{\left(\alpha {\theta }^5+120\right)}\right){\text{e}}^{-\theta {\mu }^{\beta }}\right]\\ -\frac{2\left[6\mu \left(\alpha {\theta }^5+120\right){\theta }^{\beta }-\left(\alpha {\theta }^5\Gamma \left(\beta +1,\theta {\mu }^{\beta }\right)+\Gamma \left(\beta +7,\theta {\mu }^{\beta }\right)\right)\right]}{6\left(\alpha {\theta }^5+120\right){\theta }^{\beta }}\end{array}$ (42)

$M_d=\frac{2\left[\alpha {\theta }^5\Gamma \left(\beta +1,\theta M^{\beta }\right)+\Gamma \left(\beta +7,\theta M^{\beta }\right)-3\mu \left(\alpha {\theta }^5+120\right){\theta }^{\beta }\right]}{3\left(\alpha {\theta }^5+120\right){\theta }^{\beta }}$ (43)

3.9. Bonferroni and Lorenz Curves

It has been found that the Bonferroni curve proposed by [16] and Lorenz curve proposed by [17] have applications in the fields of economics, reliability, demography, insurance, medicine, among others. So, if $X~PH\left(\theta ,\alpha ,\beta \right)$, the Bonferroni and Lorenz curves are respectively given by

$B\left(p\right)=\frac{1}{p\mu }{\displaystyle \underset{0}{\overset{q}{\int }}x{f}_{PH}\left(x\right)\text{d}x}$ (44)

and

$L\left(p\right)=\frac{1}{\mu }{\displaystyle \underset{0}{\overset{q}{\int }}x{f}_{PH}\left(x\right)\text{d}x}$ (45)

where $\mu =E\left(X\right)$ and $q=F^{-1}\left(p\right)$, $0<p<1$. Hence, for the power Hamza pdf (3), one gets

$\begin{array}{c}{\displaystyle \underset{0}{\overset{q}{\int }}x{f}_{PH}\left(x\right)\text{d}x}=\frac{\beta {\theta }^6}{\alpha {\theta }^5+120}{\displaystyle \underset{0}{\overset{q}{\int }}x\left(\alpha +\frac{\theta }{6}{x}^{6\beta }\right)x^{\beta -1}{\text{e}}^{-\theta x^{\beta }}\text{d}x}\\ =\frac{\beta {\theta }^6}{\alpha {\theta }^5+120}\left[\alpha {\displaystyle \underset{0}{\overset{q}{\int }}{x}^{\beta }{\text{e}}^{-\theta x^{\beta }}\text{d}x}+\frac{\theta }{6}{\displaystyle \underset{0}{\overset{q}{\int }}{x}^{7\beta }{\text{e}}^{-\theta x^{\beta }}\text{d}x}\right]\end{array}$ (46)

Letting $w=\theta x^{\beta }$, $x={\left(w/\theta \right)}^{1/\beta }$, $\text{d}x=\left(1/\beta \right){\left(1/\theta \right)}^{1/\beta }{w}^{1/\beta -1}\text{d}w$ and noting that $x\in \left(0,q\right)$ implies $w\in \left(0,\theta q^{\beta }\right)$, (1) becomes

$\begin{array}{l}{\displaystyle \underset{0}{\overset{q}{\int }}x{f}_{PH}\left(x\right)\text{d}x}\\ =\frac{{\theta }^6}{\alpha {\theta }^5+120}\left[\alpha {\left(\frac{1}{\theta }\right)}^{\left(\beta +\frac{1}{\beta }\right)}{\displaystyle \underset{0}{\overset{\theta q^{\beta }}{\int }}{w}^{\left(\beta +\frac{1}{\beta }\right)-1}{\text{e}}^{-w}\text{d}w}+\frac{\theta }{6}{\left(\frac{1}{\theta }\right)}^{\left(7\beta +\frac{1}{\beta }\right)}{\displaystyle \underset{0}{\overset{\theta q^{\beta }}{\int }}{w}^{\left(7\beta +\frac{1}{\beta }\right)-1}{\text{e}}^{-w}\text{d}w}\right]\\ =\frac{{\theta }^6}{\alpha {\theta }^5+120}\left[\alpha {\left(\frac{1}{\theta }\right)}^{\left(\beta +\frac{1}{\beta }\right)}\Gamma \left(\beta +\frac{1}{\beta },\theta q^{\beta }\right)+\frac{\theta }{6}{\left(\frac{1}{\theta }\right)}^{\left(7\beta +\frac{1}{\beta }\right)}\Gamma \left(7\beta +\frac{1}{\beta },\theta q^{\beta }\right)\right]\\ =\frac{{\theta }^6{\left(\frac{1}{\theta }\right)}^{\frac{{\beta }^2+1}{\beta }}\left[6\alpha {\theta }^{6\beta }\Gamma \left(\frac{{\beta }^2+1}{\beta },\theta q^{\beta }\right)+\theta \Gamma \left(\frac{7{\beta }^2+1}{\beta },\theta q^{\beta }\right)\right]}{6{\theta }^{6\beta }\left(\alpha {\theta }^5+120\right)}\end{array}$ (47)

Substituting (46) and (47) into (44) and (45), we obtain the Bonferroni and Lorenz curves, respectively, for the power Hamza distribution as

$B\left(p\right)=\frac{{\theta }^6{\left(\frac{1}{\theta }\right)}^{\frac{{\beta }^2+1}{\beta }}\left[6\alpha {\theta }^{6\beta }\Gamma \left(\frac{{\beta }^2+1}{\beta },\theta q^{\beta }\right)+\theta \Gamma \left(\frac{7{\beta }^2+1}{\beta },\theta q^{\beta }\right)\right]}{6p{\theta }^{6\beta }\mu \left(\alpha {\theta }^5+120\right)}$ (48)

$L\left(p\right)=\frac{{\theta }^6{\left(\frac{1}{\theta }\right)}^{\frac{{\beta }^2+1}{\beta }}\left[6\alpha {\theta }^{6\beta }\Gamma \left(\frac{{\beta }^2+1}{\beta },\theta q^{\beta }\right)+\theta \Gamma \left(\frac{7{\beta }^2+1}{\beta },\theta q^{\beta }\right)\right]}{6{\theta }^{6\beta }M\left(\alpha {\theta }^5+120\right)}$ (49)

3.10. Stochastic Ordering

In this section, we discuss the comparative behaviour of the power Hamza random variable using the stochastic ordering. In line with [18], a power Hamza random variable X is said to be smaller than another random variable Y in the 1) stochastic order $\left(X\le {}_{st}Y\right)$ if $F_X\left(x\right)\ge F_Y\left(x\right)\forall x$, 2) hazard rate order $\left(X\le {}_{hr}Y\right)$ if $h_X\left(x\right)\ge h_Y\left(x\right)\forall x$, 3) mean residual life order $\left(X\le {}_{mrl}Y\right)$ if $m_X\left(x\right)\ge m_Y\left(x\right)\forall x$ and 4) likelihood ratio order $\left(X\le {}_{lr}Y\right)$ if $f_X\left(x\right)/f_Y\left(x\right)$ a decreasing function of x. To show the flexibility of the power Hamza distribution, we present the following Proposition.

Proposition 8. Let $X~\left({\theta }_1,{\alpha }_1,{\beta }_1\right)$ and $Y~\left({\theta }_2,{\alpha }_2,{\beta }_2\right)$ be two independent random variables. If 1) ${\theta }_1>{\theta }_2$, ${\alpha }_1={\alpha }_2$ and ${\beta }_1={\beta }_2$ ; 2) ${\alpha }_1>{\alpha }_2$, ${\theta }_1={\theta }_2$, ${\beta }_1={\beta }_2$ ; 3) ${\beta }_1<{\beta }_2$, ${\theta }_1={\theta }_2$, ${\alpha }_1={\alpha }_2$ and 4) ${\theta }_1>{\theta }_2$, ${\alpha }_2<{\alpha }_1$, ${\beta }_1<{\beta }_2$, then $X{\le }_{lr}Y$, $X{\le }_{hr}Y$, $X{\le }_{mlr}Y$ and $X{\le }_{st}Y$.

Proof. The likelihood ratio is

$\frac{f_X\left(x;{\theta }_1,{\alpha }_1,{\beta }_1\right)}{f_Y\left(x;{\theta }_2,{\alpha }_2,{\beta }_2\right)}=\frac{{\beta }_1{\theta }_1^6\left({\alpha }_2{\theta }_2^5+120\right)}{{\beta }_2{\theta }_2^6\left({\alpha }_1{\theta }_1^5+120\right)}{x}^{{\beta }_1-{\beta }_2}\left(\frac{6{\alpha }_1+{\theta }_1{x}^{6{\beta }_1}}{6{\alpha }_2+{\theta }_2{x}^{6{\beta }_2}}\right){\text{e}}^{-\left({\theta }_1{x}^{{\beta }_1}-{\theta }_2{x}^{{\beta }_2}\right)}$ (50)

The log-likelihood ratio is

$\begin{array}{l}\ln \left\{\frac{f_X\left(x;{\theta }_1,{\alpha }_1,{\beta }_1\right)}{f_Y\left(x;{\theta }_2,{\alpha }_2,{\beta }_2\right)}\right\}=\ln \left[\frac{{\beta }_1{\theta }_1^6\left({\alpha }_2{\theta }_2^5+120\right)}{{\beta }_2{\theta }_2^6\left({\alpha }_1{\theta }_1^5+120\right)}\right]+\left({\beta }_1-{\beta }_2\right)\ln x\\ +\ln \left(\frac{6{\alpha }_1+{\theta }_1{x}^{6{\beta }_1}}{6{\alpha }_2+{\theta }_2{x}^{6{\beta }_2}}\right)-\left({\theta }_1{x}^{{\beta }_1}-{\theta }_2{x}^{{\beta }_2}\right)\end{array}$ (51)

Differentiating the log-likelihood with respect to x, we get

$\begin{array}{l}\frac{\text{d}}{\text{d}x}\ln \left\{\frac{f_X\left(x;{\theta }_1,{\alpha }_1,{\beta }_1\right)}{f_Y\left(x;{\theta }_2,{\alpha }_2,{\beta }_2\right)}\right\}\\ =\frac{36\left(6{\theta }_1{\alpha }_2{\beta }_1{x}^{6{\beta }_1-1}+{\theta }_2{\alpha }_1{\beta }_2{x}^{6{\beta }_2-1}\right)+6{\theta }_1{\theta }_2\left({\beta }_1-{\beta }_2\right)x^{6\left({\beta }_1-{\beta }_2\right)-1}}{\left(6{\alpha }_1+{\theta }_1{x}^{6{\beta }_1}\right)\left(6{\alpha }_2+{\theta }_2{x}^{6{\beta }_2}\right)}\\ +\frac{\left({\beta }_1-{\beta }_2\right)}{x}-\left({\theta }_1{\beta }_1{x}^{{\beta }_1-1}-{\theta }_2{\beta }_2{x}^{{\beta }_2-1}\right)\end{array}$ (52)

Since $\frac{\text{d}}{\text{d}x}\ln \left\{\frac{f_X\left(x;{\theta }_1,{\alpha }_1,{\beta }_1\right)}{f_X\left(x;{\theta }_1,{\alpha }_1,{\beta }_1\right)}\right\}<0$ for conditions 1), 2), 3) and 4), then $X{\le }_{lr}Y$ and hence, $X{\le }_{hr}Y$, $X{\le }_{mlr}Y$ and $X{\le }_{st}Y$, which completes the proof of Proposition 8.

3.11. Rényi Entropy

To quantify the amount of information (such as the diversity, uncertainty, or randomness) contained in a random sample drawn from a population, the entropy is utilized. It may be noted that a large value of entropy indicates that the data contains greater uncertainty. Several studies have applied entropy in the fields of physics, probability and statistics, communication theory, economics, among others. Therefore, we derive one of the commonly used entropies, namely the R nyi entropy. For $\rho >0$ and $\rho \ne 1$, the Rényi entropy due to [19] is defined for a continuous random variable as

$E_{\rho }=\frac{1}{1-\rho }\log \left({\displaystyle \underset{0}{\overset{\infty }{\int }}{f}_{PH}^{{}^{\rho }}\left(x\right)\text{d}x}\right)$ (53)

Using the pdf (3) in (53), we obtain

$\begin{array}{c}{E}_{\rho }=\frac{1}{1-\rho }\log \left[{\displaystyle \underset{0}{\overset{\infty }{\int }}{\left(\frac{\beta {\theta }^6}{\alpha {\theta }^5+120}\left(\alpha +\frac{\theta }{6}{x}_i^{6\beta }\right)x_i^{\beta -1}{\text{e}}^{-\theta x_i^{\beta }}\right)}^{\rho }\text{d}x}\right]\\ =\frac{1}{1-\rho }\log \left[{\left(\frac{\beta {\theta }^6}{\alpha {\theta }^5+120}\right)}^{\rho }{\displaystyle \underset{0}{\overset{\infty }{\int }}{\left(\alpha +\frac{\theta }{6}{x}^{6\beta }\right)}^{\rho }{x}^{\rho \left(\beta -1\right)}{\text{e}}^{-\theta \rho x^{\beta }}\text{d}x}\right]\\ =\frac{1}{1-\rho }\log \left[{\displaystyle \underset{j=1}{\overset{\rho }{\sum }}\left(\begin{array}{c}\rho \\ j\end{array}\right)\frac{{\theta }^{6\rho +j}{\alpha }^{\rho -j}{\beta }^{\rho }}{6^j{\left(\alpha {\theta }^5+120\right)}^{\rho }}}{\displaystyle \underset{0}{\overset{\infty }{\int }}{x}^{6\beta +\rho \left(\beta -1\right)}{\text{e}}^{-\theta \rho x^{\beta }}\text{d}x}\right]\end{array}$

$=\frac{1}{1-\rho }\log \left[{\displaystyle \underset{j=1}{\overset{\rho }{\sum }}\left(\begin{array}{c}\rho \\ j\end{array}\right)\frac{6^{-j}{\theta }^{6\rho +j}{\alpha }^{\rho -j}{\beta }^{\rho -1}}{{\left(\alpha {\theta }^5+120\right)}^{\rho }}}{\displaystyle \underset{0}{\overset{\infty }{\int }}{y}^{\left(6j+\rho -1/\rho +1/\beta \right)-1}{\text{e}}^{-\theta \rho y}\text{d}y}\right]$ $\because y=x^{\beta }$

$E_{\rho }=\frac{1}{1-\rho }\log \left[{\displaystyle \underset{j=1}{\overset{\rho }{\sum }}\left(\begin{array}{c}\rho \\ j\end{array}\right)\frac{6^{-j}{\theta }^{6\rho +j}{\alpha }^{\rho -j}{\beta }^{\rho -1}}{{\left(\alpha {\theta }^5+120\right)}^{\rho }}}\frac{\Gamma \left(6j+\rho -1/\rho +1/\beta \right)}{{\left(\theta \rho \right)}^{6j+\rho -1/\rho +1/\beta }}\right]$ (54)

3.12. Distribution of Order Statistics for the PH Distribution

Suppose $X_{\left(1\right)}<X_{\left(2\right)}<\cdots <X_{\left(n\right)}$ constitutes the order statistics for a random sample $X_1,X_2,\cdots ,X_n$ drawn from the power Hamza distribution with pdf (3) and cdf (7). Then the pdf of the rth order statistic $X_{\left(r\right)}$ can be written as

$\begin{array}{l}{f}_{X_{\left(r\right)}}\left(x\right)=\frac{n!}{\left(r-1\right)!\left(n-r\right)!}{f}_{PH}\left(x\right){\left(F_{PH}\left(x\right)\right)}^{r-1}{\left(1-F_{PH}\left(x\right)\right)}^{n-r};r=1,2,\cdots ,n\\ =\frac{\beta {\theta }^{6}n!}{\left(\alpha {\theta }^5+120\right)\left(r-1\right)!\left(n-r\right)!}{\displaystyle \underset{u=0}{\overset{n-r}{\sum }}{\left(-1\right)}^u\left(\begin{array}{c}n-r\\ u\end{array}\right)x^{\beta -1}{\text{e}}^{-\theta x^{\beta }}}\left(\alpha +\frac{\theta }{6}{x}^{6\beta }\right)\\ \times {\left[1-\left(1+\frac{\theta x^{\beta }\left({\theta }^5{x}^{5\beta }+6{\theta }^4{x}^{4\beta }+30{\theta }^3{x}^{3\beta }+120{\theta }^2{x}^{2\beta }+360\theta x^{\beta }+720\right)}{\left(\alpha {\theta }^5+120\right)}\right){\text{e}}^{-\theta x^{\beta }}\right]}^{r+u-1}\end{array}$ (55)

Putting $r=1$ into (55) gives the pdf of the first order statistic $X_{\left(1\right)}$ as

$\begin{array}{l}{f}_{X_{\left(1\right)}}\left(x\right)=\frac{n\beta {\theta }^6}{\alpha {\theta }^5+120}{x}^{\beta -1}{\text{e}}^{-\theta x^{\beta }}\left(\alpha +\frac{\theta }{6}{x}^{6\beta }\right)\\ \times {\left[\left(1+\frac{\theta x^{\beta }\left({\theta }^5{x}^{5\beta }+6{\theta }^4{x}^{4\beta }+30{\theta }^3{x}^{3\beta }+120{\theta }^2{x}^{2\beta }+360\theta x^{\beta }+720\right)}{\left(\alpha {\theta }^5+120\right)}\right){\text{e}}^{-\theta x^{\beta }}\right]}^{n-1}\end{array}$ (56)

Putting $r=n$ into (55) gives the pdf of the nth order statistic $X_{\left(n\right)}$ as

$\begin{array}{l}{f}_{X_{\left(n\right)}}\left(x\right)=\frac{n\beta {\theta }^6}{\alpha {\theta }^5+120}{x}^{\beta -1}{\text{e}}^{-\theta x^{\beta }}\left(\alpha +\frac{\theta }{6}{x}^{6\beta }\right)\\ \times {\left(1-\left[1+\frac{\theta x^{\beta }\left({\theta }^5{x}^{5\beta }+6{\theta }^4{x}^{4\beta }+30{\theta }^3{x}^{3\beta }+120{\theta }^2{x}^{2\beta }+360\theta x^{\beta }+720\right)}{\left(\alpha {\theta }^5+120\right)}\right]{\text{e}}^{-\theta x^{\beta }}\right)}^{n-1}\end{array}$ (57)

4. Maximum Likelihood Estimators of the Power Hamza Distribution

Let $X_1,X_2,\cdots ,X_n$ denote a random sample of size n from the PH distribution having parameters $\theta$, $\alpha$ and $\beta$. To estimate the parameters $\theta$, $\alpha$ and $\beta$ using the maximum likelihood method, we define the likelihood function of the random sample form the PH distribution as

$L\left(\theta ,\alpha ,\beta \right)={\displaystyle \underset{i=1}{\overset{n}{\prod }}f\left(x_i;\theta ,\alpha ,\beta \right)}={\left(\frac{\beta {\theta }^6}{\alpha {\theta }^5+120}\right)}^n{\text{e}}^{-\theta {\displaystyle \underset{i=1}{\overset{n}{\sum }}{x}_i^{\beta }}}{\displaystyle \underset{i=1}{\overset{n}{\prod }}\left(\alpha +\frac{\theta }{6}{x}_i^{6\beta }\right)x_i^{\beta -1}}$ (58)

Taking the natural log of (58), we obtain the log-likelihood function as

$\begin{array}{l}\ln L\left(\theta ,\alpha ,\beta \right)=n\left[\ln \left(\beta \right)+6\ln \left(\theta \right)-\ln \left(\alpha {\theta }^5+120\right)\right]+{\displaystyle \underset{i=1}{\overset{n}{\sum }}\ln \left(\alpha +\frac{\theta }{6}{x}_i^{6\beta }\right)}\\ +\left(\beta -1\right){\displaystyle \underset{i=1}{\overset{n}{\sum }}{x}_i}-\theta {\displaystyle \underset{i=1}{\overset{n}{\sum }}{x}_i^{\beta }}\end{array}$ (59)

Differentiating (59) with respect to $\theta$, $\alpha$ and $\beta$ respectively and equating the resulting derivatives to zero, one obtains

$\frac{\partial \ln L\left(\theta ,\alpha ,\beta \right)}{\partial \theta }=\frac{6n}{\theta }-\frac{5n\alpha {\theta }^4}{\alpha {\theta }^5+120}+{\displaystyle \underset{i=1}{\overset{n}{\sum }}\frac{x_i^{6\beta }}{6\alpha +\theta x_i^{6\beta }}}-{\displaystyle \underset{i=1}{\overset{n}{\sum }}{x}_i^{\beta }}=0$ (60)

$\frac{\partial \ln L\left(\theta ,\alpha ,\beta \right)}{\partial \alpha }=\frac{-n{\theta }^5}{\alpha {\theta }^5+120}+{\displaystyle \underset{i=1}{\overset{n}{\sum }}\frac{6}{6\alpha +\theta x_i^{6\beta }}}=0$ (61)

$\frac{\partial \ln L\left(\theta ,\alpha ,\beta \right)}{\partial \beta }=\frac{n}{\beta }+6\theta {\displaystyle \underset{i=1}{\overset{n}{\sum }}\frac{x_i^{6\beta }\ln \left(x_i\right)}{6\alpha +\theta x_i^{6\beta }}}+{\displaystyle \underset{i=1}{\overset{n}{\sum }}\ln \left(x_i\right)}-\theta {\displaystyle \underset{i=1}{\overset{n}{\sum }}{x}_i^{\beta }\ln \left(x_i\right)}=0$ (62)

The above non-linear systems of equations are solved by numerical iteration technique and maximum likelihood estimates are obtained. Since the maximum likelihood estimates for $\theta$, $\alpha$ and $\beta$ are not in closed form we use the large sample behaviour of maximum likelihood estimators to obtain the confidence intervals for model parameters.

5. Asymptotic Confidence Intervals of the Power Hamza Distribution

In this section, we present the asymptotic confidence intervals for the parameters of the PH distribution. Let $\hat{\psi }={\left(\hat{\theta },\hat{\alpha },\hat{\beta }\right)}^{\text{T}}$ be the maximum likelihood estimate of $\psi ={\left(\theta ,\alpha ,\beta \right)}^{\text{T}}$. Under the conditions that the parameters are in the interior of the parameter space, but not on the boundary, the asymptotic distribution of $\sqrt{n}\left(\hat{\psi }-\psi \right)$ is $N_3\left(0,I^{-1}\left(\psi \right)\right)$, where $I\left(\psi \right)$ is the expected Fisher information matrix. The asymptotic behaviour of the expected information matrix can be approximated by the observed information matrix, denoted by $I_n\left(\hat{\psi }\right)$. The observed information matrix of the power Hamza distribution is given by

$I_n\left(\hat{\psi }\right)=\left(\begin{array}{ccc}\frac{{\partial }^2\ln L\left(\theta ,\alpha ,\beta \right)}{\partial {\theta }^2}& \frac{{\partial }^2\ln L\left(\theta ,\alpha ,\beta \right)}{\partial \hat{\theta }\partial \hat{\alpha }}& \frac{{\partial }^2\ln L\left(\theta ,\alpha ,\beta \right)}{\partial \hat{\theta }\partial \hat{\beta }}\\ \frac{{\partial }^2\ln L\left(\theta ,\alpha ,\beta \right)}{\partial \hat{\alpha }\partial \hat{\theta }}& \frac{{\partial }^2\ln L\left(\theta ,\alpha ,\beta \right)}{\partial {\hat{\alpha }}^2}& \frac{{\partial }^2\ln L\left(\theta ,\alpha ,\beta \right)}{\partial \hat{\alpha }\partial \hat{\beta }}\\ \frac{{\partial }^2\ln L\left(\theta ,\alpha ,\beta \right)}{\partial \hat{\beta }\partial \hat{\theta }}& \frac{{\partial }^2\ln L\left(\theta ,\alpha ,\beta \right)}{\partial \hat{\beta }\partial \hat{\alpha }}& \frac{{\partial }^2\ln L\left(\theta ,\alpha ,\beta \right)}{\partial {\hat{\beta }}^2}\end{array}\right)$ (63)

Thus,

$I_n^{-1}\left(\hat{\psi }\right)={\left(nI\left(\hat{\psi }\right)\right)}^{-1}=\left(\begin{array}{ccc}\mathrm{var}\left(\hat{\theta }\right)& \mathrm{cov}\left(\hat{\theta },\hat{\alpha }\right)& \mathrm{cov}\left(\hat{\theta },\hat{\beta }\right)\\ \mathrm{cov}\left(\hat{\alpha },\hat{\theta }\right)& \mathrm{var}\left(\hat{\alpha }\right)& \mathrm{cov}\left(\hat{\alpha },\hat{\beta }\right)\\ \mathrm{cov}\left(\hat{\beta },\hat{\theta }\right)& \mathrm{cov}\left(\hat{\beta },\hat{\alpha }\right)& \mathrm{var}\left(\hat{\beta }\right)\end{array}\right)$ (64)

Taking the second order derivatives of (59) with respect to $\theta$, $\alpha$ and $\beta$ are, respectively, we obtain the entries of (63) as follows

$\frac{{\partial }^2\ln L\left(\theta ,\alpha ,\beta \right)}{\partial {\theta }^2}=-\frac{6n}{\theta }+\frac{5n\alpha {\theta }^3\left(\alpha {\theta }^5-480\right)}{{\left(\alpha {\theta }^5+120\right)}^2}-{\displaystyle \underset{i=1}{\overset{n}{\sum }}\frac{x_i^{12\beta }}{{\left(6\alpha +\theta x_i^{6\beta }\right)}^2}}$ (65)

$\frac{{\partial }^2\ln L\left(\theta ,\alpha ,\beta \right)}{\partial {\alpha }^2}=\frac{n{\theta }^{10}}{{\left(\alpha {\theta }^5+120\right)}^2}-{\displaystyle \underset{i=1}{\overset{n}{\sum }}\frac{36}{{\left(6\alpha +\theta x_i^{6\beta }\right)}^2}}$ (66)

$\frac{{\partial }^2\ln L\left(\theta ,\alpha ,\beta \right)}{\partial {\beta }^2}=-\frac{n}{{\beta }^2}+216\alpha \theta {\displaystyle \underset{i=1}{\overset{n}{\sum }}\frac{x_i^{6\beta }{\ln }^2\left(x_i\right)}{{\left(6\alpha +\theta x_i^{6\beta }\right)}^2}}-\theta {\displaystyle \underset{i=1}{\overset{n}{\sum }}{x}_i^{\beta }{\ln }^2\left(x_i\right)}$ (67)

$\frac{{\partial }^2\ln L\left(\theta ,\alpha ,\beta \right)}{\partial \alpha \partial \theta }=\frac{-600n{\theta }^4}{{\left(\alpha {\theta }^5+120\right)}^2}-6{\displaystyle \underset{i=1}{\overset{n}{\sum }}\frac{x_i^{6\beta }}{{\left(6\alpha +\theta x_i^{6\beta }\right)}^2}}$ (68)

$\frac{{\partial }^2\ln L\left(\theta ,\alpha ,\beta \right)}{\partial \beta \partial \theta }=36\alpha {\displaystyle \underset{i=1}{\overset{n}{\sum }}\frac{x_i^{6\beta }\ln \left(x_i\right)}{{\left(6\alpha +\theta x_i^{6\beta }\right)}^2}}-{\displaystyle \underset{i=1}{\overset{n}{\sum }}{x}_i^{\beta }\ln \left(x_i\right)}$ (69)

$\frac{{\partial }^2\ln L\left(\theta ,\alpha ,\beta \right)}{\partial \alpha \partial \beta }=-36{\displaystyle \underset{i=1}{\overset{n}{\sum }}\frac{x_i^{6\beta }\ln \left(x_i\right)}{{\left(6\alpha +\theta x_i^{6\beta }\right)}^2}}$ (70)

The expectations in the Fisher information matrix can be obtained numerically. The multivariate normal distribution with mean vector ${\left(0,0,0\right)}^{\text{T}}$ and covariance matrix $I^{-1}\left(\psi \right)$ can be used to construct confidence intervals for the model parameters. The approximate $100\left(1-\eta \right)\%$ two-sided confidence intervals for $\theta$, $\alpha$ and $\beta$ are determined by

$\hat{\theta }\pm Z_{\eta /2}\sqrt{\mathrm{var}\left(\hat{\theta }\right)}$, $\hat{\alpha }\pm Z_{\eta /2}\sqrt{\mathrm{var}\left(\hat{\alpha }\right)}$ and $\hat{\beta }\pm Z_{\eta /2}\sqrt{\mathrm{var}\left(\hat{\beta }\right)}$ (71)

respectively, where $Z_{\eta /2}$ is the upper $\left(\eta /2\right)\text{th}$ percentile of a standard normal distribution.

6. Applications

In this section, we provide an application to real data set to demonstrate the importance and flexibility of the PH distribution.

The data set is on the breaking strength of carbon fibres of 50 mm length (GPa). The data has been previously used by [20] and [21]. The data is as follows:

0.39, 0.85, 1.08, 1.25, 1.47, 1.57, 1.61, 1.61, 1.69, 1.80, 1.84, 1.87, 1.89, 2.03, 2.03, 2.05, 2.12, 2.35, 2.41, 2.43, 2.48, 2.50, 2.53, 2.55, 2.55, 2.56, 2.59, 2.67, 2.73, 2.74, 2.79, 2.81, 2.82, 2.85, 2.87, 2.88, 2.93, 2.95, 2.96, 2.97, 3.09, 3.11, 3.11, 3.15, 3.15, 3.19, 3.22, 3.22, 3.27, 3.28, 3.31, 3.31, 3.33, 3.39, 3.39, 3.56, 3.60, 3.65, 3.68, 3.70, 3.75, 4.20, 4.38, 4.42, 4.70, 4.90.

We fitted the PH distribution to the data set by using the method of maximum likelihood and the results are compared with five other competitive lifetime distributions namely,

1) Hamza distribution (HD) defined in Equation (2),

2) Weighted Weibull distribution (WWD):

$f\left(x;\theta ,\beta \right)=\frac{\theta {\beta }^{c+1}{x}^{\beta \left(c+1\right)-1}{\text{e}}^{-\theta x^{\beta }}}{\Gamma \left(c+1\right)};x>0,\theta >0,\beta >0,c>-1$,

3) Two-Parameter Weibull distribution (TPWD):

$f\left(x;\theta ,\beta \right)=\theta \beta x^{\beta -1}{\text{e}}^{-\theta x^{\beta }};x>0,\theta >0,\beta >0$,

4) Pareto distribution (PD):

$f\left(x;\alpha \right)=\frac{\alpha x_m^{\alpha }}{x_m^{\alpha +1}};x>0,\alpha >0$,

5) Exponential distribution (ED):

$f\left(x;\theta \right)=\theta {\text{e}}^{-\theta x};x>0,\theta >0$.

We used the goodness-of-fit test based on the Kolmogorov-Smirnov test due to ( [22] [23] [24] [25] ) with its corresponding p-value to verify that the data set under consideration actually follow the proposed distribution. The computational formula for this goodness-of-fit test is given by

$\text{KS}=\max \left\{\frac{i}{n}-\hat{F}\left(x_{\left(i\right)}\right),\hat{F}\left(x_{\left(i\right)}\right)-\frac{i-1}{n}\right\}$ (72)

where $\hat{F}\left(x_i\right)$ is the estimated distribution function under the ordered data. Since there is more than one distribution to be compared, the distribution with the largest KS p-value will be more appropriate to fit the given sample.

We shall also determine the appropriate model from among all models compared for the real data set by considering three discrimination criteria, based on the log-likelihood function evaluated at the maximum likelihood estimates, the Akaike information criterion (AIC) due to [26] and the Bayesian information criterion (BIC) due to [27], respectively. To compute the AIC and BIC, the following formulae are used

Table 1. Parameter estimates, standard errors, log-likelihood values and goodness-of-fit measures.

$\text{AIC}=-2l+2k$ (73)

$\text{BIC}=-2l+k\log \left(n\right)$ (74)

where l denotes the log-likelihood function evaluated at the maximum likelihood estimates, k is the number of parameters in the statistical model and n is the sample size of the fitted data respectively. All the computations for (72)-(74) were performed using R software. Generally, for the given data-sets, we consider a distribution to be best among all competing distributions if it has smallest AIC value, the smallest BIC value, the smallest log-likelihood value and the largest p-value.

As shown in Table 1, the PH distribution has the largest KS p-value and smallest AIC, BIC and log-likelihood values as compared to other fitted distributions. Hence, the PH distribution is better than the other distributions in Table 1 for fitting the data under consideration.

7. Conclusion

This study introduced a new distribution, called the power Hamza distribution using power transformation method. The contribution of this paper has to do with addition of skewness to the Hamza distribution, which depends only on one parameter. The density function of the power Hamza distribution can take various forms depending on its shape parameter. The hazard rate function of the power Hamza distribution exhibits heavy-tailed shape, upside-down bathtub shape and J-shape, which implies that the distribution can be used for analyzing lifetime and survival time datasets. A detailed discussion of the properties of the proposed distribution has been given. Estimates of the three unknown parameters of the PH distribution are obtained using the maximum likelihood estimation method. The PH distribution was fitted to a real dataset and compared to five distributions and the results showed that the proposed distribution outperformed all of them in modelling the data under consideration.

Acknowledgements

The authors are thankful to the editor and reviewers of this article for providing very useful comments which led to the improvement of this paper.

Authors’ Contribution

All authors contributed equally in developing the article.

Funding Statement

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1]	Aijaz, A., Jallal, M., Ain Ul, S.Q. and Tripathi, R. (2020) The Hamza Distribution with Statistical Properties and Applications. Asian Journal of Probability and Statistics, 8, 28-42. https://doi.org/10.9734/ajpas/2020/v8i130198
[2]	Lindley, D.V. (1958) Fiducial Distributions and Bayes. Journal of the Royal Statistical Society, Series B, 20, 102-107. https://doi.org/10.1111/j.2517-6161.1958.tb00278.x
[3]	Shanker, R. and Shukla, K.K. (2017) Ishita Distribution and Its Applications. Biometrics & Biostatistics International Journal, 5, 1-9. https://doi.org/10.15406/bbij.2017.05.00126
[4]	Shukla, K.K. (2018) Pranav Distribution with Properties and Its Applications. Biometrics and Biostatistics International Journal, 7, 244-254. https://doi.org/10.15406/bbij.2018.07.00215
[5]	Ghitany, M.E., Al-Mutairi, D.K., Balakrishnan, N. and Al-Enezi, L.J. (2013) Power Lindley Distribution and Associated Inference. Computational Statistics and Data Analysis, 64, 20-33. https://doi.org/10.1016/j.csda.2013.02.026
[6]	Shanker, R. and Shukla, K.K. (2017) Power Akash Distribution and Its Applications. Journal of Applied Quantitative Methods, 12, 1-1012.
[7]	Shanker, R. and Shukla, K.K. (2017) Power Shanker Distribution and Its Application. Turkiye Klinikleri Journal of Biostatistics, 9, 175-187. https://doi.org/10.5336/biostatic.2017-56601
[8]	Shukla, K.K. and Shanker, R. (2018) Power Ishita Distribution and Its Applications to Model Lifetime Data. Statistics in Transition New Series, 19, 135-148. https://doi.org/10.21307/stattrans-2018-008
[9]	Shanker, R. and Shukla, K.K. (2018) A Two-Parameter Power Aradhana Distribution with Properties and Application. Indian Journal of Industrial and Applied Mathematics, 9, 210-220. https://doi.org/10.5958/1945-919X.2018.00016.6
[10]	Abebe, B., Tesfay, M., Eyob, T. and Shanker, R. (2019) A Two-Parameter Power Rama Distribution with Properties and Applications. Biometrics and Biostatistics International Journal, 8, 6-11. https://doi.org/10.15406/bbij.2019.08.00262
[11]	Abebe, B., Tesfay, M., Eyob, T. and Shanker, R. (2019) A Two-Parameter Power Garima Distribution with Properties and Applications. Annals of Biostatistics and Biometric Applications, 1, 1-7. https://doi.org/10.33552/ABBA.2019.01.000515
[12]	Shukla, K.K. (2019) Power Pranav Distribution and Its Applications. Journal of Applied Quantitative Methods, 14, 1-13.
[13]	Shanker, R. and Shukla, K.K. (2019) A Two-Parameter Power Sujatha Distribution with Properties and Application. International Journal of Mathematics and Statistics, 20, 11-22.
[14]	Shukla, K.K. and Shanker, R. (2020) Power Prakaamy Distribution and Its Applications. International Journal of Computations and Theoretical Statistics, 7, 26-36.
[15]	Hogg, R.V., McKean, J.W. and Craig, A.T. (2019) Introduction to Mathematical Statistics. 8th Edition, Pearson Education, Inc., London.
[16]	Bonferroni, C.E. (1930) Elementi di Statistica Generale, Seeber, Firenze.
[17]	Lorenz, M.O. (1905) Methods of Measuring the Concentration of Wealth. Publications of the American Statistical Association, 9, 209-219. https://doi.org/10.2307/2276207
[18]	Shaked, M. and Shanthikumar, J.G. (1994) Stochastic Orders and Their Applications. Academic Press, Boston.
[19]	Rényi, A. (1961) On Measure of Entropy and Information. In: Proceedings of the 4th Berkeley Symposium on Mathematical Statistics and Probability, Vol. 1, University of California Press, Berkeley, 547-561.
[20]	Cordeiro, G.M. and Lemonte, A.J. (2011) The β-Birnbaum-Saunders Distribution: An Improved Distribution for Fatigue Life Modelling. Computational Statistics and Data Analysis, 55, 1445-1461. https://doi.org/10.1016/j.csda.2010.10.007
[21]	Al-Aqtash, R., Lee, C. and Famoye, F. (2014) Gumbel-Weibull Distribution: Properties and Applications. Journal of Modern Applied Statistical Methods, 13, 201-225. https://doi.org/10.22237/jmasm/1414815000
[22]	Kolmogorov, A. (1933) Sulla determinazione empirica di una lgge di distribuzione. Giornale dell’Istituto Italiano degli Attuari, 4, 83-91.
[23]	Smirnoff, N. (1939) Sur les écarts de la courbe de distribution empirique. Matematicheskii Sbornik, 48, 3-26.
[24]	Scheffé, H. (1943) Statistical Inference in the Non-Parametric Case. The Annals of Mathematical Statistics, 14, 305-332. https://doi.org/10.1214/aoms/1177731355
[25]	Wolfowitz, J. (1949) Non-Parametric Statistical Inference. In: Proceedings of the Berkeley Symposium on Mathematical Statistics and Probability, University of California Press, Berkeley, 93-113.
[26]	Akaike, H. (1974) A New Look at the Statistical Model Identification. IEEE Transactions on Automatic Control, 19, 716-723. https://doi.org/10.1109/TAC.1974.1100705
[27]	Schwarz, G. (1978) Estimating the Dimension of a Model. Annals of Statistics, 6, 461-464. https://doi.org/10.1214/aos/1176344136