Exponentiated Kumaraswamy-Dagum distribution with applications to income and lifetime data

A new family of distributions called exponentiated Kumaraswamy-Dagum (EKD) distribution is proposed and studied. This family includes several well known sub-models, such as Dagum (D), Burr III (BIII), Fisk or Log-logistic (F or LLog), and new sub-models, namely, Kumaraswamy-Dagum (KD), Kumaraswamy-Burr III (KBIII), Kumaraswamy-Fisk or Kumaraswamy-Log-logistic (KF or KLLog), exponentiated Kumaraswamy-Burr III (EKBIII), and exponentiated Kumaraswamy-Fisk or exponentiated Kumaraswamy-Log-logistic (EKF or EKLLog) distributions. Statistical properties including series representation of the probability density function, hazard and reverse hazard functions, moments, mean and median deviations, reliability, Bonferroni and Lorenz curves, as well as entropy measures for this class of distributions and the sub-models are presented. Maximum likelihood estimates of the model parameters are obtained. Simulation studies are conducted. Examples and applications as well as comparisons of the EKD and its sub-distributions with other distributions are given.

http://www.jsdajournal.com/content/1/1/8 statistical properties of this distribution. Domma et al. (2011) obtained the maximum likelihood estimates of the parameters of Dagum distribution for censored data. Domma and Condino (2013) presented the beta-Dagum distribution. Cordeiro et al. (2013) proposed the beta exponentiated Weibull distribution. Cordeiro et al. (2010) introduced and studied some mathematical properties of the Kumaraswamy Weibull distribution. Oluyede and Rajasooriya (2013) developed the Mc-Dagum distribution and presented its statistical properties. See references therein for additional results.
The pdf and cumulative distribution function (cdf ) of Dagum distribution are given by: and for x > 0, where λ is a scale parameter, δ and β are shape parameters. Dagum (1977) refers to his model as the generalized logistic-Burr distribution. The k th raw or non central moments are given by for k < δ, and λ, β > 0, where B(·, ·) is the beta function. The q th percentile is In this paper, we present generalizations of the Dagum distribution via Kumaraswamy distribution and its exponentiated version. This leads to the exponentiated Kumaraswamy Dagum distribution.
The motivation for the development of this distribution is the modeling of size distribution of personal income and lifetime data with a diverse model that takes into consideration not only shape, and scale but also skewness, kurtosis and tail variation. Also, the EKD distribution and its sub-models has desirable features of exhibiting a non-monotone failure rate, thereby accommodating different shapes for the hazard rate function and should be an attractive choice for survival and reliability data analysis. This paper is organized as follows. In section 3, we present the exponentiated Kumaraswamy-Dagum distribution and its sub models, as well as series expansion, hazard and reverse hazard functions. Moments, moment generating function, Lorenz and Bonferroni curves, mean and median deviations, and reliability are obtained in section 4. Section 5 contains results on the distribution of the order statistics and Renyi entropy. Estimation of model parameters via the method of maximum likelihood is presented in section 6. In section 7, various simulations are conducted for different sample sizes. Section 8 contains examples and applications of the EKD distribution and its sub-models, followed by concluding remarks.

Methods, results and discussions
Methods, results and discussions for the class of EKD distributions are presented in sections 3 to 8. These sections include the sub-models, series expansion of the pdf, closed form expressions for the hazard and reverse hazard functions, moments, moment generating function, Bonferroni and Lorenz curves, reliability, mean and median deviations, http://www.jsdajournal.com/content/1/1/8 distribution of order statistics and entropy, as well as estimation of model parameters and applications.

The exponentiated Kumaraswamy-Dagum distribution
In this section, we present the proposed distribution and its sub-models. Series expansion, hazard and reverse hazard functions are also studied in this section.

The EKD distribution
In general, the EKD distribution is For large values of x, and for θ > 1(< 1), the multiplicative factor θ[F KD (x)] θ−1 > 1(< 1), respectively. The reverse statement holds for smaller values of x.
Consequently, this implies that the ordinary moments of g EKD (x) are larger (smaller) than those of f KD (x) when θ > 1(< 1).
Replacing the dependent parameter βψ by α, the cdf and pdf of the EKD distribution are given by G EKD (x; α, λ, δ, φ, θ) and for α, λ, δ, φ, θ > 0, and x > 0, respectively. The quantile function of the EKD distribution is in closed form, Plots of the pdf for some combinations of values of the model parameters are given in Figure 1. The plots indicate that the EKD pdf can be decreasing or right skewed. The EKD distribution has a positive asymmetry. http://www.jsdajournal.com/content/1/1/8 Figure 1 Graph of pdfs.

Series expansion of EKD distribution
We apply the series expansion for b > 0 and |z| < 1, to obtain the series expansion of the EKD distribution. By using equation (6), where Dagum(α, δ, λ) distribution, α and δ are shape parameters, and λ is a scale parameter. In the Exponentiated-Kumaraswamy(ψ, φ, θ) distribution, ψ is a skewness parameter, φ is a tail variation parameter, and the parameter θ characterizes the skewness, kurtosis, and tail of the distribution.
Consequently, for the EKD(α, λ, δ, φ, θ) distribution, α is shape and skewness parameter, δ is shape parameter, λ is a scale parameter, φ is a tail variation parameter, and the parameter θ characterizes the skewness, kurtosis, and tail of the distribution.

Hazard and reverse hazard function
The hazard function of the EKD distribution is Plots of the hazard function are presented in Figure 2. The plots show various shapes including monotonically decreasing, unimodal, and bathtub followed by upside down bathtub shapes with five combinations of the values of the parameters. This attractive flexibility makes the EKD hazard rate function useful and suitable for non-monotone empirical hazard behaviors which are more likely to be encountered or observed in real life situations. Unfortunately, the analytical analysis of the shape of both the density (except for zero modal when αδ ≤ 1, and unimodal if αδ > 1, both for φ = θ = 1,) and http://www.jsdajournal.com/content/1/1/8

Figure 2 Graphs of hazard functions.
hazard rate function seems to be very complicated. We could not determine any specific rules for the shapes of the hazard rate function.
The reverse hazard function of the EKD distribution is

Moments, moment generating function, Bonferroni and Lorenz curves, mean and median deviations, and reliability
In this section, we present the moments, moment generating function, Bonferroni and Lorenz curves, mean and median deviations as well as the reliability of the EKD distribution. The moments of the sub-models can be readily obtained from the general results. (7), then the s th raw moment of the EKD distribution is given by

Moments and moment generating function
where , and s < δ. http://www.jsdajournal.com/content/1/1/8 The moment generating function of the EKD distribution is given by for r < δ.

Bonferroni and Lorenz curves
Bonferroni and Lorenz curves are widely used tool for analyzing and visualizing income inequality. Lorenz curve, L(p) can be regarded as the proportion of total income volume accumulated by those units with income lower than or equal to the volume a, and Bonferroni curve, B(p) is the scaled conditional mean curve, that is, ratio of group mean income of the population. Plots of Bonferroni and Lorenz curves are given in Figure 3.
, then Bonferroni and Lorenz curves are given by respectively, for 0 ≤ p ≤ 1, and q = G −1 EKD (p). The mean of the EKD distribution is obtained from equation (10) with s = 1 and the quantile function is given in equation (5). Consequently,

Mean and median deviations
If X has the EKD distribution, we can derive the mean deviation about the mean μ = E(X) and the median deviation about the median M from respectively. The mean μ is obtained from equation (10) with s = 1, and the median M is given by equation (5) when q = 1 2 .

Order statistics and entropy
In this section, the distribution of the k th order statistic and Renyi entropy (Renyi 1960) for the EKD distribution are presented. The entropy of a random variable is a measure of variation of the uncertainty.

Order statistics
The pdf of the k th order statistics from a pdf f (x) is Using equation (6), the pdf of the k th order statistic from EKD distribution is given by

Estimation of model parameters
In this section, we present estimates of the parameters of the EKD distribution via method of maximum likelihood estimation. The elements of the score function are presented. There are no closed form solutions to the nonlinear equations obtained by setting the elements of the score function to zero. Thus, the estimates of the model parameters must be obtained via numerical methods.

Asymptotic confidence intervals
In this section, we present the asymptotic confidence intervals for the parameters of the EKD distribution. The expectations in the Fisher Information Matrix (FIM) can be obtained numerically. Letˆ = (α,λ,δ,φ,θ) be the maximum likelihood estimate of = (α, λ, δ, φ, θ). Under the usual regularity conditions and that the parameters are in the interior of the parameter space, but not on the boundary, we have: The asymptotic behavior is still valid if I( ) is replaced by the observed information matrix evaluated atˆ , that is J(ˆ ). The multivariate normal distribution where the mean vector 0 = (0, 0, 0, 0, 0) T , can be used to construct confidence intervals and confidence regions for the individual model parameters and for the survival and hazard rate functions. The approximate 100(1 − η)% two-sided confidence intervals for α, λ, δ, φ and θ are given by: is the upper η 2 th percentile of a standard normal distribution.
We can use the likelihood ratio (LR) test to compare the fit of the EKD distribution with its sub-models for a given data set. For example, to test θ = 1, the LR statistic is ω = 2 ln L α,λ,δ,φ,θ − ln L α,λ,δ,φ, 1 , whereα,λ,δ,φ andθ are the unrestricted estimates, andα,λ,δ andφ are the restricted estimates. The LR test rejects the null hypothesis if ω > χ 2 d , where χ 2 d denote the upper 100d% point of the χ 2 distribution with 1 degrees of freedom.

Simulation study
In this section, we examine the performance of the EKD distribution by conducting various simulations for different sizes (n=200, 400, 800, 1200) via the subroutine NLP in SAS. We simulate 2000 samples for the true parameters values I : α = 2, λ = 1, δ = 3, φ = 2, θ = 2 and II : α = 1, λ = 1, δ = 1, φ = 1, θ = 1. Table 1 lists the means MLEs of the five model parameters along with the respective root mean squared errors (RMSE). From the results, we can verify that as the sample size n increases, the mean estimates of the parameters tend to be closer to the true parameter values, since RMSEs decay toward zero.

Application: EKD and sub-distributions
In this section, applications based on real data, as well as comparison of the EKD distribution with its sub-models are presented. We provide examples to illustrate the flexibility of the EKD distribution in contrast to other models, including the exponentiated http://www.jsdajournal.com/content/1/1/8 Kumaraswamy-Weibull (EKW), and beta-Kumaraswamy-Weibull (BKW) distributions for data modeling. The pdfs of EKW and BKW distributions are respectively. The first data set consists of the number of successive failures for the air conditioning system of each member in a fleet of 13 Boeing 720 jet airplanes (Proschan 1963). The data is presented in Table 2. The second data set consists of the salaries of 818 professional baseball players for the year 2009 (USA TODAY).
The third data set represents the poverty rate of 533 districts with more than 15,000 students in 2009 (Digest of Education Statistics "http://nces.ed.gov/programs/digest/ d11/tables/dt11_096.asp"). These data sets are modeled by the EKD distribution and compared with the corresponding sub-models, the Kumaraswamy-Dagum and Dagum http://www.jsdajournal.com/content/ distributions, and as well as EKW, BKW distributions. Table 3 gives a descriptive summary of each sample. The air conditioning system sample has far more variability and the baseball player salary sample has the lowest variability. The maximum likelihood estimates (MLEs) of the parameters are computed by maximizing the objective function via the subroutine NLMIXED in SAS. The estimated values of the parameters (standard error in parenthesis), -2 Log-likelihood statistic, Akaike Information Criterion, AIC = 2p − 2 ln(L), Bayesian Information Criterion, BIC = p ln(n) − 2 ln(L), and Consistent Akaike Information Criterion, AICC = AIC + 2 p(p+1) n−p−1 , where L = L(ˆ ) is the value of the likelihood function evaluated at the parameter estimates, n is the number of observations, and p is the number of estimated parameters for the EKD distribution and its sub-distributions are tabulated. See Table 4, Table 5 and  Table 6.
Plots of the fitted EKD, KD, D and the histogram of the data are given in Figure 4. The probability plots (Chambers et al. 1983) consists of plots of the observed probabilities, against the probabilities predicted by the fitted model are also presented in Figure 5. For the EKD distribution, we plotted for example,     was calculated for each plot. The plot with the smallest SS corresponds to the model with points that are closer to the diagonal line. Plots of the empirical and estimated survival functions for the models are also presented in Figure 6. For the air conditioning system data, initial values α = 1, λ = 2, δ = 0.6, φ = 3, θ = 1 are used in SAS code for EKD model. The LR statistics for the test of the hypothesis H 0 : KD against H a : EKD and H 0 : D against H a : EKD are 1.9 (p-value = 0.17) and 13.4 (p-value = 0.0012). Consequently, KD distribution is the best distribution based on the LR statistic. The KD distribution gives smaller SS value than Dagum distribution and slightly bigger than EKD. For the non nested models, the values of AIC and AICC for KD and EKW models are very close, however the BIC value for KD distribution is slightly smaller than the corresponding value for the EKW distribution. We conclude that KD model compares favorably with the EKW distribution and thus provides a good fit for the air conditioning system data.
For the baseball player salary data set, initial values for EKD model in SAS code are α = 70, λ = 0.01, δ = 1.026, φ = 0.1, θ = 1. The EKD distribution is a better fit than KD and Dagum distributions for this data, as well as the other distributions. The values of the statistics AIC, AICC and BIC for KD distribution are smaller compared to the non nested distributions. The LR statistics for the test of the hypotheses H 0 : KD against H a : EKD and H 0 : D against H a : EKD are 93.1 (p-value< 0.0001) and 361.5 http://www.jsdajournal.com/content/1/1/8 Figure 5 Observed probability vs predicted probability for data sets.
(p-value < 0.0001). Consequently, we reject the null hypothesis in favor of the EKD distribution and conclude that the EKD distribution is significantly better than the KD and Dagum distributions based on the LR statistic. The value of AIC, AICC and BIC statistics are lower for the EKD distribution when compared to those for the EKW and BKW distributions.
For poverty rate data, initial values for EKD model are α = 73, λ = 0.1, δ = 0.15, φ = 60, θ = 0.33. The LR statistic for the test of the hypotheses H 0 : KD against H a : EKD and H 0 : D against H a : EKD are 8.2 (p-value = 0.0042) and 81.1 (p-value < 0.0001), respectively. The values of AIC, AICC and BIC statistics shows EKD distributions is a better model and the SS value of EKD model is comparatively smaller than the corresponding values for the KD and D distributions. Consequently, we conclude that EKD distribution is the best fit for the poverty rate data.