The odd generalized exponential family of distributions with applications

We propose a new family of continuous distributions called the odd generalized exponential family, whose hazard rate could be increasing, decreasing, J, reversed-J, bathtub and upside-down bathtub. It includes as a special case the widely known exponentiated-Weibull distribution. We present and discuss three special models in the family. Its density function can be expressed as a mixture of exponentiated densities based on the same baseline distribution. We derive explicit expressions for the ordinary and incomplete moments, quantile and generating functions, Bonferroni and Lorenz curves, Shannon and Rényi entropies and order statistics. For the first time, we obtain the generating function of the Fréchet distribution. Two useful characterizations of the family are also proposed. The parameters of the new family are estimated by the method of maximum likelihood. Its usefulness is illustrated by means of two real lifetime data sets.


Introduction
The art of proposing generalized classes of distributions has attracted theoretical and applied statisticians due to their flexible properties. Most of the generalizations are developed for one or more of the following reasons: a physical or statistical theoretical argument to explain the mechanism of the generated data, an appropriate model that has previously been used successfully, and a model whose empirical fit is good to the data. The last decade is full on new classes of distributions that become precious for applied statisticians. There are two main approaches for adding new shape parameter(s) to a baseline distribution.

Approach 1: Adding one shape parameter
The first approach of generalization was suggested by Marshall and Olkin (1997) by adding one parameter to the survival function G(x) = 1 − G(x), where G(x) is the cumulative distribution function (cdf ) of the baseline distribution. Gupta et al. (1998) added one parameter to the cdf, G(x), of the baseline distribution to define the exponentiated-G ("exp-G" for short) class of distributions based on Lehmann-type alternatives (see Lehmann 1953). Following Gupta et al. 's class, Gupta and Kundu (1999) studied the twoparameter generalized-exponential (GE) distribution as an extension of the exponential distribution based on Lehmann type I alternative. The GE distribution is also known as the exponentiated exponential (EE) distribution. Since it is the most attractive generalization of the exponential distribution, the GE model has received increased attention and many authors have studied its various properties and also proposed comparisons with other distributions. Some significant references are: Gupta and Kundu (2001a;2001b;2003;2004;2007;2008;, Kundu et al. (2005), Nadarajah and Kotz (2006), Dey and Kundu (2009), Pakyari (2010) and . In fact, the GE model has been proven to be a good alternative to the gamma, Weibull and log-normal distributions, all of them with two-parameters. The GE distribution can be used quite effectively for analyzing lifetime data which have monotonic hazard rate function (hrf ) but unfortunately it cannot be used if the hrf is upside-down, J or reversed-J shapes. Differently, the new family can have increasing, decreasing, J, reversed-J, bathtub and upside-down bathtub and, therefore, it can be used effectively for analyzing lifetime data of various types. The generalization of the exp-G class to other distributions is beyond the scope of the paper.
We propose a new class of distributions called the odd generalized exponential ("OGE" for short) family, which is flexible because of the hazard rate shapes: increasing, decreasing, J, reversed-J, bathtub and upsidedown bathtub. In Section 2, we define the OGE family of distributions. Two special cases of this family are presented in Section 3. The density and hazard rate functions are described analytically in Section 4. A useful mixture representation for the pdf of the new family is derived in Section 5. In Section 6, we obtain explicit expressions for the moments, generating function, mean deviations and entropies. In Section 7, we determine a power series for the quantile function (qf ). In Section 8, we investigate the order statistics. Section 9 refers to some characterizations of the OGE family. In Section 10, the parameters of the new family are estimated by the method of maximum likelihood. In Section 11, we illustrate its performance by means of two applications to real life data sets. The paper is concluded in Section 12.

The new family
Consider a parent distribution depending on a parameter vector ξ with cdf G(x; ξ ), survival function G(x; ξ ) = 1 − G(x; ξ ) and probability density function (pdf ) g(x; ξ ). The cdf of the GE model with positive parameters λ and α is given by (x) = 1 − e −λ x α (for x > 0). For α = 1, it reduces to the exponential distribution with mean λ −1 . In fact, the density and hazard functions of the GE and gamma distributions are quite similar. The GE density function is always right-skewed and can be used quite effectively to analyze skewed data. Its hrf can be increasing, decreasing and constant depending on the shape parameter in a similar manner of the gamma distribution. Its applications have been wide-spread as a model to power system equipment, rainfall data, software reliability and analysis of animal behavior, among others.
We define the cdf of the OGE family by replacing x by G(x; ξ )/G(x; ξ ) in the cdf (x) leading to where α > 0 and λ > 0 are two additional parameters.
The pdf corresponding to (1) is given by where g(x; ξ ) is the baseline pdf. We can omit the dependence on the vector of parameters ξ and write simply G(x) = G(x; ξ ). Equation 2 will be most tractable when the cdf G(x) and pdf g(x) have explicit expressions. Hereafter, a random variable X with density function (2) is denoted by X ∼ OGE(α, λ, ξ ). The main motivations for using the OGE family are to make the kurtosis more flexible (compared to the baseline model) and possible to construct heavy-tailed distributions that are not long-tailed for modeling real data. The family (1) contains some special models as those listed in Table 1. In particular, it includes as a special case the widely known exponentiated Weibull (EW) distribution pioneered by Mudholkar and Srivastava (1993) and Mudholkar et al. (1995). For a detail survey on the EW model, the reader is referred to Nadarajah et al. (2013). It has been well-established in the literature that the EW distribution gives significantly better fits than the exponential, gamma, Weibull and lognormal distributions. We offer a physical interpretation of X when α is an integer. Consider a system formed by α independent components following the odd exponential-G class (Bourguignon et al. 2014) given by Suppose the system fails if all α components fail and let X denote the lifetime of the entire system. Then, the cdf of X is F(x; α, λ, ξ ) = H(x; λ, ξ ) α , which is identical to (1). Generalized exponential distribution (Gupta and Kundu 1999) 3 1 - Exponentiated Weibull distribution (Mudholkar and Srivastava 1993) 4 --1− e −a x Generalized Gompertz distribution (El-Gohary et al. 2013) In Table 1, we provide some members of the OGE family. The hrf of X is given by Theorem 1 provides some relations of the OGE family with other distributions.
The OGE family is easily simulated by inverting (1) as follows: if U has a uniform U(0, 1) distribution, then has the density function (2), where Q G (·) = G −1 (·) is the baseline qf.

Special OGE distributions
Lifetime distributions play a fundamental role in survival analysis, biomedical science, engineering and social sciences. Typically, lifetime refers to human life length, the life span of a device before it fails or the survival time of a patient with serious disease from the diagnosis to the death. Here, we present three special OGE distributions that can be useful in applied survival analysis.
The cdf and pdf of the OGE-W distribution are given by (for x > 0) respectively, where α > 0 and β > 0 are shape parameters and λ > 0 and θ > 0 are scale parameters.
The applications of the OGE-W and EW distributions can be directed to model extreme value observations in floods, software reliability, insurance data, tree diameters, carbon fibrous composites, firmware system failure, reliability prediction and fracture toughness, among others.

The OGE-Fréchet (OGE-Fr) distribution
The OGE-Fr distribution is defined from (2) by taking G(x; ξ ) = e −(b/x) a and g(x; ξ ) = a b a x −(a+1) e −(b/x) a to be the cdf and pdf of the Fréchet distribution with parameters a and b, respectively, and ξ = (a, b).
The cdf and pdf of the OGE-Fr distribution are given by (for x > 0) respectively, where α > 0 and a > 0 are shape parameters and λ > 0 and b > 0 are scale parameters.
The OGE-Fr distribution has applications ranging from accelerated life testing, rainfall and floods, queues in supermarkets, sea currents, wind speeds, and track race records, among others.
The cdf and pdf of the OGE-N distribution are given by (for x ∈ ) respectively, where α > 0 is a shape parameter, μ ∈ is a location parameter, and λ > 0 and σ > 0 are scale parameters. In Figures 1, 2, 3, 4 and 5, we display some plots of the pdf and hrf of the OGE-W, OGE-Fr and OGE-N distributions for selected parameter values. Figures 1, 2 and 3 indicate that the OGE family generates distributions with various shapes such as symmetrical, left-skewed, right-skewed and reversed-J. Further, Figures 4 and 5 reveal that this family can produce hazard rate shapes such as increasing, decreasing, J, reversed-J, bathtub and upside-down bathtub. This fact implies that the OGE family can be very useful to fit data sets under various shapes.
Proposition 2. The asymptotics of equations (1), (2) and (3) as x → ∞ are given by  The shapes of the density and hazard rate functions can be described analytically. The critical points of the OGE-G density functions are the roots of the equation:  There may be more than one root to (8).
If x = x 0 is a root of (8) then it corresponds to a local maximum if λ(x) > 0 for all x < x 0 and λ(x) < 0 for all x > x 0 . It corresponds to a local minimum if λ(x) < 0 for all x < x 0 and λ(x) > 0 for all x > x 0 . It refers to a point of inflexion if either λ(x) > 0 for all The critical point of h(x) are obtained from the equation There may be more than one root to (9). Let τ (

Useful expansions
First, we consider the quantity G(x;ξ ) . By expanding the exponential function in power series and using the generalized binomial expansion, we obtain The quantity M can be expressed as where, for k ≥ 0, b k = a k+1 , and, for k ≥ 1, Second, we can expand z α in Taylor series as where Here and from now on, we use the notation (α) l = α(α − 1) . . . (α − l + 1) for the descending factorial with (α) 0 = 1. We use throughout the paper a result of Gradshteyn and Ryzhik (2000, Section 0.314) for a power series raised to a positive integer n (for n ≥ 1) where the coefficients c n,i (for n = 1, 2, . . .) can be determined from the recurrence equation , and c n,0 = b n 0 . Third, based on Eqs. (10)-(12), the cdf (1) can be expressed as where c m,k can be obtained from the quantities b 0 , . . . , b k as in equation (12).
Then, we have where Finally, the density function of X can be expressed in the mixture form where h α+k (x) = (α + k) G(x; ξ ) α+k−1 g(x; ξ ) is the exp-G density function with power parameter α + k. Equation (15) reveals that the OGE density function is a linear combination of exp-G densities. Thus, some mathematical properties of the new model such as the ordinary and incomplete moments, and moment generating function (mgf ) can be derived from those properties of the exp-G distribution. Some exp-G structural properties are studied by Mudholkar and Srivastava (1993), Mudholkar et al. (1995), Mudholkar and Hutson (1996), Gupta et al. (1998), Gupta and Kundu (2001a;2001b), Nadarajah and Kotz (2006) and  6 Mathematical properties

Moments
The need for necessity and the importance of moments in Statistics especially in applications is obvious. Some of the most important features and characteristics of a distribution can be studied through moments. Let Y k be a random variable having the exp-G pdf h α+k (x) with power parameter α + k. A first formula for the nth moment of X follows from (15) as Expressions for moments of several exp-G distributions are given by Nadarajah and Kotz (2006), which can be used to obtain μ n . We provide two applications of (16) for the OGE-W and OGE-Fr models discussed in Section 3. First, the moments of the OGE-W model can be derived from closed-forms moments of the EW distribution given by Choudhury (2005). We obtain where the quantities A j,k (for j, k ≥ 0) are given by The double infinite series on the right-hand side converges for all parameter values. Second, the moments of the OGE-Fr model (for n < a) follow immediately from the moments of the Fréchet distribution as A second alternative formula for μ n is obtained from (16) in terms of the baseline qf as where (for a > 0) τ (n, a) = 1 0 Q G (u) n u a du. Cordeiro and Nadarajah (2011) obtained τ (n, a) for some well-known distribution such as the normal, beta, gamma and Weibull distributions, which can be applied to obtain raw moments of the corresponding OGE distributions.
For empirical purposes, the shapes of many distributions can be usefully described by what we call the first incomplete moment which plays an important role for measuring inequality, for example, income quantiles and Lorenz and Bonferroni curves.
The nth incomplete moment of X can be expressed as The last integral can be computed for most G distributions.

Generating function
The generating function provides the basis of an alternative route to analytical results compared with working directly with the pdf and cdf and it is widely used in the characterization of the distributions, and the application of goodness-of-fit tests. Let M(t) = E(e t X ) be the mgf of X. Then, we can write from (15) where M k (t) is the mgf of Y k . Hence, M(t) can be determined from the exp-G generating function. Now, we provide two applications of (20) for the OGE-W and OGE-Fr models. Nadarajah et al. (2015) derived a formula for the mgf of Y k which holds for the OGE-W model depending on the complex parameter Wright generalized hypergeometric function with p numerator and q denominator parameters (Kilbas et al. 2006, Equation (1.9)) given by Nadarajah et al. (2015) demonstrated that the mgf of Y k can be expressed as Combining (20) and the last equation gives the mgf of the OGE-W distribution. The mgf of the Fréchet distribution can be obtained by The calculations of the integral in (21) involve the generalized hypergeometric function defined in equation (2.3.1.14) (Prudnikov et al. 1986, p. 322).
For a > 0 and s > 0, if b = p/q (with p ≥ 1 and q ≥ 1 are co-primes integers), we can write where z = (−1) p+q s p a q / (p p q q ), (k, a) denote the sequence (k, a) = (a/k, (a + 1)/k, . . . , (a + k − 1)/k), m F n denotes the generalized hypergeometric function defined by and (c) k = c(c + 1) . . . (c + k − 1) denotes the ascending factorial. Numerical routines for computing the generalized hypergeometric function are available in most mathematical packages, e.g., Maple, Mathematica and Matlab. Nadarajah and Kotz (2005), and Nadarajah (2007) used also this result to obtain the properties of the distribution of the difference between two independent Gumbel variates, and to Iacbellis and Fiorentino (2000), and Fiorentino et al. (2006)'s model for peak streamflow.
By using (22), the mgf of the Fréchet distribution follows as Using (15), the mgf of the OGE-Fr model can be expressed as and then By using (23), M(t) reduces to A second alternative formula for M(t) can be derived from (15) as where (for We can obtain the mgf 's of several OGE distributions directly from equation (24). Equations (20) and (24) are the main results of this section.

Mean deviations
The mean deviations about the mean (δ 1 = E |X − μ 1 | ) and about the median (δ 2 = E(|X − M|)) of X can be expressed as respectively, where μ 1 = E(X), M = Median(X) = Q(0.5) is the median which comes from (4), F μ 1 is easily calculated from the cdf (1) and m 1 (z) is the first incomplete moment determined from (19) with n = 1. Now, we provide two alternative ways to compute δ 1 and δ 2 . A general formula for m 1 (z) can be derived from equation (15) as where Hence, the mean deviations in (25) depend only on quantity J k (z).
We give two applications of J k (z) for the the OGE-W and OGE-Fr distributions. For the first model, the EW pdf (for x > 0) with power parameter α + k is given by The last integral is equal to the incomplete gamma function and then the mean deviation for the OGE-W distribution is given by where γ (a, z) = z 0 w a−1 e −w dw is the incomplete gamma function. The first incomplete moment of the Fréchet distribution is given by So, the first incomplete moment of the OGE-Fr model can be derived from the last equation and (15) as where (a, z) = ∞ z w a−1 e −w dw is the complementary incomplete gamma function. A second general formula for m 1 (z) can be derived by setting u = G(x) in (15) where Applications of equations (26) and (27) can be directed to the Bonferroni and Lorenz curves defined for a given probability π by B(π) = m 1 (q)/(πμ 1 ) and L(π) = m 1 (q)/μ 1 , respectively, where μ 1 = E(X) and q = Q(π) is the qf of X at π.

Entropies
An entropy is a measure of variation or uncertainty of a random variable X. Two popular entropy measures are the Rényi and Shannon entropies (Rényi 1961;Shannon 1948). The Rényi entropy of a random variable with pdf f (x) is defined by for γ > 0 and γ = 1. The Shannon entropy of a random variable X is defined by E − log f (X) . It is the special case of the Rényi entropy when γ ↑ 1. Direct calculation yields Using Taylor and generalized binomial expansions, and after some algebraic manipulations, we obtain So, we can write: Proposition 3. Let X be a random variable with pdf (2). Then, The simplest formula for the entropy of X can be reduced to After some algebraic developments, we obtain an expression for I R (γ ) where Y j,k has a beta distribution with parameters j + k + 1 and one, and

Quantile power series
Power series methods are at the heart of many aspects of applied mathematics and statistics. Here, we derive a power series for the qf x = Q(u) = F −1 (u) of X by expanding (4). If the baseline qf Q G (u) does not have a closed-form expression, it can usually be expressed in terms of a power series as where the coefficients a i 's are suitably chosen real numbers which depend on the parameters of the G distribution. For several important distributions, such as the normal, Student t, gamma and beta distributions, Q G (u) does not have explicit expressions but it can be expanded as in equation (29). As a simple example, for the normal N(0, 1) distribution, a i = 0 for i = 0, 2, 4, . . . and a 1 = 1, a 3 = 1/6, a 5 = 7/120 and a 7 = 127/7560, . . . Next, we derive an expansion for the argument of Q G (·) in equation (4) .
For |u| < 1, we have − log(1 − u) = ∞ i=0 u i+1 /(i + 1). Then, using the generalized binomial expansion, we obtain Using the generalized binomial expansion again, we can write where b * 0 = λ + α * 0 , b * k = α * k for k ≥ 1 and the coefficient δ k (for k ≥ 0) can be determined from the recurrence equation (for k ≥ 1 and δ 0 = α * 0 /b * 0 ) Then, the qf of X can be expressed from (4) as Combining (29) and (30), we obtain Let W (·) be any integrable function in the positive real line. We can write Equations (31) and (32) are the main results of this section since from them we can obtain various OGE mathematical quantities. In fact, various of these properties can follow by using the second integral for special W (·) functions, which are usually more simple than if they are based on the first integral. Established algebraic expansions to determine mathematical quantities of the OGE distributions based on these equations can be more efficient then using numerical integration of the pdf (2), which can be prone to rounding off errors among others. For the great majority of these quantities, we can adopt twenty terms in the power series (31).
The effects of the shape parameter α on the skewness and kurtosis of X can be studied using quantile measures. The shortcomings of the classical kurtosis measure are wellknown. The Bowley skewness is given by Since only the middle two quartiles are considered and the outer two quartiles are ignored, this adds robustness to the measure. The Moors kurtosis is based on octiles These measures are less sensitive to outliers and they exist even for distributions without moments. In Figures 6 and 7, we plot the measures B and M for the OGE-W and OGE-Fr distributions. These plots indicate how both measures B and M vary on the parameter α.

Order statistics
Order statistics make their appearance in many areas of statistical theory and practice. Suppose X 1 , . . . , X n is a random sample from the OGE-G distribution. Let X i:n denote the ith order statistic. Using (14) and (15), the pdf of X i:n can be expressed as Figure 6 Skewness (a) and kurtosis (b) plots for the OGE-W distribution.

Figure 7 Skewness (c) and kurtosis (d) plots for the OGE-Fr distribution.
Hence, where Equation (33) is the main result of this section. It reveals that the pdf of the OGE order statistics is a triple linear combination of exp-G distributions. So, several mathematical quantities of these order statistics like ordinary and incomplete moments, factorial moments, mgf, mean deviations and several others come from these quantities of the OGE distributions.

Characterization of the OGE family
Characterizations of distributions are important to many researchers in the applied fields. An investigator will be vitally interested to know if their model fits the requirements of a particular distribution. To this end, one will depend on the characterizations of this distribution which provide conditions under which the underlying distribution is indeed that particular distribution. Various characterizations of distributions have been established in many different directions. In this section, two main characterizations of the OGE family are presented. These characterizations are based on: (i) a simple relationship between two truncated moments; (ii) a single function of the random variable.

Characterizations based on truncated moments
In this subsection we present characterizations of the OGE family in terms of a simple relationship between two truncated moments. Our characterization results presented here will employ an interesting result due to Glänzel (1987) (Theorem 2, below). The advantage of these characterizations is that, the cdf F(x) does not require to have a closedform and are given in terms of an integral whose integrand depends on the solution of a first order differential equation, which can serve as a bridge between probability and differential equation.
Theorem 2. Let ( , , P) be a given probability space and let H = [a, b] be an interval for some a < b a = −∞, b = ∞ might as well be allowed . Let X : → H be a continuous random variable with the distribution function F and let q 1 and q 2 be two real functions defined on H such that E q 1 (X) |X ≥ x = E q 2 (X) |X ≥ x η (x) , x ∈ H, is defined with some real function η. Consider that q 1 , q 2 ∈ C 1 (H), η ∈ C 2 (H) and F(x) is twice continuously differentiable and strictly monotone function on the set H. Finally, assume that the equation q 2 η = q 1 has no real solution in the interior of H. Then, F is uniquely determined by the functions q 1 , q 2 and η, particularly where the function s is a solution of the differential equation s = η q 2 ηq 2 −q 1 and C is a constant chosen to make H dF = 1.
Clearly, Theorem 2 can be stated in terms of two functions q 1 and η by taking q 2 (x) ≡ 1, which will reduce the condition given in Theorem 2 to E q 1 (X) |X ≥ x = η (x). However, adding an extra function will give a lot more flexibility, as far as its application is concerned.
Proof. Let X has density (2), then Finally, Conversely, if η is given as above, then Now, in view of Theorem 2, X has density (2).

Corollary 1. Let X :
→ R be a continuous random variable and let q 2 (x) be as in Proposition 4. The pdf of X is (2) if and only if there exist functions q 1 and η defined in Theorem 2 satisfying the differential equation , x > 0.

Remark 1. (a)
The general solution of the differential equation in Corollary 1 is where D is a constant. One set of appropriate functions is given in Proposition 4 with D = 1 2 .

(b)
Clearly there are other triplets of functions (q 2 , q 1 , η) satisfying the conditions of Theorem 2. We presented one such triplet in Proposition 4.

Characterizations based on single function of the random variable
In this subsection we employ a single function ψ of X and state characterization results in terms of ψ (X).

Estimation
Inference can be carried out in three different ways: point estimation, interval estimation and hypothesis testing. Several approaches for parameter point estimation were proposed in the literature but the maximum likelihood method is the most commonly employed. The maximum likelihood estimates (MLEs) enjoy desirable properties and can be used when constructing confidence intervals and also in test-statistics. Large sample theory for these estimates delivers simple approximations that work well in finite samples. Statisticians often seek to approximate quantities such as the density of a test-statistic that depend on the sample size in order to obtain better approximate distributions. The resulting approximation for the MLEs in distribution theory is easily handled either analytically or numerically.
Here, consider the estimation of the unknown parameters of the OGE family by the method of maximum likelihood from complete samples only. Let x 1 , . . . , x n be observed values from this family with parameters α, λ and ξ . Let = (α, λ, ξ ) be the (r × 1) parameter vector. The total log-likelihood function for is given by where The components of the score function U n ( ) = U α , U λ , U ξ k are where v (ξ ) (·) means the derivative of the function v with respect to ξ . Setting these equations to zero and solving them simultaneously yields the MLEs = ( α, λ, ξ ) of = (α, λ, ξ ) . These equations cannot be solved analytically, and analytical softwares are required to solve them numerically.

Simulation study
We evaluate the performance of the maximum likelihood method for estimating the OGE-W parameters using Monte Carlo simulation for a total of twenty four parameter combinations and the process is repeated 200 times. Two different sample sizes n = 100 and 300 are considered. The MLEs and their standard deviations of the parameters are listed in Table 2. The estimates of α, β and λ are determined by solving the nonlinear equations U( ) = 0. Based on Table 2, we note that the ML method performs well for estimating the model parameters. Also, as the sample size increases, the biases and the standard deviations of the MLEs decrease as expected. In this section, we provide two applications to real data to illustrate the importance of the OGE family by means of the OGE-W, OGE-Fr and OGE-N models presented in Section 3. We consider θ = 1, b = 1 and μ = 0 for the OGE-W, OGE-Fr and OGE-N models, respectively, due to the fact that one scale parameter is enough for fitting these univariate models. The MLEs of the parameters for the these models are calculated and four goodness-of-fit statistics are used to compare the new family with its sub-models. The first real data set represents the survival times of 121 patients with breast cancer obtained from a large hospital in a period from 1929 to 1938 (Lee 1992). The data examined by Ramos et al. (2013) are: 0.3, 0.3, 4.0, 5.0, 5.6, 6.2, 6.3, 6.6, 6.8, 7.4, 7.5, 8.4, 8.4, 10.3, 11.0, 11.8, 12.2, 12.3, 13.5, 14.4, 14.4, 14.8, 15.5, 15.7, 16.2, 16.3, 16.5, 16.8, 17.2, 17.3, 17.5, 17.9, 19.8, 20.4, 20.9, 21.0, 21.0, 21.1, 23.0, 23.4, 23.6, 24.0, 24.0, 27.9, 28.2, 29.1, 30. The second data set consists of 63 observations of the strengths of 1.5 cm glass fibres, originally obtained by workers at the UK National Physical Laboratory. Unfortunately, the units of measurement are not given in the paper.  Smith and Naylor (1987).
The MLEs are computed using the Limited-Memory Quasi-Newton Code for Bound-Constrained Optimization (L-BFGS-B) and the log-likelihood function evaluated at the MLEs (ˆ ). The measures of goodness of fit including the Akaike information criterion (AIC), Anderson-Darling (A * ), Cramér-von Mises (W * ) and Kolmogrov-Smirnov (K-S) statistics are computed to compare the fitted models. The statistics A * and W * are described in details in Chen and Balakrishnan (1995). In general, the smaller the values of these statistics, the better the fit to the data. The required computations are carried out in the R-language.
Tables 3 and 4 list the MLEs and their corresponding standard errors (in parentheses) of the model parameters. The numerical values of the model selection statisticsˆ , AIC, A * , W * and K-S are listed in Tables 5 and 6. We note from the figures in Table 5 that both    Table 6 The statisticsˆ , AIC, A * , W * , K-S and K-S p-value for the second data set OGE-W and OGE-Fr models have the lowest values of the AIC, A * , W * and K-S statistics (for the first data set) as compared to their sub-models, suggesting that the OGE-W and OGE-Fr models provide the best fits. The histogram of the first data and the estimated pdfs and cdfs of the OGE-W and OGE-Fr models and their sub-models are displayed in Figure 8. Similarly, it is also evident from Table 6 that the OGE-N gives the lowest values for the AIC, A * and W * statistics while OGE-W models gives the lowest values for the AIC and K-S statistics (for the second data set) as compared to their sub-models,  and therefore these models can be chosen as the best ones. The histogram of the second data and estimated pdfs and cdfs of the OGE-N and OGE-W distributions and their submodels are displayed in Figure 9. It is clear from Tables 5 and 6, and Figures 8 and 9 that the OGE-W, OGE-Fr and OGE-N models provide better fits than their sub-models to these two data sets.

Concluding remarks
The generalized continuous distributions have been widely studied in the literature. We propose a new class of distributions called the odd generalized exponential family. We study some structural properties of the new family including an expansion for its density function. We obtain explicit expressions for the moments, generating function, mean deviations, quantile function and order statistics. The maximum likelihood method is employed to estimate the family parameters. We fit three special models of the proposed family to two real data sets to demonstrate the usefulness of the family. We use four goodness-of-fit statistics in order to verify which distribution provides better fit to these data. We conclude that these three special models provide consistently better fits than other competing models. We hope that the proposed family and its generated models will attract wider applications in several areas such as reliability engineering, insurance, hydrology, economics and survival analysis.