Modified Kumaraswamy distribution#Related distributions
{{Short description|Continuous probability distribution}}
{{Orphan|date=May 2025}}
{{Probability distribution |
name =Modified Kumaraswamy|
type =density|
pdf_image =Image:MKpdf4.svg|
cdf_image =Image:MKcdf2.svg|
support = |
pdf =|
cdf =|
mean =
\begin{pmatrix}
\beta -1 \\
i
\end{pmatrix} \mathrm{e}^{\alpha i} \Gamma \left[ 0, \left( i+1 \right) \alpha \right]|
variance =
\begin{pmatrix}
\beta -1 \\
i
\end{pmatrix} \mathrm{e}^{\alpha i}(i+1) \Gamma \left[ -1, \left( i+1 \right) \alpha \right] - \mu^2|
quantile =|
mgf =
\begin{pmatrix}
\beta -1 \\
i
\end{pmatrix} \mathrm{e}^{\alpha i}(\alpha + \alpha i)^{h-1} \Gamma \left[1-h, \left( i+1 \right) \alpha \right]|
}}
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 {{Cite journal|last1=Sagrillo|first1=M.|last2=Guerra|first2=R. R.|last3=Bayer|first3=F. M.|date=2021|title=Modified Kumaraswamy distributions for double bounded hydro-environmental data|journal=Journal of Hydrology|volume=603|doi=10.1016/j.jhydrol.2021.127021|bibcode=2021JHyd..60327021S }} 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} \left( x;\boldsymbol{\theta} \right) = \frac{\alpha \beta x^{\alpha - \alpha/x} (1-\mathrm{e}^{\alpha - \alpha/x })^{\beta-1}}{x^2}
where , and are shape parameters.
= Cumulative distribution function =
The cumulative distribution function of Modified Kumaraswamy is given by
:
where , and are shape parameters.
=Quantile function=
The inverse cumulative distribution function (quantile function) is
:
Properties
= Moments =
The hth statistical moment of X is given by:
:
\begin{pmatrix}
\beta -1 \\
i
\end{pmatrix} \mathrm{e}^{\alpha i}(\alpha + \alpha i)^{h-1} \Gamma \left[1-h, \left( i+1 \right) \alpha \right]
= Mean and Variance =
Measure of central tendency, the mean of X is:
:
\begin{pmatrix}
\beta -1 \\
i
\end{pmatrix} \mathrm{e}^{\alpha i} \Gamma \left[ 0, \left( i+1 \right) \alpha \right]
And its variance :
:
\begin{pmatrix}
\beta -1 \\
i
\end{pmatrix} \mathrm{e}^{\alpha i}(i+1) \Gamma \left[ -1, \left( i+1 \right) \alpha \right] - \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 , is:
:
\begin{align}
\ell (\boldsymbol{\theta}) = &\,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{align}
The components of the score vector are
:
\begin{align}
\frac{\partial \ell (\boldsymbol{\theta})}{\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{align}
and
:
\begin{align}
\frac{\partial \ell (\boldsymbol{\theta})}{\partial \beta} = \frac{n}{\beta} + \sum_{i=1}^{n} \log (1-\mathrm{e}^{\alpha-\alpha/x_i})
\end{align}
The MLEs of , denoted by , are obtained as the simultaneous solution of , where is a two-dimensional null vector.
Related distributions
- If , then (Kumaraswamy distribution)
- If , then Exponentiated exponential (EE) distribution{{Cite journal|last1=Gupta|first1=R.D.|last2=Kundu|first2=D|date=1999|title=Theory & Methods: Generalized exponential distributions|journal=Australian & New Zealand Journal of Statistics|volume=41|issue=2 |pages=173–188|doi=10.1111/1467-842X.00072}}
- If , then . (Beta distribution)
- If , then .
- If , then (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.{{Cite journal|last1=Kohansal|first1=Akram|last2=Pérez-González|first2=Carlos J|last3=Fernández|first3=Arturo J|date=2023|title=Inference on the stress-strength reliability of multi-component systems based on progressive first failure censored samples|journal=Proceedings of the Institution of Mechanical Engineers, Part O: Journal of Risk and Reliability|volume=238|issue=5 |pages=1053–1073 |doi=10.1177/1748006X231188075}}
See also
References
{{Reflist}}
External links
- [https://github.com/murilosagrillo/MKRMK Github with the distribution functions in R]