The new Burr distribution and its application

This paper derives a new family of Burr-type distributions as new Burr distribution. This particular skewed distribution that can be used quite effectively in analyzing lifetime data. It is observed that the new distribution has modified unimodal hazard function. Various properties of the new Burr distribution, such that moments, quantile functions, hazard function, and Shannon’s entropy are obtained. The exact form of the probability density function and moments of ith\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i{\mathrm{th}}$$\end{document}-order statistics in a sample of size n from new Burr distribution are derived. Estimation of parameters and change-point of hazard function by the maximum likelihood method are discussed. Change-point of hazard function is usually of great interest in medical or industrial applications. The flexibility of the new model is illustrated with an application to a real data set. In addition, a goodness-of-fit test statistic based on the Rényi Kullback–Leibler information is used.


Introduction
Burr [2] developed the system of Burr distributions. The Burr system of distributions includes 12 types of cumulative distribution functions which yield a variety of density shapes. The attractiveness of this relatively unknown family of distributions for model fitting is that it combines a simple mathematical expression for cumulative frequency function with coverage in the skewness-kurtosis plane. Many standard theoretical distributions, including the Weibull, exponential, logistic, generalized logistic, Gompertz, normal, extreme value, and uniform distributions, are special cases or limiting cases of the Burr system of distributions (see [11]). Family of Burr-type distributions is a very popular distribution family for modelling lifetime data and for modelling phenomenon with monotone and unimodal failure rates (see, for example, [13,18]).
Analogous to the Pearson system of distributions, the Burr distributions are solutions to a differential equation, which has the form: where y equal to F(x) and g(x, y) must be positive for y in the unit interval and x in the support of F(x). Different functional forms of g(x, y) result in different solutions F(x), which define the families of the Burr system. For example, Burr II distribution is obtained when gðx; yÞ ¼ gðxÞ ¼ ke Àx ð1þe Àx Þ kÀ1 ð1þe Àx Þ k À1 . In this paper, we derive a new distribution of Burr-type distributions which is more flexible by replacing g(x, y) with gðxÞ ¼ 3px 2 e Àx 3 ð1þe Àx 3 Þ pÀ1 ð1þe Àx 3 Þ p À1 , (p [ 0). We refer to this new distribution as the new Burr distribution. If g(x, y) is taken to be g(x), then the solution of the differential Eq.
The shapes of density and hazard functions of the new Burr distribution for different values of shape parameter p are illustrated in Fig. 1. New Burr distribution has unimodal and bimodal pdfs. None of the 12 types of Burr distributions has this feature. Data that exhibit bimodal behavior arises in many different disciplines. In medicine, urine mercury excretion has two peaks, see, for example, [5]. In material characterization, a study conducted by [4], grain size distribution data reveals a bimodal structure. In meteorology, [19] indicated that, water vapor in tropics, commonly have bimodal distributions. To see more applications of bimodal distributions, see [7][8][9]16].
The reminder of the paper is organized as follows: properties of the new Burr distribution, such that moments, quantile functions, hazard function, Shannon's entropy, and distribution of its order statistics are discussed in Sects. 2, 3, and 4. In Sect. 5, estimation of parameters and change-point of hazard function by the maximum likelihood method are discussed, and in Sect. 6, we establish a goodness-of-fit test statistic based on the Rényi Kullback-Leibler information for testing new Burr model. Finally, in Sect. 7, we present an illustrative example. Section 8 provides conclusions.

Properties of the new Burr distribution
New Burr distribution has unimodal and bimodal pdfs. The modes of distribution are provided by differentiating the density of new Burr distribution in 1.6 with respect to x: ð2:1Þ The derivative f 0 ðx;l;r;pÞ exists every where, hence critical point(s) satisfy equation f 0 ðx;l;r;pÞ ¼ 0. In 2.1, set l ¼ 0 and r ¼ 1, because location and scale parameters will not affect the distribution shape. Thus, equation f 0 ðx;l;r;pÞ ¼ 0 simplifies to 3x 3 ð1 À pe Àx 3 Þ À 2ð1 þ e Àx 3 Þ ¼ 0: ð2:2Þ Analytical solution of 2.2 is not possible. Numerical approximation of modes using the midpoint method is applied to study the modes. The distance between the two p=0.01,μ=4,σ=1/3 p=1/2,μ=4,σ=1/3 p=1, μ=4,σ=1/3 p=2, μ=4,σ=1/3 p=10,μ=4,σ=1/3  Table 1. From Table 1, it is observed that when p increases, the distance between two modes decreases, and for 0\p\1, when p decreases, value of pdf in the second mode decreases to zero and pdf will be almost unimodal, and for p ¼ 1, values of pdf in two modes are the same but for p [ 1, and when p increases, value of pdf in the first mode decreases to zero and pdf will be almost unimodal. Hence, the new Burr distribution can be used to analyse different kinds of lifetime data sets with unimodal and bimodal shapes of pdf.
The new Burr distribution has modified unimodal (unimodal followed by increasing) hazard function, and when p increases, hazard function will be almost increasing.
The main purpose in this paper is to describe and fit the data sets with non-monotonic hazard function, such as the bathtub, unimodal and modified unimodal hazard function. Many modifications of important lifetime distributions have achieved the above purpose, but unfortunately, the number of parameters has increased, the forms of survival and hazard functions have been complicated, and estimation problems have risen. More over some of the modifications do not have a closed form for their cdfs. However, this new distribution with one parameter and simple form of cdf achieves this purpose. Now, we discuss the reverse hazard function of the new Burr distribution. The reverse hazard function of any distribution function F(x) can be defined as rðxÞ ¼ f ðxÞ FðxÞ . Consequently, the reversed hazard function of new Burr distribution with zero location parameter and unit scale parameter is given by The reversed hazard function has recently attracted considerable interest of researchers (see, for example, [1,3]). In a reliability setting, the reversed hazard function (multiplied by dx) defines the conditional probability of a failure of an object in ðx À dx; x given that the failure had occurred in [0, x]. The reversed hazard function of new Burr distribution with zero location parameter and unit scale parameter is a linear function of p.
The rth moment about origin of the new Burr distribution is given by using the change of variable, t ¼ 1 1þe Àð xÀl r Þ 3 , 0\t\1, we obtain Now, using 1 t À 1 ¼ e u , 0\u\1, we obtain where E q ð:Þ denotes expectation for X $ q and q is the standard exponential distribution and Using the importance sampling method, the importance sampling estimate of l r is given bŷ ð2:3Þ Using n ¼ 1000, the importance sampling estimate of mean and variance of the new Burr distribution as l ¼ 0 and r ¼ 1 for different values of p is demonstrated in Table 2. From Table 2, it is observed that when p increases, mean increases and variance decreases. Mean and variancê l r q are given by To form a confidence interval for l r , we need to estimate varðÊ q ðgðXÞÞÞ. Because X k are sampled from q, the natural variance estimate is Then, an approximate 99% confidence interval for l r iŝ l r q AE 2:58v The quantile function, Q(u), 0\u\1, for the new Burr distribution can be computed using the formula: The median of a new Burr distribution occurs at rðÀ lnðð 1 2 Þ À 1 p À 1ÞÞ 1 3 þ l, and clearly, it is a decreasing function of p as p 1 but an increasing function of p as p ! 1.
Skewness and kurtosis of a parametric distribution are often measured by a 3 ¼ l 3 r 3 and a 4 ¼ l 4 r 4 , respectively. When the third or fourth moment does not exist, for example, Cauchy, Lévy, and Pareto distributions, a 3 and a 4 , cannot be computed. For the new Burr distribution, skewness and kurtosis can be approximated by approximations of l 3 and l 4 or alternative measures for skewness and kurtosis, based on quantile functions. The measure of skewness S defined by [6] and the measure of kurtosis K defined by [12] are based on quantile functions and they are defined as To investigate the effect of the shape parameter p on the new Burr density function, Eqs. 2.4 and 2.5 are used to obtain Galton's skewness and Moors' kurtosis. Figure 2 displays the Galton's skewness and Moors' kurtosis for the new Burr distribution in terms of the parameter p when l ¼ 0 and r ¼ 1.

Shannon's entropy
The entropy of a random variable X is a measure of variation of uncertainty. Shannon's entropy [17] for a random variable X with pdf f(x) is defined as EðÀ logðf ðxÞÞÞ. In recent years, Shannon's entropy has been used in many applications in fields of engineering, physics, and economics.
Denote by H sh ðXÞ the well-known Shannon's entropy. The following theorem gives the Shannon's entropy of the new Burr distribution.

ð3:4Þ
In the same way, by calculating Eðð1 þ e Àð XÀl r Þ 3 Þ r Þ and then differentiating with respect to r at r ¼ 0, we obtain By replacing 3.2, 3.4, and 3.5 in relation 3.1, the proof is completed. h

Distribution of order statistics
The pdf of X i:n ði ¼ 1; . . .; nÞ is given by f i:n ðx; l; r; pÞ ¼ n! ði À 1Þ!ðn À iÞ! f ðx; l; r; pÞF iÀ1 Â ðx; l; r; pÞð1 À Fðx; l; r; pÞÞ nÀi ; where f ðx; l; r; pÞ and Fðx; l; r; pÞ are pdf and cdf given in 1.5 and 1.6, respectively: where X has new Burr distribution with parameters l, r and pði þ jÞ and Y has q distribution, standard exponential distribution, and gðyÞ ¼ ðe y þ 1Þ ÀpðiþjÞÀ1 y k 3 e 2y . Then, the importance sampling estimate of the rth moment about origin of X i:n is given by The cdf of X i:n ð1 i nÞ is given by where I x ða; bÞ is lower incomplete gamma function. There, the 100uth percentile of X i:n can be obtained by solving F i:n ðxÞ ¼ u: ð4:3Þ The percentage points of X i:n can be evaluated from 4.3 using tables of incomplete beta function (see [15]). However, for i ¼ 1, Eq. 4.3 reduces to ð1 À ð1 þ e Àð xÀl r Þ 3 Þ Àp Þ n ¼ 1 À u. Thus, the 100u-percentage point of the smallest order statistic X 1:n is given by F À1 1:n ðu; p; l; rÞ ¼ l þ r À ln 1 À ð1 À uÞ Similarly, for i ¼ n, the 100u-percentage point of the largest order statistic is F À1 n:n ðu; p; l; rÞ ¼ l þ r À ln u À 1 np À 1 1 3 :

Hazard change-point estimation-classical approach
Hazard function plays an important role in reliability and survival analysis. New Burr distribution has modified unimodal (unimodal followed by increasing) hazard function. In some medical situations, for example, breast cancer, the hazard rate of death of breast cancer patients represents a modified unimodal shape. A modified unimodal shape has three phases: first increasing, then decreasing, and then again increasing. It can be interpreted as a description of three groups of patients, first group is represented by the first phase that contains the weak patients, so the hazard rate of this group is increasing, while the second phase represents the group with strong patients, their bodies have became familiar with the disease and they are getting better. The hazard rate of death of these patients is decreasing. In the third phase, they become weaker and their ability to cope with the disease declines, then the hazard rate of death increases.
For situations, where the hazard function is modified unimodal shaped, usually, we have interest in the estimation of lifetime change-point, that is , the point at which the hazard function reaches to a maximum (minimum) and then decreases (increase). In reliability, the change-point of a hazard function is useful in assessing the hazard in the useful life phase. One of change-points of hazard function of the new Burr distribution is location parameter. In this section, we consider maximum likelihood estimation procedure for change-points of the hazard function.
Let us assume that x 1 ; . . .; x n is a random sample of size n of lifetimes generated by a new Burr distribution with parameters l, r, and p. The log-likelihood function is given by The maximum likelihood estimates for l, r, and p denoted byl,r, andp, respectively, are obtained solving the likelihood equations, ( ol ol ¼ 0, ol or ¼ 0, and ol op ¼ 0Þ. According to the above, maximum likelihood estimator of one of change-points isl.
From the invariance property of maximum likelihood estimators, we can obtain maximum likelihood estimators for functions of l, r and p. For / ¼ gðl; r; pÞ, a one-to-one function of l, r, and p, and we have/ ¼ gðl;r;pÞ. Taking

Test statistics
Suppose that we are interested in a goodness-of-fit test for H 1 : f ðxÞ 6 ¼ f 0 ðx; l; r; pÞ; where l, r, and p are unknown. We will denote the complete samples as X 1:n \X 2:n \ Á Á Á \X n:n . For a null pdf f 0 ðxÞ, the Rényi Kullback-Leibler information from complete data is defined as ðf X1:n;...;Xn:n ðx 1:n ; . . .; x n:n ÞÞ a ðf 0 X1:n;...;Xn:n ðx 1:n ; . . .; x n:n ÞÞ aÀ1 dx 1 Á Á Á dx n ; where r [ 0 and r 6 ¼ 1. Because D r ðf ; f 0 Þ has the property that D r ðf ; f 0 Þ ! 0, and the equality holds if and only if f ¼ f 0 , the estimate of the Rényi Kullback-Leibler information can be consider as a goodness-of-fit test statistic. For that purpose, the Rényi Kullback-Leibler information can be estimated by D r ðf ; f 0 Þ ¼ ÀH r ðX 1:n ; . . .; X n:n Þ À X n j¼1 f 0 ðx j Þ: Thus, the test statistics based on D r ðf ;f 0 Þ n is given by wherel,r, andp are MLEs of l, r, and p, respectively, andĤ r ðX 1:n ; . . .; X n:n Þ is an estimate of Rényi entropy for sample X 1:n \X 2:n \ Á Á Á \X n:n . Under the null hypothesis, T r for r close to 1 will be close to 0, and therefore, large values of T r will lead to the rejection of H 0 .
In this paper, we use estimation of Rényi entropy based on generalized nearest-neighbor graphs that is introduced by [14]. The basic tool to define their estimator was the generalized nearest-neighbor graph. This graph on vertex set V is a directed graph on V. The edge set of it contains for each i 2 S (S is a finite non-empty set of positive integers), an edge from each x 2 V to its i th nearest neighbor according to the Euclidean distance to x.
For p ! 0 denote by L p ðVÞ, the sum of the p th powers of Euclidean lengths of its edges. According to proven theorem in [14] lim n!1 L p ðX 1:n ; . . .; X n:n Þ n 1À p d ¼ c [ 0 a:s:; where p ¼ dð1 À rÞ and d is dimension of sample members. Based on described graph, they estimated Rényi entropy bŷ H r ðX 1:n ; . . .; X n:n Þ ¼ 1 1 À r log L p ðX 1:n ; . . .; X n:n Þ cn 1À p d :

Application
In this section, we consider an uncensored data set corresponding to remission times (in months) of a random sample of 128 bladder cancer patients. These data were previously reported in [10]. TTT plot for considered data is concave then convex indicating an increasing then decreasing hazard function, and is properly accommodated by new Burr distribution. Because in the system of Burr distributions, only Burr X and Burr XII distributions have unimodal hazard functions, and because of the similarity of cdf of the new Burr distribution with the Burr II distribution compared to the rest of distributions in Burr family, compare the fits of the new Burr distribution and those of Burr X, Burr XII, and Burr II and generalized Burr II. Plot of the estimated cdfs of models fitted to the data set is given in Fig. 3. Figure 3 and also values of defined test statistics in the previous section that are shown in Table 3 confirm that the new Burr distribution provides a significantly better fit than Burr X, Burr XII, Burr II, and generalized Burr II distributions. The required numerical evaluations are implemented using Matlab (version 2013) and R software (version 3.3.1).

Conclusions
We introduced a new family of Burr-type distributions as new Burr distribution. Various properties of the distribution are investigated. The distribution is found to be unimodal and bimodal. This new distribution with one parameter and simple form of cdf has modified unimodal (unimodal followed by increasing) hazard function. Hence, this new distribution can be used quite effectively in analyzing lifetime data with non-monotonic hazard function.  Fig. 3 cdfs of the new Burr, Burr X, Burr XII, Burr II, and generalized Burr II models for the remission times of bladder cancer data Table 3 Values of test statistics for the remission times of bladder cancer data