Research Article  Open Access
Ramadan A. ZeinEldin, Muhammad Ahsan ul Haq, Sharqa Hashmi, Mahmoud Elsehety, M. Elgarhy, "Type II Half Logistic Kumaraswamy Distribution with Applications", Journal of Function Spaces, vol. 2020, Article ID 1343596, 15 pages, 2020. https://doi.org/10.1155/2020/1343596
Type II Half Logistic Kumaraswamy Distribution with Applications
Abstract
In this paper, a new distribution with a unit interval named type II half logistic Kumaraswamy (TIIHLKw) distribution is proposed. Its density and distribution functions are presented using alternate expressions. This distribution is obtained by adding an extra parameter in the existing model to rise its ability fitting complex data sets. Some important statistical properties of TIIHLKw distribution are derived. The estimation of the parameters is obtained by numerous wellrecognized approaches and simulation study confirmed the efficiencies of estimates such obtained. We apply the related model to practical datasets, and it is concluded that the proposed model is the best by model selection criteria than other competitive models.
1. Introduction
Over the last few years, inspired by the increasing demand of probability distributions in many fields, many generalized distributions have been studied. Most of them are proposed by the addition of more parameters to wellknown probability distributions that exist in literature to make them flexible. For instance, Haq et al. [1] proposed and studied the generalized odd Burr III (GOBIII) family of distributions along with its important characterizations. Ahmed [2] proposed a new model derived by the transformation of the baseline model. Different shapes of failure function can be formed such as increasing and bathtub. Some mathematical properties are obtained. So many research works have been done for proposing more flexible generalized probability distributions such as Alzaatreh et al. [3], Haq et al. [4], Hashmi et al. [5], Elgarhy et al. [6], and ZeinEldin et al. [7]. Alshenawy [8] suggested a new oneparameter distribution. Several statistical properties and characteristics of the proposed distribution are derived along with estimation under Type II censoring. The research is concluded on the basis of a simulation study and real data analysis.
The Kumaraswamy (Kw) distribution was introduced and studied by Kumaraswamy [9] with unit interval, denoted by 𝐾w(a, b), with cumulative distribution function (cdf) is and its related probability density function (pdf) is
The Kw density has one of these shapes depending upon parameter values, unimodal (), bathtub (), increasing ( and ), decreasing (D) ( and ), or steady (both a and b equal to 1).
The behavior of Kw distribution is analogous to the beta distribution but simpler due to closedform of both its pdf and cdf. Boundary behavior and the major special models are also the same in both Beta and Kw distributions. This distribution could be a good substitute in situations wherein reality the bounds are finite i.e., (0, 1).
The Kw is originally developed as a lifetime distribution. Many authors studied and developed the generalizations of Kw distribution such as exponentiated Kw distribution studied by Lemonte et al. [10], ElSherpieny, and Ahmed [11] proposed Kw Kw distribution, transmuted Kw distribution studied by Khan et al. [12], Sharma and Chakrabarty [13] studied size biased Kw distribution, George and Thobias [14] introduced MarshallOlkin Kw distribution, exponentiated generalized Kw distribution studied by Elgarhy et al. [6], type II Topp Leone inverted Kw by ZeinEldin et al. [15], type I half logistic inverted Kw by ZeinEldin et al. [16], truncated inverted Kw by Bantan et al. [17], and Ghosh [18] introduced bivariate and multivariate weighted Kw distributions.
Hassan et al. [19] proposed Type II half logisticG (TIIHLG). The TIIHLG distribution is expressed by its cdf given by where is cdf of baseline model with parameter vector and is cdf derived by the TX generator proposed by Alzaatreh [3]. The pdf of the TIIHLG family is given as
This article is dedicated to both its mathematical and application features. A significant portion is kept for the estimation of the parameters through various methods including maximum likelihood estimation (MLE), leastsquare estimation (LSE), weighted least square estimation (WLSE), percentiles estimation (PCE), and Cramervon Mises estimation (CVE).
The core purpose of this research is to suggest a simpler and more flexible model called Type II Half Logistic Kw (TIIHLKw) distribution. This article is organized in the following manner: Section 2 is dedicated for the proposition of type II half logisticKw (TIIHLKw) distribution. Section 3 deals with the leading statistical properties of this model. Important binomial expansions of density and distribution functions are presented which involve binomial expansions. In Section 4, an extensive study of five different methods of estimation is carried out, with all derivations and detailed discussions. A comprehensive simulation study is conducted to compute the biases and efficiency for parameters and compare the performances of five estimation approaches stated above in the next section. Section 6 is devoted to the real data application of TIIHLKwD to show the importance of TIIHLKw distribution. Lastly, the conclusion is given in Section 7.
2. The TIIHLKw Distribution
In this section, we examine the usefulness, and flexibility of a new associate of type II halflogisticG family having (0, 1), using Kw distribution as the baseline. The pdf and cdf of Kw distribution (2 shape parameters: ) are given as follows
The random variable (r.v.) X follows TIIHLKw distribution if its cdf is obtained by inserting (6) in (3)
The pdf of TIIHLKw distribution is as follows
The survival function of TIIHLKw distribution is
The failure rate function of TIIHLKwD is
The flexibility of TIIHLKw distribution can be illustrated in Figures 1, 2, and 3. The pdf plots for the TIIHLKw distribution are given in Figures 1 and 2, and plots of hrf are given in Figure 3.
3. Some Mathematical Properties
The properties of TIIHLKwD are derived here. After this, we consider a r.v. X follows the pdf (8) and cdf (7).
3.1. Quantile Function
The quantile function of is denoted by , defined as the inverse function of is , given in (6).
Also,
Simulated values from TIIHLKwD can be utilized in the simulation study. The r. v. follows rectangular distribution on the interval 0 to 1, i.e., follows TIIHLKw distribution. In particular, the first three quartiles are obtained by putting , and 3/4, respectively, in (11) and (12).
On differentiation, we can have the density of quantile function
3.2. Alternate Representation
Here, we present the expansion of the pdf and cdf for TIIHLKw distribution for further mathematical manipulation.
The binomial theorem, for and , can be expressed as
Then, by applying (13) in (8), the pdf of TIIHLKwD becomes for .
Another form of pdf (8) can be obtained by means of the following expansion for and . By applying (16) in pdf (8), we get where . Now considering (13), we have
Inserting this expansion in (17), we have pdf of the following form where .
Another formula can be formed from pdf (19), which is given in (20) using an infinite linear combination where
Proposition 1. Let “h” be a positive integer. The expansion of cdf can be written in following form: where .
Proof. The cdf is obtained, for an integer by using (16).
Once more binomial expansion is applied to and it can be written as
The required result is obtained by combing some expression together, completing the proof.
3.3. Probability Weighted Moments:
The probabilityweighted moments (PWMs) are used to study some more characteristics of the probability distribution. Under the specified setting discussed above, PWMs are denoted by , can be defined as
The PWMs of TIIHLKw are obtained by substituting (21) and (22) into (25), as follows
Then, where Also, is beta function.
3.4. Moments
The moments have a vital role in the study of the distribution and real data applications. Now, we get the moment for the TIIHLKwD. Under the specified assumptions, the moment is obtained as where
The mean () and variance (var) of TIIHLKw distribution can be derived as
Also, the coefficients of skewness and kurtosis of TIIHLKw distribution are given by
The summary measures, mean, variance, skewness () and kurtosis () values are presented in Table 1. The plots of mean and var for the TIIHLKwD are given in Figure 4 and graphs of skewness and kurtosis are presented in Figure 5 for different parameter ranges.

We see a monotonic variation in these measures caused by variation in parameters , and .
3.5. MomentGenerating Function
The momentgenerating function of , using moments about the origin (29), is obtained as
3.6. Incomplete Moments
The incomplete moment of TIIHLKwD can be obtained by using (19) is where is incomplete beta function.
3.7. The Mean Deviation
Following are expressions used to get mean deviation about mean and median, respectively where = Median of and is the initial incomplete moment. Now
3.8. Bonferroni and Lorenz Curves
Bonferroni and Lorenz curves are important applications of the first incomplete moment. These curves are mostly used in different fields of life such as economics, reliability analysis, demographic studies, life testing, life insurance, and medical technology. The Lorenz and Bonferroni curves are obtained, respectively, as follows and,
3.9. Order Statistics
In statistical theory, order statistics is widely applied and practiced. Let be r.vs. with their corresponding cdfs . Let be the related ordered r. sample of size n, then the density of order statistic is given as where and β(.,.) is the beta function. The pdf of the order statistic of is obtained by putting (21) and (22) in (38), changing with where
More, the moment of order statistics for TIIHLKwD is given by
By substituting (39) in (41), we have
Then,
3.10. Rényi Entropy
Renyi (1961) proposed and used this measure. It can be obtained by
Applying binomial expansion (16) in (8) then can be written as
So, the Rényi entropy of TIIHLKw distribution is as follows where .
3.11. StressStrength Reliability
This subsection deals with the stressstrength parameter of TIIHLKw distribution. Let be the strength of a structure with a stress , and if follows and follows , provided and are statistically independent r.vs.,
Proposition 2. Under the assumption discussed above, we have
Proof. The reliability is defined by
Then, we can write where .which completes the proof.
4. Inference
This section is dedicated to estimation aspects of TIIHLKw distribution, assuming that population parameters () are unknown and can be estimated using different methods of estimation including ML, LS, WLS, PC, and CV.
4.1. Maximum Likelihood Estimation
For a random sample of from the distribution, the loglikelihood function for is
The members of are given below equating (52), (53), and (54) to zero and solve them simultaneously give the ML estimates (MLEs) of The iterative algorithm is used to obtain the numerical solution of these nonlinear equations such as the NewtonRaphson method.
4.2. Ordinary and Weighted LS Estimation (LSEs and WLSEs):
The LSEs of , and can be obtained by minimizing the sum of squares of errors with respect to parameters. Suppose is a random sample of size from TIIHLKwD and suppose be the related ordered sample.
The sum is independent of the unknown parameters.
Equivalently,
The following function is minimized with respect to and to get WLS estimators of model parameters.
i.e., , and .
4.3. PC Estimators (PCEs)
Under the specification defined above including order statistics having relationship . In the PC method of estimation, the estimators of , and are derived by minimizing the following expression with respect to , and .
4.4. The Cramervon Mises Minimum Distance Estimation
The CVE is another estimation method based on minimum distance. The CV estimators are obtained by minimizing, with respect to , and .
CV minimum distance estimators provide empirical evidence that the bias of this estimator is smaller than the other minimum distance estimators.
5. Simulation Study
In this section, we present a Monte Carlo (MC) simulation study in order to illustrate the behavior of different estimates. We consider four random sample sizes: , 100, 200, and 500 from the TIIHLKw distribution and the samples are drawn 10,000 times. Four specifications of the parameters are used in this simulation study, given as, , , , . For sample generated, the MLE, LSE, WLSE, CVE, PCE, and MPSE of estimators are computed numerically. Then, the estimates of all methods and their mean square errors (MSEs) are documented in Tables 1–4.
Estimates of parameters and their corresponding MSEs are calculated using the same sample under different methods of estimation for four sets of parameters mentioned above for small sample size () to sufficient large sample size, that is, .
The entries of Table 2 to Table 5 show that estimates are reliable and consistent. MSEs reduce as sample size () increases under each method of estimation. The summary of these four tables is presented in Table 5.
 
10^{^3}. 

 
10^{^3}. 
 
10^{^3}. 



Using the entries of Table 6 for different parametric combinations, we can conclude that the MLE method outperforms than all other estimation methods (with an overall score of 18.5). Therefore, depending on the simulation study, the MLE method performs best for TIIHLKwD.
6. Applications
The TIIHLKw distribution aims at providing an alternative distribution to fit data on the unit interval to other distributions available in the literature. Here, we used the following probability distributions as competitor models; (i)The Kumaraswamy (Kw) distribution (ii)Transmuted Kumaraswamy (TKw) distribution (iii)Size biased Kumaraswamy (SBKw) distribution (iv)Beta distribution
We use the following accuracy measures for model comparison: Akaike Information Criterion (AIC), Bayesian Information Criterion (BIC), loglikelihood, AndersonDarling (A), and Cramer–von Mises (W). RLanguage is used for numerical computations.
1^{st} data set. The data set was taken from Dasgupta [20], and considered on burr (in millimeters), with hole diameter and a sheet thickness of 12 mm and 3.15 mm, respectively.
2^{nd} data set. This data set was taken from [21], and considered measurements on petroleum rock samples from a petroleum reservoir.
The total test time (TTT) for both datasets are presented in Figure 6. We can observe that the shape of TTT plots is concave for both datasets, which demonstrates an increasing failure rate.
The MLEs for TIIHLKw distribution along with some adequacy measures are presented in Tables 7 and 8.
Hence, it is concluded that the new model provides the better fit. Figures 7 and 8 show the estimated densities, cdfs, estimated survival functions, and PP plots for the considered distributions of both data sets, respectively. We note that the proposed model is more appropriated to fit the data than the other competing models.
7. Conclusion
In this article, a new TypeII Half Logistic Kumaraswamy distribution is proposed. Some characteristics of the TIIHLKw distribution including linear combination expressions for the density function, probability weighted moments, moments, incomplete moments, quantile function, mean deviation about mean and about median, Bonferroni and Lorenz curves, order statistics, stressstrength reliability, and Rényi entropy are derived. The ML method is used to estimate model parameters. An extensive simulation study is conducted to compare several wellknown estimation methods, including the method of maximum likelihood estimation, methods of least squares and weighted least squares estimation, and method of Cramervon Mises minimum distance estimation. The simulation study showed the reliability and efficiency of the estimates. Finally, by considering the method of maximum likelihood estimation, the new model is fitted to two practical data sets. The applications on real data sets validated the significance of the new distribution.
Appendix
A.1. Data I
0.04, 0.02, 0.06, 0.12, 0.14, 0.08, 0.22, 0.12, 0.08, 0.26, 0.24, 0.04, 0.14, 0.16, 0.08, 0.26, 0.32, 0.28, 0.14, 0.16, 0.24, 0.22, 0.12, 0.18, 0.24, 0.32, 0.16, 0.14, 0.08, 0.16, 0.24, 0.16, 0.32, 0.18, 0.24, 0.22, 0.16, 0.12, 0.24, 0.06, 0.02, 0.18, 0.22, 0.14, 0.06, 0.04, 0.14, 0.26, 0.18, 0.16.
A.2. Data II
0.0903296, 0.2036540, 0.2043140, 0.2808870, 0.1976530, 0.3286410, 0.1486220, 0.1623940, 0.2627270, 0.1794550, 0.3266350, 0.2300810, 0.1833120, 0.1509440, 0.2000710, 0.1918020, 0.1541920, 0.4641250, 0.1170630, 0.1481410, 0.1448100, 0.1330830, 0.2760160, 0.4204770, 0.1224170, 0.2285950, 0.1138520, 0.2252140, 0.1769690, 0.2007440, 0.1670450, 0.2316230, 0.2910290, 0.3412730, 0.4387120, 0.2626510, 0.1896510, 0.1725670, 0.2400770, 0.3116460, 0.1635860, 0.1824530, 0.1641270, 0.1534810, 0.1618650, 0.2760160, 0.2538320, 0.2004470.
Data Availability
Data is present in appendix.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was funded by the Deanship of Scientific Research (DSR), King Abdul Aziz University, Jeddah, under grant No. (DF2873051441). The authors gratefully acknowledge the DSR technical and financial support.
References
 M. A. u. Haq, M. Elgarhy, and S. Hashmi, “The generalized odd Burr III family of distributions: properties, applications and characterizations,” Journal of Taibah University for Science, vol. 13, no. 1, pp. 961–971, 2019. View at: Publisher Site  Google Scholar
 A. M. T. Abd ElBar and I. E. Ragab, “On weighted exponentialGompertz distribution: properties and application,” Journal of Taibah University for Science, vol. 13, no. 1, pp. 616–627, 2019. View at: Publisher Site  Google Scholar
 A. Alzaatreh, C. Lee, and F. Famoye, “A new method for generating families of continuous distributions,” Metron, vol. 71, no. 1, pp. 63–79, 2013. View at: Publisher Site  Google Scholar
 M. A. u. Haq, M. Elgarhy, S. Hashmi, G. Ozel, and Q. u. Ain, “Transmuted Weibull power function distribution: its properties and applications,” Journal of Data Science, vol. 397, p. 418, 2018. View at: Publisher Site  Google Scholar
 S. Hashmi, M. A. u. Haq, and R. M. Usman, “A generalized exponential distribution with increasing, decreasing and constant shape hazard curves,” Electronic Journal of Applied Statistical Analysis, vol. 12, no. 1, pp. 223–244, 2019. View at: Google Scholar
 M. Elgarhy, M. A. u. Haq, and Q. u. Ain, “Exponentiated generalized Kumaraswamy distribution with applications,” Annals of Data Science, vol. 5, no. 2, pp. 273–292, 2018. View at: Publisher Site  Google Scholar
 R. A. ZeinEldin, M. Ahsan Ul Haq, S. Hashmi, and M. Elsehety, “Alpha power transformed inverse Lomax distribution with different methods of estimation and applications,” Complexity, vol. 2020, Article ID 1860813, 15 pages, 2020. View at: Publisher Site  Google Scholar
 R. Alshenawy, “A new one parameter distribution: properties and estimation with applications to complete and type II censored data,” Journal of Taibah University for Science, vol. 14, no. 1, pp. 11–18, 2020. View at: Publisher Site  Google Scholar
 P. Kumaraswamy, “A generalized probability density function for doublebounded random processes,” Journal of Hydrology, vol. 46, no. 12, pp. 79–88, 1980. View at: Publisher Site  Google Scholar
 A. J. Lemonte, W. B. Souzaa, and G. M. Cordeiro, “The exponentiated Kumaraswamy distribution and its logtransform,” Brazilian Journal of Probability and Statistics, vol. 27, no. 1, pp. 31–53, 2013. View at: Publisher Site  Google Scholar
 E.S. ElSherpieny and M. A. Ahmed, “On the Kumaraswamy Kumaraswamy distribution,” International Journal of Basic and Applied Sciences, vol. 3, no. 4, 2014. View at: Publisher Site  Google Scholar
 M. S. Khan, R. King, and I. L. Hudson, “Transmuted Kumaraswamy distribution,” Statistics in Transition New Series, vol. 17, no. 2, pp. 183–210, 2016. View at: Publisher Site  Google Scholar
 D. Sharma and T. K. Chakrabarty, “On size biased Kumaraswamy distribution,” Statistics, Optimization & Information Computing, vol. 4, no. 3, 2016. View at: Publisher Site  Google Scholar
 R. George and S. Thobias, “MarshallOlkin Kumaraswamy distribution,” International Mathematical Forum, vol. 12, no. 2, pp. 47–69, 2017. View at: Publisher Site  Google Scholar
 R. A. ZeinEldin, F. Jamal, C. Chesneau, and M. Elgarhy, “Type II ToppLeone inverted Kumaraswamy distribution with statistical inference and applications,” Symmetry, vol. 11, no. 12, p. 1459, 2019. View at: Publisher Site  Google Scholar
 R. A. ZeinEldin, C. Chesneau, F. Jamal, and M. Elgarhy, “Statistical properties and different methods of estimation for type I half logistic inverted Kumaraswamy distribution,” Mathematics, vol. 7, no. 10, p. 1002, 2019. View at: Publisher Site  Google Scholar
 R. A. Bantan, F. Jamal, C. Chesneau, and M. Elgarhy, “Truncated inverted Kumaraswamy generated family of distributions with applications,” Entropy, vol. 21, no. 11, p. 1089, 2019. View at: Publisher Site  Google Scholar
 I. Ghosh, “Bivariate and multivariate weighted Kumaraswamy distributions: theory and applications,” Journal of Statistical Theory and Applications, vol. 18, no. 3, p. 198, 2019. View at: Publisher Site  Google Scholar
 A. S. Hassan, M. Elgarhy, and M. Shakil, “Type II half logistic family of distributions with applications,” Pakistan Journal of Statistics and Operation Research, vol. 13, no. 2, pp. 245–264, 2017. View at: Publisher Site  Google Scholar
 R. Dasgupta, “On the distribution of burr with applications,” Sankhya B, vol. 73, no. 1, pp. 1–19, 2011. View at: Publisher Site  Google Scholar
 G. M. Cordeiro and R. dos Santos Brito, “The beta power distribution,” Brazilian Journal of Probability and Statistics, vol. 26, no. 1, pp. 88–112, 2012. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Ramadan A. ZeinEldin et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.