The Marshall-Olkin-Kumaraswamy-G family of distributions
Laba Handique and Subrata Chakraborty*
Department of Statistics, Dibrugarh University
Dibrugarh-786004, India
*Corresponding author: Email: subrata_stats@dibru.ac.in
(August 21, 2016)
Abstract
A new family of continuous distribution is proposed by using Kumaraswamy-G (Cordeiro and de
Castro, 2011) distribution as the base line distribution in the Marshal-Olkin (Marshall and Olkin,
1997) construction. A number of known distributions are derived as particular cases. Various
properties of the proposed family like formulation of the pdf as different mixture of exponentiated
baseline distributions, order statistics, moments, moment generating function, Renyi entropy, quantile
function and random sample generation have been investigated. Asymptotes, shapes and stochastic
ordering are also investigated. The parameter estimation by methods of maximum likelihood, their
large sample standard errors and confidence intervals and method of moment are also presented. Two
members of the proposed family are compared with corresponding members of KumaraswamyMarshal-Olkin-G family (Alizadeh et al., 2015) by fitting of two real life data sets.
Key words: Kumaraswamy-G distribution, Marshall-Olkin family, Exponentiated family, AIC, BIC,
Power weighted moment.

Preprint

1

1. Introduction
One of the preferred area of research in the filed of the probability distribution is that of
generating new distributions starting with a base line distribution by adding one or more additional
parameters notable among them are Azzalini’s skewed family (Azzalini, 1985), Marshall-Olkin
extended (MOE ) family (Marshall and Olkin, 1997), exponentiated family (EF ) of distributions
(Gupta et al., 1998), or by

composite methods of combining two or more known competing

distribution through transformations like beta-generated (beta  G ) family (Eugene et al., 2002; Jones
2004), gamma-generated (G  G ) families (Zografos Balakrishnan 2009, Ristic and Balakrishnan,
2012), Kumaraswamy - G ( Kw  G ) family (Cordeiro and de Castro, 2011 ; Nadarajah et al., 2012;
Hussain, 2013), McDonald-G ( Mc  G ) (Alexander et al., 2012), beta extended-G family (Cordeiro
et al., 2012), Kumaraswamy -beta generalized family (Pescim et al., 2012), exponentiated transformed
transformer family ( ET  X ) , (Alzaghal et al. 2013), exponentiated generalized ( Exp  G ) family
(Cordeiro et al., 2013), geometric exponential-Poisson family (Nadarajah et al., 2013a), truncatedexponential skew symmetric family (Nadarajah et al., 2013b), logistic-generated ( Lo  G ) family
(Torabi and Montazari, 2014) and Kumaraswamy Marshal - Olkin family (Alizadeh et al., 2015).
While the additional parameter(s) bring in more flexibility at the same time they also complicates the
mathematical form of the resulting distribution, often considerably enough to render it not amenable to
further analytical and numerical manipulations. But with the advent of sophisticated powerful
mathematical and statistical softwares unlike in past now a days more and more complex distributions
are getting accepted as viable models of data analysis. Tahir et al. (2015) provided a detail account of
how new families of univariate continuous distributions are being generated through introduction of
additional parameters.
1.1 Formulas and notations
Here we list some formulas to be used in the subsequent sections of this article.
If T is a continuous random variable with pdf, f (t ) and cdf F (t )  P [T  t ] , then its
Survival function (sf): F (t )  P [T  t ]  1  F (t ) ,
Hazard rate function (hrf): h(t )  f (t ) / F (t ) ,
Reverse hazard rate function (rhrf): r (t )  f (t ) / F (t ) ,
Cumulative hazard rate function (chrf): H (t )   log [ F (t )] .


th

( p, q, r ) Power Weighted Moment (PWM): p, q , r   t p [ F (t ) ] q [1  F (t ) ] r f (t ) dt ,


2



Entropy: The Rényi entropy of T is given by I R ( )  (1   ) 1 log  f (t )  dt  .
 


1.2 Marshall-Olkin Extended (MOE ) family of distributions
Starting with a given baseline distribution with probability density function (pdf) f (t ) , cumulative
distribution function (cdf) F (t ) , Marshall and Olkin (1997) proposed a new flexible semi parametric
family of distributions and defined a new survival function (sf) F MO (t ) by introducing an additional
parameter   0. The sf F MO (t ) of the MOE family of distributions is defined by

F MO (t ) 

 F (t )
 F (t )
 [1  F (t )]


1   F (t ) F (t )   F (t )    F (t )

(1)

where    t  ,   0 and   1   . The parameter '  ' is known as the tilt parameter as since the
hrf of the new family is shifted below (above) for   1(0    1) the hrf of the base line distribution
(Nanda and Das, 2012). That is for all t  0, hMO (t )  h(t ) when   1 , and h MO (t )  h(t )
when 0    1 , where h MO (t ) and h (t ) are the hrf’s of the MOE and baseline distributions
respectively. Now
F MO (t )  1  F MO (t )  1 



and

f

MO

(t ) 

F (t )
F (t )
 F (t )
1  (   ) F (t )



1   F (t )
1   F (t )
1   F (t ) 1  (1   ) F (t )

F (t )
1  F (t )
F (t )
or

F (t )   F (t )
   F (t )    F (t )

(2)

d
F MO ( t )  d  F (t )    f (t ) 2 or  f (t ) 2
dt
dt    F (t )  [   F (t )]
[1   F (t )]

(3)

Where    t  ,   0 and   1   . If   1 , then we have F MO (t )  F (t ) . Other reliability
measures like the hrf, rhrf and chrf associated with (1) are
h MO (t ) 

f MO (t )
 f (t )

MO
F (t ) [1   F (t )]2

r MO (t ) 

f MO (t )
f (t )
1
 r (t )
 r (t )
or


MO
F (t ) [1   F (t )] 1   F (t )
   F (t )
F (t )

 F (t )
f (t )
1
h (t )


1   F (t ) F (t ) [1   F (t )] 1   F (t )

  F (t ) 
 [1  F (t )] 
or  log 
H MO (t )   log 


1   F (t ) 
    F (t ) 

respectively.

Where h(t ) and r (t ) is the hrf and rhrf of the baseline distribution.

3

and

It is obvious that many new families can be derived from Marshall-Olkin set up by considering
different base line distribution in the equation (1). These new families are usually termed as MarshallOlkin extended F distribution. For details see Tahir et al. (2015).
Jayakumar and Mathew (2008) generalized the Marshall-Olkin set up by exponentiating the
Marshall-Olkin survival function


  F (t ) 
F GMO (t )  
 .
1   F (t ) 
This method is called method of Lehmann alternative 1 due to Lehmann (1953)
Tahir et al. (2015) propose another generalization through Lehmann alternative 2 due to
Lehmann (1953) by exponentiating the Marshall-Olkin cdf as


F

GMO


 F (t ) 
(t )  1 
 .
 1   F (t ) 

[for more on Marshall - Olkin family check Barakat et al. (2009), Jose (2011), Krishna (2011),
Barreto-Souza et al. (2013), Cordeiro et al. (2014a), Alizadeh et al., (2015)]
1.3 Kumaraswamy-G ( Kw  G ) family of distributions
For a baseline cdf G(t ) with pdf g (t ) , Cordeiro and de Castro (2011) defined Kw  G distribution
with cdf and pdf

and

F KwG (t )  1  [1  G (t ) a ] b

(4)

f KwG (t )  a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1

(5)

Where t  0, g (t ) 

d
G (t ) and a  0, b  0 are shape parameters in addition to those in the baseline
dt

distribution which partly govern skewness and variation in tail weights. For a lifetime random
variable 't ' , the sf, hrf, rhrf and chrf for distribution in (4) are given respectively by

F KwG (t )  1  F KwG (t )  [1  G (t ) a ] b
h KwG (t )  f

KwG

(t ) F

KwG

(t )  a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1 [1  G (t ) a ] b
 a b g (t ) G (t ) a 1[1  G (t ) a ] 1

r KwG (t )  f

KwG

(t ) F

KwG

(t )  a b g (t ) G (t ) a1 [1  G (t ) a ] b1 1  [1  G (t ) a ] b



 a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1 1  [1  G (t ) a ] b
and

H KwG (t )  b log [1  G (t ) a ] .

4



1

Recently, Alizadeh et al., (2015) proposed the Kumaraswamy Marshal-Olkin family of distributions
by using the Marshal-Olkin (Marshall and Olkin 1997) cdf in that Kumaraswamy-G (Cordeiro and de
Castro 2011) family and studied its many properties.
The main motivation behind the present article is to propose another family of continuous
probability distribution that generalizes the Kumaraswamy-G (Cordeiro and de Castro 2011) family
and the Marshall-Olkin Extended family (Marshall and Olkin 1997) by integrating the former as the
base line distribution in the later. We call this new family the Marshall-Olkin Kumaraswamy-G

( MOKw  G ) family of distribution which encompasses many known families of distributions and
study some its general properties, parameter estimation and real life application in the present article.
The rest of this article is organized in seven sections. In section 2 the new family is defined
along with its physical basis. Next section is devoted to presenting some important special cases of the
family along with their shape and main reliability characteristics. In section 4 we discuss some
important general results of the proposed family, while different methods of estimation of parameters
are presented in section 5. In section 6 we present two examples of comparative data fitting. The
article ends with a concluding discussion and remark in the final section followed by an appendix to
derive asymptotic confidence bounds.
2. New Marshall-Olkin Kumaraswamy-G ( MOKw  G ) family of distributions
We now propose a new extension of the MO (Marshall and Olkin 1997) family by considering the cdf
and pdf of

Kw  G (Cordeiro and de Castro 2011) distribution in (4) and (5) as the f (t ) and F (t )

respectively in the MO formulation in (3) and call it MOKw  G distribution. The pdf of
MOKw  G is given by

f

MOKwG

(t )   a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1 [   {1  [1  G (t ) a ] b }] 2
  a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1 [(   )   [1  G (t ) a ] b ] 2
  a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1 [1   [1  G (t ) a ] b ] 2

(6)

,    t  ,   0, a  0, b  0
Similarly using equation (4) in (2) we get the cdf and hence the sf of MOKw  G respectively as





F MOKwG (t )  1  [1  G (t ) a ] b    1  [1  G (t ) a ] b  1  [1  G (t ) a ] b 1   [1  G (t ) a ]b

(7)

and F MOKwG (t )  1  F MOKwG (t )  [1  G (t )a ]b (1   ) 1   [1  G (t )a ] b

  [1  G(t )a ]b 1   [1  G (t ) a ] b
hrf: h MOKwG (t ) 

f MOKwG (t )  a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1 [1   [1  G (t ) a ] b ] 2

 [1  G (t )a ] b 1   [1  G (t ) a ] b
F MOKwG (t )

5

(8)



rhrf : r MOKwG (t ) 

a b g (t ) G (t ) a 1[1  G (t ) a ] 1
1   [1  G (t ) a ] b

(9)

f MOKwG (t )  a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1 [1   [1  G (t ) a ] b ] 2

1  [1  G (t ) a ] b 1   [1  G (t ) a ] b
F MOKwG (t )


 a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1 [1  [1  G (t ) a ] b ] 1
1   [1  G (t ) a ] b



chrf: H MOKwG (t )   log [ F MOKwG (t ) ]   log  [1  G (t ) a ] b 1   [1  G (t ) a ] b



For
(i)   1 , the pdf in (6) reduces to that of Kw  G (Cordeiro and de Castro, 2011 )
(ii) a  b  1 , the pdf in (6) reduces to that of MO (Marshall and Olkin, 1997).
2.1 Genesis of the distribution
Let Ti ( i  1, 2,, N ) be a sequence of i. i. d. random variables with survival function [1  G (t ) a ] b , and
let N be a positive integer random variable independent of the Ti ’s given by the pgf of a geometric
distribution with parameter  , say  ( x ;  )   x (1   x)1 . Then the inverse of  ( x ;  ) becomes

 1 ( x ;  )   ( x ;  1 ) .We can verify that equation (8) comes from F MOKwG (t )   [{1  G (t ) a } b ;  ] for
0    1 and F MOKwG (t )   [{1  G (t ) a } b ;  1 ] for   1 . For both cases, equation (8) represents the
survival function of the minimum and maximum of Ti ( i  1, 2,, N ) , where N has probability
generating function  (x ; .) with parameter  and  1 respectively.
2.2 Shape of the density and hazard functions
Here we have plotted the pdf and hrf of the MOKw  G for some choices of the distribution G and
parameter values to study the variety of shapes assumed by the family.

6

MOKw-E

1.5
1.0
0.0

0.0

1.0

1.5

2.0

2.5

3.0

0

1

2

3

X

X

(a)

(b)

MOKw-L

MOKw-Fr
2.0

0.5

2.0

0.0

4

1.5

a=1.2,b=140.7,=0.3,=29.3,=65.5
a=1,b=125.9,=0.5,=20,=105.4
a=0.55,b=16.9,=0.2,=40,=15.4
a=0.7,b=18.9,=0.2,=40,=10.4
a=1,b=75,=0.39,=22,=102.5
a=1,b=90,=0.399,=20,=102.5

1.0
0.0

0.0

0.5

1.0

Density

1.5

a=2.55,b=1.5,=5,=1.5,=0.99
a=5.5,b=0.5,=2,=1.5,=0.1
a=0.05,b=0.45,=2,=1.5,=0.5
a=46,b=0.5,=2.5,=0.5,=0.15
a=0.15,b=0.6,=2,=1.5,=0.5
a=3.55,b=0.5,=5,=1.5,=0.99

0.5

Density

a=2,b=0.9,=2,=1.5,=0.9
a=2.1,b=0.5,=2,=1.5,=0.1
a=40,b=3.5,=2.5,=0.5,=0.5
a=0.05,b=0.45,=2,=1.5,=0.5
a=0.15,b=0.6,=2,=1.5,=0.5
a=33,b=3.5,=2.1,=0.5,=0.5

0.5

Density

1.0

1.5

a=7.5,b=2.5,=1.5,=1.8
a=1.4,b=1.7,=2.5,=1.9
a=6.5,b=2.5,=1.5,=5.5
a=2.5,b=1.5,=1.7,=0.6
a=0.02,b=0.42,=2.1,=1.5
a=2.5,b=2.5,=1.5,=2.5

0.5

Density

MOKw-W

0

1

2

3

4

0.0

0.5

1.0

1.5

2.0

X

X

(c)

(d)

2.5

3.0

3.5

Fig 1: Density plots (a) MOKw  E , (b) MOKw  W , (c) MOKw  L and (d) MOKw  Fr distributions.

7

MOKw-E

1.5
1.0

h(x)

1.0

1.5

a=1.5,b=0.45,=1.9,=1.5
a=1.5,b=0.5,=1.9,=1.5
a=0.05,b=0.6,=0.05,=1.5
a=0.05,b=0.52,=0.05,=1.5
a=0.15,b=0.52,=2,=1.5
a=0.25,b=0.5,=2,=1.5

1.0

1.5

2.0

0.0

1.0

X

X

(a)

(b)

MOKw-L

MOKw-Fr

a=53,b=9.55,=2.65,=0.5,=52
a=5.6,b=0.5, =2.2,=1.5,=0.1
a=0.05,b=0.06,=0.001,=0.09,=0.01
a=6.5,b=0.6, =2,=1.5,=0.1
a=50,b=9.6, =2.65,=0.5,=55
a=0.06,b=0.04,=0.003,=0.09,=0.05

1.5

2.0

a=1,b=125.9,=0.5,=20,=105.4
a=1,b=124.9,=0.55,=15,=104.4
a=0.8,b=20.9,=0.01,=39.5,=30.4
a=0.7,b=15.9,=0.01,=40,=25.5
a=1,b=124.9,=0.49,=22,=109.5
a=0.7,b=15.9,=0.02,=40,=25.5

0.0

0.0

0.2

0.5

0.4

0.6

h(x)

1.0

0.8

1.0

0.5

1.5

0.5

1.2

0.0

h(x)

a=50,b=9.5,=2.55,=0.5,=5.5
a=45,b=10,=2.5,=0.5,=5.5
a=0.06,b=0.5,=0.15,=4.5,=0.4
a=0.05,b=0.45,=0.15,=4.7,=0.9
a=0.01,b=0.5,=0.15,=4.5,=0.7
a=0.01,b=0.19,=0.15,=6.2,=0.9

0.0

0.0

0.5

0.5

h(x)

MOKw-W

0.0

0.5

1.0

1.5

0.0

0.5

1.0

X

X

(c)

(d)

1.5

2.0

Fig 2: Hazard plots (a) MOKw  E , (b) MOKw  W , (c) MOKw  L and (d) MOKw  Fr distributions.
From the plots in figure 1 and 2 it can be seen that the family is very flexible and can offer
many different types of shapes of density and hazard rate function including the bath tub shaped for
hazard.
3. Special models
In this section we provide some special cases of the MOKw  G family of distributions and list their
main distributional characteristics.

8

3.1 The MOKw  exponential ( MOKw  E ) distribution
Let the base line distribution be exponential with parameter   0, g (t )   e  t , t  0 and
G (t )  1  e  t ,

t  0 then for the MOKw  E model we get

pdf:

f MOKwE (t )   a b  e  t (1  e  t ) a 1 [1  (1  e t ) a ] b 1 [1   [1  (1  e  t ) a ] b ] 2

cdf :

F MOKwE (t )  1  [1  (1  e  t ) a ] b 1   [1  (1  e  t ) a ] b

sf :

F MOKwE (t )   [1  (1  e t ) a ] b 1   [1  (1  e t ) a ] b

hrf :

h MOKwE (t ) 

a b  e  t (1  e  t ) a 1 [1  (1  e  t ) a ] 1
1   [1  (1  e  t ) a ] b

rhrf: :

r MOKwE (t ) 

 a b  e  t (1  e  t ) a 1 [1  (1  e  t ) a ] b 1{1  [1  (1  e  t ) a ] b } 1
1   [1  (1  e  t ) a ] b

H MOKwE (t )   log [ F

chrf :

MOKwE

(t ) ]   log [  [1  (1  e  t ) a ] b 1   [1  (1  e  t ) a ] b ]

3.2 The MOKw  Lomax ( MOKw  L ) distribution
Considering the

Lomax distribution (Ghitany et al., 2007) with pdf and

cdf given

by g (t )  (   ) [ 1  (t  )]  (  1) , t  0 , and G ( t )  1  [ 1  ( t  ) ]   ,   0 and   0 the pdf of the
MOKw  L distribution is given by

pdf: f MOKwL (t ) 

 a b (   )[1  (t  ) ]  (  1) [1  [1  (t  )]   ] a 1[1  [1  [1  (t  )]   ] a ] b 1
[1   [1  [1  [1  (t  )]   ] a ] b ] 2

cdf:

F MOKwL (t ) 

1  [1  {1  [1  (t  )]   } a ] b
1   [1  {1  [1  (t  )]   } a ] b

sf:

F MOKwL (t ) 

 [1  {1  [1  (t  )]   } a ] b
1   [1  {1  [1  (t  )]   } a ] b

hrf:
rhrf:

h MOKwL (t ) 

a b (   )[1  (t  )]  (  1) [1  [1  (t  )]   ] a 1 [1  {1  [1  (t  )]   } a ] 1
1   [1  {1  [1  (t  )]   } a ] b

r MOKwL (t ) 

 a b (   )[1  (t  ) ]  (  1) [1  [1  (t  ) ]   ] a 1

1   [1  {1  [1  (t  )]   } a ] b

  a b 1
  a b 1 
[1  {1  [1  (t  ) ] } ] [1  [1  {1  [1  (t  )] } ] ] 
chrf:

H MOKwL (t )   log [ F

MOKwE

  [1  {1  [1  (t  ) ]   } a ] b 
(t ) ]   log 
 a b 
1   [1  {1  [1  (t  ) ] } ] 

3.3 The MOKw  Weibull ( MOKw  W ) distribution

9

Considering the Weibull distribution (Ghitany et al., 2005, Zhang and Xie, 2007, and Caroni, 2010)
with

parameters

 0

 0

and

having

pdf

and

g (t )    t  1 e   t

cdf



and



G (t )  1  e   t respectively we get the pdf of MOKw  W distribution as








f MOKwW (t ) =  a b   t  1 e   t [1  e   t ] a 1 [1  [1  e   t ] a ] b 1 [1   [1  [1  e   t ] a ] b ] 2




cdf: F MOKwW (t )  {1  [1  [ 1  e   t ] a ] b } /{1   [1  [ 1  e   t ] a ] b }
sf:





F MOKwW (t )  { [ 1  [ 1  e  t ] a ] b } /{1   [ 1  [ 1  e   t ] a ] b}








hrf: h MOKwW (t )  {a b   t  1 e   t [1  e  t ]a 1 [1  [ 1  e   t ] a ] 1} /{1   [1  [ 1  e   t ] a ] b }


rhrf: r MOKwW (t ) 







 a b   t  1 e   t [1  e   t ] a 1[1  [1  e   t ] a ] b 1[1  [1  [1  e   t ] a ] b ] 1

1   [1  [1  e   t ] a ] b


chrf: H

MOKwW

(t )   log [ F

MOKwW

  [ 1  [ 1  e  t ] a ] b 
(t ) ]   log 

 t  a b
] ] 
1   [1  [1  e

3.4 The MOKw  Frechet ( MOKw  Fr ) distribution
Suppose the base line distribution is the Frechet distribution (Krishna et al., 2013) with pdf and cdf




given by g (t )     t  (  1) e  ( t ) and G (t )  e  ( t ) , t  0 respectively, and then the corresponding
pdf of MOKw  Fr distribution becomes


pdf: f

MOKwFr



cdf: F

MOKwFr

hrf:

MOKwFr

(t ) 

1  [1  [ e  (

t) 

] a ]b


1   [1  [ e  (  t ) ] a ] b



sf:



h

r

F

MOKwFr

 [1  (e  (  t ) ) a ] b
(t ) 

1   [1  [ e  (  t ) ] a ] b



a b    t  ( 1) e  ( t ) ( e  ( t ) ) a 1[1  (e  (
(t ) 

1   [1  [ e  ( t ) ] a ] b


rhrf:



 a b    t  (  1) e  ( t ) [ e  ( t ) ] a 1 [1  [ e  ( t ) ] a ] b 1
(t ) 

[1   [1  [e  ( t ) ] a ] b ] 2

MOKwFr

t )



) a ] 1





 a b    t  (  1) e ( t ) (e  ( t ) ) a 1[1  (e  ( t ) ) a ] b 1[1  [1  [ e  ( t ) ] a ] b ] 1
(t ) 

1   [1  [ e  (  t ) ] a ] b


chrf: H MOKwFr (t )   log [ F

MOKwFr

  [1  (e ( t ) ) a ] b 
(t ) ]   log 

 ( t )  a b
] ] 
1   [1  [ e

3.5 The MOKw  Gompertz ( MOKw  Go ) distribution
Next by taking the Gompertz distribution (Gieser et al., 1998) with pdf and cdf g (t )   e t e
and G (t )  1  e



 t
( e 1)




 t
( e 1)


,   0,   0, t  0 respectively, we get the pdf of MOKw  Go distribution as

10

pdf:

f MOKwGo (t ) 

 a b  e t e



 t
( e 1 )


[1  e



 t
( e 1)


[1   [1  {1  e

cdf: F MOKwGo (t ) 

1  [1  {1  e



1   [1  {1  e
hrf: h MOKwGo (t ) 

a b  e t e



 t
( e 1)



}a ]b

 t
( e 1)


 t
( e 1)


[1  e



] a 1 [1  {1  e

 t
( e 1)


sf: F MOKwGo (t ) 

} a ] b 1

 [1  {1  e



 t
( e 1)


1   [1  {1  e

 t
( e 1)
a 1


1   [1  {1  e

 t
( e 1)


}a ]b ]2

}a ]b




]



[1  [1  e

 t
( e 1)






}a ]b

 t
( e 1)


}a ]b

 t
( e 1)
a


] ] 1

}a ]b

rhrf: r MOKwGo (t )



 a b  e t e



 t
( e  1)


[1  e



 t
( e 1 )


] a 1[1  {1  e

1   [1  {1  e

chrf: H

MOKwGo

(t )   log [ F

MOKwGo



 t
( e 1 )



 ( e t 1)


} a ] b 1[1  {1  [1  e



 t
( e  1)


] a } b ] 1

}a ]b


 ( e t 1)


 [1  {1  e 
}a ]b 

(t )]   log



 ( e t 1)
} a ] b 
1   [1  {1  e 

3.6 The MOKw  EW distribution
The

extended

Weibull

(EW)

G (t )  1  exp[ E (t :  )] ,

distributions

t  D  R ,   0

of

Gurvich

et

al.

(1997)

has

the

cdf

where E (t :  ) is a non-negative monotonically

increasing function which depends on the parameter vector  . The corresponding pdf is

g (t )   exp[ E (t :  )] e(t :  ) where e (t :  ) is the derivative of E (t :  ) .
Some important particular cases of EW can be seen as follows:
(i) E (t :  ) = t : exponential distribution.
(ii) E (t :  ) = t 2 : Rayleigh (Burr type-X) distribution.
(iii) E (t :  ) = log (t k ) : Pareto distribution
(iv) E (t :  ) =  1[exp (  t )  1] : Gompertz distribution.
Here we derive MOKw  EW by considering EW as the base line distribution as
pdf: f MOKwEW (t )


 a b exp[ E (t :  )]e (t :  ) [1  exp [ E (t :  )] ] a 1[1  {1  exp [ E (t :  ) ]} a ] b 1
[1   [1  { 1  exp [ E (t :  ) ]} a ] b ] 2

cdf:

F MOKwEW (t ) 

1  [1  {1  exp [ E (t :  ) ]} a ] b
1   [1  {1  exp [ E (t :  ) ]} a ] b

11

F MOKwEW (t ) 

sf:

hrf: h MOKwEW (t ) 

 [1  {1  exp [ E (t :  ) ]} a ] b
1   [1  {1  exp [ E (t :  ) ]} a ] b

a b  exp[ E (t :  )] e (t :  )[1  exp { E (t :  ) } ] a 1 [1  { 1  exp [ E (t :  ) ]} a ] 1
1   [1  {1  exp [ E (t :  ) ]} a ] b

rhrf: r MOKwEW (t )


 a b  exp[ E (t :  )] e (t :  )[1  exp { E (t :  ) } ] a1 [1  {1  exp [ E (t :  ) ]} a ] b1
1   [1  {1  exp [ E (t :  ) ]} a ] b
 [1  [1  {1  exp [ E (t :  ) ]} a ] b ] 1

chrf:

  [1  {1  exp [  E (t :  ) ]} a ] b 
H MOKwEW (t )   log 
a b 
1   [1  {1  exp [  E (t :  ) ]} ] 

3.7 The MOKw  Extended modified Weibull ( MOKw  EMW ) distribution
The cdf and pdf of the modified Weibull (MW) distribution (Sarhan and Zaindin, 2013) is given
respectively by
G (t )  1  exp [  t   t  ] , t  0,   0,  ,   0,     0 and

g (t )  (     t  1 ) exp [   t   t  ] .
The corresponding pdf of MOKw  EMW is given by

f MOKwEMW (t )


 a b (     t  1 ) exp ( t   t  ) [1  exp ( t   t  ) ] a 1[1  [1  exp ( t   t  ) ] a ] b 1
cdf:
[1   {1  [1  exp ( t   t  )] a }b ] 2

F MOKwEMW (t ) 

1  {1  [1  exp ( t   t  )] a }b
1   {1  [1  exp ( t   t  )] a }b

sf:

F MOKwEMW (t ) 

hrf:

h MOKwEMW (t )



 {1  [1  exp ( t   t  )] a }b
1   {1  [1  exp ( t   t  )] a }b

a b (    t  1 ) exp ( t   t  ) [1  exp ( t   t  ) ] a 1[1  [1  exp ( t   t  ) ]a ] 1
1   {1  [1  exp ( t   t  )] a }b

rhrf: r MOKwEMW (t )


 a b (    t  1 ) exp ( t   t  ) [1  exp ( t   t  ) ] a 1[1  [1  exp ( t   t  ) ]a ] b 1
1   {1  [1  exp ( t   t  )] a }b
 [1  {1  [1  exp (  t   t  )] a }b ] 1

12

  {1  [1  exp ( t   t  )] a }b 
chrf: H MOKwEMW (t )   log [ F MOKwEMW (t ) ]   log 

a b
1   {1  [1  exp ( t   t )] } 
3.8 The MOKw  Power Log-normal ( MOKw  PLn ) distribution
Power Log-normal distribution introduced by Nelson and Dognanksoy (1992) by extending the Lognormal distribution and its density and cdf are respectively given by,

p    ln(t )      ln(t ) 
g (t ) 

  

t  
  


p 1

p

    ln(t ) 
and G (t )  1   
 ,

  

     ,   0, p  0, t  0 .
For p  1 , it reduces to the usual Log-normal distribution.
With this as the base line distribution we get the pdf of MOKw  PLn distribution as

f MOKwPLn (t ) 
p    ln(t )      ln(t )  
 ab

  

t  
  

a
 
p 
1  1      ln(t )   

    
  

 


b 1

p 1 

   ln(t ) 
 1   



  




1  



p

a 1





a
 
p 
1  1      ln(t )   

    
  

 


 
a b 
p



    ln(t ) 
 
cdf: F MOKwPLn (t )  1  1  1   
   
 
     
  
 




1  



b

2






 
p a




ln(
t
)

  
 

1  1    
  

 
 

b

sf: F MOKwPLn (t )



 



a
 
p 



ln(
t
)

  
1  1   

    
  

 
 

b



1  








a
 
p 



ln(
t
)

  
1  1   

    
  

 
 

b






hrf: h MOKwPLn (t )

ab


p    ln(t )      ln(t )  

  

t  
  


p  a 1 

p 1 

   ln(t ) 
1   
 
   
 



 
p a
    ln(t )   

1   1  1   
  
   
 
 

13

p

1  1      ln(t )  

    
 

  
b

a





1






rhrf:

r MOKwPLn (t )

 a b p    ln(t )      ln(t ) 

  

t
 
  



p 1

     ln(t )  p 
1   
 
 
   

a 1

p a
 

1  1      ln(t )   
    
  
 


p a
 




ln(
t
)



1   1  1   
  
    
   
 


b 1

b

b
 
a 
p



    ln(t )  
 
 1  1  1   
   
 
     
  
 


1

chrf: H MOKwPLn (t )   log [ F MOKwPLn (t ) ]


  log 


p a
 

1  1      ln(t )   
    
  
 


b

p a
 
    ln(t )   

1   1  1   
 
    
  
 


b






for p  1 , it reduces to the MOKw -Log-normal distribution.
3.9 The MOKw  Extended Exponentiated Pareto ( MOKw  EEP ) distribution
The pdf and cdf of the exponentiated Pareto distribution, of Nadarajah (2005) are given respectively
by

g (t )   k  k t ( k 1) [1  ( t ) k ]  1 and G (t )  [1  ( t ) k ]  , x   ;  , k ,   0 .
Thus the pdf of MOKw  EEP distribution is given by

f MOKwEEP (t ) 

 a b  k  k t  ( k 1) {1  ( t ) k }  1[{1  ( t ) k }  ] a 1[1  [{1  ( t )k }  ] a ] b 1
[1   [1  {[1  ( t )k ]  } a ] b ] 2

cdf: F MOKwEEP (t ) 

1  [1  {[1  ( t ) k ]  }a ] b
 [1  {[1  ( t ) k ]  }a ] b
MOKwEEP
sf:
F
(t
)

1   [1  {[1  ( t ) k ]  }a ] b
1   [1  {[1  ( t ) k ]  }a ] b

hrf: h MOKwEEP (t ) 

a b  k  k t ( k 1) [1  ( t ) k ]  1 [1  ( t ) k ]  ( a1) [1  [1  ( t ) k ]  a ]1
1   [1  [1  ( t ) k ]  a ] b

rhrf: r MOKwEEP (t ) 

 a b  k  k t  ( k 1) [1  ( t ) k ]  1 [1  ( t ) k ]  ( a 1) [1  [1  ( t ) k ]  a ] b 1
1   [1 [1  ( t ) k ]  a ] b
 [1  [1 {1  ( t ) k }  a ] b ] 1

chrf:

  [1  {[1  ( t ) k ]  }a ] b 
H MOKwEEP (t ) =  log 
k  a b
1   [1  {[1  ( t ) ] } ] 

3.10 The MOKw  Extended Power ( MOKw  EP ) distribution

14

The cdf and pdf of the extended power distribution are G (t )  ( t ) k and g (t )  k  k t k 1 , t  (0,1  )
and   0 . The corresponding pdf of MOKw  EP distribution is then given by
pdf:

f MOKwEP (t ) 

cdf: F MOKwEP (t ) 

 a b k  k t k 1[( t ) k ] a1`[1  [ ( t ) k ] a ] b1
[1   {1  [ ( t ) k ] a } b ] 2

1  [1  [( t ) k ] a ] b
 [1  [( t ) k ] a ] b
MOKwEP
,
sf
:
F
(t
)

1   [1  [( t ) k ] a ] b
1   [1  [( t ) k ] a ] b

rhf:

h MOKwEP (t ) 

a b k k t k 1[( t )k ] a 1`[1  [( t ) k ] a ] 1
1   [1  [( t ) k ] a ] b

rhrf:

r MOKwEP (t ) 

 a b k k t k 1[( t )k ] a 1`[1  [( t ) k ] a ] b 1 [1  {1  [( t ) k ] a } b ] 1
1   [1  [( t ) k ] a ] b

chrf:

  [1  [( t ) k ] a ] b 
H MOKwEP (t )   log 
k a b
1   [1  [( t ) ] ] 

4. General results for the Marshall-Olkin’s Kumaraswamy-G ( MOKw  G ) family of
distributions
In this section we derive certain general results for the proposed MOKw  G family following the
methods described in Barreto-Souza et al. (2013), Cordeiro et al. (2014b), and Alizadeh et al. (2015).
4.1 Series expansions


( k  j ) j
z .
j  0 ( k ) j!

Consider the series representation (1  z ) k  

(10)

This is valid for z  1 and k  0, where  (.) is the gamma function. If   (0,1) using (10) in (6), we
obtain

f

MOKwG

(t )   a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1[1   [1  G (t ) a ] b ]  2


  a b g (t ) G (t ) a 1[1  G (t ) a ] b 1  ( j  1)[ [1  G (t ) a ] b ] j
j 0


  a b g (t ) G (t ) a 1[1  G (t ) a ] b 1  ( j  1) (1   ) j {[1  G (t ) a ] b } j
j 0



 f KwG (t )  j [ F KwG (t ) ] j

(11)

j 0






j0

/
j

d
[F
dt

KwG

(t ) ]

j 1

15

Where  j   j ( )  ( j  1) (1   ) j ;  /j   /j ( )   (1   ) j and
F KwG (t )  1  F KwG (t )  [1  G (t ) a ] b is the sf of Kw  G distribution.

MOKwG

Again, f

(t )   a b g (t ) G (t ) a 1[1  G (t ) a ] b 1  ( j  1) (1   ) j [1  [1  {[1  G (t ) a ] b } j ]]
j 0

j

 j
  a b g (t ) G (t ) a 1[1  G (t ) a ] b 1  ( j  1) (1   ) j    [  {1  [1  G (t ) a ] b } ] j  k
j 0
k 0  k 
j

 j
  a b g (t ) G (t ) a 1[1  G (t ) a ] b 1  ( j  1) (1   ) j (1) j  k    [ {1  [1  G (t ) a ] b } ] j  k
j 0
k 0  k 
j

 j
 f KwG (t )    (1   ) j (1) j  k ( j  1)  [ F KwG (t ) ] j  k
j  0 k 0
k 
j



 f KwG (t )    j , k [ F KwG (t ) ] j  k

(12)

j  0 k 0





j



j 0 k 0




 j , k d KwG
[ F (t ) ] j  k 1
j  k  1 dt

j



/
j ,k

j 0 k 0

d
[ F KwG (t ) ] j  k 1
dt

(13)

j
 j  1

Where  j , k   j , k ( )   (1   ) j (1) j  k ( j  1)  ;  /j , k   /j , k ( )   (1   ) j (1) j  k 
 k
 k 
and F KwG (t )  1  [1  G (t ) a ] b is the cdf of Kw  G (Cordeiro and de Castro 2011) distribution in
equation (4). Similarly an expansion for the survival function of MOKw  G [for   (0,1) ] can be
derives as

F MOKwG (t )   [1  G (t )a ]b 1   [1  G(t ) a ] b
  [1  G (t ) a ] b {1   [1  G (t ) a ] b } 1


  [1  G (t ) a ] b 
j 0

( j  1)
{ [1  G (t ) a ] b } j
(1) j !



  F KwG (t ; a, b)  (1   ) j {F KwG (t; a, b)} j
j 0



   j [ F KwG (t ; a, b)] j 1
j 0


    j F KwG (t; a, b( j  1))
j 0

16

The density function (6) can also be expressed as
f

MOKwG

(t ) 

a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1
 (  1) [1  [1  G (t ) a ] b ] 
 1 





2

Hence for   1 using (10) we get

a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1  (  1)[1  [1  G (t ) a ] b ] 

1 






  1 
 a b g (t ) G (t ) a 1 [ 1  G (t ) a ] b 1  ( j  1) 

  
j 0

j

2

1  [ 1  G(t ) ] 
b

j


  1
a 1
a b 1
a b
 
 ( j  1)a b g (t ) G (t ) [1  G (t ) ] 1  [1  G (t ) ]


j 0 



j



a







 f KwG (t )  j F KwG (t )

j

(14)

j 0



   /j
j0

d
F KwG (t )
dt





j 1

(15)

j

where  j   j ( ) 

j 1 1 
1
1
/
/
1   and  j   j ( )  1  
 
  

j

For   1 , the survival function of MOKw  G can be expressed as

F MOKwG (t ) 

[1  G (t ) a ] b
1  (  1)[1  [1  G (t ) a ] b ] / 

On using (10) we get

 (  1)[1  [1  G (t ) a ] b ] 
 [1  G (t ) ] 1 




a

1

b



 [1  G (t ) a ] b


j 0

( j  1) (  1) j
[1  [1  G (t ) a ] b ] j
(1) j !  j
j


 
 F KwG (t ; a, b )  ( 1) j   [ F KwG (t; a, b) ] j
 
j 0



 F KwG (t; a, b)  C j [ F KwG (t ; a, b) ] j where C j  ( j  1)  j .
j 0

4.2 Order statistics

17

j

Suppose T1, T 2 ,...T n is a random sample from any MOKw  G distribution. Let T i:n denote the i th order
statistics. The pdf of T i :n can be expressed as

f i:n (t ) 

n!
f MOKwG (t )[1  F MOKwG (t )] i 1 F MOKwG (t ) n  i
(i  1)! (n  i )!

i 1
 i  1
n!
 [ F MOKwG (t )] p
f MOKwG (t ) F MOKwG (t ) ni  
p
(i  1)! (n  i )!
p 0 

i 1
i

1


n!
 F MOKwG (t ) n  p i

f MOKwG (t ) (1) p 
p
(i  1)! (n  i )!
p 0





Now using the general expansion of the MOKw  G distribution pdf and sf we get the pdf of the
i th order statistics for of the MOKw  G for   (0,1) as
f i:n (t ) 

 KwG

n!
f
(
t
)
 j [ F KwG (t ) ] j 


(i  1)! (n  i)! 
j 0


i 1



p 0


 i  1 
  q [ F KwG (t ; a, b)] q  n  p 1 i 
(1) p 1 
 p  q  0


where  j and  q are defined earlier.


i1
n!

(i  1)! ( n  i)! p o

 i  1

 ( 1) p1 f
p




KwG



(t; a, b)   j  q [ F KwG (t ; a, b)] j q n  p 1 i
j  0 q 0



 f KwG (t; a, b)  X j , q [ F KwG (t ; a, b )] j  q  n  p 1  i

(16)

j,q 0



   [ X j , q /( j  q  n  p  i  2) ]
j,q 0





 X

j, q

j,q 0




 X

j, q

j,q 0

d
[ F KwG (t ; a, b)] j  q  n  p  i  2
dt

d
[ F KwG (t; a, b)] j  q  n  p  i  2)
dt
d
[ F KwG (t ; a , b ( j  q  n  p  i  2))]
dt





 X

j, q

f KwG (t ; a , b ( j  q  n  p  i  2))

(17)

j,q 0

 n  1 i 1
 
Where X j , q  n  j  q 
 i  1  p0

 i  1

 (1) p 1 and X j , q   X j , q /( j  q  n  p  i  2) .
 p

Again using the general expansion of the pdf and sf of MOKw  G distribution we get the pdf of the i th
order statistics for of the MOKw  G for   1 as

18

f i:n (t ) 


i 1
 KwG
j
n!
(t ; a, b )  j F KwG (t ; a, b )  (1) p
f
(i  1)! (n  i )! 
j 0
 p 0





 i  1


 p 


 KwG

[
F
(
t
;
a
,
b
)]
Ck [ F KwG (t ; a, b) ] k 


k 0






i 1
 KwG
j
n!
KwG
p
f
(
t
;
a
,
b
)

F
(
t
;
a
,
b
)

 (1)

j
(i  1)! (n  i )! 
j 0
 p0



[F

KwG



(t ; a, b )]

( n  p i )

n  p i

 i  1


 p 

 
KwG
k
  Ck [ F (t ; a , b) ] 
 k 0


n p i

Where  j and Ck defined in section 4.1

 
KwG
k
  Ck [ F (t ; a , b) ] 
 k 0


Now,

Where

d n  p i , k 

n pi



  d n  p  i , k [ F KwG (t ; a , b) ] k

(Nadarajah et. al 2015)

k0

1 k
 [h (n  p  i  1)  k ]Ch d n  p i , k h
k C0 h 1

The density function of the i th order statistics of MOKw  G distribution can be expressed as
f i:n (t ) 


i 1
 KwG
j
n!
(t ; a , b)  j F KwG (t ; a , b)  (1) p
f
(i  1)! (n  i)! 
j0
 p0





 i  1


 p 



[ F KwG (t ; a, b)] ( n  p i )  d n  p  i , k [ F KwG (t ; a, b) ] k
k 0



i 1
n!
 (1) p
(i  1)! (n  i)! p  0

 i  1 KwG

 f
(t ; a, b )[ F KwG (t ; a , b)] ( n  p  i )
p






j





d n  p i , k F KwG (t; a, b)

j k

j,k 0


 f KwG (t; a, b) [ F KwG (t; a, b)] ( n pi )

  F

KwG

j ,k



(t ; a, b)

jk

(18)

j ,k 0

 n  1 i 1
  (1) p
Where  j , k  n  j d n  p  i , k 
i

1

 p 0

 i  1


 p

Remark I: Equations (11)-(18) reveal that the density functions of the MOKw  G distribution and
that of its order statistics can be expressed as a product of the baseline density f (x) with an infinite
power series of G (.) and also as a mixture of exponentiated-G distributions under Lehman alternatives.
These results play important role and may be used to obtain explicit expressions for the moments and

19

moment generating function (mgf) of the MOKw  G distribution and of its order statistics in a general

framework and also for special models using the corresponding results of exponentiated - G distributions.
4.3 Probability weighted moments
The probability weighted moments (PWMs), first proposed by Greenwood et al. (1979), are
expectations of certain functions of a random variable whose mean exists. The ( p, q, r )th PWM of T is
having a base line cdf F (t ) defined by


p , q , r   t p [ F (t ) ] q [1  F (t ) ] r f (t ) dt


From equations (11), (12) and (14) the s th moment of T for   (0,1) and   1, can be written either
as




s
a 1
a b 1
a bj
 t a b g (t ) G(t ) [1  G(t ) ]  j [1  G(t ) ] dt

E (T s ) 

j 0







   j  t s [1  G (t ) a ] b j a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1 dt
j 0





   j s , 0, j or
j 0




j

s
a 1
a b 1
a b jk
 t a b g (t ) G(t ) [1  G (t ) ]   j , k [1  [1  G(t ) ] ] dt

E (T s ) 

j 0 k  0







j





  j, k

j 0 k 0


s
a b jk
a 1
a b 1
 t [1  [1  G (t ) ] ] a b g (t ) G (t ) [1  G (t ) ] dt



j

    j ,k s , jk , 0 or as
j 0 k 0



E (T s ) 





t s a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1   j [1  [1  G (t ) a ] b ] j dt
j 0




 j
j 0



t

s

[1  [1  G (t ) a ] b ] j a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1 dt





   j s , j , 0
j 0



Where p, q , r   t p {1  [1  G (t ) a ] b } q {[1  G (t ) a ] b } r [ a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1 ] dt


is the PWM of MOKw  G the for the baseline distribution G .

20

The PWM can generally be used for estimating parameters quantiles of generalized
distributions. These moments have low variance and don’t possess severe biases, and they compare
favourably with estimators obtained by maximum likelihood (Alizadeh et al., 2015).
Proceeding as above we can derive s th moment of the i th order statistic Ti:n , in a random
sample of size n from MOKw  G for   (0,1) and   1, on using equations (16) and (18), as


E (T s i ,n ) 



 X j , q  s , 0,

j  q n  p i

and E (T s i, n ) 



j,q 0

 s,

j,k

j  k , n pi

respectively. Where the quantities

j ,k 0

 j ,k ,  j ,  j , X j , q and  j ,k are defined in section 4.1 and 4.2.
4.4 Moment generating function
The moment generating function of MOKw  G family can be easily expressed in terms of those of the
exponentiated Kw  G (Cordeiro and de Castro, 2011) distribution using the results of section 4.1. For
example using equation (15) it can be seen that






M T ( s )  E[e st ]   e st f (t ) dt   e st   /j





/

   j  est
j 0





d
F KwG (t )
dt





j 1

j 0

d
F KwG (t )
dt





j 1

dt



dt    /j M X ( s )
j 0

Where X follows exponentiated Kw  G (Cordeiro and de Castro, 2011) distribution.
4.5 Renyi Entropy
The entropy of a random variable is a measure of uncertainty. The Rényi entropy is defined as



I R ( )  (1   ) log  f (t )  dt 
 

1

where   0 and   1 For furthers details, see Song (2001). For   (0,1) using expansion (10) in (6)
we can write



f (t )   a b g (t ) G (t ) a 1[1  G (t ) a ] b 1 [1   {1  G (t ) a } b ] 2



  [ a b g (t ) G (t )a 1[1  G (t )a ]b 1 ] 
 (2 )







 (1   ) j (2  j )
j 0

[{1  G (t ) a }b ] j
j!

Thus for   (0,1) the Rényi entropy of T can be obtained as

 


I R ( )  (1   ) 1 log   q j  a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1 [{1  G (t ) a } b ] j dt 
 j 0 




Where q j  q j ( ) 



  (1   ) j (2  j )
 (2 ) j!

21

The density function (6) can be expressed as

f (t )

MOKwG

 a b g (t ) G (t )

a 1

 (  1)[1  [1  G (t ) a ] b ] 
 1 





a b 1

[1  G (t ) ]

2

Therefore
2

 (  1)[1  {1  G (t )a } b ]  
f (t )  a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1  1 
 



 





For   1, using expansion (10) we get,
f (t ) 

[a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1 ] 
   j (2 )



 (  1) j (2  j )
j 0

{1  [1  G (t ) a ] b } j
j!

Thus for   1, the Rényi entropy of T also can be obtained as

 

I R ( )  (1   ) 1 log   r j  [ a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1 ]  [1  [1  G (t ) a ] b ] j dt 
 j 0 


Where rj  rj ( ) 

(  1) j  (2  j )
   j (2 ) j!

4.6 Quantile function and random sample generation
We shall now present a formula for generating MOKw  G random variable by using inversion
method by inverting the cdf or the survival function.

F MO (t ) 

F (t )
 F (t ) [1   F MO (t )]   F MO (t )
   F (t )

 F (t ) 

MO
 F MO (t )
(t ) 
1   F

t

F


MO
MO
1   F (t )
1   F ( t ) 

(19)

Now for the Kw  G we have

F (t )  [1  G (t ) a ] b . Therefore
log F (t )  b log [1  G (t ) a ]  1  G (t ) a  F (t )
1

 G (t )  [1  {1  F (t ) } b ]

1

a

1

b

1
 a log G (t )  log 1  F (t ) b 



1

 t  G 1 [1  {1  F (t ) } b ]

1

a

(20)

Now combining (19) and (20) we get for MOKw  G
1
 
 u  b 
1

t  G 1  1 

  1 u  



1

a

(21)

Where u ~ U (0,1)

22

The pth Quantile t p for MOKw  G can be easily obtained from (21) as
1
 
 p  b 
1

t p  G 1  1 

  1 p 



1

a

For example, let the G be exponential distribution with parameter   0, having pdf and cdf as
g (t :  )   e  t , t  0 and G (t :  )  1  e t , respectively. Then the p th quantile is obtained as

 (1 /  ) log [1  p ] . Therefore, the p th quantile t p , of MOKw  E is given by
1
 
1
a
b


1

p
 
 
 
t p   log 1  1  1 

1   p   


 
 


4.7 Asymptotes and shapes
Here we investigate the asymptotic shapes of the proposed family following the method followed in
Alizadeh et al., (2015).
Proposition 1. The asymptotes of equations (6), (7) and (9) as t  0 are given by
f (t ) ~

a b g (t ) G (t ) a 1


F (t ) ~ 0

h (t ) ~

a b g (t ) G (t ) a 1


as G (t )  0
as G (t )  0
as G (t )  0

Proposition 2. The asymptotes of equations (6), (7) and (9) as t   are given by
f (t ) ~ a b g (t )[1  G (t ) a ] b 1 as t  
F (t ) ~ 1  [1  G (t ) a ] b

as t  

h(t ) ~ a b g (t )[1  G (t ) a ] 1

as t  

The shapes of the density and hazard rate functions can be described analytically. The critical points of
the MOKw  G density function are the roots of the equation:

d log [ f (t )]
dt


g (t )
g (t )
g (t )G (t ) a 1
 a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1
 (a  1)
 a (1  b)

2
g (t )
G (t )
1  G (t ) a
1   [1  G (t )a ] b

23

(22)

The number of roots of (22) may be more than one. In particular, if t  t0 is a root of (22) then it
corresponds to a local maximum, a local minimum or a point of inflexion depending on whether

 (t0 )  0,  (t0 )  0 or  (t0 )  0 where
 (t ) 

d2
log [ f (t )]
dt 2

 (t ) 

g (t ) g (t )  [ g (t )] 2
G (t ) g (t )  g (t ) 2

(
a

1
)
g (t ) 2
G (t )

 a ( 1  b) [



g (t ) G (t ) a 1 (a  1) g (t ) 2 G(t ) a  2 a g (t ) 2 G (t ) 2 a  2


]
1  G (t )a
1  G (t )a
[1  G (t )a ] 2

2  a b [ g (t ) G (t ) a 1 [1  G (t ) a ] b 1 ] 2  a b [ (a  1) g (t ) 2 G (t ) a  2 [ 1  G (t ) a ] b 1 ]

1   [1  G (t ) a ] b
1   [1  G (t ) a ] b

2  a b [ a (b  1) g (t ) 2 G (t ) 2 a  2 [1  G (t ) a ] b  2 ]   a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1 



1   [1  G (t ) a ] b
1   [1  G (t ) a ] b




2

g (t ) g (t )  [ g (t )] 2
G (t ) g (t )  g (t ) 2

(
a

1
)
g (t ) 2
G (t )

 a ( 1  b) [

g (t ) G (t ) a 1 (a  1) g (t ) 2 G(t ) a  2 a g (t ) 2 G (t ) 2 a  2
2 g (t )
[1  G (t ) a ] b h MOKwG (t )


]
a
a
a 2
g (t )
1  G (t )
1  G (t )
[1  G (t ) ]

 2  ( a  1) g (t ) G (t ) 1 [ 1  G (t ) a ] b h MOKwG (t )  2  a (b  1) g (t ) G(t ) a 1 (t ) {1  G (t ) a } b 1 h MOKwG (t )
 {[1  G (t ) a ] b  h MOKwG (t )} 2

The critical points of h(t ) are the roots of the equation
g (t ) a g (t ) G (t )a 1  a b g (t ) G (t )a 1 [1  G (t ) a ] b 1
d log [h (t )] g (t )

 (a  1)


dt
g (t )
G (t )
1  G (t ) a
1   [1  G(t )a ] b

(23)

The number of roots of (23) may be more than one. In particular, ff t  t0 is a root of (23) then it
corresponds to a local maximum, a local minimum or a point of inflexion depending on whether

 (t0 )  0,  (t0 )  0 or  (t0 )  0 where  (t ) 

 (t ) 

g (t ) g (t )  [ g (t )] 2
G (t ) g (t )  g (t ) 2

(
a

1
)
g (t ) 2
G (t )

a [



d2
log [h(t )]
dt 2

g (t ) G (t ) a 1 (a  1) g (t ) 2 G (t ) a  2 a g (t ) 2 G (t ) 2 a  2


]
1  G (t )a
1  G (t ) a
[1  G (t ) a ] 2

 a b [ g (t ) G (t ) a 1 [1  G (t ) a ] b 1 ]  a b [ (a  1) g (t ) 2 G (t ) a  2 [1  G (t ) a ] b 1 ]

1   [1  G (t ) a ] b
1   [1  G (t ) a ] b

24

 a b [ a (b  1) g (t ) 2 G (t ) 2 a  2 [ 1  G (t ) a ] b  2 ]  a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1 



1   [1  G (t ) a ] b
1   [1  G (t ) a ] b


g (t ) g (t )  [ g (t )] 2
G (t ) g (t )  g (t ) 2

(
a

1
)
g (t ) 2
G (t )



a [

2

g (t ) G (t ) a 1 (a  1) g (t ) 2 G (t ) a  2 a g (t ) 2 G(t ) 2 a  2
 g (t )


]
[1  G (t ) a ] b h MOKwG (t )
a
a
a 2
1  G (t )
1  G (t )
[1  G (t ) ]
g (t )

  ( a  1) g (t ) G (t ) a 1 (t ) [ 1  G (t ) a ] b h MOKwG (t )   a (b  1) g (t ) G (t ) a 1 (t ) {1  G (t ) a }b 1 h MOKwG (t )

 {[1  G (t ) a ] b  hMOKwG (t )} 2
4.8 Stochastic orderings
In

this

section

we

study

the

reliability

properties

and

stochastic

ordering

of

the

MOKw  G distributions Stochastic ordering properties have applications in diverse fields such as

economics, reliability, survival analysis, insurance, finance, actuarial and management sciences
(Shaked and Shanthikumar, 2007).
Let X and Y be two random variables with cfds F and G, respectively, survival functions F  1  F
and G  1  G , and corresponding pdf’s f, g. Then X is said to be smaller than Y in the likelihood ratio
order ( X  lr Y ) if f (t ) g (t ) is decreasing in t  0 ; stochastic order ( X  st Y ) if F (t )  G (t ) for
all t  0 ; hazard rate order ( X  hr Y ) if F (t ) G (t ) is decreasing in t  0 ; reversed hazard rate order
( X  rhr Y ) if F (t ) G (t ) is decreasing in t  0 . These four stochastic orders are related to each other,
as

X  rhr Y  X lr Y  X  hr Y  X  st Y

(24)

Let X ~ MOKwG ( 1 , a, b) and Y ~ MOKwG ( 2 , a, b) .

Theorem1:

If 1   2 ,

then X lr Y

( X  hr Y , X  rhr Y , X  st Y ) .

Proof:

f (t ) 1 a b g (t ) G (t ) a 1 [ 1  G (t ) a ] b 1 [1  1 [1  G(t ) a ] b ] 2

g (t )  2 a b g (t ) G (t ) a 1 [ 1  G (t ) a ] b 1 [1   2 [1  G (t ) a ] b ] 2


 1
2

1   2 [ 1  G (t ) a ] b 

a b 
 1  1 [1  G(t ) ] 

2

Since 1   2
d
dt



 f (t ) 
1 1   2 [1  G (t ) a ] b

2


 2 1  1 [ 1  G (t ) a ] b
 g (t ) 

{1  1 [ 1  G (t ) a ] b }[ 2 a b [ 1  G (t ) a ] b 1 G (t ) a 1 g (t )] {1   2 [1  G(t ) a ] b }[1a b [ 1  G (t ) a ] b 1 G(t ) a 1 g (t )]
{1  1 [ 1  G(t ) a ] b } 2
21{1   2 [ 1  G (t ) a ] b } a b [ 1  G (t ) a ] b 1 G (t ) a 1 g (t )[{1  1 [ 1  G (t ) a ] b } 2  {1   2 [ 1  G (t ) a ] b }1 ]

 2 {1  1 [ 1  G (t ) a ] b } 3

25



21{1   2 [ 1  G (t ) a ] b } a b [ 1  G (t ) a ] b 1 G (t ) a 1 g (t ) ( 2  1 )
 2 {1  1 [ 1  G (t ) a ] b } 3



21 (1   2 ){1   2 [ 1  G (t ) a ] b } a b [ 1  G (t ) a ] b 1 G (t ) a 1 g (t ) . Which is always less than 0.
 2 {1  1 [ 1  G (t ) a ] b } 3

Hence, f (t ) / g (t ) is decreasing in t. That is X lr Y .The remaining statements follow from the
implications (24).
5. Estimation
5.1 Maximum likelihood estimation for MOKw  G
The model parameters of the MOKw  G distribution can be estimated by maximum likelihood. Let

t  (t1 , t 2 ,...t n ) T be

a

random

sample

of

size

from

n

MOKw  G

with

parameter

vector θ  ( , a, b, βT ) T , where β  ( 1 ,  2 ,...  q ) T corresponds to the parameter vector of the baseline
distribution G. Then the log-likelihood function for θ is given by

   (θ )
n

n

 n log    log [ a b g (ti , β) G (ti , β) a 1[1  G (ti , β) a ] b 1 ]  2  log [1   [1  G (ti , β) a ] b ]
i 1

i 1

n

n

n

 n log   n log (a b )   log [ g (ti , β) ]  (a  1)  log [G (ti , β)]  (b  1)  log [1  G(ti , β) a ]
i 0

i 0

i 0

n

 2 log [1   [1  G (ti , β) a ] b ]

(25)

i 1

This log-likelihood function can not be solved analytically because of its complex form but it can be
maximized numerically by employing global optimization methods available with softwares like R,
Mathematica.
By taking the partial derivatives of the log-likelihood function with respect to  , a, b and β
components of the score vector U θ  (U  , U a , U b , U  T ) T can be obtained as follows:
U 

n

n
[1  G (ti , β) a ] b
  2
a b
 
i  0 1   [1  G (ti , β ) ]

Ua 

n
 n n
G (ti , β)a log [G (ti , β)]
   log [G (ti , β)]  (1  b) 
a a i  0
1  G (ti , β)a
i 0

b  [1  G (ti , β) a ] b 1 G (ti , β) a log [G (ti , β)]
1   [1  G (ti , β) a ] b
i 0
n

 2

Ub 

 n
 
b b

n

n

 log [1  G (ti , β)a ]  2 
i 0

i 0

 [1  G (ti , β) a ] b log [1  G (ti , β) a ]
1   [1  G (ti , β) a ] b

26

U 

n
n
n

g (β ) (ti , β)
G ( β ) (ti , β)
a G (ti , β) a 1 G ( β ) (ti , β)

 (a  1) 
 (1  b) 
β i  0 g (ti , β)
1  G (ti , β) a
i  0 G (ti , β )
i 0
n

 2
i 0

b  [1  G (ti , β) a ] b1 a G (ti , β) a1 G ( β ) (ti , β)
1   [1  G (ti , β) a ] b

Solving the equations U  (U  ,U a , U b , U  T ) T  0 simultaneously gives the maximum likelihood
estimate (MLE) ˆ  (ˆ , aˆ , bˆ, βˆ T ) T of   ( , a, b, β T ) T .
5.2 Asymptotic standard error and confidence interval for the mles
The asymptotic variance-covariance matrix of the MLEs of parameters can obtained by inverting the
Fisher information matrix I (θ) which can be derived using the second partial derivatives of the loglikelihood function with respect to each parameter. The i j th elements of I n (θ) are given by

  2 l (θ) 
 ,
I i j  E 
   i  j 



i, j  1, 2, , 3  q

The exact evaluation of the above expectations may be cumbersome. In practice one can estimate

I n (θ) by the observed Fisher’s information matrix Î n (θˆ ) is defined as:
  2 l (θ) 

Î i j   
    
i
j



,

i, j  1, 2,, 3  q

θ  θˆ

Using the general theory of MLEs under some regularity conditions on the parameters as n   the
asymptotic distribution of

n (θˆ  θ) is N k (0,Vn ) where Vn  (v jj )  I n 1 (θ) . The asymptotic

behaviour remains valid if Vn is replaced by Vˆn  Î 1 (θˆ ) . This result can be used to provide large
sample standard errors and also construct confidence intervals for the model parameters. Thus an
approximate standard error and (1   / 2) 100% confidence interval for the mle of jth parameter  j are
respectively given by

v̂ j j and ˆ j  Z  / 2 vˆ j j , where Z  / 2 is the  / 2 point of standard normal

distribution.
As an illustration on the MLE method its large sample standard errors, confidence interval in
the case of MOKw  E ( , a , b,  ) is discussed in an appendix.
5.3 Estimation by method of moments
Here an alternative method to estimation of the parameters is discussed. Since the moments are not in
closed form, the estimation by the usual method of moments may not be tractable. Therefore following
(Barreto-Souzai et al., 2013) we get

27



a b

E [1   {1   [1  G (t ) ] }] 

a b
 1   [1  G (t ) ]

 a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1

1   [1  G(t ) ] 

a b 2



dt




Let, u  1   [1  G (t ) a ]b
1

 a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1
dt

1   [1  G (t ) a ]b


then

du   a b g (t ) G (t ) a 1 [1  G (t ) a ] b 1 dt

 1



du   log u  1   log   
log 
u


(1   )


 



E [ [ 1   { 1   [1  G ( t ) a ] b } ] v ] 

a 1

 a b g (t ) G (t ) [1  G (t )
 1   [1  G (t ) ] 
[1   [1  G (t ) ] ]
a b v

a b

a b 1

]

2

dt






 {1   [1  G (t )

] } v  2  a b g (t ) G (t ) a 1[1  G (t ) a ] b 1 dt

a b



1





1

 v 2
  u v 1 
 (1   v 1 )
u du 




  v  1 
 (v  1)

  log ( ) 1   , v  1
E [ [1   {1   [1  G (t )a ]b } ]v ]  
v 1
 (1   )  (v  1) , v  (2,3,...)

(28)

One can use (28) to give a new method of estimation i.e. for a random sample t1 , t 2 ,...t n from a
population with survival function (8) the model parameters can be estimated from the equations
  log ( ) 1   , v  1
1 n
[[1   1   [1  G (ti ) a ] b ] v ]  

v 1
n i 1
 (1   )  (v  1) , v  2,3,...q  1





6. Real life applications
In this subsection, we consider fitting of two real data sets to show that the distributions from the
proposed MOKw  G distribution can be a better model than the corresponding distributions from
KwMO  G (Alizadeh et al., 2015) by considering the Frechet and the Exponential distribution as our
G . Here the parameters are estimated by numerical maximization of log-likelihood function and their

asymptotic standard errors and 95% confidence intervals are computed using large sample approach.
In order to compare the distributions, we have considered known criteria like AIC (Akaike
Information Criterion), BIC (Bayesian Information Criterion), CAIC (Consistent Akaike Information
Criterion) and HQIC (Hannan-Quinn Information Criterion)
It may be noted that AIC  2 k  2 l ; BIC  k log (n)  2 l ; CAIC  AIC  (2 k (k  1)) /(n  k  1) ;
and HQIC  2 k log [log (n )]  2l , where k is the number of parameters in the statistical model, n the

28

sample size and l is the maximized value of the log-likelihood function under the considered model.
We have used the maximum likelihood method for estimating the model parameters.
Example I:
Here we work with the following data about 346 nicotine measurements made from several brands of
cigarettes in 1998. The data have been collected by the Federal Trade Commission which is an
independent agency of the US government, whose main mission is the promotion of consumer
protection. [http: //www.ftc.gov/ reports/tobacco or http: // pw1.netcom.com/ rdavis2/ smoke. html.]

{1.3, 1.0, 1.2, 0.9, 1.1, 0.8, 0.5, 1.0, 0.7, 0.5, 1.7, 1.1, 0.8, 0.5, 1.2, 0.8, 1.1, 0.9, 1.2, 0.9, 0.8, 0.6, 0.3,
0.8, 0.6, 0.4, 1.1, 1.1, 0.2, 0.8, 0.5, 1.1, 0.1, 0.8, 1.7, 1.0, 0.8, 1.0, 0.8, 1.0, 0.2, 0.8, 0.4, 1.0, 0.2, 0.8,
1.4, 0.8, 0.5, 1.1, 0.9, 1.3, 0.9, 0.4, 1.4, 0.9, 0.5, 1.7, 0.9, 0.8, 0.8, 1.2, 0.9, 0.8, 0.5, 1.0, 0.6, 0.1, 0.2,
0.5, 0.1, 0.1, 0.9, 0.6, 0.9, 0.6, 1.2, 1.5, 1.1, 1.4, 1.2, 1.7, 1.4, 1.0, 0.7, 0.4, 0.9, 0.7, 0.8, 0.7, 0.4, 0.9,
0.6, 0.4, 1.2, 2.0, 0.7, 0.5, 0.9, 0.5, 0.9, 0.7, 0.9, 0.7, 0.4, 1.0, 0.7, 0.9, 0.7, 0.5, 1.3, 0.9, 0.8, 1.0, 0.7,
0.7, 0.6, 0.8, 1.1, 0.9, 0.9, 0.8, 0.8, 0.7, 0.7, 0.4, 0.5, 0.4, 0.9, 0.9, 0.7, 1.0, 1.0, 0.7, 1.3, 1.0, 1.1, 1.1,
0.9, 1.1, 0.8, 1.0, 0.7, 1.6, 0.8, 0.6, 0.8, 0.6, 1.2, 0.9, 0.6, 0.8, 1.0, 0.5, 0.8, 1.0, 1.1, 0.8, 0.8, 0.5, 1.1,
0.8, 0.9, 1.1, 0.8, 1.2, 1.1, 1.2, 1.1, 1.2, 0.2, 0.5, 0.7, 0.2, 0.5, 0.6, 0.1, 0.4, 0.6, 0.2, 0.5, 1.1, 0.8, 0.6,
1.1, 0.9, 0.6, 0.3, 0.9, 0.8, 0.8, 0.6, 0.4, 1.2, 1.3, 1.0, 0.6, 1.2, 0.9, 1.2, 0.9, 0.5, 0.8, 1.0, 0.7, 0.9, 1.0,
0.1, 0.2, 0.1, 0.1, 1.1, 1.0, 1.1, 0.7, 1.1, 0.7, 1.8, 1.2, 0.9, 1.7, 1.2, 1.3, 1.2, 0.9, 0.7, 0.7, 1.2, 1.0, 0.9,
1.6, 0.8, 0.8, 1.1, 1.1, 0.8, 0.6, 1.0, 0.8, 1.1, 0.8, 0.5, 1.5, 1.1, 0.8, 0.6, 1.1, 0.8, 1.1, 0.8, 1.5, 1.1, 0.8,
0.4, 1.0, 0.8, 1.4, 0.9, 0.9, 1.0, 0.9, 1.3, 0.8, 1.0, 0.5, 1.0, 0.7, 0.5, 1.4, 1.2, 0.9, 1.1, 0.9, 1.1, 1.0, 0.9,
1.2, 0.9, 1.2, 0.9, 0.5, 0.9, 0.7, 0.3, 1.0, 0.6, 1.0, 0.9, 1.0, 1.1, 0.8, 0.5, 1.1, 0.8, 1.2, 0.8, 0.5, 1.5, 1.5,
1.0, 0.8, 1.0, 0.5, 1.7, 0.3, 0.6, 0.6, 0.4, 0.5, 0.5, 0.7, 0.4, 0.5, 0.8, 0.5, 1.3, 0.9, 1.3, 0.9, 0.5, 1.2, 0.9,
1.1, 0.9, 0.5, 0.7, 0.5, 1.1, 1.1, 0.5, 0.8, 0.6, 1.2, 0.8, 0.4, 1.3, 0.8, 0.5, 1.2, 0.7, 0.5, 0.9, 1.3, 0.8, 1.2,
0.9}
In Table 1, the MLEs and standard errors (SEs) and 95% confidence intervals (in parentheses)
of the parameters from all the fitted distributions along with AIC, BIC, CAIC and HQIC values are
presented. According to the lowest values of the AIC, BIC, CAIC and HQIC, the MOKw  Fr
( MOKw  E ) distribution could be chosen as the better model than KwMO  Fr ( KwMO  E )
distribution respectively.
More information is provided for a visual comparison is presented in the form of a histogram
of the data with the fitted densities in Figures 3. The fitted cdf’s are also displayed in Figures 4. These
plots indicate that the proposed distributions provide a good fit to these data.

29

Table 1: MLEs, standard errors and 95% confidence intervals (in parentheses) and
the AIC, BIC, CAIC and HQIC values for the nicotine measurements data.
Parameters
â

b̂

̂
ˆ

̂
log-likelihood
AIC
BIC
CAIC
HQIC

KwMO  Fr

MOKw  Fr

KwMO  E

MOKw  E

10.462
(2.669)
(5.23, 15.69)
75.085
(19.512)
(36.84,113.33)
0.251
(0.014)
(0.22, 0.28)
16.270
(4.834)
(6.79, 25.75)

1.176
(0.497)
(0.20, 2.15)
126.659
(100.668)
(-70.65, 323.97)
0.303
(0.063)
(0.17, 0.43)
29.261
(23.497)
(-16.79,75.32)

1.434
(0.236)
(0.97, 1.80)
1.959
(0.779)
(0.43,3.49)
3.557
(0.677)
(2.23, 4.86)

1.511
(0.515)
(0.50, 2.52)
11.063
(35.759)
(-59.02, 81.15)
0.586
(1.572)
(-2.49, 3.67)

---

---

0.075
(0.013)
(0.04, 0.10)
-155.49
320.98
340.21
321.16
328.64

60.442
(85.106)
(-107.37, 227.25)
-110.27
230.54
249.77
230.72
238.19

21.214
(10.377)
(0.88, 41.56)
-107.89
223.78
239.17
223.89
229.91

20.96
(19.978)
(-18.19, 60.12)
-106.61
221.22
236.61
221.34
227.35

Estimated pdf's

Estimated pdf's

MOKw-Fr
KwMO-Fr

1.5
0.0

0.0

0.5

1.0

Density

1.0
0.5

Density

1.5

MOKw-E
KwMO-E

0.0

0.5

1.0

1.5

2.0

0.0

0.5

1.0

X

X

(a)

(b)

1.5

2.0

Fig: 3 Plots of the observed histogram and estimated pdf’s for the (a) MOKw  Fr , KwMO  Fr and
(b) MOKw  E , KwMO  E distributions for example I.

30

Estimated cdf's
1.0

1.0

Estimated cdf's
MOKw-E
KwMO-E

0.0

0.0

0.2

0.2

0.4

0.4

cdf

cdf

0.6

0.6

0.8

0.8

MOKw-Fr
KwMO-Fr

0.0

0.5

1.0

1.5

2.0

0.0

X

0.5

1.0

1.5

2.0

X

(a)

(b)

Fig: 4 Plots of the observed ogive and estimated cdf’s for the (a) MOKw  Fr , KwMO  Fr and (b)
MOKw  E , KwMO  E distributions for example I.

Example II:
This data set consists of 100 observations of breaking stress of carbon fibres (in Gba) given by Nichols
and Padgett (2006).

{3.70, 2.74, 2.73, 2.50, 3.60, 3.11, 3.27, 2.87, 1.47, 3.11, 4.42, 2.40, 3.15, 2.67,3.31, 2.81, 0.98, 5.56,
5.08, 0.39, 1.57, 3.19, 4.90, 2.93, 2.85, 2.77, 2.76, 1.73, 2.48, 3.68, 1.08, 3.22, 3.75, 3.22, 2.56, 2.17,
4.91, 1.59, 1.18, 2.48, 2.03, 1.69, 2.43, 3.39, 3.56, 2.83, 3.68, 2.00, 3.51, 0.85, 1.61, 3.28, 2.95, 2.81,
3.15, 1.92, 1.84, 1.22, 2.17, 1.61, 2.12, 3.09, 2.97, 4.20, 2.35, 1.41, 1.59, 1.12, 1.69, 2.79, 1.89, 1.87,
3.39, 3.33, 2.55, 3.68, 3.19, 1.71, 1.25, 4.70, 2.88, 2.96, 2.55, 2.59, 2.97, 1.57, 2.17, 4.38, 2.03, 2.82,
2.53, 3.31, 2.38, 1.36, 0.81, 1.17, 1.84, 1.80, 2.05, 3.65}.

In Table 2, the MLEs and standard errors (SEs) and 95% confidence intervals (in parentheses)
of the parameters from all the fitted distributions along with AIC, BIC, CAIC and HQIC values are
presented. According to the lowest values of the AIC, BIC, CAIC and HQIC, the MOKw  Fr
( MOKw  E ) distribution could be chosen as the better model than KwMO  Fr ( KwMO  E )
distribution respectively.

31

Table 2: MLEs, standard errors and 95% confidence intervals (in parentheses) and
the AIC, BIC, CAIC and HQIC values for the breaking stress of carbon fibres data
Parameters
â

b̂

̂
ˆ
̂
log-likelihood
AIC
BIC
CAIC
HQIC

KwMO  Fr
9.243
(2.958)
(3.45, 15.04)
19.203
(14.763)
(-9.73, 48.14)
0.953
(0.172)
(0.62, 1.29)
0.051
(0.026)
(0.00004, 0.10)
17.222
(7.881)
(1.78, 32.67)
-144.83
299.66
312.69
300.29
304.94

MOKw  Fr
6.777
(4.936)
(-2.89, 16.45)
38.279
(67.729)
(-94.47, 171.03)
0.522
(0.259)
(0.01, 1.03)
0.416
(0.500)
(-0.56, 1.39)
16.978
(31.897)
(-45.54, 79.49)
-141.63
293.26
306.29
293.89
298.54

KwMO  E
2.647
(2.066)
(-1.40, 6.69)
4.571
(14.461)
(-23.77, 32.91)
0.591
(1.103)
(-1.57, 0.75)

MOKw  E
3.226
(1.409)
(0.46, 5.99)
9.065
(42.865)
(-74.95, 93.08)
0.295
(0.834)
(-1.34, 0.93)

---

---

3.899
(10.290)
(-16.27, 24.07)
-141.25
290.50
300.92
290.92
294.73

2.566
(5.709)
(-8.62, 13.76)
-141.09
290.18
300.60
290.60
294.40

Like for the last example here also we have presented the histogram of the observed data with
the fitted densities in Figures 3 and fitted cdf’s with empirical cdf in Figures 4 for visual inspection of
the closeness of the fittings. These plots indicate that the distributions from the proposed family
provide a closer fit to this data.

32

0.5

Estimated pdf's

0.5

Estimated pdf's

0.3

0.4

MOKw-E
KwMO-E

0.0

0.1

0.2

Density

0.3
0.2
0.0

0.1

Density

0.4

MOKw-Fr
KwMO-Fr

0

1

2

3

4

5

6

0

1

2

X

3

4

5

6

X

(a)

(b)

Fig: 5 Plots of the observed histogram and estimated pdf’s for the (a) MOKw  Fr , KwMO  Fr and
(b) MOKw  E , KwMO  E distributions for example II.
Estimated cdf's
1.0

1.0

Estimated cdf's
MOKw-E
KwMO-E

0.0

0.0

0.2

0.2

0.4

0.4

cdf

cdf

0.6

0.6

0.8

0.8

MOKw-Fr
KwMO-Fr

0

1

2

3

4

5

6

0

1

2

3

X

X

(a)

(b)

4

5

6

Fig: 6 Plots of the observed ogive and estimated cdf’s for the (a) MOKw  Fr , KwMO  Fr and (b)
MOKw  E , KwMO  E distributions for example II.

7. Conclusion
In this paper, the new family of Marshall-Olkin extended Kumaraswamy generalized distributions is
introduced and some of its distributional and mathematical properties are investigated. The maximum
likelihood and moment method for estimating the parameters are also discussed. Comparative data
33

modeling by two important members of the proposed MOKw  G family with corresponding members
of KwMO  G family revealed formers superiority in all examples considered here.

Appendix: Maximum likelihood estimation for MOKw  E
The pdf of the MOKw  E distribution is given by
f

MOKwGE

(t )   a b  e  t (1  e  t ) a 1 [1  (1  e t ) a ] b 1 [1   [1  (1  e  t ) a ] b ] 2

for   0,  0, a  0, b  0, t  0
For a random sample of size n from this distribution, the log-likelihood function for the parameter
vector θ  ( , a, b,  )T is given by
n

n

 (θ)  n log   n log (a b)  n log     ti  (a  1)  log (1  e  ti )
i 0

n

i 0

n

 (b  1)  log [1  (1  e  ti ) a ]  2  log [1   [1  (1  e  ti ) a ] b ]
i0

i 0

The components of the score vector θ  ( , a , b,  )T are
n
  (θ) n
[1  (1  e ti ) a ] b
  2
 t i a b


) ]
i  0 1   [1  (1  e
n
(1  e  ti ) a log (1  e  ti )
  (θ) n n
 t i

log
(
1

e
)

(
1

b
)



a i 0
1  (1  e  t i ) a
a
i 0

 b [1  (1  e ti ) a ] b 1 (1  e  t i ) a log [(1  e  ti ) ]
1   [1  (1  e ti ) a ] b
i 0
n

 2

n
  (θ) n n
 [1  (1  e ti ) a ] b log [1  (1  e  ti ) a ]
   log [1  (1  e ti ) a ]  2 
b
b i0
1   [1  (1  e ti ) a ] b
i0

  (θ) n n
   ti  (a  1)

 i0

n
 e  ti
a  (1  e  ti ) a 1 e ti

(
1

b
)


 t i
1  (1  e ti ) a
i 0 1  e
i 0
n

 a b  [1  (1  e ti ) a ] b 1 (1  e  t i ) a 1 e  ti
.
1   [1  (1  e  t i ) a ] b
i 0
n

 2

The asymptotic variance covariance matrix for mles of the parameters of MOKw  E ( , a, b,  )
distribution is estimated by

34

cov (ˆ , bˆ) cov (ˆ , ˆ ) 
var (aˆ )
cov (aˆ bˆ) cov (aˆ , ˆ ) 

cov (bˆ, aˆ ) var (bˆ)
cov (bˆ , ˆ ) 

cov (ˆ , aˆ ) cov (ˆ, bˆ) var (ˆ ) 
   2 l (θ )  

Where the elements of the information matrix Iˆ n (θˆ )    
can be derived using the
   i   j  
  θ  θˆ
 
following second partial derivatives:
 var (ˆ )

 cov (aˆ ˆ )
Î n 1 (θˆ )  
 cov (bˆ,ˆ )

ˆ
 cov ( ,ˆ )

cov (ˆ , aˆ )

n
 2
n
[1  (1  e ti ) a ] 2 b



2

 t i a b 2
 2
2
) ] ]
i  0 [1   [1  (1  e

n 
n
(1  e ti ) 2 a log (1  e  ti ) 2 (1  e  ti ) a log (1  e t i ) 2 
 2



(
1

b
)





a2
a2
{1  (1  e  ti ) a } 2
i 0 
1  (1  e  t i ) a


 2 b 2 (1  e  ti ) 2 a [1  (1  e  ti ) a ] 2 (b 1) log (1  e ti ) 2
{1   [1  (1  e  t i ) a ] b } 2
i 0
n

 2

b (b  1)  (1  e  t i ) 2 a [1  (1  e  t i ) a ] ( b  2) log (1  e t i ) 2
1   [1  (1  e  t i ) a ] b
i 0
n

 2

b  (1  e  t i ) a [1  (1  e t i ) a ] ( b 1) log (1  e  ti ) 2
1   [1  (1  e  t i ) a ] b
i 0
n

 2

n
 2
n
 2 [1  (1  e ti ) a ] 2 b log [1  (1  e  ti ) a ] 2



2

b 2
b2
{1   [1  (1  e  t i ) a ] b } 2
i 0

 [1  (1  e  t i ) a ] b log [1  (1  e  ti ) a ] 2
1   [1  (1  e ti ) a ] b
i 0
n

 2

2
2
2
n 
n
 2
n
e  2  ti ti
e   t i ti 
a (1  e  ti ) 2 ( a 1) e 2  ti ti


  2  (a  1)   

 (1  b) 
 2

(1  e ti ) 2 1  e  t i 
{1  (1  e ti ) a } 2
i 0 
i 0
n

 (1  b) 
i0

2

n
a (a  1) (1  e  ti ) a  2 e2  ti ti
a (1  e  t i ) a 1 e   ti ti

(
1

b
)

1  (1  e  ti ) a
1  (1  e  t i ) a
i 0

 2 a 2 b 2 e 2  ti (1  e  ti ) 2 ( a 1) [1  (1  e ti ) a ] 2 ( b 1) ti 2
 2
{1   [1  (1  e  ti ) a ] b } 2
i 0
n

n

 2
i 0
n

 2
i 0

a 2  b (b  1) e2  ti (1  e  t i ) 2 ( a 1) [1  (1  e ti ) a ] b  2 ti
1   [1  (1  e  ti ) a ] b
a (a  1) b e2  ti (1  e  ti ) a  2 [1  (1  e  ti ) a ] b 1 ti
1   [1  (1  e  t i ) a ] b

35

2

2

2

a b  e  ti (1  e ti ) a 1 [1  (1  e ti ) a ] b 1 ti
 2
1   [1  (1  e  ti ) a ] b
i 0
n

2

n
 2
b  (1  e  ti ) a [1  (1  e ti ) a ] 2 b 1 log (1  e  t i )
 2
  a
[1   [1  (1  e  t i )a ] b ] 2
i 0

b (1  e  ti ) a [1  (1  e  ti ) a ] b 1 log (1  e t i )
1   [1  (1  e  t i )a ] b
i 0
n

 2

n
n
 2
 [1  (1  e ti ) a ] 2 b log [1  (1  e  t i ) a ]
[1  (1  e ti ) a ] b log [1  (1  e  t i ) a ]
 2 

2

  b
{1   [1  (1  e ti ) a ] b } 2
1   [1  (1  e  ti ) a ] b
i 0
i 0
n
 2
a b  e  t i [1  (1  e  t i ) a ] 2 b 1 (1  e ti ) a 1 ti
 2
  
[1   [1  (1  e  ti ) a ] b ] 2
i 0
n

 2
i 0

a b e  ti [1  (1  e  ti ) a ] b1 (1  e ti ) a1 ti
1   [1  (1  e ti ) a ] b

n
n
 2
(1  e  ti ) a log (1  e  ti )
 (1  e  t i ) a [1  (1  e  ti ) a ] b 1 log (1  e  t i )
 

2

 ab
1  (1  e  t i ) a
1   [1  (1  e ti ) a ] b
i 0
i0
n

 2
i 0

b  2 (1  e ti ) a [1  (1  e  ti ) a ] 2 b1 log (1  e  ti ) log [1  (1  e  ti ) a ]
{1   [1  (1  e ti ) a ] b } 2

b  (1  e ti ) a [1  (1  e ti ) a ] b 1 log (1  e ti ) log [1  (1  e  ti ) a ]
1   [1  (1  e ti ) a ] b
i 0
n

 2
 2

a

n
e   ti ti
a (1  e  ti ) 2 a 1 e   ti ti log (1  e  t i )

(
1

b
)


 t i
{1  (1  e  ti ) a } 2
i0 1  e
i 0
n

n

 (1  b) 
i0

(1  e  ti ) a 1 e  ti ti
1  (1  e  t i ) a

n

 (1  b) 
i0

a (1  e  ti ) a 1 e   ti ti log (1  e  ti )
1  (1  e  ti ) a

b  e   ti (1  e  t i ) a 1 [1  (1  e ti ) a ] b 1 ti
1   [1  (1  e  ti ) a ] b
i 0
n

 2

a b 2  2 e   t i (1  e  t i ) 2 a 1 [1  (1  e  t i ) a ] 2 (b 1) ti log (1  e  t i )
{1   [1  (1  e t i ) a ] b } 2
i 0
n

 2

a b (b  1) e   t i (1  e ti ) 2 a 1 [1  (1  e t i ) a ] b  2 ti log (1  e  t i )
 2
1   [1  (1  e  t i ) a ] b
i 0
n

a b  e   t i (1  e  t i ) a 1 [1  (1  e ti ) a ] b 1 ti log (1  e  t i )
 2
1   [1  (1  e  ti ) a ] b
i 0
n

n
n
 2
a e   t i (1  e  ti ) a 1 ti
a  e   t i (1  e t i ) a 1 [1  (1  e  t i ) a ] b 1 ti
 
 2
b
1  (1  e  t i ) a
1   [1  (1  e  t i ) a ] b
i 0
i 0

36

n

 2
i 0
n

 2
i 0

a b  2 e   ti (1  e  ti ) a 1 [1  (1  e  ti ) a ] 2 b 1 ti log [1  (1  e  ti ) a ]
{1   [1  (1  e ti ) a ] b } 2
a b  e   t i (1  e  t i ) a 1 [1  (1  e  t i ) a ] b 1 ti log [1  (1  e ti ) a ]
1   [1  (1  e ti ) a ] b

Where  (.) is the derivative of the digamma function.

References
1. Alexander C., Cordeiro G.M., Ortega E.M.M. and Sarabia J.M. (2012). Generalized betagenerated distributions. Comput Stat Data Anal 56: 1880-1897.
2. Alizadeh M., Tahir M.H., Cordeiro G.M., Zubair M. and Hamedani G.G. (2015). The
Kumaraswamy Marshal-Olkin family of distributions. Journal of the Egyptian Mathematical
Society, In press.
3. Alzaghal A., Lee C. and Famoye F. (2013). Exponentiated T-X family of distributions with
some applications. Int J Statist Probab 2: 31-49.
4. Azzalini A. (1985). A class of distributions which includes the normal ones. Scand J Statist 12:
171-178.
5. Barakat H.M., Ghitany M.E. and AL-Hussaini E.K. (2009). Asymptotic distributions of order
statistics and record values under the Marshall-Olkin parameterization operation. Commun
Statist Theor Meth 38: 2267-2273.
6. Barreto-Souza W., Lemonte A.J. and Cordeiro G.M. (2013). General results for Marshall and
Olkin's family of distributions. An Acad Bras Cienc 85: 3-21.
7. Caroni C. (2010). Testing for the Marshall-Olkin extended form of the Weibull distribution.
Stat Pap 51: 325-336.
8. Cordeiro G.M., Lemonte A.J. and Ortega E.E.M. (2014a). The Marshall-Olkin family of
distributions: Mathematical properties and new models. J Stat Theor Pract 8: 343-366.
9. Cordeiro G.M. and De Castro M. (2011). A new family of generalized distributions. J Stat
Comput Simul 81: 883-893.
10. Cordeiro G.M., Silva G.O. and Ortega E.M.M. (2012). The beta extended Weibull distribution.
J Probab Stat Sci 10: 15-40.
11. Cordeiro G.M., Ortega E.M.M. and DA Cunha D.C.C. (2013). The exponentiated generalized
class of distributions. J Data Sci 11: 1-27.
12. Cordeiro G.M., Lemonte A.J. and Ortega E.E.M. (2014b). The Marshall-Olkin family of
distributions: Mathematical properties and new models. J Stat Theor Pract 8: 343-366.

37

13. Eugene N., Lee C. and Famoye F. (2002). Beta-normal distribution and its applications.
Commun Statist Theor Meth 31: 497-512.
14. Ghitany M.E., Al-Awadhi F.A. and Alkhalfan L.A. (2007). Marshall-Olkin extended Lomax
distribution and its application to censored data. Commun Stat Theory Method 36: 1855-1866.
15. Ghitany M.E., Al-Hussaini E.K. and Al-Jarallah R.A. (2005). Marshall-Olkin extended
Weibull distribution and its application to censored data. J Appl Stat 32: 1025-1034.
16. Gieser P.W., Chang M.N., Rao P.V., Shuster J.J. and Pullen J. (1998). Modelling cure rates
using the gompertz model with covariate information. Stat Med 17(8):831–839.
17. Greenwood J.A., Landwehr J.M., Matalas N.C. and Wallis J.R. (1979). Probability weighted
moments: definition and relation to parameters of several distributions expressible in inverse
form. Water Resour Res 15: 1049-1054.
18. Gupta R.C., Gupta P.I. and Gupta R.D. (1998). Modeling failure time data by Lehmann
alternatives. Commun Statist Theor Meth 27: 887-904.
19. Gurvich M., DiBenedetto A. and Ranade S. (1997). A new statistical distribution for
characterizing the random strength of brittle materials. J. Mater. Sci. 32: 2559–2564.
20. Hussain M.A. (2013). Estimation of P(Y<X) for the class of Kumaraswamy G distribution.
Aust J Basic Appl Sci 7: 158-169.
21. Jayakumar K. and Mathew T. (2008). On a generalization to Marshall-Olkin scheme and its
application to Burr type XII distribution. Stat Pap 49: 421-439.
22. Jones M.C. (2004). Families of distributions arising from the distributions of order statistics.
Test 13: 1-43.
23. Jose K.K. (2011). Marshall-Olkin family of distributions and their applications in reliability
theory, time series modeling and stress strength analysis. Proc. ISI 58th World Statist. Congr
Int Stat Inst, 21st-26th August 2011, Dublin, p. 3918-3923.
24. Jose K.K. and Krishna E. (2011). Marshall-Olkin extended uniform distribution. ProbStat
Forum 4: 78-88.
25. Krishna E., Jose K.K. and Risti´c M. (2013). Applications of Marshal-Olkin Fréchet
distribution. Comm. Stat. Simulat. Comput. 42: 76–89.
26. Krishna E. (2011). Marshall-Olkin Generalization of Some Distributions and Their
Applications. Ph. D. Thesis, Mahatama Ghandi University, Kottayam, Kerala, India.
(Unpublished).
27. Marshall A. and Olkin I. (1997). A new method for adding a parameter to a family of
distributions with applications to the exponential and Weibull families. Biometrika 84: 641652.
38

28. Morad Alizadeh, Gauss M. Cordeiro and Clarice Garcia B. Demétrio (2015). The beta
Marshall-Olkin family of distributions. Journal of Statistical Distributions and Applications 2:4
doi:10.1186/s40488-015-0027-7.
29. Nadarajah S. (2005). Exponentiated Pareto distributions. Statistics 39: 255-260.
30. Nadarajah S., Cordeiro G.M. and Ortega E.M.M. (2012). General results for the
Kumaraswamy G distribution. J Stat Comput Simul 87: 951-979.
31. Nadarajah S., Cancho V.G. and Ortega E.M.M. (2013a). The geometric exponential Poisson
distribution. Stat Methods Appl 22: 355-380.
32. Nadarajah S., Jayakumar K. and Ristic M.M. (2013b). A new family of lifetime models. J Stat
Comput Simul 83: 1389-1404.
33. Nanda A.K. and Das S. (2012). Stochastic orders of the Marshall-Olkin extended distribution.
Statist Probab Lett 82: 295-302.
34. Nelson W. and Dognanksoy N. (1992). Computer program pownor for fitting the powernormal and-lognormal models to life or strength data from specimens of various sizes. NASA
STI/Recon Technical Report N, 92: 29813.
35. Nichols M.D. and Padgett W.J.A. (2006). A bootstrap control chart for Weibull percentiles.
Quality and Reliability Engineering International, v. 22: 141-151.
36. Pescim R.R., Cordeiro G.M., Demetrio C.G.B., Ortega E.M.M. and Nadarajah S. (2012). The
new class of Kummer beta generalized distributions. SORT 36: 153-180.
37. Ristic M.M. and Balakrishnan (2012). The gamma-exponentiated exponential distribution. J
Stat Comput Simul 82: 1191-1206.
38. Sarhan A.M. and Apaloo J. (2013). Exponentiated modified Weibull extension distribution.
Reliab Eng Syst Safety 112: 137-144.
39. Shaked M. and Shanthikumar J.G. (2007). Stochastic Orders, Springer, New York.
40. Song KS (2001). Rényi information, loglikelihood and an intrinsic distribution measure. EM J
Stat Plan Infer 93: 51-69.
41. Tahir M.H. and Nadarajah S. (2015). Parameter induction in continuous univariate
distributions: Well-established G families. Anais da Academia Brasileira de Ciências (Annals
of the Brazilian Academy of Sciences) 87(2): 539-568.
42. Torabi H. and Montazari N.H. (2014). The logistic-uniform distribution and its application.
Commun Statist Simul Comput 43: 2551-2569.
43. Zhang T. and Xie M. (2007). Failure data analysis with extended Weibull distribution.
Communication in Statistics – Simulation and Computation 36: 579–592.

39

44. Zografos K. and Balakrishnan (2009). On families of beta- and generalized gamma generated
distributions and associated inference. Stat Methodol 6: 344-362.

40

