Next Article in Journal
A Structural Credit Risk Model Driven by the Lévy Process under Knightian Uncertainty
Previous Article in Journal
Control Analysis of Stochastic Lagging Discrete Ecosystems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Flexible Univariate and Bivariate Family of Distributions for Unit Interval (0, 1)

by
Mahmoud El-Morshedy
1,2,*,
Muhammad H. Tahir
3,
Muhammad Adnan Hussain
3,
Afrah Al-Bossly
1 and
Mohamed S. Eliwa
4,5
1
Department of Mathematics, College of Science and Humanities in Al-Kharj, Prince Sattam Bin Abdulaziz University, Al-Kharj 11942, Saudi Arabia
2
Department of Statistics and Computer Science, Faculty of Science, Mansoura University, Mansoura 35516, Egypt
3
Department of Statistics, Faculty of Computing, The Islamia University of Bahawalpur, Bahawalpur 63100, Pakistan
4
Department of Statistics and Operation Research, College of Science, Qassim University, P.O. Box 6644, Buraydah 51482, Saudi Arabia
5
Department of Mathematics, Faculty of Science, Mansoura University, Mansoura 35516, Egypt
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(5), 1040; https://doi.org/10.3390/sym14051040
Submission received: 2 March 2022 / Revised: 11 May 2022 / Accepted: 13 May 2022 / Published: 19 May 2022

Abstract

:
We propose a new generator for unit interval which is used to establish univariate and bivariate families of distributions. The univariate family can serve as an alternate to the Kumaraswamy-G univariate family proposed earlier by Cordeiro and de-Castro in 2011. Further, the new generator can also be used to develop more alternate univariate and bivariate G-classes such as beta-G, McDonald-G, Topp-Leone-G, Marshall-Olkin-G and Transmuted-G for support (0, 1). Some structural properties of the univariate family are derived and the estimation of parameters is dealt. The properties of a special model of this new univariate family called a New Kumaraswamy-Weibull (NKwW) distribution are obtained and parameter estimation is considered. A Monte Carlo simulation is reported to assess NKwW model parameters. The bivariate extension of the family is proposed and the estimation of parameters is described. The simulation study is also conducted for bivariate model. Finally, the usefulness of the univariate NKwW model is illustrated empirically by means of three real-life data sets on Air Conditioned Failures, Flood and Breaking Strength of Fibers, and one real-life data on UEFA Champion’s League for bivariate model.

1. Introduction

Although, the art of proposing new model(s) and families is an old practice but in the past a few decades ago, this activity has received increased attention. This activity has been termed as “parameter addition” and “parameter induction”. The main objective of such proposals is to develop model(s) or to explore real-life phenomenon(s) through data sets available in different fields of life. Consequently, the concentration gets diverted in proposing flexible models which may be through modification, extension or generalized class (abbreviated as G-class) of distributions etc. among other methodologies.
The flexibility of a model can mainly be considered, assessed or investigated through  (i) the flexible density and hazard rate shapes of a model; which have the ability to accommodate different types of data sets and to explore lifetime and reliability phenomenons,  (ii) some nice attractive mathematical properties,  (iii) performance evaluation of estimators through simulation study, and (iv) yielding better goodness-of-fit statistics results (numerically and graphically), among other floating criterions. It has been observed in recent literature that the induction of shape parameter(s) in a model or a family, the number of useful parameters and/or functional form of the model play an important role in achieving flexible models.
In the past, the objectives of proposing flexible models were obtained through families such as Pearson, Pareto, Burr, Johnson, Tukey’s g- and -h, Feller’s Pareto, Arnold’s Pareto etc. among others. These pioneer works motivated the researchers and practitioners to move forward and as a result many new G-classes were proposed in literature; some peculiar ones are beta-G [1], Kumarasawmy-G [2], gamma-G [3], transmuted-G [4], Weibull-G [5], Lomax-G [6], Gumbel-G [7], Burr-G [8] etc. among others. It has been noted that some G-classes were developed through the logit of well-established parent models chosen from the literature with the intensions to add influence or role of proposed G-classes parameters in the special, extended, modified or generalized models. In literature, a very few G-classes have been proposed directly from the function, formula or transformation, for example Marshall-Olkin-G, transmuted-G, exponentiated-generalized-G (extended Topp-Leone model) and alpha-power transformed-G. More detail about G-classes, their developments and range of random variable T for support { ( 0 , 1 ) , ( 0 , ) , ( , ) } is available in [5,9,10,11]. Therefore, for submitting a new proposal about a G-class, there is need to select a new parent model or alternately to propose a new generator, say G ( x ) , which is a function that defines a generator in terms of the cumulative distribution function of a baseline (or parent) model denoted by G ( · ) .

2. Brief Review and Methodology

In the literature, a very few generators have been proposed for a unit interval (0, 1) with the objective to propose a new G-class from them. Also some G-classes of distributions have been reported in literature for support (0, 1) through the only available generator, G ( x ) and by using the logit of a parent model, viz. beta-G family (Eugene et al. [1]), Kumaraswamy-G (Kw-G) (Cordeiro and de-Castro [2]), McDonald-G (Mc-G) (Alexander et al. [12]) and Topp-Leone-G (TL-G) (Rezaei et al. [13]). For more G-classes, the reader is referred to [14,15,16,17].
The main objective of this article is to introduce a new generator for support (0, 1) and to proposed a new G-class from this newly proposed generator. For this we consider proposing a new Kw-G class since the existing Kw-G class [2] has received increased attention having nice mathematical properties including mainly the closed-form expression of its distribution function. To fulfill our objective, here a new Kw-G family has been proposed from a newly proposed generator 1 [ 1 G ( x ) ] exp { G ( x ) }   (as an alternative to G ( x ) ). To the best of our knowledge, there exists only two generators { G ( x ) and 1 [ 1 G ( x ) ] G ( x ) } for support (0, 1) including the recently reported one 1 [ 1 G ( x ) ] G ( x ) by Tahir et al. [18].
Alzaatreh et al. [5] pioneered a general method for proposing G-classes which they referred to as “transformed-transformer” (T-X) approach; that has received increased attention in recent years for developing generators and G-classes. Let r ( t ) be the cumulative distribution function (cdf) and R ( t ) be the probability density function (pdf) of a random variable (rv) T for interval [ a , b ] such that < a < b < . Further, let us suppose that [ G ( x ) ] is a function that defines a generator in terms of G ( x ) or survival function (sf) G ¯ ( x ) = 1 G ( x ) related to any baseline rv which satisfies the following conditions:  
(i)
[ G ( x ) ] [ a , b ] ,
(ii)
[ G ( x ) ] is differentiable and monotonically non-decreasing, and
(iii)
lim x [ G ( x ) ] = a and lim x [ G ( x ) ] = b .
The cdf for T–X family is given by
F T X ( x ) = a [ G ( x ) ] r ( t ) d t = R G ( x ) ,
where the function [ G ( x ) ] satisfy the conditions (i)–(iii). The pdf corresponding to Equation (1) is given by
f T X ( x ) = r G ( x ) d d x G ( x ) .
Kumaraswamy [19] introduced a two-parameter model for a unit interval (0, 1). The cdf and pdf for a rv T are, respectively, given by
R ( t ) = 1 ( 1 t a ) b , 0 < t < 1 ,
and
r ( t ) = a b t a 1 ( 1 t a ) b 1 , 0 < t < 1 ,
where a > 0 and b > 0 are both shape parameters.
Cordeiro and de-Castro [2] introduced the cdf and pdf of Kw-G family as
F K w G ( x ; a , b , ξ ) = 1 1 G ( x ; ξ ) a b , 0 < x < 1 ,
and
f K w G ( x ; a , b , ξ ) = a b g ( x ; ξ ) G ( x ; ξ ) a 1 1 G ( x ; ξ ) a b 1 , 0 < x < 1 ,
where a > 0 and b > 0 are two additional shape parameters of this family, and ξ is the vector for baseline parameters.
The purpose of our article is to introduce a new Kumaraswamy-G class from a new generator 1 [ 1 G ( x ) ] exp { G ( x ) } other than G ( x ) and 1 [ 1 G ( x ) ] G ( x ) for support (0, 1). This new generator can also serves for developing new G-classes for any support. For example, the support of T is (0, 1) then the new beta-G, new Mc-G and new TL-G classes can also be introduced from this new generator. The detailed properties of the new generator 1 [ 1 G ( x ) ] exp { G ( x ) } are reported in another communication.
The paper is organized as follows: In Section 3, we define the univariate New Kumaraswamy-G (NKw-G) family of distributions. In Section 4, some mathematical properties of this new univariate Kw-G class are reported. A special model of New Kumaraswamy Weibull (NKwW) distribution of the new Kw-G class is investigated in detail in Section 5. A Monte Carlo simulation is also conducted in Section 6 to assess the performance of maximum likelihood estimators of a special model of the univariate G-class. In Section 7, we define the Bivariate New Kumaraswamy-G family where some properties and simulation study are reported. In Section 8 and Section 9, the usefulness of the new univariate and bivariate families are illustrated through four real-life data sets. We conclude the paper in Section 10 with some final remarks.

3. The NKw-G Family

From Equations (1) and (4), and then inducting the new generator 1 G ¯ ( x ; ξ ) exp [ G ( x ; ξ ) ] , the cdf of the NKw-G family is defined by
F ( x ) = a b 0 1 G ¯ ( x ; ξ ) exp [ G ( x ; ξ ) ] t a 1 1 t a b 1 d t = 1 1 1 G ¯ ( x ; ξ ) exp [ G ( x ; ξ ) ] a b , 0 < x < 1 .
where the parameters a > 0 and b > 0 refine the shape and ξ is the vector of the baseline parameters.
The pdf corresponding to Equation (7) is
f ( x ) = a b g ( x ; ξ ) G ( x ; ξ ) 1 G ¯ ( x ; ξ ) exp [ G ( x ; ξ ) ] a 1 × exp [ G ( x ; ξ ) ] 1 1 G ¯ ( x ; ξ ) exp [ G ( x ; ξ ) ] a b 1 0 < x < 1 .
Henceforth, a rv having density Equation (8) is denoted by X N K w G ( a , b , ξ ) . The sf S ( x ) and hazard rate function (hrf) h ( x ) of X are, respectively, given by
S ( x ) = 1 1 G ¯ ( x ; ξ ) exp [ G ( x ; ξ ) ] a b 0 < x < 1 .
and
h ( x ) = a b g ( x ; ξ ) G ( x ; ξ ) 1 G ¯ ( x ; ξ ) exp [ G ( x ; ξ ) ] a 1 × exp [ G ( x ; ξ ) ] 1 1 G ¯ ( x ; ξ ) exp [ G ( x ; ξ ) ] a 1 0 < x < 1 .

4. Properties of the NKw-G Family

In this section, we obtain some mathematical properties of the NKw-G family.

4.1. Quantile Function

Random variates are commonly generated through quantile function (qf), Q ( u ) = F 1 ( u ) = min { x ; F ( x ) u } , where F 1 ( . ) is inverse of cdf. Thus, the qf of the NKw-G class is obtained by inverting Equation (7)
Q ( u ) = G 1 1 ( 1 u ) 1 / b 1 / a 1 e + 1 , 0 < u < 1 ,
where the Lambert-W function is the inverse function of F ( W ) = W exp { W } . The Lambert-W function was introduced by Corless et al. [20] which is used to solve a variety of exponential and logarithmic equations in the form of W exp { W } . This function is also known as the product log function, and also help to solve many sorts of problems where the variable acts as both a base and as an exponent.

4.2. Linear Representation for the NKw-G Density

In this section, we give a useful expansion for Equation (7). The purpose of expanding density by using appropriate series is to reduce the complexity of density function, and converting it into solvable form to find the statistical properties for baseline model or G-class. There are many ways to get the properties of a model or a G-class. One way could be conventional, and other way by using series and functions. The expansion of the density through linear representation by utilizing series or function is recent and have been reported in many articles (as widely accepted). The main reason for the getting linear representation using series and functions is to obtain reasonable properties of a model or a G-class such as moments, incomplete moments, moment generating function, probability weighted moments, entropy, order statistics etc. among others. Simply, the properties of any special model developed from a G-class can be obtained by just utilizing the exponentiated (or generalized) properties of the special model of that class. For an arbitrary baseline cdf G ( x ) , a rv is said to have the exponentiated-G (exp-G) distribution with power parameter λ > 0 if its cdf and pdf are H λ ( x ) = G ( x ) λ and h λ ( x ) = λ g ( x ) G ( x ) λ 1 , respectively.
If ν is any real non-integer and | z | < 1 , then by using generalized binomial expansion
( 1 z ) ν = c = 0 ( 1 ) c ν c z c ,
the Equation  (7) becomes
F ( x ) = 1 i = 0 ( 1 ) i b i 1 G ¯ ( x ; ξ ) exp [ G ( x ; ξ ) ] i a .
Now expanding Equation (11) and then applying binomial expansion again, we have the following form
F ( x ) = 1 i = 1 j = 1 a i ( 1 ) i + j b i a i j G ¯ ( x ; ξ ) j exp [ j G ( x ; ξ ) ] .
Again using generalized binomial series expansion and exponential power series, and after some simplification, Equation (12) reduced to
F ( x ) = 1 k = 0 s = 0 j i = 1 j = 1 a i ( 1 ) i + j + s j k k ! b i a i j j s G ( x ; ξ ) k + s .
The above result can be expressed in reduced form as
F ( x ) = 1 + k = 0 s = 0 j ω k , s G ( x ; ξ ) k + s ,
where
ω k , s = i = 1 j = 1 a i ( 1 ) i + j + s + 1 j k k ! b i a i j j s .
Now, differentiating Equation (13), we get the NKw-G density as
f ( x ) = k = 0 s = 0 j ω k , s ( k + s ) g ( x ; ξ ) G ( x ; ξ ) k + s 1 .
Equation (14) reveals that the family density function is a linear combination of exp-G densities. Thus, some structural properties of X, for example, moments, mean deviations, Bonferroni and Lorenz curves, and generating function can follow from exp-G distribution.

4.3. Mathematical Properties

Here we reported some mathematical properties of NKw-G family by using Equation (14) and those of the properties of the exp-G distribution.
Let Y k , s be a rv relating to exp-G model having power parameter ( k + s ) , then the rth ordinary moment of X, say E ( X r ) , can be expressed from Equation (14) as
E ( X r ) = k = 0 s = 0 j ω k , s E ( Y k , s r ) = k = 0 s = 0 j ω k , s ( k + s ) τ r , k , s ,
where τ r , k , s = x r G ( x ; ξ ) k + s 1 g ( x ; ξ ) d x = 0 1 Q G ( u ; ξ ) r u k + s d u , and Q G ( u ; ξ ) is the qf of the baseline G.
The rth lower incomplete moment of X, say m r ( y ) = y x r f ( x ) d x , is
m r ( y ) = k = 0 s = 0 j ω k , s y x r h k , s ( x ) d x = k = 0 s = 0 j ω k , s 0 G ( y ; ξ ) ) Q G ( u ; ξ ) r u k + s d u .
The last two integrals can easily be evaluated numerically for most of the baseline G models.
The moment generating function (mgf)  M ( t ) = E ( e t X ) of X follows from Equation (14) is
M ( t ) = k = 0 s = 0 j ω k , s M k , s ( t ) = k = 0 s = 0 j ω k , s ( k + s ) ρ k , s ( t ) ,
where M k , s ( t ) is the mgf of Y k , s and ρ k , s ( t ) = 0 1 exp [ t Q G ( u ; ξ ) ] u k + s d u . Hence, we can obtain the mgfs of many special NKw-G distributions directly from exp-G generating function and Equation (17).
The mode of NKw-G distribution is obtained by taking log of Equation (8) by
log f x ( x ; ξ ) = log ( a b ) + log g ( x ; ξ ) + G ( x ; ξ ) + ( a 1 ) log [ 1 ( 1 G ( x ; ξ ) ) exp [ G ( x ; ξ ) ] ] + log [ G ( x ; ξ ) ] + ( b 1 ) log 1 [ 1 ( 1 G ( x ; ξ ) ) exp [ G ( x ; ξ ) ] ] a .
Taking derivative with respect to x, the mode of NKw-G satisfies
g ( x ; ξ ) g ( x ; ξ ) + g ( x ; ξ ) a exp [ G ( x ; ξ ) ] G ( x ; ξ ) 1 b [ 1 ( 1 G ( x ; ξ ) ) exp [ G ( x ; ξ ) ] ] a 1 [ 1 ( 1 G ( x ; ξ ) ) exp [ G ( x ; ξ ) ] ] a + 1 e G ( x ; ξ ) G ( x ; ξ ) + 1 1 ( 1 G ( x ; ξ ) ) exp [ G ( x ; ξ ) ] = 0 .

4.4. Estimation of Univariate G-Family Parameters

Here, we discuss the maximum likelihood estimation of NKw-G family parameters based on random samples X 1 , , X n taken from any G distribution. The log-likelihood function ( . ) for vector θ = ( a , b , ξ ) is
( θ ) = n log ( a b ) + i = 1 n log g ( x i ; ξ ) + i = 1 n log G ( x i ; ξ ) + ( a 1 ) i = 1 n log 1 G ¯ ( x i ; ξ ) exp [ G ( x i ; ξ ) ] + i = 1 n G ( x i ; ξ ) + ( b 1 ) i = 1 n log 1 1 G ¯ ( x i ; ξ ) exp [ G ( x i ; ξ ) ] a .
The MLEs θ ^ of θ can be evaluated by maximizing ( θ ) . There exists several routines in literature which can be used for numerical maximization of ( θ ) particulary in R-language (optim function) and in Ox the (sub-routine MaxBFGS).
The score components obtained after differentiating the log-likelihood with respect to a, b and ξ are:
a = n a ( b 1 ) i = 1 n 1 G ¯ ( x i ; ξ ) exp [ G ( x i ; ξ ) ] a log 1 G ¯ ( x i ; ξ ) exp [ G ( x i ; ξ ) ] 1 1 G ¯ ( x i ; ξ ) exp [ G ( x i ; ξ ) ] a + i = 1 n log 1 G ¯ ( x i ; ξ ) exp [ G ( x i ; ξ ) ] , b = n b + i = 1 n log 1 1 G ¯ ( x i ; ξ ) exp [ G ( x i ; ξ ) ] a , ξ = i = 1 n g i ξ g ( x i ; ξ ) i = 1 n G i ξ + i = 1 n G i ξ G ( x i ; ξ ) + ( a 1 ) i = 1 n G ( x i ; ξ ) exp [ G ( x i ; ξ ) ] G i ξ 1 G ¯ ( x ; ξ ) exp [ G ( x i ; ξ ) ] ( b 1 ) i = 1 n a G ( x i ; ξ ) exp [ G ( x i ; ξ ) ] G i ξ 1 G ¯ ( x i ; ξ ) exp [ G ( x i ; ξ ) ] a 1 1 1 G ¯ ( x i ; ξ ) exp [ G ( x i ; ξ ) ] a ,
where g i ξ = ξ g ( x i ; ξ ) and G i ξ = ξ G ( x i ; ξ ) are derivatives of column vectors of the same dimension of ξ .
The MLEs of the model parameters can be obtained from the above score components after equating to zero, and then solving them simultaneously. If analytical solution is impossible, then analytical results for these parameters can be obtained through iterative algorithms.
For interval estimation of the parameters, we require the 3 × 3 observed information matrix I ( Θ ) = { I q s } (for q , s = a , b , ξ ). Standard maximization routines provide only an approximation for the observed information matrix. The observed information matrix for the parameter vector Θ = ( a , b , ξ ) is given by
I ( Θ ) = 2 ( Θ ) Θ Θ = I a a I a b I a ξ I b b I b ξ I ξ ξ .
Under standard regularity conditions, the multivariate normal N 3 ( 0 , I ( Θ ^ ) 1 ) distribution can be used to construct approximate confidence intervals for the model parameters. Here, I ( Θ ^ ) is the total observed information matrix evaluated at Θ ^ . Then, the 100 ( 1 σ ) % confidence intervals for a, b and ξ are given by a ^ ± z σ * / 2 × v a r ( a ^ ) , b ^ ± z σ * / 2 × v a r ( b ^ ) and ξ ^ ± z σ * / 2 × v a r ( ξ ^ ) , respectively, where the v a r ( · ) ’s denote the diagonal elements of I ( Θ ^ ) 1 corresponding to the model parameters, and z σ * / 2 is the quantile ( 1 σ * / 2 ) of the standard normal distribution.
The standard errors of the estimated parameters are obtained by the inverse of I ( Θ ) .

5. The Univariate NKwW Distribution

We now define the NKwW distribution by taking the Weibull baseline with cdf G ( x ) = 1 e ( α x ) β and pdf g ( x ) = α β β x β 1 exp [ ( α x ) β ] . Then, the cdf and pdf of the NKwW distribution are, respectively, given by
F N K w W ( x ) = 1 1 1 exp [ ( α x ) β ] e 1 exp [ ( α x ) β ] a b , x > 0
and
f N K w W ( x ) = a b α β β x β 1 exp [ ( α x ) β ] e 1 exp [ ( α x ) β ] 1 exp [ ( α x ) β ] e 1 exp [ ( α x ) β ] a 1 × 1 exp [ ( α x ) β ] 1 1 exp [ ( α x ) β ] e 1 exp [ ( α x ) β ] a b 1 , x > 0 ,
where α > 0 is scale parameter and β > 0 is shape parameter.
Henceforth, a random variable having density in Equation (20) is denoted by X NKwW ( a , b , α , β ) . The sf and hrf of X are given by
S ( x ) = 1 1 exp [ ( α x ) β ] e 1 exp [ ( α x ) β ] a b , x > 0
and
h ( x ) = a b α β β x β 1 exp [ ( α x ) β ] e 1 exp [ ( α x ) β ] 1 exp [ ( α x ) β ] e 1 exp [ ( α x ) β ] a 1 × 1 exp [ ( α x ) β ] 1 1 exp [ ( α x ) β ] e 1 exp [ ( α x ) β ] a 1 , x > 0 .
In Figure 1a,b, we display some plots of the density function and hrf of NKwW model for selected parameter values. From Figure 1a, it can be observed that NKwW density exhibits shapes such as unimodal (right-skewed, left skewed, symmetrical), J and reversed-J. From Figure 1b, one can depict that NKwW hazard rate shapes could be monotonic (increasing and decreasing) and non-monotonic (bathtub and upside-down bathtub).

5.1. Properties of Univariate NKwW Model

First, we will deduce linear representation of NKwW density to obtain useful properties of that model.
Following Equation (14), the NKwW density will become
f N K w W ( x ) = k = 0 s = 0 j ω k , s ( k + s ) α β β x β 1 exp { ( α x ) β } 1 exp { ( α x ) β } k + s 1 ,
is the density of exp-Weibull ( k + s , α , β ) , which after using generalized binomial expansion becomes
f N K w W ( x ) = p = 0 k + s 1 t p π ( x ; α , β ) ,
where t p = ( 1 ) p k + s 1 p k = 0 s = 0 j ω k , s ( k + s ) and π ( x ; α , β ) is the Weibull density.
From Equation (22) it is evident that NKwW density is actually a linear combination of Weibull densities. Therefore, several mathematical properties of NKwW model can easily be obtained from the properties of Weibull model.
Let Z p follows a density π ( x , α , β ) , then some properties of X can obtained from the properties of Z p .
The expression for the qf of the NKwW distribution cannot be obtained explicitly. However, we can define it in term of Lambert-W function as
Q ( u ) = 1 α 1 log 1 1 ( 1 u ) 1 / b 1 / a + W 1 ( 1 u ) 1 / b 1 / a 1 e 1 / β .
For more detail about Lambert-W function, the reader is referred to Benkthehelifa [21].
The nth moment of X, say μ r = 0 x r f ( x ) d x , is
μ r = p = 0 k + s 1 t p α β β 0 x r x β 1 exp { ( p + 1 ) ( α x ) β } d x .
After solving the last integral, the nth moment expression for X will be
μ r = Γ r β + 1 p = 0 k + s 1 t p α r ( p + 1 ) r / β .
In general, some characteristics of a model can be obtained through moments especially mean, variance, skewness and kurtosis. Table 1 displays the first four moments of NKwW distribution on different parametric values using Equation (24). The skewness and kurtosis plots of NKwW model are displayed in Figure 2. These plots reveal that the parameters a and b play a significant role in modeling the skewness and kurtosis behaviors of X.
The rth incomplete moment expression can be written as
m r ( z ) = p = 0 k + s 1 v p α r ( p + 1 ) r / β γ r β + 1 , ( p + 1 ) ( α x ) β ,
where γ ( q , x ) = 0 x q 1 exp ( x ) d x is Euler gamma function.

5.2. Estimation of Univariate NKw Family Parameters

The log-likelihood function for the parameters vector θ = ( α , β , a , b ) of NKwW model given in Equation (20)
= n log ( a b α β β ) + ( β 1 ) i = 1 n log ( x i ) + i = 1 n log 1 exp [ ( α x i ) β ] i = 1 n ( α x i ) β + i = 1 n 1 exp [ ( α x i ) β ] + ( a 1 ) i = 1 n log 1 exp [ ( α x i ) β ] e 1 exp [ ( α x i ) β ] + ( b 1 ) i = 1 n log 1 1 exp [ ( α x i ) β ] e 1 exp [ ( α x i ) β ] a .
when maximizing the log-likelihood for constrained parameters, the parameters are often transformed to avoid the constraints and failure of Newton-Raphson [22]. This is the case of the proposed model since all the parameters are constrained (>0).

6. Simulation Study

Using Monte Carlo simulations, we now evaluate the asymptotic properties of the MLEs for the parameters of the NKwW model. The simulation study is repeated 1000 times with sample sizes of n = 50, 100, 200 and 500 with the following parameter scenarios: I: a = 0.9, b =1.5, α = 0.1, and β = 2.5, II: a = 2.5, b =2.8, α = 0.4, and β = 0.8, and III: a = 0.5, b =5.5, α = 1.5, and β = 1.5. The random samples are generated using Equation (23). Table 2, Table 3 and Table 4 display the average estimate (Average), average bias (Bias) of the MLEs, mean square error (MSE), average standard error (Average SE), standard deviation (SD), model-based coverage probabilities (CP), lower bound (LB) and upper bound (UB) for the parameters a, b, α and β , respectively, for these scenarios and sample sizes. We conclude from the simulation findings that MLEs perform well in estimating parameters a, b, α and β . The CPs of the confidence intervals are quite close to the 95% nominal levels. It is apparent that the bias and MSE decreases as the sample size increases. This gives an indication that maximum likelihood method has performed well while estimation of NKwW model parameters. The related bias, MSE, SD and Average SE are, respectively, defined by
B i a s ( θ ^ ) = i = 1 N θ i ^ N θ , M S E ( θ ^ ) = i = 1 N ( θ i ^ θ ) 2 N ,
S D ( θ ^ ) = i = 1 N θ i ^ θ 2 N 1 and A v e r a g e S E ( θ ^ ) = S D N = N N 1 i = 1 N θ i ^ θ 2 .
For the 95% confidence intervals of the MLEs, we calculate the empirical coverage probability (CP) which is defined by
C P ( θ ^ ) = 1 N i = 1 N I θ i ^ 1.9599 s θ i ^ < θ < θ i ^ + 1.9599 s θ i ^ ,
where I ( . ) is the indicator function and s θ i ^ = s α i ^ , s β i ^ , s a i ^ , s b i ^ are the standard errors of the MLEs.

7. Bivariate New Kumaraswamy (BvNKw) G-Family

In this Section, we introduce a bivariate extension of the NKw-G univariate family by considering Marshall and Olkin shock model (see, Marshall and Olkin [23]). This approach has been adopted by Kundu and Gupta [24], Barreto-Souza and Lemonte [25], Ghosh and Hamedani [26], Eliwa and El-Morshedy [27], El-Morshedy et al. [28,29] etc. among others, for developing bivariate extensions of the well-known parent models. The BvNKw-G family is constructed from three independent NKw-G families using a minimization process. Assuming that the three independent rvs are Z 1 NKw G ( b 1 , a , ξ ) , Z 2 NKw G ( b 2 , a , ξ ) and Z 3 NKw G ( b 3 , a , ξ ) . Define X 1 = min { Z 1 , Z 3 } and X 2 = min { Z 2 , Z 3 } . Thus, the random vector X = ( X 1 , X 2 ) is said to have the BvNKw-G family with parameters vector Ω = ( b 1 , b 2 , b 3 , a , ξ ) if its joint reliability function (jrf) is given by
F ¯ X 1 , X 2 ( x 1 , x 2 ) = S NKw - G ( x 1 ; b 1 , a , ξ ) S NKw - G ( x 2 ; b 2 + b 3 , a , ξ ) if x 1 x 2 , S NKw - G ( x 1 ; b 1 + b 3 , a , ξ ) S NKw - G ( x 2 ; b 2 , a , ξ ) if x 1 > x 2 .
The marginal rfs of the BvNKw-G family are
F ¯ X i ( x i ) = Pr min { T i , T 3 } > x i = F ¯ NKw - G ( x i ; b i + b 3 , a , ξ ) ; i = 1 , 2 .
The joint pdf (jpdf) corresponding to Equation (27) is
f X 1 , X 2 ( x 1 , x 2 ) = f NKw - G ( x 1 ; b 1 , a , ξ ) f NKw - G ( x 2 ; b 2 + b 3 , a , ξ ) if 0 < x 1 < x 2 < , f NKw - G ( x 1 ; b 1 + b 3 , a , ξ ) f NKw - G ( x 2 ; b 2 , a , ξ ) if 0 < x 2 < x 1 < , b 3 b 1 + b 2 + b 3 f NKw - G ( x ; b 1 + b 2 + b 3 , a , ξ ) if 0 < x 1 = x 2 = x < .
The jpdf over the domains 0 < x 1 < x 2 < and 0 < x 2 < x 1 < can be obtained by differentiating Equation (27) with respect to x 1 and x 2 , whereas the jpdf over the domain 0 < x 1 = x 2 = x < can be derived from a well-known formula (see Eliwa and El-Morshedy [27]). The marginal pdfs corresponding to Equation (28) can be proposed as
f X i ( x i ) = f NKw - G ( x i ; b i + b 3 , a , ξ ) ; i = 1 , 2 .
Using Equations (27) and (29), the joint hrf (jhrf) can easily be derived by using h X 1 , X 2 ( x 1 , x 2 ) = f X 1 , X 2 ( x 1 , x 2 ) F ¯ X 1 , X 2 ( x 1 , x 2 ) . In Figure 3, Figure 4 and Figure 5, some plots of the jpdf, jhrf, and jrf are presented for selected values of the BvNKw-Weibull (BvNKwW) parameters.
If X have the BvNKw-G family, then the distributions of max { X 1 , X 2 } and min { X 1 , X 2 } are
F max { X 1 , X 2 } ( t ) = i = 1 3 F N K w G ( t ; b i , a , ξ ) and F min { X 1 , X 2 } ( t ) = 1 i = 1 3 F ¯ N K w G ( t ; b i , a , ξ ) ,
respectively. If X 1 NKw G ( b 1 + b 3 , a , ξ ) and X 2 NKw G ( b 2 + b 3 , a , ξ ) , then the qf for the marginals is
Q ( u ) X i = G 1 W e 1 1 ( 1 u ) 1 / ( b i + b 3 ) 1 / a 1 + 1 ; i = 1 , 2 .
Using (31), the coefficient of correlation between X 1 and X 2 is
Q ( u ) X 1 , X 2 = 4 F NKw - G Q ( u ) X 2 ; b 2 + b 3 , a , ξ F NKw - G Q ( u ) X 1 ; b 1 , a , ξ 1 if x 1 < x 2 4 F NKw - G Q ( u ) X 1 ; b 1 + b 3 , a , ξ F NKw - G Q ( u ) X 2 ; b 2 , a , ξ 1 if x 1 > x 2 .
The coefficient of median correlation between X 1 and X 2 can be derived at u = 0.5 . The BvNKw-G family exhibits both characteristics such as absolute continuous as well as singular part likewise reported in Marshall and Olkin’s bivariate exponential model. Thus, the joint distribution function (jdf) of X 1 and X 2 exhibits a singular part along the line x 1 = x 2 having weight b 3 / ( b 1 + b 2 + b 3 ) , and an absolute continuous part for 0 < x 1 x 2 <∞ having weight ( b 1 + b 2 ) / ( b 1 + b 2 + b 3 ) . If γ i = F ¯ X i ( x i ) where X i NKw G ( b i + b 3 , a , ξ ) ; i = 1 , 2 , then the jrf by utilizing the Marshall-Olkin copula will be as follows:
F ¯ X 1 , X 2 ( x 1 , x 2 ) = γ 1 1 Ω 1 γ 2 1 Ω 2 max γ 1 Ω 1 , γ 2 Ω 2 , for 0 < Ω 1 , Ω 2 < 1 ,
where Ω i = b 3 b i + b 3 .

7.1. The MLE for the BvNKw-G Parameters

Here we dealt with the estimation of BvNKw-G family parameters through maximum likelihood method. Suppose that ( x 11 , x 21 ) , ( x 12 , x 22 ) , , ( x 1 m , x 2 m ) is a sample of size m from the BvNKw-G family where Ψ 1 = { x 1 i < x 2 i } , Ψ 2 = { x 1 i > x 2 i } , Ψ 3 = { x 1 i = x 2 i = x i } , Ψ 1 = m 1 , Ψ 2 = m 2 , Ψ 3 = m 3 and Ψ = i = 1 3 m i = m . Using Equation (29), the likelihood function ( Ω ) can be written as
( Ω ) = i = 1 m 1 f NKw - G ( x 1 i ; b 1 , a , ξ ) f NKw - G ( x 2 i ; b 2 + b 3 , a , ξ ) i = 1 m 2 f NKw - G ( x 1 i ; b 1 + b 3 , a , ξ ) × f NKw - G ( x 2 i ; b 2 , a , ξ ) i = 1 m 3 b 3 b 1 + b 2 + b 3 f NKw - G ( x i ; b 1 + b 2 + b 3 , a , ξ ) .
Differentiating after taking logarithm of the above likelihood function with respect to parameters b 1 , b 2 , b 3 , a and ξ , and then equating them equal to zero can help in obtaining estimates of the bivariate model. If the solution is not possible, then optimization methods can be adopted to obtain numerical results.

7.2. Kendall’s Rank Correlation ( τ )

Assume ( W 1 , W 2 ) and ( Z 1 , Z 2 ) are two independent and identically distributed random vectors, then the τ is defined as follows
τ = P W 1 Z 1 W 2 Z 2 > 0 P W 1 Z 1 W 2 Z 2 < 0 .
Nelsen [30] has shown that the τ is also a copula property. The BvNKwW distribution can be obtained using the copula function also. For every bivariate distribution function, with continuous marginals F X 1 ( x 1 ) and F X 2 ( x 2 ) , corresponds a unique function C : [ 0 , 1 ] 2 [ 0 , 1 ] , called a copula such that
F X 1 , X 2 ( x 1 , x 2 ) = C ( F X 1 ( x 1 ) , F X 2 ( x 2 ) ) ,
where C ( . , . ) is the Clayton copula (see Nelsen [30,31]). Thus, if ( Z 1 , Z 2 ) has the BvNKwW family, the τ between Z 1 and Z 2 is given by
τ = 0 1 0 1 C ( u , v ) 2 C ( u , v ) u v d u d v 1 .

7.3. Simulation for the BvNKw-G Family

In this section, the MLE method is used to estimate the parameters b 1 , b 2 , b 3 , α , β and a of the BvNKwW distribution. The population parameters are generated using software R package. The sampling distributions are obtained for different sample sizes n = 50, 100, 200 from N = 1000 replications. This study proposed an assessment of the properties of the MLE for the parameters. Table 5 reports the simulation results for the BvNKwW model for b 1 = 3 , b 2 = 3 , b 3 = 3 , α = 0.5 , a = 0.3 , β = 1.5 with different sample sizes.
From Table 5, it is observed that the bias and MSE are reduced as the sample size is increased for the MLE. Further, the results of simulation indicate that the BvNKwW distribution works well under the situation where no censoring occurs.

8. Empirical Illustrations of the Proposed Univariate Model

We compare the newly proposed 4-parameter model, NKwW distribution, with some well-established models viz. 4-parameter Kumaraswamy-Weibull [32], 4-parameter beta-Weibull [33], 4-parameter exponentiated-generalized Weibull [34], 5-parameter exponentiated Kumaraswamy-Weibull [35], 3-parameter gamma-Weibull [36], 3-parameter exponentiated-Weibull [37] and 2-parameter Weibull [38] distributions. Table 6 shows the abbreviation of some well–established competitive models.
Three real-life data sets, namely, Air Conditioned Failures Data [39], Flood Data [40] and Breaking Strength Data [41] (see Appendix A for data sets) are used to compare and illustrate the potentiality of NKwW model.
The MLEs ( ^ ) and the standard errors (SEs) of the estimates, and the results of goodness-of-fit (GoF) criterion viz. Akaike information criterion (AIC), Bayesian Information Criterion (BIC), Hannan-Quinn Information Criterion (HQIC), Anderson-Darling (AD), Cramér–von Mises (CvM) and Kolmogrov-Smirnov (KS) are computed using AdequacyModel package in R-Statistical Computing Environment written by Marinho et al. [42]. Then, GoF values of all models are compared to search for the best model. The lower values of GoF statistics (GoFS) and high KS p-values indicate good fits.
In Table 7, Table 8 and Table 9, the values of MLEs and their SEs are listed for the NKwW, KwW, BW, EGW, EKwW, GaW, EW and W models, when applied to three real-life data sets. In Table 10, Table 11 and Table 12, the numerical values of the GoFS are reported. It is observed from the results of Table 10, Table 11 and Table 12 that our proposed univariate model NKwW shows minimum value of GoFS in comparison to other competitive models. This superiority is valid for the three real-life data sets. The graphs presented in Figure 6, Figure 7 and Figure 8 also support that “the model NKwW is better as compared to other well-known and well-established models”.

9. Empirical Illustration of the Proposed Bivariate Model

In this section, we empirically show that our proposed bivariate NKwW (BvNKwW) model is better in performance in comparison to some selected bivariate models. The MLEs of bivariate models are computed and are then compared on the basis of GoFS viz. AIC, BIC and HQIC. The fitting of marginal of the BvNKwW distribution is assessed through KS GoFS and its p-value. The UEFA Champion’s League data (Meintanis, [43]) is considered for the empirical illustration of the proposed bivariate model. The proposed BvNKwW model and other bivariate models, bivariate exponentiated modified Weibull extension (BExMW) (El-Gohary et al. [44]), bivariate exponentiated Weibull-Gompertz (BvExWGz) (El-Bassiouny et al. [45]), bivariate Weibull (BvW) (Kundu and Dey, [24]), bivariate odd Weibull exponential (BvOWE) (Eliwa and El-Morshedy [46]), bivariate exponentiated Weibull (BvExW) (new), bivariate generalized power Weibull (BvGPW) (new) models are fitted to UEFA Champion’s League data. Firstly, we fit the marginals X 1 , X 2 and max ( X 1 , X 2 ) to the UEFA League data, and obtain the MLEs of the parameters (b, a, α , β ) of the corresponding NKwW distribution for X 1 , X 2 and max ( X 1 , X 2 ) yielding (124.9434, 16.3415, 2.8149, 0.1685), (123.7074, 3.4388, 0.0159, 0.2619) and (7.82 × 10 4 , 39.6294, 854.4094, 0.0780), respectively. The K-S and its corresponding p-values for the marginals are reported in Table 13. It is evident from the graphs in Figure 9, Figure 10 and Figure 11, that our proposed model fits well on the basis of numerical values of Table 13.
The sample τ ^ and the model-based τ are approximately equal where τ ^ τ s a m p l e = 0.505 and τ B v N K w W = 0.513 having p-value =0.879.
It is clear from the above plots that the BvNKwW distribution fits well to the UEFA Champion’s League data. In Table 14, we report the MLEs, ^ , AIC, BIC and HQIC values for the BvNKwW and some competitive bivariate models.
It is evident from the results of Table 14 that the proposed BvNKwW model provides a better goodness-of-fits results as compared to other competitive bivariate distributions.

10. Concluding Remarks

We considered Kumaraswamy distribution and proposed new univariate and bivariate Kumaraswamy-G family from a newly proposed generator 1 G ¯ ( x ) exp [ G ( x ) ] for support (0, 1) with the objective to provide an alternative family in comparison to Kumaraswamy-G family of distributions already introduced in literature by Cordeiro and de-Castro in 2011. By using this new generator, alternative beta-G, Mc-G and TL-G families of distributions can also be developed. One can easily inspect that the proposed generator 1 G ¯ ( x ) exp { G ( x ) } has a different functional form as compared the existing generator G ( x ) for (0, 1). Some mathematical properties of the new univariate family are reported. Also, some structural properties of a special model of this class called the New Kumaraswamy-Weibull distribution were also obtained. After comparing the proposed special univariate model with some well-known extended Weibull models such as Kumaraswamy-Weibull, beta-Weibull, exponentiated-generalized Weibull, exponentiated Kumaraswamy-Weibull, gamma-Weibull and exponentiated-Weibull models, we observed that our proposed univariate model is better in performance as compared to all competitive models when applied to these three real-life data sets. Similarly, the BvNKwW distribution is proposed, and then compared with other well-known bivariate models such as bivariate exponentiated modified Weibull extension, bivariate exponentiated Weibull-Gompertz, bivariate Weibull, bivariate odd Weibull exponential, bivariate exponentiated Weibull (new) and bivariate generalized power Weibull (new) distributions. The results of some popular goodness-of-fit statistics showed that the proposed bivariate Kumaraswamy model is better in comparison to other bivariate models.

Author Contributions

The main results are a part of M.A.H.’s dissertation under M.H.T.’s supervision. Another authors A.A.-B., M.E.-M. and M.S.E. contributed to the mathematical properties, and bivariate section and its application. In general, all authors have contributed to manuscript refinement, preparation, and revision. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors are grateful to the Editor, and the three anonymous referees for their useful comments and suggestions, which led to improve our earlier version of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Data Set 1. Air Conditioned Failure Data
The first data set represents the successive failures of 179 air conditioning system given by Kus [39]. The data are as follows: 50, 130, 487, 57, 102, 15, 14, 10, 57, 320, 261, 51, 44, 9, 254, 493, 33, 18, 209, 41, 58, 60, 48, 56, 87, 11, 102, 12, 5, 14, 14, 29, 37, 186, 29, 104, 7, 4, 72, 270, 283, 7, 61, 100, 61, 502, 220, 120, 141, 22, 603, 35, 98, 54, 100, 11, 181, 65, 49, 12, 239, 14, 18, 39, 3, 12, 5, 32, 9, 438, 43, 134, 184, 20, 386, 182, 71, 80, 188, 230, 152, 5, 36, 79, 59, 33, 246, 1, 79, 3, 27, 201, 84, 27, 156, 21, 16, 88, 130, 14, 118, 44, 15, 42, 106, 46, 230, 26, 59, 153, 104, 20, 206, 5, 66, 34, 29, 26, 35, 5, 82, 31, 118, 326, 12, 54, 36, 34, 18, 25, 120, 31, 22, 18, 216, 139, 67, 310, 3, 46, 210, 57, 76, 14, 111, 97, 62, 39, 30, 7, 44, 11, 63, 23, 22, 23, 14, 18, 13, 34, 16, 18, 130, 90, 163, 208, 1, 24, 70, 16, 101, 52, 208, 95, 62, 11, 191, 14, 71.
Data Set 2. Flood Data
The second data set represents the maximum annual flood discharges (in units of 1000 cubic feet per second) of the North Saskachevan River at Edmonton, over a period of 48 years given by Asgharzadeh et al. [40]. The data are as follows: 19.885, 20.940, 21.820, 23.700, 24.888, 25.460, 25.760, 26.720, 27.500, 28.100, 28.600, 30.200, 30.380, 31.500, 32.600, 32.680, 34.400, 35.347, 35.700, 38.100, 39.020, 39.200, 40.000, 40.400, 40.400, 42.250, 44.020, 44.730, 44.900, 46.300, 50.330, 51.442, 57.220, 58.700, 58.800, 61.200, 61.740, 65.440, 65.597, 66.000, 74.100, 75.800, 84.100, 106.600, 109.700, 121.970, 121.970, 185.560.
Data Set 3. Breaking Strength Data
The third data set represents breaking strengths of 100 yarn given by Duncan [41]. The data are as follows: 66, 117, 132, 111, 107, 85, 89, 79, 91, 97, 138, 103, 111, 86, 78, 96, 93, 101, 102, 110, 95, 96, 88, 122, 115, 92, 137, 91, 84, 96, 97, 100, 105, 104, 137, 80, 104, 104, 106, 84, 92, 86, 104, 132, 94, 99, 102, 101, 104, 107, 99, 85, 95, 89, 102, 100, 98, 97, 104, 114, 111, 98, 99, 102, 91, 95, 111, 104, 97, 98, 102, 109, 88, 91, 103, 94, 105, 103, 96, 100, 101, 98, 97, 97, 101, 102, 98, 94, 100, 98, 99, 92, 102, 87, 99, 62, 92, 100, 96, 98.

References

  1. Eugene, N.; Lee, C.; Famoye, F. Beta-normal distribution and its applications. Commun. Stat. Theory Methods 2002, 31, 497–512. [Google Scholar] [CrossRef]
  2. Cordeiro, G.M.; de-Castro, M. A new family of generalized distributions. J. Stat. Comput. Simul. 2011, 81, 883–898. [Google Scholar] [CrossRef]
  3. Zografos, K.; Balakrishnan, N. On families of beta- and generalized gamma generated distributions and associated inference. Stat. Methodol. 2009, 6, 344–362. [Google Scholar] [CrossRef]
  4. Shaw, W.T.; Buckley, I.R. The Alchemy of Probability Distributions: Beyond Gram-Charlier Expansions, and a Skew-Kurtotic-Normal Distribution from a Rank Transmutation Map. UCL Discovery Repository. 2009. Available online: http://discovery.ucl.ac.uk/id/eprint/643923 (accessed on 12 May 2022).
  5. Alzaatreh, A.; Lee, C.; Famoye, F. A new method for generating families of continuous distributions. Metron 2013, 71, 63–79. [Google Scholar] [CrossRef] [Green Version]
  6. Cordeiro, G.M.; Ortega, E.M.M.; Popović, B.V.; Pescim, R.R. The Lomax generator of distributions: Properties, minification process and regression model. Appl. Math. Comput. 2014, 247, 465–486. [Google Scholar] [CrossRef]
  7. Al-Aqtash, R.; Famoye, F.; Lee, C. On generating a new family of distributions using the logit function. J. Probab. Stat. Sci. 2015, 13, 135–152. [Google Scholar]
  8. Nasir, M.A.; Tahir, M.H.; Jamal, F.; Õzel, G. A new generalized Burr family of distributions for the lifetime data. J. Statist. Appl. Probab. 2017, 6, 401–417. [Google Scholar] [CrossRef]
  9. Lee, C.; Famoye, F.; Alzaatreh, A. Methods for generating families of univariate continuous distributions in the recent decades. WIREs Comput Stat. 2013, 5, 219–238. [Google Scholar] [CrossRef]
  10. Tahir, M.H.; Nadarajah, S. Parameter induction in continuous univariate distributions: Well-established G families. An. Acad. Bras. Ciênc. 2015, 87, 539–568. [Google Scholar] [CrossRef]
  11. Tahir, M.H.; Cordeiro, G.M. Compounding of distributions: A survey and new generalized classes. J. Stat. Dist. Applic. 2016, 3, 35. [Google Scholar] [CrossRef] [Green Version]
  12. Alexander, C.; Cordeiro, G.M.; Ortega, E.M.M.; Sarabia, J.M. Generalized beta-generated distributions. Comput. Stat. Data Anal. 2012, 56, 1880–1897. [Google Scholar] [CrossRef]
  13. Rezaei, S.; Sadr, B.B.; Alizadeh, M.; Nadarajah, S. Topp-Leone generated family of distributions: Properties and applications. Commun. Stat. Theory Methods 2017, 46, 2893–2909. [Google Scholar] [CrossRef]
  14. Cordeiro, G.M.; Carrasco, J.M.F.; Ortega, E.M.M. A generalized modified Weibull distribution for lifetime modeling. Comput. Stat. Data Anal. 2008, 53, 450–462. [Google Scholar]
  15. Rocha, R.; Nadarajah, S.; Tomazella, V.; Louzada, F.; Eudes, A. New defective models based on the Kumaraswamy family of distributions with application to cancer data sets. Stat. Methods Med. Res. 2017, 26, 1737–1755. [Google Scholar] [CrossRef] [PubMed]
  16. Handique, L.; Chakraborty, S. A new four-parameter extension of Burr-XII distribution: Its properties and applications. Jpn. J. Stat. Data Sci. 2018, 1, 271–296. [Google Scholar] [CrossRef]
  17. Klakattawi, H.S. The Weibull-gamma distribution: Properties and applications. Entropy 2019, 21, 438. [Google Scholar] [CrossRef] [Green Version]
  18. Tahir, M.H.; Hussain, M.A.; Cordeiro, G.M.; El-Morshedy, M.; Eliwa, M.S. A new Kumaraswamy generalized family of distributions with properties, applications, and bivariate extension. Mathematics 2020, 8, 1989. [Google Scholar] [CrossRef]
  19. Kumaraswamy, P. Generalized probability density function for double bounded random processes. J. Hydrol. 1980, 46, 79–88. [Google Scholar] [CrossRef]
  20. Corless, R.M.; Gonnet, G.H.; Hare, D.E.G.; Jeffrey, D.J.; Knuth, D.E. On the Lambert W function. Adv. Comput. Math. 1996, 5, 329–359. [Google Scholar] [CrossRef]
  21. Benkhelifa, L. The Marshall-Olkin extended generalized Lindley distribution: Properties and applications. Commun. Stat. Simul. Comput. 2017, 46, 8306–8330. [Google Scholar] [CrossRef]
  22. MacDonald, I.L. Does Newton-Raphson really fail? Stat. Methods Med. Res. 2014, 23, 308–311. [Google Scholar] [CrossRef]
  23. Marshall, A.W.; Olkin, I. A multivariate exponential distribution. J. Am. Statist. Assoc. 1967, 62, 30–44. [Google Scholar] [CrossRef]
  24. Kundu, D.; Gupta, R.D. Bivariate generalized exponential distribution. J. Multivar. Anal. 2009, 100, 581–593. [Google Scholar] [CrossRef] [Green Version]
  25. Barreto-Souza, W.; Lemonte, A.J. Bivariate Kumaraswamy distribution: Properties and a new method to generate bivariate classes. Statistics 2013, 47, 1321–1342. [Google Scholar] [CrossRef]
  26. Ghosh, I.; Hamedani, G.G. On the Ristic-Balakrishnan distribution: Bivariate extension and characterizations. J. Statist. Theory Prac. 2018, 12, 436–449. [Google Scholar] [CrossRef] [Green Version]
  27. Eliwa, M.S.; El-Morshedy, M. Bivariate Gumbel-G family of distributions: Statistical properties, Bayesian and non-Bayesian estimation with application. Ann. Data Sci. 2019, 6, 39–60. [Google Scholar] [CrossRef]
  28. El-Morshedy, M.; Alhussain, Z.A.; Atta, D.; Almetwally, E.M.; Eliwa, M.S. Bivariate Burr X generator of distributions: Properties and estimation methods with applications to complete and type-II censored samples. Mathematics 2020, 8, 264. [Google Scholar] [CrossRef] [Green Version]
  29. El-Morshedy, M.; Eliwa, M.S.; El-Gohary, A.; Khalil, A.A. Bivariate exponentiated discrete Weibull distribution: Statistical properties, estimation, simulation and applications. Math. Sci. 2020, 14, 29–42. [Google Scholar] [CrossRef] [Green Version]
  30. Nelsen, R.B. An Introduction to Copulas; Springer: New York, NY, USA, 1999. [Google Scholar]
  31. Nelsen, R.B. An Introduction to Copulas, 2nd ed.; Springer: New York, NY, USA, 2006. [Google Scholar]
  32. Cordeiro, G.M.; Ortega, E.M.M.; Nadarajah, S. The Kumaraswamy Weibull distribution with application to failure data. J. Franklin Inst. 2010, 347, 1399–1429. [Google Scholar] [CrossRef]
  33. Lee, C.; Famoye, F.; Olumolade, O. Beta-Weibull distribution: Some properties and applications to censored data. J. Mod. Appl. Stat. Methods. 2007, 6, 173–186. [Google Scholar] [CrossRef]
  34. Oguntunde, P.E.; Odetunmibi, O.A.; Adejum, A.O. On the exponentiated generalized Weibull distribution: A generalization of the Weibull distribution. Indian J. Sci. Tech. 2015, 8. [Google Scholar] [CrossRef]
  35. Eissa, F.H. The exponentiated Kumaraswamy-Weibull distribution with application to real data. Int J. Stat. Prob. 2017, 6, 167–182. [Google Scholar] [CrossRef]
  36. Cordeiro, G.M.; Aristizábal, W.D.; Suárez, D.M.; Lozano, S. The gamma modified Weibull distribution. Chilean J. Statist. 2011, 6, 37–48. [Google Scholar]
  37. Pal, M.; Ali, M.M.; Woo, J. Exponentiated Weibull distribution. Statistica 2006, 2, 139–147. [Google Scholar]
  38. Weibull, W. A statistical distribution function of wide applicability. J. Appl. Mech. 1951, 18, 293–297. [Google Scholar] [CrossRef]
  39. Kus, C. A new lifetime distribution. Comput. Stat. Data Anal. 2007, 51, 4497–4509. [Google Scholar] [CrossRef]
  40. Asgharzadeh, A.; Bakouch, H.S.; Habibi, M. A generalized binomial exponential 2 distribution: Modeling and applications to hydrologic events. J. Appl. Statist. 2017, 44, 2368–2387. [Google Scholar] [CrossRef]
  41. Duncan, A.J. Quality Control and Industrial Statistics; Irwin: Homewood, IL, USA, 1974. [Google Scholar]
  42. Marinho, P.R.D.; Bourguignon, M.; Dias, C.R.B. Adequacy Model 1.0.8: Adequacy of Probabilistic Models and Generation of Pseudo-Random Numbers. 2013. Available online: http://cran.rproject.org/web/packages/AdequacyModel/AdequacyModel.pdf (accessed on 12 December 2013).
  43. Meintanis, S.G. Test of fit for Marshall-Olkin distributions with applications. J. Stat. Plann. Infer. 2007, 137, 3954–3963. [Google Scholar] [CrossRef]
  44. El-Gohary, A.; El-Bassiouny, A.H.; El-Morshedy, M. Bivariate exponentiated modified Weibull extension distribution. J. Stat. Appl. Probab. 2016, 5, 67–78. [Google Scholar] [CrossRef]
  45. El-Bassiouny, A.H.; EL-Damcese, M.; Abdelfattah, M.; Eliwa, M.S. Bivariate exponentaited generalized Weibull-Gompertz distribution. J. Appl. Probab. Stat. 2016, 11, 25–46. [Google Scholar]
  46. Eliwa, M.S.; El-Morshedy, M. Bivariate odd Weibull-G family of distributions: Properties, Bayesian and non-Bayesian estimation with bootstrap confidence intervals and application. J. Taibah Uni. Sci. 2020, 14, 331–345. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Plots of the (a) densities, and (b) hazard rate of the NKwW model.
Figure 1. Plots of the (a) densities, and (b) hazard rate of the NKwW model.
Symmetry 14 01040 g001
Figure 2. Plots of the skewness and kurtosis of NKwW model for different parametric values.
Figure 2. Plots of the skewness and kurtosis of NKwW model for different parametric values.
Symmetry 14 01040 g002
Figure 3. The surface plots of the jpdf, jhrf and the jrf of the BvNKwW model for b 1 = 3 , b 2 = 3 , b 3 = 3 , a = 0.5 , α = 0.3 and β = 1.5 .
Figure 3. The surface plots of the jpdf, jhrf and the jrf of the BvNKwW model for b 1 = 3 , b 2 = 3 , b 3 = 3 , a = 0.5 , α = 0.3 and β = 1.5 .
Symmetry 14 01040 g003
Figure 4. The surface plots of the jpdf, jhrf and the jrf of the BvNKwW model for b 1 = 0.3 , b 2 = 0.3 , b 3 = 0.3 , a = 1.5 , α = 0.3 and β = 0.5 .
Figure 4. The surface plots of the jpdf, jhrf and the jrf of the BvNKwW model for b 1 = 0.3 , b 2 = 0.3 , b 3 = 0.3 , a = 1.5 , α = 0.3 and β = 0.5 .
Symmetry 14 01040 g004
Figure 5. The surface plots of the jpdf, jhrf and the jrf of the BvNKwW model for b 1 = 7 , b 2 = 7 , b 3 = 7 , a = 2 , α = 2.3 and β = 0.5 .
Figure 5. The surface plots of the jpdf, jhrf and the jrf of the BvNKwW model for b 1 = 7 , b 2 = 7 , b 3 = 7 , a = 2 , α = 2.3 and β = 0.5 .
Symmetry 14 01040 g005
Figure 6. The plots of estimated densities for the (a) Air Conditioned Failures Data (b) Flood Data (c) Breaking Strength Data sets.
Figure 6. The plots of estimated densities for the (a) Air Conditioned Failures Data (b) Flood Data (c) Breaking Strength Data sets.
Symmetry 14 01040 g006
Figure 7. The plots of the estimated sf for the (a) Air Conditioned Failures Data (b) Flood Data (c) Breaking Strength Data sets.
Figure 7. The plots of the estimated sf for the (a) Air Conditioned Failures Data (b) Flood Data (c) Breaking Strength Data sets.
Symmetry 14 01040 g007
Figure 8. The plots of estimated hazard rate for the (a) Air Conditioned Failures Data (b) Flood Data (c) Breaking Strength Data sets.
Figure 8. The plots of estimated hazard rate for the (a) Air Conditioned Failures Data (b) Flood Data (c) Breaking Strength Data sets.
Symmetry 14 01040 g008
Figure 9. Fitted density plots for X 1 , X 2 and max ( X 1 , X 2 ) .
Figure 9. Fitted density plots for X 1 , X 2 and max ( X 1 , X 2 ) .
Symmetry 14 01040 g009
Figure 10. Estimated CDFs plots for X 1 , X 2 and max ( X 1 , X 2 ) .
Figure 10. Estimated CDFs plots for X 1 , X 2 and max ( X 1 , X 2 ) .
Symmetry 14 01040 g010
Figure 11. pp-plots for X 1 , X 2 and max ( X 1 , X 2 ) .
Figure 11. pp-plots for X 1 , X 2 and max ( X 1 , X 2 ) .
Symmetry 14 01040 g011
Table 1. First four moments of NKwW distribution for different parametric values.
Table 1. First four moments of NKwW distribution for different parametric values.
ab α β μ 1 μ 2 μ 3 μ 4
0.52.50.50.50.93874.345145.8866863.8300
2.52.50.51.52.83038.562427.514193.4544
0.91.20.51.82.30236.133918.337960.3563
5.03.01.03.01.29971.70462.25553.0102
Table 2. Average, Bias, MSE, SD, Average SE, CP, LB and UB for scenario-I.
Table 2. Average, Bias, MSE, SD, Average SE, CP, LB and UB for scenario-I.
n = 50 n = 100
ab α β ab α β
Average0.9141.5010.12.5150.9131.5010.12.509
Bias0.0140.0010.00.0110.0130.0010.00.009
MSE0.0050.0000.00.0040.0050.0000.00.003
SD0.1030.0310.0060.0670.0670.0140.0030.046
Average SE0.0130.0040.0010.0090.0060.0010.0010.004
CP0.9280.8861.0000.9970.9210.8791.0000.995
LB0.34513.0200.2450.9131.01610.3410.3950.982
UB1.35015.5440.5762.9052.30112.2860.7923.353
n = 200 n = 500
ab α β ab α β
Average0.9051.5020.12.5070.9041.5010.12.503
Bias0.0050.0020.00.0070.0040.0010.00.003
MSE0.0030.0000.00.0030.0010.0000.00.000
SD0.0620.0230.0030.0610.0280.0080.0020.024
Average SE0.0040.0010.00010.0040.0010.00020.0000.001
CP0.930.9101.0000.990.9470.9881.0000.982
LB0.21111.1300.1660.9630.2447.9100.1081.056
UB0.98314.1640.4522.5850.79511.1500.3662.205
Table 3. Average, Bias, MSE, SD, Average SE, CP, LB and UB for scenario-II.
Table 3. Average, Bias, MSE, SD, Average SE, CP, LB and UB for scenario-II.
n = 50 n = 100
ab α β ab α β
Average2.5192.7750.4050.8312.5202.7750.4040.821
Bias0.019−0.0250.0050.0310.020−0.0250.0040.021
MSE0.1680.0540.0070.0120.1480.0510.0060.007
SD0.4240.2400.0880.1090.3440.1990.0700.081
Average SE0.0580.0330.0120.0150.0390.0220.0070.008
CP0.8830.8680.9980.9610.8730.8551.0000.949
LB13.68215.3543.0521.3796.78012.9971.2131.269
UB19.07417.7754.2163.86811.04615.6782.0963.521
n = 200 n = 500
ab α β ab α β
Average2.4942.8010.3990.8092.4902.8000.3990.806
Bias−0.0060.002−0.0010.009−0.0060.001−0.0010.006
MSE0.0420.0120.0020.0030.0400.0110.0020.002
SD0.2180.1170.0440.0530.2170.1120.0420.041
Average SE0.0160.0090.0030.0040.0100.0050.0020.002
CP0.8660.8700.9880.9400.8480.8990.9980.879
LB5.0429.6921.0560.8793.4508.7290.5470.740
UB9.14512.7391.8392.9217.24012.4221.1222.520
Table 4. Average, Bias, MSE, SD, Average SE, CP, LB and UB for scenario-III.
Table 4. Average, Bias, MSE, SD, Average SE, CP, LB and UB for scenario-III.
n = 50 n = 100
ab α β ab α β
Average0.5905.5251.7401.5660.5395.5121.6091.550
Bias0.0900.0250.2400.0660.0390.0120.1090.050
MSE0.1020.0170.4460.2600.0410.0070.1570.143
SD0.2820.1110.5750.4740.1920.0770.3620.364
Average SE0.0390.0160.0770.0700.0210.0080.0390.038
CP0.9980.8491.0000.9980.9970.8511.0001.000
LB3.63156.39358.8031.3291.04442.40411.5211.173
UB5.31960.69673.5843.4632.25946.95822.1513.020
n = 200 n = 500
ab α β ab α β
Average0.5235.5071.5671.5520.5145.5031.5311.504
Bias0.0230.0070.0670.0520.0140.0030.0310.004
MSE0.0260.0050.0820.1180.0080.0010.0230.042
SD0.1480.0710.2560.3370.0860.0320.1420.199
Average SE0.0100.00430.0170.0230.0040.0010.0060.009
CP0.9970.8910.9990.9990.9930.9031.0000.998
LB0.55635.4307.1041.0330.34025.2854.6951.009
UB1.52840.66816.1442.6741.02131.07712.40512.271
Table 5. Simulation results for the BvNKwW model.
Table 5. Simulation results for the BvNKwW model.
n = 50
b 1 b 2 b 3 α β a
Average3.2543.3253.1450.4851.5360.334
Bias0.014 0.0340.029−0.0150.0850.074
MSE0.024 0.0370.0190.0180.0350.085
    SD0.0210.0310.0140.0120.0270.079
    CP0.8570.9670.8980.9230.9910.935
    LB2.9872.8542.3340.2591.1830.227
    UP4.1174.5984.9690.743 1.8870.633
n = 100
b 1 b 2 b 3 α β a
Average3.1653.2143.0980.4891.5230.321
Bias0.008 0.0230.015−0.0110.0490.036
MSE0.008 0.0170.0120.0120.0240.056
    SD0.0040.0120.0070.0030.0170.044
    CP0.9980.8760.9640.8760.9100.962
    LB2.5972.3672.1220.3771.1010.245
    UP3.9113.6373.8760.8642.0180.886
n = 200
b 1 b 2 b 3 α β a
Average3.0723.1103.0250.4961.5130.313
Bias0.003 0.0070.008−0.0070.0130.022
MSE0.006 0.0090.0080.0070.0140.037
    SD0.0010.0040.0020.0010.0070.024
    CP0.9310.8940.9100.8830.9050.899
    LB2.1572.3672.1100.1591.0170.119
    UP4.3364.9804.7300.9692.3400.887
Table 6. Some well-established competitive models.
Table 6. Some well-established competitive models.
Some Well—Established ModelsAbbrivations
Kumaraswamy-WeibullKwW
Beta-WeibullBW
Exponentiated-generalized WeibullEGW
Exponentiated Kumaraswamy-WeibullEKwW  
Gamma-WeibullGaW
Exponentiated-WeibullEW
WeibullW
Table 7. The MLEs and their SEs (in parentheses) of all competing models for Air Conditioned Failures Data set.
Table 7. The MLEs and their SEs (in parentheses) of all competing models for Air Conditioned Failures Data set.
Modelab α β θ
NKwW1.42340.14760.15700.7035-
(0.3004)(0.0123)(0.0093)(0.0059)-
KwW6.98780.13710.44370.6404
(0.0674)(0.0104)(0.0033)(0.0026)-
BW3.86960.14360.36620.6566-
(0.7402)(0.0124)(0.0047)(0.0063)-
EGW1.36591.60630.01280.7089-
(0.8344)(0.4839)(0.0103)(0.1106)-
EKwW3.57430.15250.17560.75160.8511
(0.1978)(0.0192)(0.0165)(0.0110)(0.0888)
GaW1.2854-0.01740.7854-
(0.3523)-(0.0082)(0.1273)-
EW1.5398-0.01880.7278-
(0.4638)-(0.0068)(0.1143)-
W--0.01180.9057-
--(0.0010)(0.0512)-
Table 8. The MLEs and their standard errors (in parentheses) for Flood Data set.
Table 8. The MLEs and their standard errors (in parentheses) for Flood Data set.
Modelab α β θ
NKwW47.48530.22450.24130.8860-
(1.2203)(0.0747)(0.0365)(0.0898)-
KwW54.78250.20410.16091.0252-
(0.1358)(0.0382)(0.0153)(0.0276)-
BW23.06020.19400.13201.1080-
(8.7941)(0.0324)(0.0073)(0.0068)-
EGW5.596610.54930.00900.7774-
(2.0458)(5.6821)(0.0041)(0.1370)-
EKwW13.51030.26251.06620.668210.4554
(1.4869)(0.0374)(0.0114)(0.0093)(3.4956)
GaW14.7225-4.61440.4983-
(1.7239)-(0.1518)(0.0217)-
EW113.6840-0.76980.4651-
(142.9147)-(1.0663)(0.1165)-
W--0.01711.7719-
--(0.0015)(0.1776)-
Table 9. The MLEs and their SEs (in parentheses) of all competing models for Breaking Strength Data set.
Table 9. The MLEs and their SEs (in parentheses) of all competing models for Breaking Strength Data set.
Modelab α β θ
NKwW10.30031.99550.01792.0187-
(5.5072)(0.9568)(0.0029)(0.3505)-
KwW12.38501.43120.01432.7814-
(6.5412)(0.8195)(0.0018)(0.6625)-
BW6.85610.60980.01334.1015-
(4.8522)(0.3226)(0.0016)(1.0461)-
EGW0.40748.64590.01713.5077-
(0.2631)(5.8742)(0.0037)(0.7875)-
EKwW1.15593.16800.00963.03906.7019
(0.2743)(1.4522)(0.0013)(0.4169)(2.1035)
GaW14.4806-0.03532.1128-
(2.5218)-(0.0050)(0.1240)-
EW42.8741-0.01992.1092-
(21.3162)-(0.0026)(0.2529)-
W--0.00957.5137-
--(0.0001)( 0.5450)-
Table 10. The values of GoFS AIC, BIC, HQIC, AD, CvM, KS and KS p-value of all competing models for Air Conditioned Failures Data set.
Table 10. The values of GoFS AIC, BIC, HQIC, AD, CvM, KS and KS p-value of all competing models for Air Conditioned Failures Data set.
KS
Model ^ AICBICHQICADCvMKSp-Value
NKwW976.70811961.41601974.16601966.58600.19540.02380.03770.9615
KwW978.26751964.53501977.28501969.70500.29250.03130.03880.9501
BW977.08031962.16101974.91001967.33000.20010.02190.03910.9470
EGW978.73621965.47201978.22201970.64200.56200.08780.04200.9100
EKwW977.69411963.38801979.32501969.85100.19180.02290.03880.9505
GaW979.84451965.68901975.25101969.56600.76330.12190.05040.7530
EW978.88591963.77201974.33401967.64900.59570.09350.04400.8797
W981.14771966.29501975.67001968.88000.97550.15770.05670.6123
Table 11. The values of GoFS AIC, BIC, HQIC, AD, CvM, KS and KS p-value of all competing models for Flood Data set.
Table 11. The values of GoFS AIC, BIC, HQIC, AD, CvM, KS and KS p-value of all competing models for Flood Data set.
KS
Model ^ AICBICHQICADCvMKSp-Value
NKwW215.0072438.0144445.4992440.84290.17250.02350.06910.9761
KwW215.5195439.0389446.5238441.86750.24950.03470.08340.8924
BW216.1573440.3147447.7995443.14320.33870.04770.09730.7538
EGW218.1801444.3601451.8449447.18870.61470.09130.09730.7543
EKwW216.8837443.7674453.1234447.30310.41830.06090.08930.8387
GaW219.4700444.9401450.5537447.06150.82780.12500.11760.5203
EW216.1707438.7413445.9549440.46270.30060.04300.07630.9428
W225.7065455.4131459.1555456.82731.72860.27650.13990.3048
Table 12. The values of GoFS AIC, BIC, HQIC, AD, CvM, KS and KS p-value of all competing models for Breaking Strength Data set.
Table 12. The values of GoFS AIC, BIC, HQIC, AD, CvM, KS and KS p-value of all competing models for Breaking Strength Data set.
KS
Model ^ AICBICHQICADCvMKSp-Value
NKwW391.4242790.8485801.2692795.06592.17650.36840.13490.0525
KwW391.6541791.3081801.7288795.52562.24780.37810.13720.0463
BW392.1492792.2985802.7192796.51592.36900.39830.14460.0305
EGW391.8419791.6838802.1045795.90122.31050.38780.14170.0361
EKwW393.8688797.7376810.7635803.00942.32060.38920.15440.0170
GaW393.0037792.0073799.8228795.17042.57270.42760.15110.0208
EW393.8216793.6431801.6586796.80622.40390.42110.13630.0505
W404.7118813.4235818.6339815.53234.74880.79850.19780.0008
Table 13. The ^ , KS and KS p-values for X 1 , X 2 and max ( X 1 , X 2 ) .
Table 13. The ^ , KS and KS p-values for X 1 , X 2 and max ( X 1 , X 2 ) .
X 1 X 2 max ( X 1 , X 2 )
Model KSKS p-Value KSKS p-Value KSKS p-Value
NKwW164.70870.10370.8212163.67990.11650.6973166.18810.12880.5716
Table 14. The MLEs with its SE and GoFS for UEFA Champion’s League data.
Table 14. The MLEs with its SE and GoFS for UEFA Champion’s League data.
Model
StatisticBvNKwWBvWBvExWBvGPWBvExMWBvOWEBvExWGz
b 1 ^ 18.9630.3971.2273.2290.1670.1350.5474
SE0.0250.0630.7724.2520.281 1.4 × 10 3 0.9043
b 2 ^ 17.6460.2740.3821.9830.0610.3020.1920
SE0.6350.0660.3562.5800.1010.00010.3137
b 3 ^ 37.3920.3390.6614.0840.1390.2650.4437
SE0.9850.0670.4545.3400.227 9.0 × 10 4 0.7173
α ^ 4.6890.0830.0120.03785.9180.0250.4109
SE0.0980.0250.0330.04833.829 9.8 × 10 3 1.9960
β ^ 0.1661.2684.5050.0797
SE 0.007 0.6096.9241.2465
a ^ 14.4730.0251.0940.0048
SE0.7430.054 2.9 × 10 5 0.0325
ϑ ^ 1.3582
SE1.3611
287.666346.00298.930344.76294.135291.129294.610
AIC587.332700.00607.860697.53600.280592.259603.220
BIC596.997706.44615.914703.97609.945600.31307.082
HQIC590.739702.27610.699699.79603.687595.099607.195
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

El-Morshedy, M.; Tahir, M.H.; Hussain, M.A.; Al-Bossly, A.; Eliwa, M.S. A New Flexible Univariate and Bivariate Family of Distributions for Unit Interval (0, 1). Symmetry 2022, 14, 1040. https://doi.org/10.3390/sym14051040

AMA Style

El-Morshedy M, Tahir MH, Hussain MA, Al-Bossly A, Eliwa MS. A New Flexible Univariate and Bivariate Family of Distributions for Unit Interval (0, 1). Symmetry. 2022; 14(5):1040. https://doi.org/10.3390/sym14051040

Chicago/Turabian Style

El-Morshedy, Mahmoud, Muhammad H. Tahir, Muhammad Adnan Hussain, Afrah Al-Bossly, and Mohamed S. Eliwa. 2022. "A New Flexible Univariate and Bivariate Family of Distributions for Unit Interval (0, 1)" Symmetry 14, no. 5: 1040. https://doi.org/10.3390/sym14051040

APA Style

El-Morshedy, M., Tahir, M. H., Hussain, M. A., Al-Bossly, A., & Eliwa, M. S. (2022). A New Flexible Univariate and Bivariate Family of Distributions for Unit Interval (0, 1). Symmetry, 14(5), 1040. https://doi.org/10.3390/sym14051040

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop