Weibull Generalized Exponential Distribution
Abdelfattah Mustafa, Beih S. El-Desouky and Shamsan AL-Garash

arXiv:1606.07378v1 [] 2 Jun 2016

Department of Mathematics, Faculty of Science, Mansoura University, Mansoura 35516, Egypt.

Abstract
This paper introduces a new three-parameters model called the Weibull-G exponential
distribution (WGED) distribution which exhibits bathtub-shaped hazard rate. Some of it’s
statistical properties are obtained including quantile, moments, generating functions, reliability and order statistics. The method of maximum likelihood is used for estimating the
model parameters and the observed Fisher’s information matrix is derived. We illustrate
the usefulness of the proposed model by applications to real data.
Keywords: Weibull-G class; exponential distribution; generalized exponential; maximum likelihood
estimation.

1

Introduction

The exponential distribution (ED), [3, 12], has a wide range of applications including life testing experiments, reliability analysis, applied statistics and clinical studies. This distribution is a special case of
the two parameter Weibull distribution with the shape parameter equal to 1. The origin and other aspects
of this distribution can be found in [3, 4, 5, 6]. A random variable X is said to have the exponential
distribution (ED) with parameters λ > 0 if it’s probability density function (pdf) is given by
g(x) = λe−λx ,

x > 0,

(1.1)

while the cumulative distribution function (cdf) is given by
G(x) = 1 − e−λx ,

x > 0.

(1.2)

The survival function is given by the equation
S(x) = 1 − G(x) = e−λx ,

x > 0,

(1.3)

and the hazard function is
h(x) = λ

(1.4)

Weibull distribution introduced by [18] is a popular distribution for modeling phenomenon with monotonic failure rates. But this distribution does not provide a good fit to data sets with bathtub shaped or
upside-down bathtub shaped (unimodal) failure rates, often encountered in reliability, engineering and

1

2

Abdelfattah Mustafa, B. S. El-Desouky and Shamsan AL-Garash

biological studies. Hence a number of new distributions modeling the data in a better way have been
constructed in literature as ramifications of Weibull distribution. Bourguignon et al. [8] introduced and
studied generality a family of univariate distributions with two additional parameters, similarly as the
extended Weibull, Gurvich et al. [15] and gamma families, Zografos and Balakrshnan [17], using the
G(x)
Weibull generator applied to the odds ratio 1−G(x)
. If G(x) is the baseline cumulative distribution function (cdf) of a random variable, with probability density function (pdf) g(x) and the Weibull cumulative
distribution function is
b

F (x; a, b) = 1 − e−ax ,

x ≥ 0,

(1.5)

G(x)
with parameters a and b are positive. Based on this density, by replacing x with ratio 1−G(x)
. The cdf
of Weibull- generalized distribution, say Weibull-G distribution with two extra parameters a and b, is
defined by, Bourguignon et al. [8]

Z
F (x; a, b, λ) =

G(x;λ)
1−G(x;λ)

b

abtb−1 e−ax dt

0
G(x;λ)

= 1−e

−a[ 1−G(x;λ) ]b

,

x ≥ 0,

a, b ≥ 0,

(1.6)

where G(x; λ) is a baseline cdf, which depends on a parameter λ. The corresponding family pdf becomes
f (x; a, b, λ) = ab g(x; λ)

G(x;λ) b
[G(x; λ)]b−1 −a[ 1−G(x;λ)
]
e
b+1
[1 − G(x; λ)]

(1.7)

A random variable X with pdf Eq. (1.7) is denoted by X distributed weibll-G(a, b, λ), x ∈ R, a, b > 0.
the additional parameters induced by the weibull generator are sought as a manner to furnish a more
flexible distribution. If b = 1, it corresponds to the exponential- generator. An interpretation of the
weibull-G family of distributions can by given as follows ( Coorary, [9]) is a similar context.
Let Y be a lifetime random variable having a certain continuous G distribution. The odds ratio that an
G(x)
individual (or component) following the lifetime Y will die (failure) at time x is 1−G(x)
. Consider that
the variability of this odds of death is represented by the random variable X and assume that it follows
the Weibull model with scale a and shape b. We can write


G(x)
P r(Y ≤ x) = P r X ≤
= F (x; a, b, λ).
1 − G(x)
Which is given by Eq. (1.6). The survival function of the Weibull-G family is given by
G(x)

−a[ 1−G(x) ]b

R(x; a, b, λ) = 1 − F (x; a, b, λ) = e

,

(1.8)

and hazard rate function of the Weibull-G family is given by
f (x; a, b, λ)
ab g(x; λ)[G(x; λ)]b−1
=
1 − F (x; a, b, λ)
[1 − G(x; λ)]b+1
[G(x; λ)]b−1
= ab h(x; λ)
,
[1 − G(x; λ)]b

h(x; a, b, λ) =

(1.9)

3

Abdelfattah Mustafa, B. S. El-Desouky and Shamsan AL-Garash

b−1

g(x;λ)
where h(x; λ) = 1−G(x;λ)
. The multiplying quantity ab g(x;λ)[G(x;λ)]
works as a corrected factor for
[1−G(x;λ)]b
the hazard rate function of the baseline model Eq.(1.6) can deal with general situation in modeling survival data with various shapes of the hazard rate function. By using the power series for the exponential
function, we obtain
G(x)

e

−a[ 1−G(x) ]b


ib
∞
X
G(x; λ)
(−1)i ai
=
,
i!
1 − G(x; λ)

(1.10)

i=0

substituting from Eq.(1.10) into Eq. (1.7), we get
f (x; a, b, λ) = ab g(x; λ)

∞
X
(−1)i ai
i=0

i!

[G(x; λ)]b(i+1)−1
.
[1 − G(x; λ)]b(i+1)+1

(1.11)

Using the generalized binomial theorem we have
[1 − G(x; λ)]−(b(i+1)+1) =

∞
X
Γ(b(i + 1) + j + 1)
j=0

j!Γ(b(i + 1) + 1)

[G(x; λ)]j .

(1.12)

Inserting Eq. (1.12) in Eq. (1.11), the Weibull-G family density function is

f (x; a, b, λ) =

∞ X
∞
X
(−1)i ai+1 b Γ(b(i + 1) + j + 1)
i=0 j=0

i!j!Γ(b(i + 1) + 1)

g(x; λ)[G(x; λ)]b(i+1)+j−1 .

(1.13)

In section 2, we define the cumulative, density and hazard functions of the Weibull-G Exponential distribution (WGED) . In section 3 and 4, we introduced the statistical properties include, quantile function
skewness and kurtosis, rth moments and moment generating function. The distribution of the order
statistics is expressed in section 5. Finally, maximum likelihood estimation of the parameters is determined in section 6. Real data sets are analyzed in Section 7 and the results are compared with existing
distributions. Finally we introduce the conclusions in Section 8.

2

The Weibull Generalized Exponential Distribution

In this section, we study the three parameters Weibull-G exponential distribution (WGED). Using G(x)
and g(x) in Eq. (1.13) to be the cdf and pdf of Eq. (1.6) and Eq. (1.7). The cumulative distribution
function (cdf) of the Weibull-G exponential distribution (WGED) is given by
F (x; a, b, λ) = 1 − e−a[e

λx −1]b

,

x > 0,

a, b, λ > 0,

(2.1)

The pdf corresponding to Eq. (2.1) is given by
f (x; a, b, λ) = abλ eλx [eλx − 1]b−1 e−a[e

λx −1]b

,

x > 0,

where a, b > 0 and λ > 0 are two additional shape parameters.
Plots of the cdf, Eq. (2.1), of the WGED for some parameter values are displayed in Figure 1,

(2.2)

4

Abdelfattah Mustafa, B. S. El-Desouky and Shamsan AL-Garash

Figure 1: Cumulative probability function of the WGED.
We denote by X ∼ W GED(a, b, λ) a random variable having the pdf Eq. (2.1). The survival function,
S(x), hazard rate function, h(x), reversed hazard rate function, r(x) and cumulative hazard rate function
H(x) of X are given by
S(x; a, b, λ) = 1 − F (x; a, b, λ) = e−a[e

λx −1]b

h(x; a, b, λ) = abλ eλx [eλx − 1]b−1 ,

,

x > 0,

x > 0,

(2.3)

(2.4)

λx −1]b

r(x; a, b, λ) =

abλ eλx [eλx − 1]b−1 · e−a[e
1 − e−a[eλx −1]b

and

Z∞
H(x; a, b, λ) =

,

x > 0.

h λx
ib
h(x; a, b, λ)dx = a ee − 1 ,

(2.5)

(2.6)

0

respectively. Plots of S(x; a, b, λ), h(x; a, b, λ), r(x; a, b, λ) and H(x; a, b, λ) of the WGED for some
parameters values are displayed in Figure 2–6.

Figure 2: Probability density function of the WGED.

Abdelfattah Mustafa, B. S. El-Desouky and Shamsan AL-Garash

Figure 3: Survival function of the WGED.

Figure 4: Hazard rate function of the WGED.

Figure 5: Reversed hazard rate function of the WGED.

5

6

Abdelfattah Mustafa, B. S. El-Desouky and Shamsan AL-Garash

Figure 6: Cumulative hazard rate function of the WGED.
It is clear that the PDF and the hazard function have many different shapes, which allows this distribution
to fit different types of lifetime data.

3

Statistical Properties

In this section, we study the statistical properties for the WGED, specially quantile function and simulation median, skewness, kurtosis and moments.

3.1

Quantile and median

In this subsection, we determine the explicit formulas of the quantile and the median of WGED distribution. The quantile xq of the WGED is given by
F (xq ) = q,

0 < q < 1.

From Eq. (2.1), xq can be obtained as follows.
"

1 #
b
1
1
xq = ln 1 +
ln(1 − q)
.
λ
a
Setting q = 0.5 in Eq. (3.2), we get the median of WGED as follows.
"

1 #
1
ln(2) b
xq = ln 1 + −
.
λ
a

3.2

(3.1)

(3.2)

(3.3)

The mode

In this subsection, we will derive the mode of the WGED by derivation its pdf with respect to x and
equate it to zero. The mode is the solution the following equation with respect to x.
0

f (x) = 0.
By substitution PDF from Eq. (2.2) in Eq.(3.4), we have

(3.4)

7

Abdelfattah Mustafa, B. S. El-Desouky and Shamsan AL-Garash



∂
λxi
λxi
b−1
−a[eλxi −1]b
= 0,
abλ · e · [e − 1]
·e
∂x


∂
h(x; a, b, λ) · S(x; a, b, λ) = 0,
∂x


0
0
h (x; a, b, λ) · S(x; a, b, λ) + h(x; a, b, λ) · S (x; a, b, λ) = 0,
then



0

2



h (x; a, b, λ) − (h(x; a, b, λ)) S(x; a, b, λ) = 0,

(3.5)

where h(x; a, b, λ), S(x; a, b, λ) are hazard function and survival function of WGED respectively.
It is not possible to get an analytic solution in x to Eq. (3.5) in the general case. It has to be obtained
numerically by using methods such as fixed-point or bisection method.

3.3

Skewness and kurtosis

The analysis of the variability Skewness and Kurtosis on the shape parameters b, λ can be investigated
based on quantile measures. The short comings of the classical Kurtosis measure are well-known. The
Bowely’s skewness [10] based on quartiles is given by
Sk =

q(0.75) − 2q(0.5) + q(0.25)
,
q(0.75) − q(0.25)

(3.6)

and the Moors’ Kurtosis [11] is based on octiles
Ku =

q(0.875) − q(0.625) − q(0.375) + q(0.125)
,
q(0.75) − q(0.25)

(3.7)

where q(.) represents quantile function.

3.4

Moments

In this subsection, we discuss the rth moment for WGED. Moments are important in any statistical
analysis, especially in applications. It can be used to study the most important features and characteristics
of a distribution (e.g. tendency, dispersion, skewness and kurtosis).
Theorem 3.1. If X has WGED (a, b, λ), then The rth moments of random variable X, is given by the
following
0

µr =



∞ X
∞ X
∞
X
(−1)i+k ai+1 Γ(b(i + 1) + j + 1)Γ(r + 1) b(i + 1) + j + 1
.
i!j!λr (k + 1)r+1 Γ((b + i) + 1)
k

(3.8)

i=0 j=0 k=0

Proof. We start with the well known distribution of the rth moment of the random variable X with
probability density function f (x) given by
Z ∞
0
µr =
xr f (x; a, b, λ)dx.
(3.9)
0

Substituting from Eq. (1.1) and Eq. (1.2) into Eq. (1.13) we get

8

Abdelfattah Mustafa, B. S. El-Desouky and Shamsan AL-Garash

∞

0

µr =

Z
∞ X
∞
X
(−1)i ai+1 b · Γ(b(i + 1) + j + 1)
i!j!Γ((b + i) + 1)

i=0 j=0

h
ib(i+1)+j−1
dx,
xr λe−λx 1 − e−λx

0

since 0 < 1 − e−λx < 1 for x > 0, the binomial series expansion of
h
yields

1 − e−λx

ib(i+1)+j−1



∞
h
ib(i+1)+j−1 X
b(i + 1) + j − 1 −kλx
−λx
k
1−e
=
(−1) ·
e
,
k
k=0

then we get
0

µr =


Z
∞ X
∞ X
∞
X
(−1)i+k ai+1 λΓ(b(i + 1) + j + 1) b(i + 1) + j − 1
i!j!Γ((b + i) + 1)

i=0 j=0 k=0

k

∞

xr e−(k+1)λx dx,

0

by using the definition of gamma function in the form, Zwillinger [13],
Γ(Z) = x

z

Z∞

etx tz−1 dt,

z, x, > 0.

0

Finally, we obtain the rth moment of WGED in the form


∞ X
∞ X
∞
X
(−1)i+k ai+1 Γ(b(i + 1) + j + 1)Γ(r + 1) b(i + 1) + j + 1
.
µr =
i!j!λr (k + 1)r+1 Γ((b + i) + 1)
k
0

i=0 j=0 k=0

This completes the proof.

4

The Moment Generating Function

The moment generating function (mgf), MX (t), of a random variable X provides the basis of an alternative route to analytic results compared with working directly with the pdf and cdf of X.
Theorem 4.1. The moment generating function (mgf) of WGED is given by

MX (t) =



∞ X
∞ X
∞ X
∞
X
(−1)i+k ai+1 tr Γ(b(i + 1) + j + 1)Γ(r + 1) b(i + 1) + j + 1
r=0 i=0 j=0 k=0

i!j!r!λr (k + 1)r+1 Γ((b + i) + 1)

k

.

(4.1)

Proof. The moment generating function, MX (t), of the random variable X with probability density
function, f (x) is given by
Z ∞
MX (t) =
etx f (x; a, b, λ)dx,
0

using series expansion of

ext ,

we obtain

9

Abdelfattah Mustafa, B. S. El-Desouky and Shamsan AL-Garash

MX (t) =

∞ r Z
X
t

r!

r=0

∞
r

x f (x; a, b, λ)dx =

0

∞ r
X
t
r=0

r!

0

µr .

(4.2)

Substituting from Eq. 3.8 into Eq. 4.2, we obtain the moment generating function of WGED in the form


∞ X
∞ X
∞ X
∞
X
(−1)i+k ai+1 tr Γ(b(i + 1) + j + 1)Γ(r + 1) b(i + 1) + j + 1
.
MX (t) =
i!j!r!λr (k + 1)r+1 Γ((b + i) + 1)
k
r=0 i=0 j=0 k=0

This completes the proof.

5

Order Statistics

In this section, we derive closed form expressions for the pdf of the rth order statistic of the WGED. Let
X1:n , X2:n , · · · , Xn:n denote the order statistics obtained from a random sample X1 , X2 , · · · , Xn which
taken from a continuous population with cumulative distribution function (cdf), F (x; ϕ) and probability
density function (pdf), f (x; ϕ), then the probability density function of Xr:n is given
fr:n (x; ϕ) =

1
[F (x; ϕ)]r−1 [1 − F (x; ϕ)]n−r f (x; ϕ),
B(r, n − r + 1)

(5.1)

where f (x; ϕ), F (x; ϕ) are the pdf and cdf of W GED(ϕ) given by Eq. (2.2) and Eq.(2.1) respectively,
ϕ = (a, b, λ) and B(., .) is the beta function, also we define first order statistics X1:n = min(X1 , X2 , · · · , Xn ),
and the last order statistics as Xn:n = max(X1 , X2 , · · · , Xn ). Since 0 < F (x; ϕ) < 1 for x > 0, we
can use the binomial expansion of [1 − F (x; ϕ)]n−r given as follows
n−r

[1 − F (x; ϕ)]

=

n−r
X
i=0


n−r
(−1)i [F (x; ϕ)]i .
i

(5.2)

Substituting from Eq. (5.2) into Eq. (5.1), we obtain
n−r
X n − r
f (x; ϕ)
fr:n (x; ϕ) =
(−1)i [F (x; ϕ)]i+r−1 .
B(r, n − r + 1)
i

(5.3)

i=0

Substituting from Eq. (2.1) and Eq. (2.2) into Eq. (5.3), we obtain
fr:n (x; a, b, λ) =

n−r
X i+r−1
X
i=0 j=0



(−1)i+j n!
i+r−1
f (x; (j + 1)a, b, λ).
i!(r − 1)!(n − r − 1)!
j

(5.4)

Relation (5.4) show that fr:n (x; ϕ) is the weighted average of the Weibull-G exponential distribution
with different shape parameters.

6

Parameters Estimation

In this section, point and interval estimation of the unknown parameters of the WGED are derived by
using the method of maximum likelihood based on a complete sample data.

10

Abdelfattah Mustafa, B. S. El-Desouky and Shamsan AL-Garash

6.1

Maximum likelihood estimation:

Let x1 , x2 , · · · , xn denote a random sample of complete data from the WGED. The likelihood function
is given as
L=

n
Y

f (xi , a, b, λ),

(6.1)

i=1

substituting from (2.2) into (6.1), we have
L=

n
Y

λxi −1]b

abλ eλxi [eλxi − 1]b−1 e−a[e

.

i=1

The log-likelihood function is
L = n ln(abλ) + λ

n
X

xi + (b − 1)

i=1

n
X

λxi

ln(e

− 1) − a

i=1

n
X

[eλxi − 1]b .

(6.2)

i=1

The maximum likelihood estimation of the parameters are obtained by differentiating the log-likelihood
function L with respect to the parameters a, b and λ and setting the result equal to zero, we have the
following normal equations.
∂L
∂a
∂L
∂b
∂L
∂λ

=

n
ib
n X h λxi
−
e − 1 = 0,
a

(6.3)

i=1

=
=

n
+
b
n
+
λ

n
X
i=1
n
X



ln e

λxi

n h
ib 


X
−1 −a
eλxi − 1 ln eλxi − 1 = 0,

xi + (b − 1)

i=1

n
X
i=1

i=1
xi eλxi
eλxi − 1

− ab

n
X

ib−1
h
= 0.
xi eλxi eλxi − 1

(6.4)
(6.5)

i=1

The MLEs can be obtained by solving the nonlinear equations previous, (6.3)–(6.5), numerically for a, b
and λ.

6.2

Asymptotic confidence bounds

In this section, we derive the asymptotic confidence intervals of these parameters when a, b > 0 and
λ > 0 as the MLEs of the unknown parameters a, b > 0 and λ > 0 can not be obtained in closed forms,
by using variance covariance matrix I −1 see Lawless [14], where I −1 is the inverse of the observed
information matrix which defined as follows
−1
2
∂2L
∂2L
− ∂∂aL2 − ∂a∂b
− ∂a∂λ
2

∂2L
∂2L 
=  − ∂b∂a
− ∂∂bL2 − ∂b∂λ

∂2L
∂2L
∂2L
− ∂λ∂a − ∂λ∂b − ∂λ2


var(â) cov(â, b̂) cov(â, λ̂)
=  cov(b̂, â) var(b̂) cov(b̂, λ̂)  .
cov(λ̂, â) cov(λ̂, b̂) var(λ̂)


I−1

(6.6)

11

Abdelfattah Mustafa, B. S. El-Desouky and Shamsan AL-Garash

The second partial derivatives included in I are given as follows.
∂2L
∂a2
∂2L
∂a∂b
∂2L
∂a∂λ
∂2L
∂b2
∂2L
∂b∂λ
∂2L
∂λ2

= −
= −

n
,
a2
n h
X
i=1
n
X

= −b

(6.7)
eλxi − 1

ib



ln eλxi − 1 ,

(6.8)

h
ib−1
,
xi eλxi eλxi − 1

(6.9)

i=1
n

ib h 
i2
Xh
n
= − 2 −a
,
eλxi − 1 ln eλxi − 1
b
=

(6.10)

i=1
λx
e i

n
n
h
i
X
X
xi
λxi λxi
b−1
λxi
−
a
x
e
[e
−
1]
b
ln(e
−
1)
+
1
,
i
eλxi − 1
i=1

= −

n
− (b − 1)
λ2

i=1
n
X
i=1

(6.11)

n

ib−2
h

X
x2i eλxi
λxi
λxi
2 λxi
−
1
−
1
e
. (6.12)
be
−
ab
x
e
i
(eλxi − 1)2
i=1

We can derive the (1 − δ)100% confidence intervals of the parameters a, b and λ, by using variance
matrix as in the following forms
q
q
p
â ± Z δ var(â), b̂ ± Z δ var(b̂), λ̂ ± Z δ var(λ̂),
2

2

2

where Z δ is the upper ( 2δ )-th percentile of the standard normal distribution.
2

7

Application

In this section, we present the analysis of a real data set using the WGED (a, b, λ) model and compare it
with the other fitted models such as exponential distributions (ED), generalized exponential distribution
(GED), [12], beta exponential distribution (BED), [7] and the beta generalized exponential distribution
(BGED), [1] using KolmogorovSmirnov (K-S) statistic, as well as Akaike information criterion (AIC),
Akaike , [2], Bayesian information criterion (BIC) and Hannan-Quinn information criterion (HQIC)
values, Schwarz [16].
The data set is obtained from Smith and Naylor [19]. The data are the strengths of 1.5 cm glass fibres,
measured at the National Physical Laboratory, England. Unfortunately, the units of measurement are not
given in the paper. This data set is in Table 1.
Table 1: The data are the strengths of 1.5 cm glass fibres, [19].
0.55 0.93 1.25 1.36 1.49 1.52 1.58 1.61 1.64
1.68 1.73 1.81
2
0.74 1.04 1.27 1.39 1.49
1.53 1.59 1.61 1.66 1.68 1.76 1.82 2.01 0.77
1.11 1.28 1.42 1.5 1.54 1.6 1.62 1.66 1.69
1.76 1.84 2.24 0.81 1.13 1.29 1.48 1.5 1.55
1.61 1.62 1.66 1.7 1.77 1.84 0.84 1.24 1.3
1.48 1.51 1.55 1.61 1.63 1.67 1.7 1.78 1.89

12

Abdelfattah Mustafa, B. S. El-Desouky and Shamsan AL-Garash

Table 2 gives MLs of parameters of the WGED and sub-models and goodness of fit statistics are in Table
3.

Model
ED
GED
BED
BGED
WGED

Table 2: MLEs of parameters, Log-likelihood.
MLEs of parameters
K-S
λ̂=0.664
0.402316
λ̂=2.6105, α̂=31.3032
0.213118
â=17.7786, b̂=22.7222, λ̂=0.3898
0.159819
â=0.4125, b̂=93.4655, λ̂=0.92271, α̂=22.6124 0.150611
â=56.881, b̂=4.893, λ̂=0.222
0.127366

p-value
1.44529 × 10−09
0.005444
0.07220
0.10470
0.24259

Table 3: Log-likelihood, AIC, AICC, BIC and HQIC values of models fitted.
Model
L
−2L
AIC
AICC
BIC
HQIC
ED
-88.8300 -177.6600 179.6600 179.7256 181.8031 180.5029
GED
-31.3834 -62.7668
66.7668
66.9668
71.0531
68.4526
BED
-24.1270 -48.2540
54.2540
54.6608
60.6834
56.7827
BGED -15.5995 -31.1990
39.1990
39.8887
47.7715
42.5706
WGED -14.828
-29.6560
35.6560
36.0628
42.0854
38.1847
We find that the WGED distribution with the three-number of parameters provides a better fit than the
previous new modified exponential distribution which was the best in [1, 3, 7, 12]. It has the largest
likelihood, and the smallest K-S, AIC, BIC and HQIC values among those considered in this paper.
Substituting the MLE’s of the unknown parameters a, b, λ into (6.6), we get estimation of the variance
covariance matrix as the following


3.655 × 103
7.228
−2.205
7.228
0.213
1.141 × 10−3 
I0−1 = 
−2.205
1.141 × 10−3 1.505 × 10−3
The approximate 95% two sided confidence intervals of the unknown parameters a, b and λ are
[0, 175.11], [3.989, 57.785] and [0.146, 0.298], respectively.
To show that the likelihood equation have unique solution, we plot the profiles of the log-likelihood
function of a, b and λ in Figures 7 and 8.

Figure 7: The profile of the log-likelihood function of a, b.

Abdelfattah Mustafa, B. S. El-Desouky and Shamsan AL-Garash

13

Figure 8: The profile of the log-likelihood function of λ.
The nonparametric estimate of the survival function using the Kaplan-Meier method and its fitted parametric estimations when the distribution is assumed to be ED, GE, BED, BGED and W GED are
computed and plotted in Figure 9.

Figure 9: The Kaplan-Meier estimate of the survival function for the data.
Figure 10 gives the form of the CDF for the ED, GE, BED, BGED and W GED which are used to fit
the data after replacing the unknown parameters included in each distribution by their MLEs.

Figure 10: The Fitted cumulative distribution function for the data.

Abdelfattah Mustafa, B. S. El-Desouky and Shamsan AL-Garash

8

14

Conclusions

A new distribution, based on Weibull- G family distributions, has been proposed and its properties studied. The idea is to add parameter to exponential distribution, so that the hazard function is either increasing or more importantly, bathtub shaped. Using Weibull generator component, the distribution has
flexibility to model the second peak in a distribution. We have shown that the Weibull-G exponential distribution fits certain well-known data sets better than existing modifications of the generalized families
of exponential probability distribution.

References
[1] Wagner Barreto-Souza, Alessandro H. S. and Gauss M. Cordeiro (2009). The Beta Generalized
Exponential Distribution. Journal of Statistical Computation and Simulation, Vol. 80, pp. 159–172.
[2] Akaike, H. (1974). A new look at the statistical model identification. IEEE Transactions on Automatic Control, AC-19, pp. 716-23.
[3] Gupta, R. D. and Kundu, D. (2001). Exponentiated exponential family; an alternative to gamma
and Weibull. Biometrical Journal, Vol. 43, pp. 117–130.
[4] Gupta, R. D. and Kundu, D. (2001). Generalized exponential distributions: different methods of
estimation. Journal of Statistical Computation and Simulation, Vol. 69, pp. 315–338.
[5] Gupta, R. D. and Kundu, D. (2002). Generalized exponential distributions: statistical inferences.
Journal of Statistical Theory and Applications, Vol. 1, pp. 101–118.
[6] Gupta, R. D. and Kundu, D. (2003). Discriminating between the Weibull and the GE distributions.
Computational Statistics and Data Analysis, Vol. 43, pp. 179–196.
[7] Nadarajah, S. and Kotz, S., (2006). The beta exponential distribution. Reliability Engineering and
System Safety, Vol. 91, pp. 689-697
[8] Marcelo, B. Silva, R. and Cordeiro, G., (2014). The Weibull - G Family Probability Distributions.
Journal of Data Science, Vol. 12, pp. 53–68.
[9] Cooray, K. (2006). Generalization of the Weibull distribution: the odd Weibull family. Statistical
Modelling, Vol. 6, pp. 265–277.
[10] Kenney, J. and Keeping, E. Mathematics of Statistics. Vol. Princeton, (1962).
[11] Moors, J.J.A. (1998). A quantile alternative for kurtosis. The Statistician, Vol. 37, pp. 25-32.
[12] Gupta, R.D. and Kundu, D. (1999). Generalized exponential distribution. Australian and New
Zealand Journal of Statistics, Vol. 41, no. 2, pp. 173-88.
[13] Zwillinger, D. Table of integrals, series, and products. Elsevier, (2014).
[14] Lawless, J. F. Statistical Models and Methods for Lifetime Data. John Wiley and Sons, New York,
Vol. 20, pp. 1108-1113, (2003).

Abdelfattah Mustafa, B. S. El-Desouky and Shamsan AL-Garash

15

[15] Gurvich, M. R. DiBenedetto, A. T. and Ranade, S. V. (1998). A new statistical distribution for
characterizing the random strength of brittle materials. Journal of Materials Science, Vol. 32, pp.
2559–2564.
[16] Schwarz, G. (1978). Estimating the dimension of a model. Annals of Statistics, Vol. 6, pp. 461-4.
[17] Zografos, K. and Balakrishnan, N. (2009). On families of beta- and generalized gamma-generated
distributions and associated inference. Statistical Methodology, Vol. 6, pp. 344-362.
[18] Weibull, W. (1951). Wide applicability. Journal of applied mechanics, Vol. 40, pp. 203–210.
[19] Smith, R. L. and Naylor, J. C. (1987). A comparison of maximum likelihood and Bayesian estimators for the three-parameter Weibull distribution. Applied Statistics, Vol. 36, pp. 358–369

