Abstract
In this article, a structural modification of the Kumaraswamy distribution yields a new two-parameter distribution defined on (0,1), called the modified Kumaraswamy distribution. It has the advantages of being (i) original in its definition, mixing logarithmic, power and ratio functions, (ii) flexible from the modeling viewpoint, with rare functional capabilities for a bounded distribution-in particular, N-shapes are observed for both the probability density and hazard rate functions-and (iii) a solid alternative to its parental Kumaraswamy distribution in the first-order stochastic sense. Some statistical features, such as the moments and quantile function, are represented in closed form. The Lambert function and incomplete beta function are involved in this regard. The distributions of order statistics are also explored. Then, emphasis is put on the practice of the modified Kumaraswamy model in the context of data fitting. The well-known maximum likelihood approach is used to estimate the parameters, and a simulation study is conducted to examine the performance of this approach. In order to demonstrate the applicability of the suggested model, two real data sets are considered. As a notable result, for the considered data sets, statistical benchmarks indicate that the new modeling strategy outperforms the Kumaraswamy model. The transmuted Kumaraswamy, beta, unit Rayleigh, Topp-Leone and power models are also outperformed.