Abstract

INTRODUCTION

In many applied sciences such as medicine, engineering and finance, among others, modeling and analyzing lifetime data are crucial. Several lifetime models have been adopted to model different types of survival data. The quality of the procedures used in a statistical analysis depends heavily on the generated family of distributions and considerable effort has been directed to define new statistical models. However, there still remain many important problems involving real data, which do not follow any of the popular statistical models. So, the procedure of expanding a family of distributions by adding new shape parameters is well-known in the statistical literature.

The role of the extra shape parameters is to introduce skewness and to vary tail weights. Further, several models have been constructed by extending some useful lifetime distributions and investigated them with respect to different characteristics. The generalized distributions provide more flexibility to model both monotonic and non-monotonic failure rates even though the baseline failure rate may be monotonic.

A large number of models has been proposed to model lifetime data. Lindley (1958)LINDLEY DV. 1958. Fiducial distributions and Bayes’ theorem. J R Stat Soc Series B 20: 102-107. proposed the Lindley (Li) distribution as a mixture of exponential and gamma distributions to analyze failure time data. The Li distribution is motivated by its ability to model failure time data with increasing, decreasing, unimodal and bathtub shaped hazard rates.

The properties, estimation and inference of the Li distribution are described in the literature as follows. Hussain (2006)HUSSAIN E. 2006. The non-linear functions of order statistics and their properties in selected probability models. Ph.D. thesis. Department of Statistics. University of Karachi, Pakistan. used this distribution for studying stress–strength reliability modeling. Ghitany et al. (2008)GHITANY M, ATIEH B and NADARAJAH S. 2008. Lindley distribution and its application. Math Comp Simul 78: 493-506. provided a comprehensive treatment of its mathematical properties showing that it may provide a better fitting than the exponential distribution.

Many authors developed generalizations of the Li distribution. For example, Sankaran (1970)SANKARAN M. 1970. The discrete Poisson–Lindley distribution. Biometrics 26: 145-149. introduced the discrete Poisson–Li, Zakerzadeh and Dolati (2009)ZAKERZADEH H and DOLATI A. 2009. Generalized Lindley Distribution. J Math Ext 3: 13-25. defined the generalized Li, Zamani and Ismail (2010)ZAMANI H and ISMAIL N. 2010. Negative binomial-Lindley distribution and its application. Journal of Mathematics and Statistics 6: 4-9. proposed the negative binomial-Li, Mahmoudi and Zakerzadeh (2010)MAHMOUDI E and ZAKERZADEH H. 2010. Generalized Poisson–Lindley distribution. Commun Stat - Theory Methods 39: 1785-1798. investigated the generalized Poisson–Li, Nadarajah et al. (2011)NADARAJAH S, BAKOUCH HS and TAHMASBI R. 2011. A generalized Lindley distribution. Sankhya B 73: 331-359. introduced a two-parameter generalized Li as an alternative to the gamma, lognormal, Weibull and exponentiated exponential distributions, Bakouch et al. (2012)BAKOUCH HS, AL-ZAHRANI BM, AL-SHOMRANI AA, MARCHI VAA and LOUZADA F. 2012. An extended Lindley distribution. J Korean Stat Soc 41: 75-85. proposed the extended Li motivated by its ability to model lifetime data with different shapes for the failure rate, Ghitany et al. (2013)GHITANY ME, AL-MUTAIRI DK, BALAKRISHNAN N and AL-ENEZI L. 2013. Power Lindley distribution and associated inference. Comput Stat Data Anal 64: 20-33. pioneered the power Li and Nedjar and Zeghdoudi (2016)NEDJAR S and ZEGHDOUDI H. 2016. Gamma Lindley distribution and its application. J Appl Probab Stat 11: 129-138. studied the gamma Li distribution.

In this paper, based on the T-X family pioneered by Alzaatreh et al. (2013)ALZAATREH A, LEE C and FAMOYE F. 2013. A new method for generating families of distributions. Metron 71: 63-79. and the Li distribution, we construct the generalized odd Lindley-G (GOLi-G for short) family, which extends the OLi-G class (Gomes-Silva et al. 2017)GOMES-SILVA FS, PERCONTINI A, BRITO E, RAMOS MW, VENÂNCIO R and CORDEIRO GM. 2017. The Odd Lindley-G Family of Distributions. Austrian Journal of Statistics 46: 65-87., and provide some of its mathematical properties. In fact, the new generator of distributions is motivated by its ability to model lifetime data with increasing, decreasing, constant, unimodal and bathtub shaped failure rates. Further, the special models of this family are shown to provide better fits than other competitive models generated by other well-known families in the literature.

The cumulative distribution function (cdf) and probability density function (pdf) of the GOLi-G family with two additional shape parameters $\alpha >0$ and $\lambda >0$ are defined from a baseline cdf $G(x;𝛗)$ by

and

respectively, where $g(x;𝛗)=dG(x;𝛗)/dx$ and $𝛗$ is the baseline vector. Henceforth, a random variable with density (2) is denoted by $X\sim$ GOLi-G( $\alpha ,\lambda ,𝛗$ ). For $\alpha =1$ , we obtain as a special case the OLi-G class.

An interpretation of the GOLi-G family (1) can be given as follows. Let $Y$ be a random variable describing a stochastic system by the cdf $G(x{)}^{\alpha }$ (for $\alpha >0$ ). If the random variable $X$ represents the odds ratio, the risk that the system following the lifetime $Y$ will be not working at time $x$ is given by $G(x{)}^{\alpha }/[1-G(x{)}^{\alpha }]$ . If we are interested in modeling the randomness of the odds ratio by the Li pdf, $r(t)={\lambda }^2(1+t)e^{-\lambda t}/(1+\lambda )$ (for $t>0$ ), then the cdf of $X$ is given by

which is exactly the cdf (1) of the GOLi-G family.

Henceforth, we can omit the dependence on the model parameters. The hazard rate function (hrf) of $X$ is given by

The rest of the paper is organized as follows. In Section 2, we present four special models and plots of their pdfs and hazard rate functions (hrfs). In Section 3, we provide a very useful linear representation for the density function of $X$ . In Section 4, we derive some of its general mathematical properties including quantile and generating functions, asymptotics, ordinary and incomplete moments, order statistics and entropies. Maximum likelihood estimation of the model parameters is addressed in Section 5. Section 6 is devoted to simulation results to assess the performance of the maximum likelihood estimation of the unknown parameters of the generalized odd Lindley Weibull (GOLiW) distribution. In Section 7, we provide two applications to real data to illustrate the flexibility of the special models of the new family. Finally, we offer some concluding remarks in Section 8.

FOUR SPECIAL GOLi-G MODELS

In this section, we provide four special models of the GOLi-G family. The pdf (2) will be most tractable when the cdf $G(x)$ and the pdf $g(x)$ have simple analytic expressions.

THE GOLiW DISTRIBUTION

Consider the cdf (for $x>0$ ) $G(x)=1-exp(-a{x}^b)$ of the Weibull distribution with positive parameters $a$ and $b$ . Then, the pdf of the GOLiW model is given by

For $b=1$ , we obtain the GOLi-exponential (GOLiE) distribution. For $b=2$ , we obtain the GOLi-Rayleigh (GOLiR) distribution. The GOLiW reduces to the OLiW when $\alpha =1$ . We have the OLiE and OLiR when $\alpha =b=1$ and $\alpha =b=1$ , respectively. Some possible shapes of the GOLiW density and hazard functions are displayed in Figure 1.

Figure 1
Plots of the GOLiW pdf (a) and hrf (b) for some parameter values.

THE GOLi-BURR XII (GOLiBXII) DISTRIBUTION

The cdf (for $x>0$ ) of the Burr XII (BXII) distribution with positive parameters $a$ and $b$ is $G(x)=1-{(1+x^a)}^{-b}$ . Then, the GOLiBXII density reduces to

For $a=1$ and $b=1$ , we obtain the one-parameter GOLi Lomax (GOLiLx) and one-parameter GOLi log-logistic (GOLiLL) distributions, respectively. The cases $\alpha =a=1$ and $\alpha =b=1$ refer to the OLiLx and OLiLL distributions, respectively. Plots of the density and hazard functions of the GOLiBXII distribution are displayed in Figure 2.

Figure 2
Plots of the GOLiBXII pdf (a) and hrf (b) for some parameter values.

The GOLi-Lomax (GOLiLx) distribution

The cdf (for $x>0$ ) of the Lomax (Lx) distribution with positive parameters $a$ and $b$ is $G(x)=1-[1+(x/a){]}^{-b}$ . Then, the GOLiLx pdf becomes

Plots of the density and hazard functions of the GOLiLx distribution are displayed in Figure Figure 3.

Figure 3
Plots of the GOLiLx pdf (a) and hrf (b) for some parameter values.

THE GOLi-LOG-LOGISTIC (GOLiLL) DISTRIBUTION

Consider the cdf (for $x>0$ ) $G(x)=1-{[1+{\left(\frac{x}{a}\right)}^b]}^{-1}$ of the log-logistic (LL) distribution with positive parameters $a$ and $b$ . The GOLiLL density is given by

Plots of the density and hazard functions for the GOLiLL distribution for selected parameter values are displayed in Figure 4.

Figure 4
Plots of the GOLiLL pdf (a) and hrf (b) for some parameter values.

LINEAR REPRESENTATION

In this section, we provide a useful linear representation for the GOLi-G density. It can be expressed as

Using the exponential series, we can write

Consider the power series

Applying this series to the quantity $A$ , we obtain

Substituting (6) in equation (4), we can write

Then, the GOLi-G density can be expressed as a linear combination of exponentiated-G (Exp-G) densities

where

and $h_{\delta }\left(x\right)=\delta g\left(x\right)G{\left(x\right)}^{\delta -1}$ is the Exp-G density with power parameter $\delta >0$ . Thus, several mathematical properties of the GOLi-G family can be determined from those properties of the Exp-G family. Equation (8) is the main result of this section.

The cdf of the GOLi-G family can also be expressed as a linear combination of Exp-G cdfs. By integrating (8), we have

where $H_{\delta }\left(x\right)$ is the cdf of the Exp-G family with power parameter $\delta$ .

STRUCTURAL PROPERTIES

In this section, we obtain some structural properties of the GOLi-G family including quantile and generating functions, asymptotics, ordinary and incomplete moments, order statistics and entropies.

QUANTILE FUNCTION

Let $Q_G(\cdot )=G^{-1}(\cdot )$ be the quantile function (qf) of the parent G. In this section, we provide two algorithms for simulating the GOLi-G model.

The first algorithm is based on generating random data from the Li distribution using the exponential-gamma mixture.

• Algorithm 1 (Mixture form of the Li distribution)

1. Generate $U_i\sim Uniform(0,1),i=1,\dots ,n;$

2. Generate $V_i\sim Exponential(\lambda ),i=1,\dots ,n;$

3. Generate $W_i\sim Gamma(2,\lambda ),i=1,\dots ,n;$

If $U_i\le \frac{\lambda }{\lambda +1}$ , set $X_i=G^{-1}([\frac{V_i}{1+V_i}{]}^{1/\alpha })$ ;otherwise, set$X_i=G^{-1}([\frac{W_i}{1+W_i}{]}^{1/\alpha })$

The second algorithm is based on generating random data by inverting (1).

• Algorithm 2 (Inverse cdf)

1. Generate $U_i\sim Uniform(0,1),i=1,\dots ,n;$

Set (for $i=1,\dots ,n$ )

where $W(x)$ is the negative branch of the Lambert function. The branches of this function are defined by $x=W(x{e}^x),x\in ℂ$ . It is a two-valued function on the interval $[-1/e,0)$ . For $W(x)\le -1$ , the function is denoted $W_{-1}(x)$ and is called the negative branch. For $W(x)>-1$ , the function is called the principal branch of the W function. The Lambert function cannot be expressed in terms of elementary functions.

ASYMPTOTICS

Let $a=inf\{x|G(x)\}>0$ . Then, the asymptotics of equations (1), (2) and (3) when $x\to a$ are given by

The asymptotics of equations (1), (2) and (3) when $x\to \infty$ are given by

GENERATING FUNCTION

Here, we obtain the moment generating function (mgf) $M(t)=E(e^{tX})$ of $X$ . Henceforth, let $T_{(k+j+1)\alpha }$ denote the random variable having the Exp-G class with power parameter $(k+j+1)\alpha$ .

We can write from equation (8)

where ${\tau }_{k+j}(t)=(k+j+1)\alpha {\int }_0^{1}exp\left[t{Q}_G\left(u\right)\right]u^{(k+j+1)\alpha -1}du$ and $M_{(k+j+1)\alpha }(t)$ is the mgf of $T_{(k+j+1)\alpha }$ . Hence, $M(t)$ can be determined from the Exp-G generating function. If $M_{(k+j+1)\alpha }(t)$ has an explicit expression, $M(t)$ will also have a closed-form. For some baselines such as the exponential, Lomax, normal, gamma, among others, this is possible.

We define the cumulant generating function (cgf) of $X$ by $K(t)=log[M(t)]$ . The saddle–point approximation is of the main applications of the cgf in Statistics and provides highly accurate approximation formulae for the densities of the sum and mean of independent identically distributed (iid) random variables. Let $X_1,\cdots ,X_n$ be iid random variables having common GOLi-G cgf $K(t)$ . We define $S_n={\sum }_{j=1}^n{X}_j$ and obtain $\hat{\lambda }$ from the (usual nonlinear) equation $K^{\prime }(\hat{\lambda })=x/n$ and let $y=[x-n{K}^{\prime }(\lambda )]/\sqrt{n{K}^{\prime \prime }(\lambda )}$ . The density function of ${\overline{X}}_n=S_n/n$ is followed from Daniels’ saddle–point approximation as

MOMENTS

The $r$ th ordinary moment of $X$ , say ${\mu }_r^{\prime }$ , follows from (8) as

The central moment of $X$ are easily obtained from (10).

The cumulants ( ${\chi }_n$ ) of $X$ follow recursively from

where ${\chi }_1={\mu }_1^{\mathrm{\prime }}{\chi }_2={\mu }_2^{\mathrm{\prime }}-{\mu }_1^{\mathrm{\prime }2},{\chi }_3={\mu }_3^{\mathrm{\prime }}-3{\mu }_2^{\mathrm{\prime }}{\mu }_1^{\mathrm{\prime }}+{\mu }_1^{\mathrm{\prime }3}$ , etc. The measures of skewness and kurtosis can be calculated from the ordinary moments using well-known relationships.

The $s$ th incomplete moment of $X$ , say ${\phi }_s(t)={\int }_{-\infty }^t{x}^{s}f(x)dx$ , can be expressed from (8) as

where $v_{(k+j+1)\alpha }^{(s)}\left(t\right)=(k+j+1)\alpha {\int }_0^{G\left(t\right)}{Q}_G{\left(u\right)}^s{u}^{(k+j+1)\alpha -1}du$ can be computed numerically from the baseline qf $Q_G(u)$ .

For $s=1$ , we obtain the first incomplete moment ${\phi }_1(t)$ from (11). It can be applied to construct Bonferroni and Lorenz curves defined for a given probability $\pi$ by $B(\pi )={\phi }_1\left(q\right)/(\pi {\mu }_1^{\prime })$ and $L(\pi )={\phi }_1\left(q\right)/{\mu }_1^{\prime }$ , respectively, where ${\mu }_1^{\prime }$ is given by (10) with $r=1$ and $q=Q(\pi )$ is the qf of $X$ at $\pi$ obtained from (9). These curves are very useful in economics, reliability, demography, insurance and medicine.

ORDER STATISTICS

Order statistics make their appearance in many areas of statistical theory and practice. Let $X_1,\cdots ,X_n$ be a random sample from the GOLi-G family. The pdf of $X_{i:n}$ can be written as

where $B(\cdot ,\cdot )$ is the beta function.

Based on equation (1), we can write

After applying the generalized binomial series, we have

Using (2) and (13), we can write

Using the exponential series, we obtain

The Taylor series for $z^{\beta }$ defined by

holds, where $(\beta {)}_k=\beta (\beta -1)\dots (\beta -k+1)$ is the descending factorial.

Applying (14) to ${[1+\frac{\lambda G(x{)}^{\alpha }}{1-G(x{)}^{\alpha }}]}^l$ gives

Then, the last equation reduces to

Applying the power series (5) to ${[1-G(x{)}^{\alpha }]}^{-(s+h+3)}$ gives

Inserting the last equation in (12), the pdf of $X_{i:n}$ can be expressed as

where, as before, $h_{\delta }\left(x\right)$ is the Exp-G density with power parameter $\delta$ and

Then, the density function of the GOLi-G order statistics is a linear combination of Exp-G densities. Based on equation (15), we note that the mathematical properties of $X_{i:n}$ follow from those of $T_{(k+s+h+1)\alpha }$ . For example, the $q$ th moment of $X_{i:n}$ is given by

Based upon the moments in equation (16), we can derive explicit expressions for the L-moments of $X$ as infinite weighted linear combinations of the means of suitable GOLi-G order statistics.

ENTROPIES

The Rényi entropy of a random variable $X$ represents a measure of variation of the uncertainty. The Rényi entropy is given by

Using the pdf (2), we can write

Using the exponential series, we have

Applying the power series (5) to the last term, we have

Then, the Rényi entropy of the GOLi-G family reduces to

where

The Shannon entropy of $X$ , say $SI=E\{-[logf(X)]\}$ , is given by

First, we define and compute

By using the exponential series, we obtain

Using the binomial expansion, we have

The following proposition is used to determine the Shannon entropy of $X$ .Proposition 1. Let $X$ be a random variable with pdf given in (2). Then,

MAXIMUM LIKELIHOOD ESTIMATION

Several approaches for parameter estimation were proposed in the literature but the maximum likelihood method is the most commonly employed. The maximum likelihood estimators (MLEs) enjoy desirable properties and can be used to obtain confidence intervals for the model parameters. The normal approximation for these estimators in large samples can be easily handled either analytically or numerically. Here, we consider the estimation of the unknown parameters of the new family from complete samples only by maximum likelihood.

Let $X_1,\cdots ,X_n$ be a random sample from the GOLi-G family with parameters α,λ and $𝛗$ . Let $𝛉=$ ( $\alpha ,\lambda ,{𝛗}^{⊺}$ ) ${}^{⊺}$ be the $p\times 1$ parameter vector. The log-likelihood function for $𝛉$ is given by

Then, the score vector components, $𝐔\left(𝛉\right)=\frac{\partial \ell }{\partial 𝛉}={(U_{\alpha },U_{\lambda },U_{{\phi }_k})}^{⊺}$ , are

and

where $g_k^{\prime }(x_i;𝛗)=\partial g(x_i;𝛗)/\partial {\phi }_k$ and $G_k^{\prime }(x_i;𝛗)=\partial G(x_i;𝛗)/\partial {\phi }_k.$

Setting the nonlinear system of equations $U_{\alpha }=U_{\lambda }=0$ and $U_{𝛗}=0$ and solving them simultaneously yields the MLE $\widehat{𝛉}=(\hat{\alpha },\hat{\lambda },{\widehat{𝛗}}^{⊺}{)}^{⊺}$ . For doing this, it is usually more convenient to adopt nonlinear optimization methods such as the quasi-Newton algorithm to maximize $\ell$ numerically. For interval estimation of the parameters, we obtain the $p\times p$ observed information matrix $J(𝛉)=\{-\frac{{\partial }^2\ell }{\partial r\partial s}\}$ (for $r,s=\alpha ,\lambda ,{\phi }_k$ ), which can be evaluated numerically.

Under standard regularity conditions, the distribution of $\widehat{𝛉}$ can be approximated by a multivariate normal $N_p(0,J(\widehat{𝛉}{)}^{-1})$ distribution when $n\to \infty$ to obtain confidence intervals for the parameters. Here, $J(\widehat{𝛉})$ is the total observed information matrix evaluated at $\widehat{𝛉}$ . The method of the re-sampling bootstrap can be used for correcting the biases of the MLEs of the model parameters. Good interval estimates may also be obtained using the bootstrap percentile method.

SIMULATION STUDY

In this section, we present some simulation results that investigate the behavior of the MLEs in terms of the sample size $n$ . All simulations are performed using the Rprogramming language (R Core Team 2017R CORE TEAM. 2017. R: A language and environment for statistical computing. R Foundation for Statistical Computing.).

The qf of the GOLiW distribution is given by

We generate $2,000$ random samples from this distribution using the last expression for three different sample sizes $n=350$ , $n=550$ and $n=750$ . The true values of the parameters are taken as: $\alpha =(0.2,1,2),\lambda =(0.5,2),a=(0.7,2.5,4.3)$ and $b=(0.9,3,5)$ . For each sample size and parameter combination, the average MLEs and the mean square errors (MSEs) are computed. In order to save space, Table I gives only the results for $\alpha =2$ , and are not reported for $\alpha =(0.2,1)$ . It can be verified that the estimates are stable and quite close the true parameter values for all sample sizes. Further, the MSEs decrease when the sample size increases in all cases in agreement with the first-order asymptotic theory.

APPLICATIONS

In this section, we provide two applications to real data to illustrate the flexibility of the GOLiLx, GOLiLL and GOLiBXII distributions presented in Section 2. The goodness-of-fit statistics for these models are compared with other competitive models and the MLEs of the model parameters are determined.

The model selection is carried out using some goodness-of-fit measures including the Akaike information criterion ( $AIC$ ), consistent Akaike information criterion ( $CAIC$ ), Hannan-Quinn information criterion ( $HQIC$ ), Bayesian information criterion ( $BIC$ ), Anderson-Darling $(A^{*})$ and Cramér-von Mises ( $W^{*}$ ). The smaller these statistics are, the better the fit.

APPLICATION 1: TIME-TO-FAILURE OF TURBOCHARGER DATA

The first data set consists of 40 times to failures ( ${10}^3$ h) of turbocharger of one type of engine given in Xu et al. (2003)XU K, XIE M, TANG LC and HO SL. 2003. Application of neural networks in forecasting engine systems reliability. Appl Soft Comput 2: 255-268.. For these data, we shall compare the fit of the GOLiLx and GOLiLL distributions with those of other competitive models, namely: the Kumaraswamy transmuted log-logistic (KwTLL) (Afify et al. 2016AFIFY AZ, CORDEIRO GM, YOUSOF HM, ALZAATREH A and NOFAL ZM. 2016. The Kumaraswamy transmuted-G family of distributions: properties and applications. Journal of Data Science 14: 245-270.), Kumaraswamy log-logistic (KwLL) (de Santana et al. 2012DE SANTANA TVF, ORTEGA EMM, CORDEIRO GM and SILVA GO. 2012. The Kumaraswamy log-logistic distribution. Stat Theory Appl 3: 265-291.), transmuted log-logistic (TLL) (Granzotto and Louzada 2015GRANZOTTO DCT and LOUZADA F. 2015. The transmuted log-logistic distribution: modeling, inference and na application to a polled tabapua race time up to first calving data. Commun Stat - Theory Methods 44: 3387-3402.), McDonald log-logistic (McLL) (Tahir et al. 2014TAHIR MH, MANSOOR M, ZUBAIR M and HAMEDANI G. 2014. McDonald log-logistic distribution with an application to breast cancer data. Journal of Statistical Theory and Applications 13: 65-82.), beta log-logistic (BLL) (Lemonte 2014LEMONTE AJ. 2014. The beta log-logistic distribution. Braz J Probab Stat 28: 313-332.), generalized transmuted log-logistic (GTLL) (Nofal et al. 2017NOFAL ZM, AFIFY AZ, YOUSOF HM and CORDEIRO GM. 2017. The generalized transmuted-G family of distributions. Commun Stat - Theory Methods 46: 4119-4136.) and Li-log-logistic (LiLL) (Cakmakyapan and Ozel 2016CAKMAKYAPAN S and OZEL G. 2016. The Lindley family of distributions: properties and applications. Hacet J Math Stat 46: 1-27.) distributions, whose densities (for $x>0$ ) are given, respectively, by: KwTLL: $f\left(x\right)=ab\beta {\alpha }^{-\beta }{x}^{\beta -1}{[1+{\left(\frac{x}{\alpha }\right)}^{\beta }]}^{-2}{\{1-{[1+{\left(\frac{x}{\alpha }\right)}^{\beta }]}^{-1}\}}^{a-1}$

$\times \{1-\lambda +2\lambda {[1+{\left(\frac{x}{\alpha }\right)}^{\beta }]}^{-1}\}{\{1+\lambda {[1+{\left(\frac{x}{\alpha }\right)}^{\beta }]}^{-1}\}}^{a-1}$

$\times {\{1-{[1+\lambda {[1+{\left(\frac{x}{\alpha }\right)}^{\beta }]}^{-1}]}^a{\{1-{[1+{\left(\frac{x}{\alpha }\right)}^{\beta }]}^{-1}\}}^a\}}^{b-1};$

KwLL: $f(x)=\frac{ab\beta }{{\alpha }^{a\beta }}{x}^{a\beta -1}{[1+{\left(\frac{x}{\alpha }\right)}^{\beta }]}^{-a-1}{\{1-{[1-\frac{1}{1+{\left(\frac{x}{\alpha }\right)}^{\beta }}]}^a\}}^{b-1};$

TLL: $f(x)=\beta {\alpha }^{-\beta }{x}^{\beta -1}{[1+{\left(\frac{x}{\alpha }\right)}^{\beta }]}^{-2}\{1-\lambda +2\lambda {[1+{\left(\frac{x}{\alpha }\right)}^{\beta }]}^{-1}\};$

McLL: $f(x)=\frac{\beta c}{\alpha B(a/c,b)}{\left(\frac{x}{\alpha }\right)}^{\beta a-1}{[1+{\left(\frac{x}{\alpha }\right)}^{\beta }]}^{-a-1}$

$\times {\{1-{[1-{[1+{\left(\frac{x}{\alpha }\right)}^{\beta }]}^{-1}]}^c\}}^{b-1};$

BLL: $f(x)=\frac{\beta {\alpha }^{-a\beta }}{B(a,b)}{x}^{a\beta -1}{[1+{\left(\frac{x}{\alpha }\right)}^{\beta }]}^{-a-b};$

GTLL: $f(x)=\beta {\alpha }^{-\beta }{x}^{\beta -1}{[1+{\left(\frac{x}{\alpha }\right)}^{\beta }]}^{-2}{\{1-{[1+{\left(\frac{x}{\alpha }\right)}^{\beta }]}^{-1}\}}^{a-1}$

$\times \{a(1+\lambda )-\lambda (a+b){\{1-{[1+{\left(\frac{x}{\alpha }\right)}^{\beta }]}^{-1}\}}^b\};$

LiLL: $f(x)=\frac{\beta {\lambda }^2{\alpha }^{-\beta }}{(1+\lambda )}{x}^{\beta -1}{[1+{\left(\frac{x}{\alpha }\right)}^{\beta }]}^{-\lambda -1}(1-log\left\{{[1+{\left(\frac{x}{\alpha }\right)}^{\beta }]}^{-1}\right\}).$ The parameters of the above densities are all positive real numbers except for the KwTLL, TLL and GTLL models for which the parameter $\lambda$ is $\left|\lambda \right|\le 1$ .

APPLICATION 2: CANCER PATIENTS DATA

The second data set on the remission times (in months) of a random sample of 128 bladder cancer patients (Lee and Wang 2003LEE ET and WANG JW. 2003. Statistical methods for survival data analysis. 3rd ed., Wiley, New York.). For these data, we compare the fit of the GOLiBXII distribution with those of the OLiBXII, Weibull BXII (WBXII) (Afify et al. 2018AFIFY AZ, CORDEIRO GM, ORTEGA EMM, YOUSOF HM and BUTT NS. 2018. The four-parameter Burr XII distributions: properties, regression model and applications. Commun Stat - Theory Methods 47: 2605-2624.), Kumaraswamy exponentiated Burr XII (KwEBXII) (Mead and Afify 2017MEAD ME and AFIFY AZ. 2017. On five-parameter Burr XII distribution: properties and applications. South African Statistical Journal 51: 67-80.), McDonald Weibull (McW) (Cordeiro et al. 2014), beta BXII (BBXII) (Paranaíba et al. 2011PARANAÍBA PF, ORTEGA EMM, CORDEIRO GM and PESCIM RR. 2011. The beta Burr XII distribution with application to lifetime data. Comput Stat Data Anal 55: 1118-1136.), beta exponentiated BXII (BEBXII) (Mead 2014MEAD ME. 2014. The beta exponentiated Burr XII distribution. Journal of Statistics: Advances in Theory and Applications 12: 53-73.), LiBXII (Cakmakyapan and Ozel 2016) and BXII models with densities (for $x>0$ ), respectively, given by:

WBXII: $f(x)=\alpha \lambda ab{x}^{\alpha -1}{(1+x^{\alpha })}^{\lambda b-1}{[1-{(1+x^{\alpha })}^{-\lambda }]}^{b-1}$

$\times exp\{-a{[{(1+x^{\alpha })}^{\lambda }-1]}^b\};$

KwEBXII: $f(x)=\frac{abck\beta x^{c-1}}{{(1+x^c)}^{k+1}}{[1-{(1+x^c)}^{-k}]}^{a\beta -1}$