Next Article in Journal
On Sub Convexlike Optimization Problems
Next Article in Special Issue
Unit Exponential Probability Distribution: Characterization and Applications in Environmental and Engineering Data Modeling
Previous Article in Journal
Prescribed Fixed-Time Adaptive Neural Control for Manipulators with Uncertain Dynamics and Actuator Failures
Previous Article in Special Issue
Unit Distributions: A General Framework, Some Special Cases, and the Regression Unit-Dagum Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Univariate Probability-G Classes for Scattered Samples under Different Forms of Hazard: Continuous and Discrete Version with Their Inferences Tests

by
Mohamed S. Eliwa
1,2,
Muhammad H. Tahir
3,
Muhammad A. Hussain
3,
Bader Almohaimeed
4,*,
Afrah Al-Bossly
5 and
Mahmoud El-Morshedy
5,6
1
Department of Statistics and Operation Research, College of Science, Qassim University, Buraydah 51482, Saudi Arabia
2
Department of Mathematics, 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 Mathematics, College of Science, Qassim University, Buraydah 51482, Saudi Arabia
5
Department of Mathematics, College of Science and Humanities in Al-Kharj, Prince Sattam bin Abdulaziz University, Al-Kharj 11942, Saudi Arabia
6
Department of Statistics and Computer Science, Faculty of Science, Mansoura University, Mansoura 35516, Egypt
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(13), 2929; https://doi.org/10.3390/math11132929
Submission received: 5 June 2023 / Revised: 25 June 2023 / Accepted: 27 June 2023 / Published: 29 June 2023
(This article belongs to the Special Issue New Advances in Distribution Theory and Its Applications)

Abstract

:
In this paper, we define a new generator to propose continuous as well as discrete families (or classes) of distributions. This generator is used for the DAL model (acronym of the last names of the authors, Dimitrakopoulou, Adamidis, and Loukas). This newly proposed family may be called the new odd DAL (NODAL) G-class or alternate odd DAL G-class of distributions. We developed both a continuous as well as discrete version of this new odd DAL G-class. Some mathematical and statistical properties of these new G-classes are listed. The estimation of the parameters is discussed. Some structural properties of two special models of these classes are described. The introduced generators can be effectively applied to discuss and analyze the different forms of failure rates including decreasing, increasing, bathtub, and J-shaped, among others. Moreover, the two generators can be used to discuss asymmetric and symmetric data under different forms of kurtosis. A Monte Carlo simulation study is reported to assess the performance of the maximum likelihood estimators of these new models. Some real-life data sets (air conditioning, flood discharges, kidney cysts) are analyzed to show that these newly proposed models perform better as compared to well-established competitive models.

1. Introduction

There is an increasing trend in modern distribution theory by which new flexible models are being tested in different fields through modifications, extended versions, and, most preferably, through generalized classes (G-classes). Suppose that W [ G ( · ) ] is a mathematical function that helps in developing G-classes. In modern distribution theory, this function is described as a generator (i.e., a function which generates a G-class after fulfilling the desired criterion). Let T be a random variable (rv); then, the generator is basically a function of a baseline (or parent) cumulative distribution function (cdf) or survival function (sf) G ¯ ( t ) = 1 G ( t ) . F ( t ) and f ( t ) are the cdf and probability density function (pdf) of a new model or a G-class, and very few generators have been reported in the literature for developing new G-classes for an rv T [ ( 0 , 1 ) , ( 0 , ) , ( , ) ] . In the literature the following probability classes/generators have been listed so far for any rv T:
(i)
G ( t ) for range T ( 0 , 1 ) ;
(ii)
log G ¯ ( t ) , log G ( t ) , G ( t ) / G ¯ ( t ) (odds) and [ log G ¯ ( t ) ] / G ¯ ( t ) for range T ( 0 , ) ;
(iii)
log [ log G ¯ ( t ) ] and log { G ( t ) / G ¯ ( t ) } (log-odds) for range T ( , ) ] .
The main objective of this article is to present a new G-class of distributions through some odds ratio (or function).
For a lifetime rv T, let G ( t ) / G ¯ ( t ) be the basic odds function criterion, which naturally turns into ratio λ ( t ) / h ( t ) (reversed hazard rate function divided by the hazard rate function), as a useful measure for lifetime assessment of component(s) (or human organ(s)). Moreover, the probabilities of a uni-variate continuous rv spread over the range of the cdf and sf, that is, G ( t ) + G ¯ ( t ) = g ( t ) d t = 1 . So, the ratio between these two alternatives ( G ( t ) and G ¯ ( t ) ) is very useful in investigating changes occurring within a model or a phenomenon. Furthermore, these two alternatives are also the key elements for the order statistics, entropies, and records (upper and lower) density functions. The other odds function measures can be chosen as the ratio of the identities in upper and lower records [ log G ¯ ( t ) ] / [ log G ( t ) ] (cumulative hazard rate function divided by the cumulative reversed hazard function), for instance:
(i)
Ratio of Lehmann alternatives ( G α ( t ) / [ 1 G ¯ α ( t ) ] ) , where α > 0 is the power parameter (see Gupta et al. [1]);
(ii)
Log odds function ( log G ( t ) / G ¯ ( t ) ) (see Al-Aqtash et al. [2]);
(iii)
Logit function [ log ( t / ( 1 t ) ] (see Torabi and Montazeri [3] and Zubair et al. [4]).
Based on the difference of the two log-odds functions, a well-established tool in survival analysis is the proportional odds model, say log G ( t ) / G ¯ ( t ) = h 0 ( t ) + z i β , where h 0 ( t ) = log G 0 ( t ) / G 0 ¯ ( t ) is the baseline log odds function and G 0 ( t ) is the probability of failure by time t for an individual with z = 0 .
Furthermore, Cooray [5] pioneered the concept of the odd function while dealing with probability models, and then he established the odd Weibull model. Gleaton and Lynch [6], while modeling the “strength distribution of an inhomogeneous bundle of brittle elastic fibers under equal load sharing and for checking implementation of the maximum entropy principal (MEP)”, proposed the generalized log-logistic transformation, which led to the odd log-logistic G-class. These two pioneering works motivated researchers and practitioners to develop odds-based G-classes, and to investigate special models from them. Some G-classes based on odd ratio G ( t ) / G ¯ ( t ) , presented in the statistical literature, are included in Table 1. For more details about G-class, see Alzaatreh et al. [7].

2. Background

Several authors have suggested modifications and enhancements to both the exponential and Weibull models in the recent past, with the aim of enhancing their empirical performance and increasing their flexibility. The most peculiar ones are the Lomax exponentiated Weibull model (Ansari and Nofal [34]), extended exponential (ExtE) or generalized exponential(GE) (see Gupta and Kundu [35]), and Nadarajah–Haghighi (NH) (see, Nadarajah and Haghighi [36]). The cdfs of the GE and NH models are
F G E ( x ) = 1 exp λ x α , x > 0
and
F N H ( x ) = 1 exp 1 1 + λ x α , x > 0 ,
respectively, where λ > 0 is a scale parameter and α > 0 is a power (or shape) parameter. Clearly, these two models reduce to the exponential model when α = 1 .
Dimitrakopoulou, Adamidis, and Loukas [37] presented an extension of the Weibull model, in the so-called DAL. The cdf and its corresponding pdf of the DAL distribution can be formulated as
F D A L ( x ) = 1 exp 1 1 + λ x β α , x > 0
and
f D A L ( x ) = λ α β x β 1 1 + λ x β α 1 exp 1 1 + λ x β α ,
where the scale parameter is denoted by λ > 0 , while the shape (or power) parameters are indicated by α > 0 and β > 0 . For α = 1 , the DAL model reduces to the Weibull, and for α = β = 1 , the DAL distribution becomes the exponential model. Nowadays, the DAL model has also been reported as the power generalized Weibull (PGW) distribution (Nikulin and Haghighi [38,39]). However, there is some difference in parametrization of the DAL and PGW models, which is apparent from the cdf of the PGW model, as follows:
F P G W ( x ) = 1 exp 1 1 + x / λ β 1 / γ , x > 0 ,
where λ > 0 is a scale parameter and α > 0 and β > 0 are shape (or power) parameters. Some generalizations/modifications of the DAL model were derived and discussed in the literature. See, for example, exponentiated-DAL (Peña-Ramírez et al. [40]), half-logistic-DAL (Anwar and Bibi [41]), DAL-Logarithmic (Tafakori et al. [42]), MO-DAL (Afify et al. [43]), and transmuted-DAL (Khan [44]).
In the literature, some odd-based G-classes have been discussed as extensions of the exponential or Weibull models, for instance, Bourguignon et al. [9], Tahir et al. [10], Nascimento et al. [28], El-Morshedy et al. [26], El-Morshedy and Eliwa [29], and Ahmad et al. [31] proposed the OW-G, odd generalized Weibull-G (OGW-G), OGE-G, ONH-G, OFW-H, OChen-G, and ODAL-G classes of distributions.
Recently, Hussain et al. [45] defined two new generators (i) H 1 ( x ) = G ( x ) exp { G ¯ ( x ) } and (ii) H 2 ( x ) = G ¯ ( x ) exp { G ( x ) } for bounded unit interval ( 0 , 1 ) and then introduced two new Kumaraswamy G-classes of distributions from them. In this paper, we develop a new generator W [ G ( x ) ] = [ 1 H 2 ( x ) ] / H 2 ( x ) for T [ 0 , ) , which seems less complicated in comparison to earlier published generators, but it performs better when compared with other models.
On the other hand, from the last two decades, discretizing continuous probability models has received wider attention in distribution theory. The phenomenon of discretization occurs when measuring the lifespan of a product or device becomes impractical or impossible on a continuous scale. In such cases, it may be necessary to record lifetimes on a discrete scale rather than a continuous one. This has led to the study of several discrete distributions in the literature. See, for example, Roy [46], Krishna and Pundir [47], Gómez-Déniz [48], Jazi et al. [49], Gómez-Déniz and Calderín-Ojeda [50], Hussain and Ahmad [51], Hussain et al. [52], Para and Jan [53,54], El-Morshedy et al. [55], Eliwa at al. [56], Eliwa and El-Morshedy [57], Eliwa et al. [58], among others. Despite the existence of several discrete probability models in the literature, there is still space for deriving new discretized probability distributions that are appropriate for various areas. To address this, our paper presents a flexible generator of discrete distributions, known as the discrete new odd DAL-G (DNODAL-G) family, which can cater to various conditions. Our proposal for introducing new G-classes is as follows:
  • Generate probability models (ProM) with asymmetric “negatively-skewed, positively-skewed” or symmetric shapes;
  • Define special ProM with all kinds of risk/failure rate functions;
  • Propose ProM suitable for analyzing and discussing both over- and under-dispersed data;
  • Develop ProM for modeling/analyzing both lifetime and counting data sets;
  • Provide ProM that consistently produces a better fit than other ProM built using the same underlying model, in addition to other ProM known in the literature.
The article is organized as follows. A new odd G-class of distributions is introduced in Section 3. Some mathematical properties of a new G-class such as a linear representation for the density, moments, generating function, and estimation of the model parameters are addressed in Section 4. A new model (a special case of the newly proposed G-class for continuous rv) is studied in Section 5 along with a Monte Carlo simulation study. The new discrete odd G-class along with a sub-model is defined, and a Monte Carlo simulation study is investigated in Section 6. Empirical investigation of the proposed models is reported in Section 7 by means of real-life data sets. In Section 8, we conclude our paper with some remarks.

3. The New Odd DAL G-Class

Let T be an rv representing the lifetime of a stochastic system having a baseline G ( x ) distribution. If the rv X represents the odd ratio, then the risk that a system will not be working at time x is given by [ 1 H 2 ( x ) ] / H 2 ( x ) . Therefore, the randomness of X can be modelled by the cdf
F ( x ; λ , α , β , ξ ) = P r X x = Π 1 H 2 ( x ; ξ ) H 2 ( x ; ξ ) ,
where Π is the cdf of T, H 2 ( x ; ξ ) = G ¯ ( x ; ξ ) exp { G ( x ; ξ ) } and then W [ G ( x ; ξ ) ] = [ 1 H 2 ( x ; ξ ) ] / H 2 ( x ; ξ ) . The cdf of the odd DAL-G (NODALG) class is defined as
F ( x ; λ , α , β , ξ ) = 0 [ 1 H 2 ( x ) ] / H 2 ( x ) λ α β t β 1 ( 1 + λ t β ) α 1 exp { 1 ( 1 + λ t β ) α } d t = 1 exp 1 1 + λ 1 G ¯ ( x ; ξ ) exp { G ( x ; ξ ) } G ¯ ( x ; ξ ) exp { G ( x ; ξ ) } β α , x > 0 ,
where λ > 0 is a scale parameter, α > 0 and β > 0 are shape parameters, and ξ is the vector of the baseline parameters. The pdf corresponding to Equation (7) can be expressed as
f ( x ; λ , α , β , ξ ) = α β λ g ( x ; ξ ) G ( x ; ξ ) exp { G ( x , ξ ) } 1 G ¯ ( x ; ξ ) exp { G ( x ; ξ ) } β 1 × G ¯ ( x ; ξ ) exp { G ( x ; ξ ) } β 1 1 + λ 1 G ¯ ( x ; ξ ) exp { G ( x ; ξ ) } G ¯ ( x ; ξ ) exp { G ( x ; ξ ) } β α 1 × exp 1 1 + λ 1 G ¯ ( x ; ξ ) exp { G ( x ; ξ ) } G ¯ ( x ; ξ ) exp { G ( x ; ξ ) } β α .
Henceforth, the rv X with density (7) is denoted by X N O D A L G ( λ , α , β ) . The hazard rate function (hrf) h ( x ) of X has the form
h ( x ) = α β λ g ( x ; ξ ) G ( x ; ξ ) exp { G ( x , ξ ) } 1 G ¯ ( x ; ξ ) exp { G ( x ; ξ ) } β 1 × G ¯ ( x ; ξ ) exp { G ( x ; ξ ) } β 1 1 + λ 1 G ¯ ( x ; ξ ) exp { G ( x ; ξ ) } G ¯ ( x ; ξ ) exp { G ( x ; ξ ) } β α 1 .
Here, we let G ( x ; ξ ) = G ( x ) , G ¯ ( x ; ξ ) = G ¯ ( x ) and g ( x ; ξ ) = g ( x ) to omit the dependence of the parameters.
Proposition. 
Following [37], if X N O D A L G ( λ , α , β ) , then the subsets of our proposed G-class are
(i)
If Y = G ( x ) , then F Y ( y ) = 1 e 1 1 + λ 1 ( 1 y ) exp ( y ) ( 1 y ) exp ( y ) β α ,  for  0 < y < 1 ;
(ii)
If Y = log G ( x ) 1 G ( x ) , then X O D A L ( λ , α , β ) ;
(iii)
If Y = log G ( x ) 1 G ( x ) β , then X O N H ( λ , α ) ;
(iv)
If Y = 1 G ¯ ( x ) exp [ G ( x ) ] G ¯ ( x ) exp [ G ( x ) ] , then X D A L ( λ , α , β ) ;
(v)
If Y = 1 G ¯ ( x ) exp [ G ( x ) ] G ¯ ( x ) exp [ G ( x ) ] β , then X N H ( λ , α ) ;
(vi)
If Y = 1 + 1 G ¯ ( x ) exp [ G ( x ) ] G ¯ ( x ) exp [ G ( x ) ] α 1 , then X W e i ( λ , β ) ;
(vii)
If Y = 1 + 1 G ¯ ( x ) exp [ G ( x ) ] G ¯ ( x ) exp [ G ( x ) ] β α 1 , then X E x p ( λ ) ;
(viii)
If Y = log 1 + 1 G ¯ ( x ) exp [ G ( x ) ] G ¯ ( x ) exp [ G ( x ) ] β , then X Modified   Extreme   Value ( α ) ;
(ix)
If Y = log 1 + 1 G ¯ ( x ) exp [ G ( x ) ] G ¯ ( x ) exp [ G ( x ) ] β 1 / β , then X Power   exponential ( α ) ;
where Y is a random variable that can take different forms of probability generators.

4. Properties of the NKw-G Family

A G-class or a model is known from some important characteristics which they exhibit mathematically or graphically. In this segment, some mathematical and statistical features of the NKw-G class are derived, which will be useful for the readers.

4.1. Quantile Function

The quantile function (qf) is a useful statistical measure that is helpful in obtaining some useful properties, including simulation study. The qf of the NODAL-G class can be expressed as
Q ( u ) = G 1 1 + W e 1 1 + 1 λ { 1 log ( 1 u ) } 1 α 1 1 β 1 ,
where the u follows a uniform distribution over an interval ( 0 , 1 ) , G 1 is the inverse function of base line cdf, and the Lambert-W function is the inverse function of F ( w ) = w e w . The power series expansion for F ( z ) = P r o d u c t L o g [ z ] using the software Mathematica 12 yields the principal solution for w in F ( w ) = w e w = z
W ( z ) = z z 2 + 3 z 3 2 8 z 4 3 + 125 z 5 24 54 z 6 5 + 16,807 z 7 720 16,384 z 8 315 + 531,441 z 9 4480 156,250 z 10 567 ,
where P r o d u c t L o g [ z ] is used as Lambert-W function in the software Mathematica.

4.2. Linear Representation

Here, a useful expansion for Equation (8) is derived. By utilizing the exponential power series in Equation (7), we can write
F ( x ) = 1 e i = 0 ( 1 ) i i ! 1 + λ 1 G ¯ ( x ; ξ ) exp [ G ( x ; ξ ) ] G ¯ ( x ; ξ ) exp [ G ( x ; ξ ) ] β α i .
For a real non-integer, the generalized binomial expansion holds ( 1 + z ) a = j = 0 a j z j , and then applying in Equation (9) gives
F ( x ) = e i = 0 j = 1 λ j ( 1 ) i + 1 i ! α i j G ¯ ( x ; ξ ) exp { G ( x ; ξ ) } β j 1 G ¯ ( x ; ξ ) exp { G ( x ; ξ ) } β j .
Using the previous expansion and the exponential power series, we obtain after some algebra
F ( x ) = p = 0 r = 0 V p , r G ( x ; ξ ) p + r ,
where
V p , r = i , m = 0 j = 1 e λ j ( 1 ) i + m + p + r + 1 i ! p ! α i j β j m β j + m r .
By differentiating Equation (10), the G-class density follows as
f ( x ) = p = 0 r = 0 V p , r ( p + r ) g ( x ; ξ ) G ( x ; ξ ) p + r 1 .
Equation (11) reveals that the NODAL-G family density is a linear combination of exponentiated-G (exp-G) densities. Then, some statistical properties of X can be obtained from Equation (11) and well-established properties of the exp-G distributions.

4.3. Moments

In this segment, the ordinary moment (om), lower incomplete moment (lincm), and moment generating function (mgf) are derived. Let Y p , r be an rv having the exp-G family with power parameter ( p + r ) . First, the sth om of X, say E ( X s ) = x s f ( x ) d x , can be expressed from Equation (11) as
E ( X s ) = p , r = 0 V p , r E ( Y p , r s ) = p , r = 0 ( p + r ) V p , r τ s , p + r ,
where τ s , p + r = x s G ( x ; ξ ) p + r 1 g ( x ; ξ ) d x = 0 1 Q G ( u ; ξ ) s u p + r d u , and Q G ( u ; ξ ) is qf of baseline G. The well-known relationships can be used to derive the central moments and cumulants of X from Equation (12). Second, the sth lincm of X, say m s ( y ) = y x s f ( x ) d x , is
m s ( y ) = p , r = 0 V p , r y x s h p + r ( x ) d x = p , r = 0 ( p + r ) V p , r 0 G ( y ; ξ ) ) Q G ( u ; ξ ) s u p + r d u .
For most G distributions, it is possible to numerically evaluate the last two integrals in the equation. The first lincm, m 1 ( y ) , is useful in constructing popular measures such as Bonferroni and Lorenz curves in various fields such as demography, economics, reliability, medicine, and insurance. In addition, it can also be applied to determine the sum of the deviations from the mean and median of X. Furthermore, the mgf M X ( t ) = E ( e t , X ) of X can be derived from (11).
M X ( t ) = p , r = 0 V p , r M p , r ( t ) = p , r = 0 ( p + r ) V p , r ρ p , r ( t ) ,
where M p , r ( t ) is the mgf of Y p , r and ρ p , r ( t ) = 0 1 exp [ t Q G ( u ; ξ ) ] u p + r d u . Therefore, we can derive the mgfs of several particular NODAL-G models directly from Equation (14) and exp-G generating functions.

4.4. Maximum Likelihood Estimation

Uncensored maximum likelihood estimation, in which all of the data are observed without any censoring, is a technique for estimating the parameters of a probability distribution based on a sample of data. This segment deals with the estimation of the unknown NODAL-G class parameters via the maximum likelihood (ML) approach. The ML estimation refers to a method of estimating unknown parameters by selecting values that maximize the likelihood of observing a given set of data. This technique is often used in various types of statistical modeling, such as regression and classification. It is a popular approach due to its simplicity and the fact that it is easy to implement. At its core, ML estimation is a mathematical approach for finding the probability distribution of a unknown variable based on a given sample of data. This is achieved by finding the maximum value of the likelihood function, which is based on the probability distribution of the given data. The likelihood function is computed by taking the product of the probability of the observations in the data set. Max likelihood estimation is used in many areas of study, including economics, biology, engineering, and computer science. As an example, in economics, it is used to make predictions about the probability of future events based on past observations. In biology, it is used to estimate gene frequencies and relationships between genes. For more details about the ML approach and its statistical properties, see Casella and Berger [59]. The log-likelihood function ( Θ ) for the parameter vector Θ = ( λ , α , β , ξ ) can be derived from n observations/notes x 1 , , x n .
n = n log α + n log β + n log λ + i = 1 n log g ( x i ; ξ ) + ( β 1 ) i = 1 n log 1 G ¯ ( x i ; ξ ) exp { G ( x i ; ξ ) } + i = 1 n log G ( x i ; ξ ) exp { G ( x i , ξ ) } + ( α 1 ) i = 1 n log 1 + λ 1 G ¯ ( x i ; ξ ) exp { G ( x i ; ξ ) } G ¯ ( x i ; ξ ) exp { G ( x i ; ξ ) } β ( β + 1 ) i = 1 n log G ¯ ( x i ; ξ ) exp { G ( x i ; ξ ) } + i = 1 n 1 1 + λ 1 G ¯ ( x i ; ξ ) exp { G ( x i ; ξ ) } G ¯ ( x i ; ξ ) exp { G ( x i ; ξ ) } β α .
One way to find the maximum likelihood estimate (MLE) Θ ^ of Θ is to maximize the likelihood function ( Θ ) . There are several numerical optimization routines available in different programming languages such as R, SAS, and Ox that can be used to maximize ( Θ ) . For instance, the optim function in R, the PROC NLMIXED in SAS, and the (sub-routine MaxBFGS) Ox can be used for this purpose.

5. The NODAL-Weibull Distribution

In this Section, we consider a special model of the NODAL-G class, the NODAL-Weibull (NODALW) distribution, by taking the Weibull as a baseline model. The cdf and pdf of Weibull distributions are G ( x ) = 1 e k x c and g ( x ) = c k x c 1 e k x c , respectively, where k > 0 is a scale parameter and c > 0 is a shape parameter. By setting λ = 1 , the cdf of the NODALW distribution reduces to
F ( x ) = 1 exp 1 1 + 1 e k x c exp 1 e k x c e k x c exp 1 e k x c β α .
The pdf corresponding to Equation (15) is
f ( x ) = α β c k x c 1 e k x c 1 e k x c exp 1 e k x c 1 e k x c exp 1 e k x c β 1 × e k x c exp 1 e k x c β 1 1 + 1 e k x c exp 1 e k x c e k x c exp 1 e k x c β α 1 × exp 1 1 + 1 e k x c exp 1 e k x c e k x c exp 1 e k x c β α .
Henceforth, let X N O D A L W ( α , β , k , c ) be the rv with density (16). The hrf of X is
h ( x ) = α β c k x c 1 e k x c 1 e k x c exp 1 e k x c 1 e k x c exp 1 e k x c β 1 × e k x c exp 1 e k x c β 1 1 + 1 e k x c exp 1 e k x c e k x c exp 1 e k x c β α 1 .
The sketches of the density and failure rate of X are plotted in Figure 1 for some parameter values. Figure 1a displays the uni-modal (right-skewed and left-skewed) and reversed-J shapes of the density of X. Figure 1b exhibits the failure rate shapes of X such as increasing, decreasing, and bathtub.

5.1. Linear Representation of NODALW Model

The NODALW density follows from Equation (11) as
f ( x ) = p = 0 r = 0 V p , r k ( p + r ) c x c 1 e k x c 1 e k x c p + r 1 .
Using the generalized binomial expansion in Equation (17), we obtain
f ( x ) = τ = 0 ω τ k ( τ + 1 ) c x c 1 e k ( τ + 1 ) x c ,
where ω τ = ( 1 ) τ τ + 1 p + r 1 τ p , r = 0 ( p + r ) V p , r . Equation (18) reveals that the NODALW density has a linear representation in terms of Weibull densities. So, several of its structural properties can be obtained from the Weibull density. The qf of the NODALW distribution is
Q ( u ) = 1 k log W e 1 1 + 1 log ( 1 u ) 1 α 1 1 β 1 1 c , 0 < u < 1 .

5.2. Properties of NODALW Model

Let Z p be an rv with Weibull density π ( x ; p , c ) . Then, some quantities of X can follow from those of Z p . First, the sth om of X can be expressed as
μ s = Γ s c + 1 τ = 0 ω τ [ k ( τ + 1 ) ] s / c ,
where Γ ( a ) = 0 t a 1 e t d t . Using Equation (20), we can recursively calculate the cumulants ( κ s ) of X. Specifically, the sth cumulant is determined by subtracting the sum of products of previous cumulants ( κ k ) and raw moments ( μ s k ), where the products are taken over all k from 1 to s 1 , inclusive. The formula is expressed as
κ s = μ s k = 1 s 1 s 1 k 1 , κ k , μ s k
The first cumulant, κ 1 , is equal to the first raw moment, μ 1 . The skewness, γ 1 , and kurtosis, γ 2 , of X can be obtained by dividing the third and fourth standardized cumulants by the square of the second standardized cumulant, respectively. Some plots of skewness (sk) and kurtosis (ku) for X are presented in Figure 2. It can be seen that the proposed class can be used to discuss the different forms of kurtosis.
The sth incomplete moment of X, denoted by m s ( y ) = E ( X s X y ) = 0 y x s f ( x ) d x , is easily obtained by changing variables from the lower incomplete gamma function γ ( s , x ) = 0 x s 1 e x d x when calculating the corresponding moment of Z p . Then, we obtain
m s ( z ) = τ = 0 ω τ γ s c + 1 , k ( τ + 1 ) x c ( k ( τ + 1 ) ) s / c .

5.3. Simulation Study: NODALW Model

Simulation studies are a popular and effective method of testing estimator performance in a variety of scenarios. By running several simulations, it is possible to approximate real-world performance and gain insight into how well the estimator will function in the field. The first step of a simulation study is to establish the parameters of the experiment. This involves setting up criteria for the data to be sampled, including number of samples, sample size, and the sampling process. Once the parameters have been established, the next step is to generate a simulated dataset that matches the specified criteria. This dataset should have features that are as close as possible to the features of the real-world data. Once the dataset has been generated, the estimator can be applied to the data. Depending on the type of estimator, different metrics may be used to measure the performance of the estimator. These metrics may include root mean squared error, mean absolute error, and R-squared. In this Section, we conduct a Monte Carlo simulation study to assess the performance of the MLEs of the parameters α , β , c, and k. The random numbers of size 50, 100, 200, and 500 are generated by the inversion method and are repeated N = 1000 times for each sample size. The sample average biases (Bias), coverage probabilities (CPs), and mean-squared errors (MSEs) of the estimates are calculated. The following formulas are used:
B i a s ( θ ^ ) = i = 1 N θ i ^ N θ and M S E ( θ ^ ) = i = 1 N ( θ i ^ θ ) 2 N
and
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 c i ^ , s k i ^ are the standard errors of the MLEs. These quantities for some values of α , β , c, and k are reported in Table 2, Table 3 and Table 4. The figures in these tables reveal that the MLEs perform well for estimating the parameters of the NODALW distribution.

6. Discrete NODAL-G Family

According to a survival discretization approach, the rv X is said to have the discrete NODAL-G (dNODAL-G) class if its cdf can be formulated as
V ( x ; λ , α , β , ξ ) = 1 exp 1 1 + λ 1 G ¯ ( x + 1 ; ξ ) exp G ( x + 1 ; ξ ) G ¯ ( x + 1 ; ξ ) exp G ( x + 1 ; ξ ) β α , x N 0 ,
where λ > 0 is a scale parameter, α > 0 and β > 0 are shape parameters, ξ is the vector of the baseline parameters, and N 0 = 0 , 1 , 2 , . The probability mass function (pmf) and hrf corresponding to Equation (22) are
v ( x ; λ , α , β , ξ ) = exp 1 1 + λ 1 G ¯ ( x ; ξ ) exp G ( x ; ξ ) G ¯ ( x ; ξ ) exp G ( x ; ξ ) β α exp 1 1 + λ 1 G ¯ ( x + 1 ; ξ ) exp G ( x + 1 ; ξ ) G ¯ ( x + 1 ; ξ ) exp G ( x + 1 ; ξ ) β α ,
and h ( x ; λ , α , β , ξ ) = v ( x ; λ , α , β , ξ ) 1 V ( x 1 ; λ , α , β , ξ ) , respectively, where x N 0 .

6.1. The DNODAL-Geometric (DNODALGeo) Distribution

Consider the cdf of the geometric (Geo) model. Then, the pmf of the DNODALGeo distribution is
V ( x ; λ , α , β , p ) = exp 1 1 + λ 1 p x + 1 exp 1 p x + 1 p x + 1 exp 1 p x + 1 β α , x N 0 ,
where λ > 0 , α > 0 , β > 0 , and 0 < p < 1 . For convenience, let λ = 1 in Equation (24). Figure 3 and Figure 4 display the pmf and hrf of the DNODALGeo distribution for some parameter values.
PMF can be either monomodal or bimodal and can be used to analyze different types of data (positively skewed, negatively skewed, as well as symmetric). Moreover, hrf can be either an increment, a fixed increment, a J-shape, or a bathtub.

6.2. Simulation

We generate a random variable X from the DNODALGeo distribution by generating the value Z from the continuous model and then discretize X = [ Z ] , which is the largest integer less than or equal to Z. We generate 1 , 000 samples of size n = 20 , 22 , 24 , , 150 from the DNODALGeo ( 0.5 , 0.9 , 0.8 ) and DNODALGeo ( 1.5 , 1.5 , 0.5 ) models, respectively. The empirical results are given in Figure 5 and Figure 6, respectively.
It is clear from Figure 5 and Figure 6 that the biases and MSEs always tend to zero when n increases, which shows the consistency of the estimators.

7. Applications

In this Section, we illustrate the empirical importance of the NODALW and DNODALGeo models by means of three real-life data sets. First, two data sets are utilized to illustrate the flexibility of the NODALW distribution. Second, a third data set is used to test the usefulness of the DNODALGeo distribution.

7.1. Empirical Illustration of the NODAL-G Family

Here, we compare the NODALW distribution with some well-established competitive models: Kumaraswamy–Weibull (KwW) (see Cordeiro et al. [60]), beta-Weibull (BW) (see Lee et al. [61]), exponentiated-generalized Weibull (EGW) (see Oguntunde et al. [62]), McDonald-Weibull (McW) (see Cordeiro et al. [63]), gamma-Weibull (GaW) (see Cordeiro et al. [64]), and Weibull (W) to prove the flexibility of the new family. The cdf and pdf of NODALW distribution are, respectively, given as
F ( x ) = 1 exp 1 1 + 1 e k x c exp 1 e k x c e k x c exp 1 e k x c β α .
and
f ( x ) = α β c k x c 1 e k x c 1 e k x c exp 1 e k x c 1 e k x c exp 1 e k x c β 1 × e k x c exp 1 e k x c β 1 1 + 1 e k x c exp 1 e k x c e k x c exp 1 e k x c β α 1 × exp 1 1 + 1 e k x c exp 1 e k x c e k x c exp 1 e k x c β α .
  • Data Set 1. Air Conditioning Data.  The data are taken from Kus [65] representing the numbers of the successive failures for an air conditioning system. The shape of the data can be discussed through Figure 7. The data were found to be asymmetric and some extreme observations were reported.
  • Data Set 2. Precipitation Data.  The data are taken from Katz et al. [66] and Asgharzadeh et al. [67] representing 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. The shape of the data can be displayed in Figure 8.
The NODALW model and other competitive models are fitted to these two data sets using the AdequacyModel package for the R statistical computing environment written by Marinho et al. [68].
The MLEs ( ^ ) are used to evaluate the log-likelihood function, while various goodness-of-fit statistics (“GoFS”), such as 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 employed to compare models. A good fit is indicated by lower values of these statistics and higher P-values of the KS statistic.
The values of the GoFS in Table 5 and Table 7 show that the NODALW model gives small values for these statistics and then it provides the best fit as compared to other fitted distributions (KwW, BW, EGW, McW, GaW, and W) to the two data sets. Table 6 and Table 8 report the MLEs and their standard errors (SEs) for the NODALW model and other competitive models. The plots in Figure 9 and Figure 10 also support our claim. To establish the unique property of the maximum likelihood estimators, the profile log-maximum likelihood function (pllf) was plotted for each parameter for the first and second datasets. It can be noted that the values of the estimators gave the maximum likelihood function the largest value; see Figure 11 and Figure 12. As can be seen, the form of the maximum likelihood function is unimodal for each parameter.
Figure 7. Non-parametric visualization plots for data set I.
Figure 7. Non-parametric visualization plots for data set I.
Mathematics 11 02929 g007
Figure 8. Non-parametric visualization plots for data set II.
Figure 8. Non-parametric visualization plots for data set II.
Mathematics 11 02929 g008
Figure 9. Plots of estimated density, cdf, Kaplan–Meier (K-M), and hrf plots for data set 1.
Figure 9. Plots of estimated density, cdf, Kaplan–Meier (K-M), and hrf plots for data set 1.
Mathematics 11 02929 g009
Figure 10. Plots of estimated density, cdf, Kaplan–Meier (K-M), and hrf plots for data set 2.
Figure 10. Plots of estimated density, cdf, Kaplan–Meier (K-M), and hrf plots for data set 2.
Mathematics 11 02929 g010
Figure 11. The pllf of data set 1.
Figure 11. The pllf of data set 1.
Mathematics 11 02929 g011
Figure 12. The pllf of data set 2.
Figure 12. The pllf of data set 2.
Mathematics 11 02929 g012

7.2. Empirical Illustration of the DNODAL-G Family

Here, we illustrate the usefulness of the DNODALGeo model by means of an application to real count data. The data set represents the count of kidney cysts using steroids (see Chan et al. [69]). The shape of the data can be seen in Figure 13. The fitted distributions are compared using the AIC, CAIC, HQIC, and Chi-square ( χ 2 ), having a degree of freedom (df) and its p-value. The competitive fitted models are reported in Table 9.
The MLEs and their corresponding SEs are listed in Table 10, whereas Table 11 and Table 12 give the GoFS, expected frequencies (ExFr), and observed frequencies (ObFr), respectively.
The DNODALGeo model performs better than all other tested models based on the numbers in Table 11 and Table 12. Figure 14 supports the claims from these tables, and it is noted that data set 3 is explained by this model. Figure 15 shows the contour plot of log-Likelihood function of the DNODALGeo for the third data set.

8. Concluding Remarks and Future Work

In this article, a new odd DAL-G family of models is presented from a new class/generator W [ G ( x ) ] = ( 1 G ¯ ( x ; ξ ) exp { G ( x ; ξ ) } ) / [ G ¯ ( x ; ξ ) exp { G ( x ; ξ ) } ] for T ( 0 , ) . The new probability family W [ G ( · ) ] involves a different function of the cdf instead of existing generators. We obtain some structural properties of this new continuous and discussed discrete odd DAL-G family, and also studied some properties of the special models called the new odd DAL-Weibull (NODALW) and discrete new odd DAL-geometric (DNODALGeo) distributions. Both of two sub-models can be used to discuss asymmetric and symmetric data under different kinds of kurtosis. Furthermore, the two sub-models can be applied to discuss several shapes of risk/hazard rates. We compared the NODALW distribution with the well-known extended Weibull models (KwW, BW, EGW, McW, GaW, W) via six popular test statistics. Similarly, we compare the DNODALGeo distribution with the well-known extended models’ (Geo, GGeo, DR, DIR, DIW, DLi-I, DLi-II, DLi-III, NeBi, Poi, DPa, DB-XII, DLogL, DFx-I, DLo) distributions using these test statistics. We found that the new generated distributions provide better estimates and minimum values of the test statistics. The new NODALW and DNODALGeo models outperform the above-described competitive models on the basis of numerical and graphical analysis. We foresee that the new family/class will be able to attract readers and applied statisticians. As a future work, the bivariate extension of the proposed generators with its applications will be discussed. Furthermore, some prediction models will be analyzed based on these generators.

Author Contributions

Conceptualization, M.S.E. and M.A.H.; Methodology, M.E.-M.; Software, M.A.H., B.A. and M.E.-M.; Validation, M.A.H. and M.E.-M.; Formal analysis, M.S.E., M.H.T. and A.A.-B.; Resources, B.A and M.E.-M.; Data curation, M.S.E., B.A. and A.A.-B.; Writing—original draft, M.H.T., M.S.E. and M.A.H.; Writing—review & editing, M.S.E., M.A.H. and A.A.-B.; Visualization, B.A. and M.E.-M.; Supervision, M.H.T. All authors have read and agreed to the published version of the manuscript.

Funding

The paper did not receive any fund.

Data Availability Statement

The data sets are available in the paper.

Acknowledgments

The researchers would like to thank the Deanship of Scientific Research, Qassim University for funding publication of this project.

Conflicts of Interest

The authors declare no conflict of interests.

References

  1. Gupta, R.C.; Gupta, P.L.; Gupta, R.D. Modeling failure time data by Lehman alternatives. Commun. Stat. Theory Methods 1998, 27, 887–904. [Google Scholar] [CrossRef]
  2. Al-Aqtash, R.; Lee, C.; Famoye, F. Gumbel-Weibull distribution: Properties and applications. J. Mod. Appl. Stat. Methods 2014, 13, 201–225. [Google Scholar] [CrossRef]
  3. Torabi, H.; Montazeri, N.H. The logistic-uniform distribution and its application. Commun. Stat. Simul. Comput. 2014, 43, 2551–2569. [Google Scholar] [CrossRef]
  4. Zubair, M.; Pogany, T.K.; Cordeiro, G.M.; Tahir, M.H. The log-odd normal generalized family of distributions with application. An. Acad. Bras. Ciênc. 2019, 91, e20180207. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Cooray, K. Generalization of the Weibull distribution: The odd Weibull family. Stat. Methodol. 2006, 6, 265–277. [Google Scholar] [CrossRef]
  6. Gleaton, J.U.; Lynch, J.D. Properties of generalized log-logistic families of lifetime distributions. J. Probab. Stat. Sci. 2006, 4, 51–64. [Google Scholar]
  7. Alzaatreh, A.; Famoye, F.; Lee, C. A new method for generating families of continuous distributions. Metron 2013, 71, 63–79. [Google Scholar] [CrossRef] [Green Version]
  8. Torabi, H.; Montazeri, N.H. The gamma-uniform distribution and its application. Kybernetika 2012, 48, 16–30. [Google Scholar]
  9. Bourguignon, M.; Silva, R.B.; Cordeiro, G.M. The Weibull-G family of probability distributions. J. Data Sci. 2014, 12, 53–68. [Google Scholar] [CrossRef]
  10. Tahir, M.H.; Cordeiro, G.M.; Alizadeh, M.; Mansoor, M.; Zubair, M.; Hamedani, G.G. The odd generalized exponential family of distributions with applications. J. Stat. Distrib. Appl. 2015, 2, 1. [Google Scholar] [CrossRef] [Green Version]
  11. Hassan, A.S.; Hemeda, S.E. A new family of additive Weibull-generated distributions. Int. J. Math. Appl. 2016, 4, 151–164. [Google Scholar]
  12. Gomes-Silva, F.S.; Percontini, A.; de Brito, E.; Ramos, M.W.; Venancio, R.; Cordeiro, G.M. The odd Lindley-G family of distribution. Austrian J. Stat. 2017, 46, 65–87. [Google Scholar] [CrossRef] [Green Version]
  13. Cordeiro, G.M.; Alizadeh, M.; Ramires, T.G.; Ortega, E.M.M. The generalized odd half-Cauchy family of distributions: Properties and applications. Commun. Stat. Theory Methods 2017, 46, 5685–5705. [Google Scholar] [CrossRef]
  14. Afify, A.Z.; Altun, E.; Alizadeh, M.; Ozel, G.; Hamedani, G.G. The odd exponentiated half-logistic-G Family: Properties, characterizations and applications. Chil. J. Stat. 2017, 8, 65–91. [Google Scholar]
  15. Jamal, F.; Nasir, M.A.; Tahir, M.H.; Montazeri, N.H. The odd Burr-III family of distributions. J. Stat. Appl. Probab. 2017, 6, 105–122. [Google Scholar] [CrossRef]
  16. Yousof, H.M.; Afify, A.Z.; Hamedani, G.G. The Burr X generator of distributions for lifetime data. J. Stat. Theory Appl. 2017, 16, 288–305. [Google Scholar] [CrossRef] [Green Version]
  17. Cordeiro, G.M.; Yousof, H.M.; Ramires, T.G.; Ortega, E.M.M. The Burr XII system of densities: Properties, regression model and applications. J. Stat. Comput. Simul. 2018, 88, 432–456. [Google Scholar] [CrossRef]
  18. Haq, M.A.U.; Elgarhy, M. The odd Fréchet-G family of probability distributions. J. Stat. Appl. Probab. 2018, 7, 189–203. [Google Scholar] [CrossRef]
  19. Hassan, A.S.; Nassr, S.G. The inverse Weibull generator of distributions: Properties and applications. J. Data Sci. 2018, 16, 723–742. [Google Scholar] [CrossRef]
  20. Alizadeh, M.; Altun, E.; Cordeiro, G.M.; Rasekhi, M. The odd power-Cauchy family of distributions: Properties, regression models and applications. J. Stat. Comput. Simul. 2018, 88, 785–807. [Google Scholar] [CrossRef]
  21. Maiti, S.S.; Pramanik, S. A generalized Xgamma generator family of distributions. arXiv 2018, arXiv:1805.03892. [Google Scholar]
  22. Hassan, A.S.; Nassr, S.G. Power Lindley-G family of distributions. Ann. Data Sci. 2019, 6, 189–210. [Google Scholar] [CrossRef]
  23. Korkmaz, M.C.; Altun, E.; Yousof, H.M.; Hamedani, G.G. The odd power Lindley generator of probability distributions: Properties, characterizations and regression modeling. Int. J. Stat. Probab. 2019, 8, 70–89. [Google Scholar] [CrossRef]
  24. Cordeiro, G.M.; Afify, A.Z.; Ortega, E.M.M.; Suzuki, A.K.; Mead, M.E. The odd Lomax generator of distributions: Properties, estimation and applications. J. Comput. Appl. Math. 2019, 347, 222–237. [Google Scholar] [CrossRef]
  25. Kharazmi, O.; Saadatinik, A.; Alizadeh, M.; Hamedani, G.G. Odd hyperbolic cosine-FG family of lifetime distributions. J. Stat. Theory Appl. 2019, 18, 387–401. [Google Scholar] [CrossRef] [Green Version]
  26. El-Morshedy, M.; Eliwa, M.S. The odd flexible Weibull-H family of distributions: Properties and estimation with applications to complete and upper record data. Filomat 2019, 33, 2635–2652. [Google Scholar] [CrossRef]
  27. Aldahlan, M.A.; Afify, A.Z.; Ahmed, A.-H.N. The odd inverse Pareto-G class: Properties and applications. J. Nonlinear Sci. Appl. 2019, 12, 278–290. [Google Scholar] [CrossRef] [Green Version]
  28. Nascimento, D.C.; Silva, K.F.; Cordeiro, G.M.; Alizadeh, M.; Yousof, H.M.; Hamedani, G.G. The odd Nadarajah-Haghighi family of distributions: Properties and applications. Stud. Sci. Math. Hung. 2019, 56, 185–210. [Google Scholar] [CrossRef]
  29. El-Morshedy, M.; Eliwa, M.S.; Afify, A.Z. The odd Chen generator of distributions: Properties and estimation methods with applications in medicine and engineering. J. Natl. Sci. Found. Sri Lanka 2020, 48, 113–130. [Google Scholar]
  30. Anzagra, L.; Sarpong, S.; Nasiru, S. Odd Chen-G family of distributions. Ann. Data Sci. 2020, 9, 369–391. [Google Scholar] [CrossRef]
  31. Ahmad, Z.; Elgarhy, M.; Hamedani, G.G.; Butt, N.S. Odd generalized N-H generated family of distributions with application to exponential model. Pak. J. Stat. Oper. Res. 2020, 16, 53–71. [Google Scholar] [CrossRef]
  32. Nasir, M.A.; Tahir, M.H.; Chesneau, C.; Jamal, F.; Shah, M.A.A. The odds generalized gamma-G family of distributions: Properties, regression and applications. Statistica 2020, 80, 3–38. [Google Scholar] [CrossRef]
  33. Ishaq, A.I.; Abiodun, A.A. The Maxwell-Weibull distribution in modeling lifetime data sets. Ann. Data Sci. 2020, 7, 639–662. [Google Scholar] [CrossRef]
  34. Ansari, S.I.; Nofal, Z.M. The Lomax exponentiated Weibull model. Jpn. J. Stat. Data Sci. 2021, 4, 21–39. [Google Scholar] [CrossRef]
  35. Gupta, R.D.; Kundu, D. Generalized exponential distributions. Aust. N. Z. J. Stat. 1999, 41, 173–188. [Google Scholar] [CrossRef]
  36. Nadarajah, S.; Haghighi, F. An extension of the exponential distribution. Statistics 2011, 45, 54–558. [Google Scholar] [CrossRef]
  37. Dimitrakopoulou, T.; Adamidis, K.; Loukas, S. A lifetime distribution with an upside–down bathtub-shaped hazard function. IEEE Trans. Reliab. 2007, 56, 308–311. [Google Scholar] [CrossRef]
  38. Nikulin, M.; Haghighi, F. A chi-squared test for the generalized power Weibull family for the head-and-neck cancer censored data. J. Math. Sci. 2006, 133, 1333–1341. [Google Scholar] [CrossRef]
  39. Nikulin, M.; Haghighi, F. On the power generalized Weibull family: Model for cancer censored data. Metron 2009, 67, 75–86. [Google Scholar]
  40. Peńa-Ramírez, F.A.; Guerra, R.R.; Cordeiro, G.M.; Rinho, P.R.D. The exponentiated power generalized Weibull: Properties and applications. An. Acad. Bras. Ciêc. 2018, 90, 2553–2577. [Google Scholar] [CrossRef]
  41. Anwar, M.; Bibi, A. The half-logistic generalized Weibull distribution. J. Probab. Stat. 2018, 12, 8767826. [Google Scholar] [CrossRef]
  42. Tafakori, L.; Pourkhanali, A.; Nadarajah, S. A new lifetime model with different types of failure rate. Commun. Stat. Theory Methods 2018, 47, 4006–4020. [Google Scholar] [CrossRef]
  43. Afify, A.Z.; Kumar, D.; Elbatal, I. Marshall-Olkin power generalized Weibull distribution with applications in engineering and medicine. J. Stat. Theory Appl. 2020, 19, 223–237. [Google Scholar] [CrossRef]
  44. Khan, M.S. Transmuted generalized power Weibull distribution. Thail. Stat. 2018, 16, 156–172. [Google Scholar]
  45. Hussain, M.A. Some New Generalized Kumaraswamy Families of Distributions. Ph.D. Thesis, The Islamia University of Bahawalpur, Bahawalpur, Pakistan, 2020. [Google Scholar]
  46. Roy, D. Discrete Rayleigh distribution. IEEE Trans. Reliab. 2004, 53, 255–260. [Google Scholar] [CrossRef]
  47. Krishna, H.; Pundir, P.S. Discrete Burr and discrete Pareto distributions. Stat. Methodol. 2009, 6, 177–188. [Google Scholar] [CrossRef]
  48. Gómez-Déniz, E. Another generalization of the geometric distribution. Test 2010, 19, 399–415. [Google Scholar] [CrossRef]
  49. Jazi, M.A.; Lai, C.D.; Alamatsaz, M.H. A discrete inverse Weibull distribution and estimation of its parameters. Stat. Methodol. 2010, 7, 121–132. [Google Scholar] [CrossRef]
  50. Gómez-Déniz, E.; Calderín-Ojeda, E. The discrete Lindley distribution: Properties and applications. J. Stat. Comput. Simul. 2011, 81, 1405–1416. [Google Scholar] [CrossRef]
  51. Hussain, T.; Ahmad, M. Discrete inverse Rayleigh distribution. Pak. J. Stat. 2014, 30, 203–222. [Google Scholar]
  52. Hussain, T.; Aslam, M.; Ahmad, M. A two-parameter discrete Lindley distribution. Rev. Colomb. Estad. 2016, 39, 45–61. [Google Scholar] [CrossRef]
  53. Para, B.A.; Jan, T.R. Discrete version of log-logistic distribution and its applications in genetics. International J. Mod. Math. Sci. 2016, 14, 407–422. [Google Scholar]
  54. Para, B.A.; Jan, T.R. On discrete three parameter Burr type XII and discrete Lomax distributions and their applications to model count data from medical science. Biom. Biostat. Int. J. 2016, 4, 1–15. [Google Scholar]
  55. El-Morshedy, M.; Eliwa, M.S.; Altun, E. Discrete Burr-Hatke distribution with properties, estimation methods and regression model. IEEE Access 2020, 8, 74359–74370. [Google Scholar] [CrossRef]
  56. Eliwa, M.S.; Altun, E.; El-Dawoody, M.; El-Morshedy, M. A new three-parameter discrete distribution with associated INAR (1) process and applications. IEEE 2020, 8, 91150–91162. [Google Scholar]
  57. Eliwa, M.S.; El-Morshedy, M. A one-parameter discrete distribution for over-dispersed data: Statistical and reliability properties with estimation approaches and applications. J. Appl. Stat. 2022, 49, 2467–2487. [Google Scholar] [CrossRef] [PubMed]
  58. Eliwa, M.S.; Alhussain, Z.A.; El-Morshedy, M. Discrete Gompertz-G family of distributions for over-and under-dispersed data with properties, estimation, and applications. Mathematics 2020, 8, 358. [Google Scholar] [CrossRef] [Green Version]
  59. Casella, G.; Berger, R.L. Statistical Inference; Cengage Learning: Boston, MA, USA, 2021. [Google Scholar]
  60. Cordeiro, G.M.; Ortega, E.M.M.; Nadarajah, S. The Kumaraswamy Weibull distribution with application to failure data. J. Frankl. Inst. 2010, 347, 1399–1429. [Google Scholar] [CrossRef]
  61. 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]
  62. Oguntunde, P.E.; Odetunmibi, O.A.; Adejum, A.O. On the exponentiated generalized Weibull distribution: A generalization of the Weibull distribution. Indian J. Sci. Technol. 2015, 8, 1–7. [Google Scholar] [CrossRef]
  63. Cordeiro, G.M.; Hashimoto, E.M.; Ortega, E.M.M. The McDonald Weibull model. Statistics 2012, 48, 256–278. [Google Scholar] [CrossRef]
  64. Cordeiro, G.M.; Aristizábal, W.D.; Suárez, D.M.; Lozano, S. The gamma modified Weibull distribution. Chil. J. Stat. 2011, 6, 37–48. [Google Scholar]
  65. Kus, C. A new lifetime distribution. Comput. Stat. Data Anal. 2007, 51, 4497–4509. [Google Scholar] [CrossRef]
  66. Katz, R.W.; Parlange, M.B.; Naveau, P. Statistics of extremes in hydrology. Adv. Water Resour. 2002, 25, 1287–1304. [Google Scholar] [CrossRef] [Green Version]
  67. Asgharzadeh, A.; Bakouch, H.S.; Habibi, M. A generalized binomial exponential 2 distribution: Modeling and applications to hydrologic events. J. Appl. Stat. 2017, 44, 2368–2387. [Google Scholar] [CrossRef]
  68. 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. 2015. Available online: http://cran.rproject.org/web/packages/AdequacyModel/AdequacyModel.pdf (accessed on 12 November 2015).
  69. Chan, S.K.; Riley, P.R.; Price, K.L.; McEldu, F.; Winyard, P.J.; Welham, S.J.; Woolf, A.S.; Long, D.A. Corticosteroid-induced kidney dysmorphogenesis is associated with deregulated expression of known cystogenic molecules, as well as Indian hedgehog. Am. J. Physiol.-Ren. Physiol. 2010, 298, F346–F356. [Google Scholar] [CrossRef] [Green Version]
  70. Poisson, S.D. Probabilité des Jugements en Matié re Criminelle et en Matiére Civile, Précédées des Régles Génerales du Calcul des Probabilitiés; Bachelier: Paris, France, 1837; pp. 206–207. [Google Scholar]
Figure 1. Sketches of (a) the densities, and (b) the hrf of the NODALW model.
Figure 1. Sketches of (a) the densities, and (b) the hrf of the NODALW model.
Mathematics 11 02929 g001
Figure 2. Plots of the (a) sk and (b) ku of the NODALW model.
Figure 2. Plots of the (a) sk and (b) ku of the NODALW model.
Mathematics 11 02929 g002
Figure 3. The pmf of the DNODALGeo model.
Figure 3. The pmf of the DNODALGeo model.
Mathematics 11 02929 g003
Figure 4. The hrf of the DNODALGeo model.
Figure 4. The hrf of the DNODALGeo model.
Mathematics 11 02929 g004
Figure 5. The biases and MSEs for the DNODALGeo(0.5, 0.9, 0.8) model.
Figure 5. The biases and MSEs for the DNODALGeo(0.5, 0.9, 0.8) model.
Mathematics 11 02929 g005
Figure 6. The biases and MSEs for the DNODALGeo(1.5, 1.5, 0.5) model.
Figure 6. The biases and MSEs for the DNODALGeo(1.5, 1.5, 0.5) model.
Mathematics 11 02929 g006
Figure 13. Non-parametric visualization plots for data set III.
Figure 13. Non-parametric visualization plots for data set III.
Mathematics 11 02929 g013
Figure 14. The fitted PMFs to data set 3.
Figure 14. The fitted PMFs to data set 3.
Mathematics 11 02929 g014
Figure 15. Contour plot of pllf for data set 3.
Figure 15. Contour plot of pllf for data set 3.
Mathematics 11 02929 g015
Table 1. Odd ratio for G-classes of distributions.
Table 1. Odd ratio for G-classes of distributions.
No.G-ClassYearAuthors
1Odd log-logistic-G (OLL-G)2006Gleaton and Lynch [6]
2Odd Gamma-G (OGa-G)2012Torabi and Montazeri [8]
3Odd Weibull-G (OW-G)2014Bourguignon et al. [9]
4Odd generalized-exponential-G (OGE-G)2015Tahir et al. [10]
5Odd additive Weibull-G (OAddW-G)2016Hassan and Hemeda [11]
6Odd Lindley-G (OLind-G)2017Gomes-Silva et al. [12]
7Odd half-Cauchy-G (OHCa-G)2017Cordeiro et al. [13]
8Odd half-Logistic-G (OHL-G)2017Affify et al. [14]
9Odd Burr III-G (OBr3-G)2017Jamal et al. [15]
10Odd Burr X-G (OBrX-G)2017Yousof et al. [16]
11Odd Burr-XII-G (OBr-G)2018Cordeiro et al. [17]
12Odd Frechét-G (OFr-G)2018Haq et al. [18]; Hassan and Nassr [19]
13Odd power-Cauchy-G (OPCa-G)2018Alizadeh et al. [20]
14Odd Xgamma-G (OXGa-G)2018Maiti and Pramanik [21]
15Odd power-Lindley-G (OPLind-G)2019Hassan and Nassr [22]; Korkmaz et al. [23]
16Odd Lomax-G (OLx-G)2019Cordeiro et al. [24]
17Odd Hyperbolic Cosine-G (OHC-G)2019Kharazmi et al. [25]
18Odd flexible Weibull-H (OFW-H)2019El-Morshedy and Eliwa [26]
19Odd inverse Pareto-G (OIPa-G)2019Aldahlan et al. [27]
20Odd Nadarajah-Haghighi-G (ONH-G)2019Nascimento et al. [28]
21Odd Chen-G (OChen-G)2020El-Morshedy et al. [29]; Anzagra et al. [30]
22Odd DAL-G (ODAL-G)2020Ahmad et al. [31]
23Odd Stacy’s Gamma-G (OStGa-G)2020Nasir et al. [32]
24Odd Maxwell-G (OMax-G)2020Ishaq and Abiodun [33]
Table 2. Biases, MSEs, and CPs for α = 0.2 , β = 0.5 , c = 0.5 , and k = 1 .
Table 2. Biases, MSEs, and CPs for α = 0.2 , β = 0.5 , c = 0.5 , and k = 1 .
n = 50 n = 100
BiasMSECPBiasMSECP
α 0.2460.1340.980.2270.1040.95
β 0.0780.2450.930.2220.1060.97
c0.0720.0400.990.0540.0210.99
k−0.4130.3400.73−0.4050.2960.76
n = 200 n = 500
α 0.1440.0530.980.0750.0160.99
β 0.0070.0080.960.0040.0030.98
c0.0260.0060.990.0080.0010.98
k−0.3290.1940.79−0.2160.0970.86
Table 3. Biases, MSEs, and CPs for α = 0.2 , β = 0.6 , c = 1.8 , and k = 1.5 .
Table 3. Biases, MSEs, and CPs for α = 0.2 , β = 0.6 , c = 1.8 , and k = 1.5 .
n = 50 n = 100
BiasMSECPBiasMSECP
α 0.2310.1250.980.2190.1170.98
β 0.1580.5910.890.0330.2010.92
c0.2430.6271.000.2030.3190.99
k−0.5630.6330.80−0.5550.5710.80
n = 200 n = 500
α 0.1860.0950.980.1250.0510.98
β 0.0030.0770.94−0.0070.0050.95
c0.1180.1380.990.0560.0360.98
k−0.4750.4640.82−0.3720.3140.82
Table 4. Biases, MSEs, and CPs for α = 0.4 , β = 1.5 , c = 0.9 , and k = 1 .
Table 4. Biases, MSEs, and CPs for α = 0.4 , β = 1.5 , c = 0.9 , and k = 1 .
n = 50 n = 100
BiasMSECPBiasMSECP
α 0.1280.8270.780.0950.3830.80
β 0.8533.9350.890.7583.2190.94
c0.1510.4720.950.0330.2550.93
k0.2290.2600.920.2220.1990.95
n = 200 n = 500
α 0.0300.1370.82−0.0070.0500.84
β 0.5241.7450.960.3491.2030.91
c0.0220.1630.910.0170.1070.93
k0.2180.1520.970.1570.0700.98
Table 5. Some statistics and p-value for the fitted models to data set 1.
Table 5. Some statistics and p-value for the fitted models to data set 1.
KS
Distribution ^ AICBICHQICADCvMKSp-Value
NODALW976.16181960.32401973.07301965.49400.16400.01960.03110.9951
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
McW976.24091962.48201978.41901968.94400.16280.02020.03210.9928
GaW979.84451965.68901975.25101969.56600.76330.12190.05040.753
W981.14771966.29501972.67001968.88000.97550.15770.05670.6123
Table 6. MLEs and their SEs (in parentheses) for the fitted models to data set 1.
Table 6. MLEs and their SEs (in parentheses) for the fitted models to data set 1.
Distribution α β ck θ
NODALW0.14542.90360.47620.0914
(0.0580)(0.6640)(0.0970)(0.0151)
KwW0.443700.640406.987800.13710
(0.00330)(0.00260)(0.06740)(0.01040)
BW0.366200.656603.869600.14360
(0.00470)(0.00630)(0.74020)(0.01240)
EGW0.012800.708901.365901.60630
(0.01030)(0.11060)(0.83440)(0.48390)
McW0.045100.916300.183800.250409.685
(0.03280)(0.23770)(0.1420)(0.12550)(4.51590)
GaW0.017400.785401.28540
(0.00820)(0.12730)(0.35230)
W0.011800.90570
(0.0010)(0.05120)
Table 7. Some statistics and p-value for the fitted models to data set 2.
Table 7. Some statistics and p-value for the fitted models to data set 2.
KS
Distribution ^ AICBICHQICADCvMKSp-Value
NODALW215.1079438.2157445.7005441.04420.20230.02860.08270.8981
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
McW215.7566441.5132450.8692445.04890.26990.03740.08370.8895
GaW219.4700444.9401450.5537447.06150.82780.12500.11760.5203
W225.7065455.4131459.1555456.82731.72860.27650.13990.3048
Table 8. MLEs and their SEs (in parentheses) for the fitted models to data set 2.
Table 8. MLEs and their SEs (in parentheses) for the fitted models to data set 2.
Distribution α β ck θ
NODALW0.068616.03680.49690.0568
(0.0379)(9.1204)(0.2109)(0.0363)
KwW0.160901.0252054.782500.20410
(0.01530)(0.02760)(0.13580)(0.03820)
BW0.132001.1080023.060200.19400
(0.00730)(0.00680)(8.79410)(0.03240)
EGW0.009000.777405.5966010.54930
(0.00410)(0.13700)(2.04580)(5.68210)
McW0.160801.0049014.507800.221002.5180
(0.03400)(0.04660)(9.82270)(0.07570)(0.08950)
GaW4.614400.4983014.72250
(0.15180)(0.02170)(1.72390)
W0.017101.77190
(0.00150)(0.17760)
Table 9. The competitive models.
Table 9. The competitive models.
No.DistributionYearAuthor(s)
1Geometric (Geo)--
2Generalized geometric (GGeo)2010Gómez-Déniz [48]
3Discrete Rayleigh (DR)2004Roy [46]
4Discrete inverse Rayleigh (DIR)2014Hussain and Ahmad [51]
5Discrete inverse Weibull (DIW)2010Jazi et al. [49]
6One parameter discrete Lindley (DLi-I)2011Gómez-Déniz and Calderín-Ojeda [50]
7Two parameters discrete Lindley (DLi-II)2016Hussain et al. [52]
8Three parameters discrete Lindley (DLi-III)2020Eliwa at al. [56]
9Negative Binomial (NeBi)--
10Poisson (Poi)1837Poisson [70]
11Discrete Pareto (DPa)2009Krishna and Pundir [47]
12Discrete Burr type XII (DB-XII)2016Para and Jan [54]
13Discrete log-logistic (DLogL)2016Para and Jan [53]
14One parameter discrete flexible (DFx-I)2020Eliwa and El-Morshedy [57]
15Discrete Lomax distribution (DLo)2016Para and Jan [54]
Table 10. The MLEs and their SEs for data set 3.
Table 10. The MLEs and their SEs for data set 3.
Parameter ⟶p α β
ModelMLESEMLESEMLESE
DNODALGeo 0.711 0.130 1.587 0.392 0.240 0.034
Geo 0.582 0.030
DR 0.901 0.009
DIR 0.554 0.049
DLi-I 0.436 0.026
DFx-I 0.623 0.031
Poi 1.390 0.112
DPa 0.268 0.034
GGeo 0.188 0.089 0.800 0.064
DIW 1.049 0.146 0.581 0.048
DLo 0.152 0.098 1.830 0.952
DBX-II 0.278 0.045 1.053 0.167
NeBi 0.812 0.045 0.322 0.074
DLogL 0.780 0.136 1.208 0.159
DLi-II 0.581 0.045 0.001 0.058
DLi-III 0.582 0.005 358.728 11863.37 0.001 20.698
Table 11. The GoFS for data set 3 “part I”.
Table 11. The GoFS for data set 3 “part I”.
X ObFrExFr
DNODALGeoGeoDRDIRDLi-IDFx-IPoiDPa
065 64.97 45.98 10.89 60.89 40.29 45.26 27.39 65.84
114 14.69 26.76 26.62 33.99 29.83 29.09 38.08 18.27
210 8.46 15.58 29.45 8.12 18.36 16.51 26.47 8.16
36 5.73 9.06 22.29 3.00 10.34 8.89 12.26 4.51
44 4.19 5.28 12.63 1.42 5.52 4.70 4.26 2.82
52 3.18 3.07 5.54 0.78 2.85 2.49 1.19 1.91
62 2.44 1.79 1.91 0.47 1.44 1.34 0.27 1.37
72 1.87 1.04 0.53 0.31 0.71 0.73 0.05 1.02
81 1.41 0.61 0.12 0.21 0.35 0.41 0.01 0.79
91 1.04 0.35 0.02 0.15 0.17 0.23 0.01 0.63
101 0.74 0.21 0.00 0.11 0.08 0.14 0.01 0.51
112 1.28 0.27 0.00 0.55 0.06 0.21 0.00 4.17
Total110110110110110110110110110
l 166.48 178.77 277.78 186.55 189.11 182.29 246.21 171.19
AIC 338.96 359.537 557.56 375.09 380.22 366.58 494.42 344.38
HQIC 342.25 360.63 558.65 376.19 381.32 367.67 495.52 345.48
CAIC 339.19 359.577 557.59 375.13 380.26 366.61 494.46 344.42
χ 2 0.59 19.11 306.52 40.46 34.64 31.70 89.28 3.43
df 24424434
p-value 0.750 <0.001<0.001<0.001<0.001<0.001<0.001 0.489
Table 12. The GoFS for data set 3 “part II”.
Table 12. The GoFS for data set 3 “part II”.
X ObFrExFr
DNODALGeoGGeoDIWDLoDBX-IINeBiDLogLDLi-IIDLi-III
065 64.97 62.74 63.91 61.62 64.74 56.52 63.19 46.03 46.01
114 14.69 19.67 20.69 21.02 19.18 15.89 20.10 26.77 26.77
210 8.46 9.44 8.05 9.69 8.48 9.17 8.64 15.57 15.57
36 5.73 5.44 4.23 5.28 4.63 6.20 4.66 9.06 9.06
44 4.19 3.46 2.59 3.19 2.86 4.50 2.86 5.27 5.27
52 3.18 2.35 1.75 2.09 1.92 3.40 1.92 3.06 3.07
62 2.44 1.66 1.26 1.44 1.37 2.64 1.37 1.78 1.78
72 1.87 1.21 0.95 1.04 1.01 2.08 1.02 1.04 1.08
81 1.41 0.90 0.74 0.77 0.78 1.66 0.79 0.60 0.60
91 1.04 0.69 0.59 0.59 0.61 1.34 0.62 0.35 0.35
101 0.74 0.53 0.49 0.46 0.49 1.09 0.50 0.20 0.20
112 1.28 1.91 4.75 2.81 3.93 5.51 4.33 0.27 0.24
Total110110110110110110110110110110
l 166.48 168.56 172.94 170.48 171.14 168.54 171.72 178.77 178.77
AIC 338.96 341.11 349.87 344.96 346.28 340.09 347.43 361.53 363.53
HQIC 342.25 343.30 352.06 347.15 348.47 343.28 349.62 363.72 366.82
CAIC 339.19 341.23 349.98 345.07 346.39 344.49 347.55 361.65 363.76
χ 2 0.59 2.44 6.45 3.24 2.59 4.29 4.033 19.09 19.09
df 233324332
p-value 0.750 0.485 0.092 0.356 0.274 0.369 0.258 0.0003 <0.0001
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Eliwa, M.S.; Tahir, M.H.; Hussain, M.A.; Almohaimeed, B.; Al-Bossly, A.; El-Morshedy, M. Univariate Probability-G Classes for Scattered Samples under Different Forms of Hazard: Continuous and Discrete Version with Their Inferences Tests. Mathematics 2023, 11, 2929. https://doi.org/10.3390/math11132929

AMA Style

Eliwa MS, Tahir MH, Hussain MA, Almohaimeed B, Al-Bossly A, El-Morshedy M. Univariate Probability-G Classes for Scattered Samples under Different Forms of Hazard: Continuous and Discrete Version with Their Inferences Tests. Mathematics. 2023; 11(13):2929. https://doi.org/10.3390/math11132929

Chicago/Turabian Style

Eliwa, Mohamed S., Muhammad H. Tahir, Muhammad A. Hussain, Bader Almohaimeed, Afrah Al-Bossly, and Mahmoud El-Morshedy. 2023. "Univariate Probability-G Classes for Scattered Samples under Different Forms of Hazard: Continuous and Discrete Version with Their Inferences Tests" Mathematics 11, no. 13: 2929. https://doi.org/10.3390/math11132929

APA Style

Eliwa, M. S., Tahir, M. H., Hussain, M. A., Almohaimeed, B., Al-Bossly, A., & El-Morshedy, M. (2023). Univariate Probability-G Classes for Scattered Samples under Different Forms of Hazard: Continuous and Discrete Version with Their Inferences Tests. Mathematics, 11(13), 2929. https://doi.org/10.3390/math11132929

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