Modified Kumaraswamy distribution

Template:Probability distribution In probability theory, the Modified Kumaraswamy (MK) distribution is a two-parameter continuous probability distribution defined on the interval (0,1). It serves as an alternative to the beta and Kumaraswamy distributions for modeling double-bounded random variables. The MK distribution was originally proposed by Sagrillo, Guerra, and Bayer through a transformation of the Kumaraswamy distribution. Its density exhibits an increasing-decreasing-increasing shape, which is not characteristic of the beta or Kumaraswamy distributions. The motivation for this proposal stemmed from applications in hydro-environmental problems.

Definitions

Probability density function

The probability density function of the Modified Kumaraswamy distribution is

$f_X(x;𝜽)=\frac{\alpha \beta x^{\alpha -\alpha /x}(1-{\mathrm{e}}^{\alpha -\alpha /x}{)}^{\beta -1}}{x^2}$

where $𝜽=(\alpha ,\beta {)}^{\mathrm{\top }}$ , $\alpha >0$ and $\beta >0$ are shape parameters.

Cumulative distribution function

The cumulative distribution function of Modified Kumaraswamy is given by

$F_X(x;𝜽)=1-(1-{\mathrm{e}}^{\alpha -\alpha /x}{)}^{\beta }$

where $𝜽=(\alpha ,\beta {)}^{\mathrm{\top }}$ , $\alpha >0$ and $\beta >0$ are shape parameters.

Quantile function

The inverse cumulative distribution function (quantile function) is

$Q_X(u;𝜽)=\frac{\alpha }{\alpha -\log (1-(1-u{)}^{1/\beta })}$

Properties

Moments

The hth statistical moment of X is given by:

$\text{E}\left(X^h\right)=\alpha \beta {\mathrm{e}}^{\alpha }\sum _{i=0}^{\mathrm{\infty }}(-1{)}^i\left(\begin{array}{c}\beta -1\\ i\end{array}\right){\mathrm{e}}^{\alpha i}(\alpha +\alpha i{)}^{h-1}\mathrm{\Gamma }[1-h,(i+1)\alpha ]$

Mean and Variance

Measure of central tendency, the mean $(\mu )$ of X is:

$\mu =\text{E}(X)=\alpha \beta {\mathrm{e}}^{\alpha }\sum _{i=0}^{\mathrm{\infty }}(-1{)}^i\left(\begin{array}{c}\beta -1\\ i\end{array}\right){\mathrm{e}}^{\alpha i}\mathrm{\Gamma }[0,(i+1)\alpha ]$

And its variance $({\sigma }^2)$:

${\sigma }^2=\text{E}(X^2)={\alpha }^2\beta {\mathrm{e}}^{\alpha }\sum _{i=0}^{\mathrm{\infty }}(-1{)}^i\left(\begin{array}{c}\beta -1\\ i\end{array}\right){\mathrm{e}}^{\alpha i}(i+1)\mathrm{\Gamma }[-1,(i+1)\alpha ]-{\mu }^2$

Parameter estimation

Sagrillo, Guerra, and Bayer suggested using the maximum likelihood method for parameter estimation of the MK distribution. The log-likelihood function for the MK distribution, given a sample $x_1,\dots ,x_n$, is:

$\begin{array}{cc}\ell (𝜽)=& n\alpha +n\log \left(\alpha \right)+n\log \left(\beta \right)-\alpha \sum _{i=1}^n\frac{1}{x_i}-2\sum _{i=1}^n\log (x_i)\\ & +(\beta -1)\sum _{i=1}^n\log (1-{\mathrm{e}}^{\alpha -\alpha /x_i}).\end{array}$

The components of the score vector $U\left(𝜽\right)=[\frac{\partial \ell (𝜽)}{\partial \alpha },\frac{\partial \ell (𝜽)}{\partial \beta }]$ are

$\begin{array}{c}\frac{\partial \ell (𝜽)}{\partial \alpha }=n+\frac{n}{\alpha }+(\beta -1){\mathrm{e}}^{\alpha }\sum _{i=1}^n\frac{x_i-1}{x_i({\mathrm{e}}^{\alpha }-{\mathrm{e}}^{\alpha /x_i})}-\sum _{i=1}^n\frac{1}{x_i}\end{array}$

and

$\begin{array}{c}\frac{\partial \ell (𝜽)}{\partial \beta }=\frac{n}{\beta }+\sum _{i=1}^n\log (1-{\mathrm{e}}^{\alpha -\alpha /x_i})\end{array}$

The MLEs of $𝜽$, denoted by $\widehat{𝜽}={(\hat{\alpha },\hat{\beta })}^{\mathrm{\top }}$, are obtained as the simultaneous solution of $𝑼(𝜽)=0$, where $0$ is a two-dimensional null vector.

Related distributions

• If $X\sim \text{MK}(\alpha ,\beta )$, then $\{1-\frac{1}{X}\}\sim \text{K}(\alpha ,\beta )$ (Kumaraswamy distribution)
• If $X\sim \text{MK}(\alpha ,\beta )$, then $\frac{1}{X}-1\sim$ Exponentiated exponential (EE) distribution
• If $X\sim \text{MK}(1,\beta )$, then $\exp \{1-\frac{1}{X}\}\sim \text{Beta}(1,\beta )$. (Beta distribution)
• If $X\sim \text{MK}(\alpha ,1)$, then $\exp \{1-\frac{1}{X}\}\sim \text{Beta}(\alpha ,1)$.
• If $X\sim \text{MK}(\alpha ,\beta )$, then $\frac{1}{X}-1\sim \text{Exp}(\alpha )$ (Exponential distribution).

Applications

The Modified Kumaraswamy distribution was introduced for modeling hydro-environmental data. It has been shown to outperform the Beta and Kumaraswamy distributions for the useful volume of water reservoirs in Brazil. It was also used in the statistical estimation of the stress-strength reliability of systems.

See also

• Kumaraswamy distribution

External links

• Github with the distribution functions in R