Journal of Probability and Statistics
Volume 2018 (2018), Article ID 8767826, 12 pages
https://doi.org/10.1155/2018/8767826

Research Article

The Half-Logistic Generalized Weibull Distribution

Masood Anwar and Amna Bibi

Department of Mathematics, COMSATS Institute of Information Technology, Park Road, Chak Shahzad, Islamabad, Pakistan

Correspondence should be addressed to Masood Anwar; masoodanwar@comsats.edu.pk

Received 8 August 2017; Revised 17 November 2017; Accepted 28 November 2017; Published 4 January 2018

Academic Editor: Ahmed Z. Afify

Copyright © 2018 Masood Anwar and Amna Bibi. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

A new three-parameter generalized distribution, namely, half-logistic generalized Weibull (HLGW) distribution, is proposed. The proposed distribution exhibits increasing, decreasing, bathtub-shaped, unimodal, and decreasing-increasing-decreasing hazard rates. The distribution is a compound distribution of type I half-logistic-G and Dimitrakopoulou distribution. The new model includes half-logistic Weibull distribution, half-logistic exponential distribution, and half-logistic Nadarajah-Haghighi distribution as submodels. Some distributional properties of the new model are investigated which include the density function shapes and the failure rate function, raw moments, moment generating function, order statistics, L-moments, and quantile function. The parameters involved in the model are estimated using the method of maximum likelihood estimation. The asymptotic distribution of the estimators is also investigated via Fisher’s information matrix. The likelihood ratio (LR) test is used to compare the HLGW distribution with its submodels. Some applications of the proposed distribution using real data sets are included to examine the usefulness of the distribution.

1. Introduction

Statistical distributions are the basic aspects of all parametric statistical techniques including inference, modeling, survival analysis, and reliability. For the analysis of lifetime data, it is an important task to fit the data by a statistical model. A number of lifetime distributions have been developed in literature for this purpose. The widely used lifetime models usually have a limited range of behaviors. Such type of distributions cannot give a better fit to model all the practical situations. Recently, several authors have developed a number of families of statistical models by applying different techniques. Various techniques have been introduced in the literature to derive new flexible models as discussed by Lai [1].

Marshall and Olkin [2] introduced an effective technique to extend a family of distributions by addition of another parameter. By applying this technique, they generalized the exponential and the Weibull distributions. Al-Zahrani and Sagor [3] proposed the Poisson Lomax model by compounding the Poisson Lomax distributions. Bidram and Nadarajah [4] introduced the exponentiated EG2 distribution by using the method of resilience parameter. Kus [5] considered compounding of Poisson and exponential distribution. Gurvich et al. [6] generalized the Weibull distribution offering a new distribution function elucidating a wide range of functional forms of the effect of size on the strength distribution, using a simple method of evaluation of the basic statistical parameters. Nadarajah and Kotz [7], Lai et al. [8], Lai et al. [9], and Xie et al. [10] further discussed some modifications of the Weibull model.

In this paper, another extension of the extended Weibull distribution is introduced using half-logistic-G generator. So we propose the half-logistic generalized Weibull (HLGW) distribution without adding any extra parameter to the baseline model. The new model is the compound distribution of two previously known distributions, one of which follows the class proposed by Gurvich et al. [6] and the other is type I half-logistic-G model. The proposed model is able to depict more complex hazard rates and provides a good alternate to the Weibull distribution that does not exhibit upside-down bathtub-shaped or unimodal failure rates.

Dimitrakopoulou et al. [11] established a three-parameter lifetime model with PDF

where and are the shape parameters and is a scale parameter. The CDF corresponding to (1) is

In this paper, a three-parameter lifetime model is presented which is the compound model of the previously known models introduced by Dimitrakopoulou et al. [11] and half-logistic-G (HL-G) distribution called half-logistic generalized Weibull (HLGW) distribution. The half-logistic-G distribution is presented by Cordeiro et al. [12] with the CDF

where is the CDF of the baseline distribution and is the shape parameter. As a special case, for , the TIHL-G is the half-logistic-G (HL-G) model with cumulative distribution function

The corresponding PDF to (4) is given by

where and .

The rest of the paper is unfolded as follows. Section 2 contains the introduction of the half-logistic generalized Weibull (HLGW) distribution and provides the plots of its density function. Section 3 explores the distributional properties of the HLGW model. In Section 4, the method of maximum likelihood estimation is used to obtain the estimators of unknown parameters. The asymptotic distribution of the estimators is also investigated in this section via Fisher’s information matrix. A simulation study is discussed in Section 5 to check out the accuracy of point and interval estimates of the HLGW parameters. Section 6 involves some applications of the HLGW distribution using lifetime data sets to examine the fitness of the proposed model. Section 7 provides concluding remarks about the paper.

2. The Half-Logistic Generalized Weibull Distribution

Substitution of (1) and (2) in (5) results the following PDF of the HLGW distribution:

The CDF associated with (6) is as follows:

The parameters are the shape parameters and is a scale parameter. From now on, a random variable with PDF (6) will be written as ~.

3. Distributional Properties

This section deals with the investigation of the distributional properties of HLGW distribution. The statistical properties include the plots of the density function, the failure rate function, raw moments, moment generating function, order statistics, L-moments, and quantile function.

3.1. Special Cases

The HLGW distribution includes the following distributions as special cases:

(i)For , the HLGW model reduces to the half-logistic Weibull (HLW) model with the PDF

 where is the shape parameter and is the scale parameter. For , the half-logistic Weibull model is also called the power half-logistic distribution (PHLD) proposed and studied by Krishnarani [13].(ii)For , the HLGW model generates a new model, the half-logistic exponential (HLE) model, with scale parameter and the PDF

(iii)For , the HLGW model gives another new model, the half-logistic Nadarajah-Haghighi (HLNH) model, with the PDF

 where and are the shape and scale parameters, respectively.

3.2. The Shapes of HLGW Distribution

The following are the conditions under which the PDF of the HLGW distribution (6) shows different behaviors:

(i)For , the PDF is monotone decreasing with , .(ii)For , the same shape is exhibited with and .(iii)For , the PDF has the value zero at the origin; then it increases to a higher value and then decreases, falling towards the value of zero at infinity.

Various behaviors of the PDF are shown in Figure 1, for some parameter values; . The mode of density (6) can be obtained from .

Figure 1: The shapes of the PDF of the HLGW distribution.

3.3. Hazard Rate Function (HRF)

For is a random variable, the HRF is given as , where represents the survival function given by

The HRF of ~ can be written as

for , demonstrating different shapes for different parameter values. By differentiating (12), it can be easily checked that

(a)for and , the value of HRF is zero at the origin; then it increases to its maximum; after that it is constant,(b)for and or and , has increasing or decreasing behavior,(c)for and ,

(i)if , is decreasing or decreasing-increasing-decreasing (DID),(ii)if , is increasing or bathtub-shaped,

(d)for and ,

(i)if , has upside-down bathtub shape (unimodal),(ii)if , is monotone increasing.

The different hazard shapes are shown in Figure 2 for different parameter values. By restricting , the failure rate function (12) can also be depicted as

for , where and . Thus the can be written as the sum of terms and hence (13) is the failure rate function of a series system of components.

Figure 2: The HRF of HLGW distribution for different parameter values.

3.4. Moments

In this section, the th moment of the HLGW model is presented as an infinite sum representation. The first four moments for have been calculated accordingly.

Theorem 1. Let ~ be a random variable, where , and the th moment of about the origin is as follows:

where and .

Corollary 2. Let ~ be a random variable, where and . The first four moments of the random variable are

3.5. Moment Generating Function (MGF)

The MGF of is retrieved from .

Theorem 3. let ~ be a random variable, and the moment generating function (MGF) of is given by

where .

3.6. Order Statistics

Let be a random sample of size from a distribution with PDF and CDF and are the analogous order statistics. The PDF and CDF of , , are

where is the beta function.

Theorem 4. Let and be the PDF and CDF of random variable ~; then the PDF of is

where .
The CDF corresponding to (18) is

3.7. L-Moments

Suppose that we have a random sample collected from ~. The th population L-moments are as follows:

We use the substitution , so and , where . Therefore, we get

By working out the integration, we arrive at the following formula:

where

The relation (22) can be used to find out the first L-moments of ; that is, for we get .

The other two moments, and , are, respectively, given by

3.8. Quantile Function

For to be a random variable with the PDF (6), the quantile function is

The above relation is used to find the quantile function of HLGW distribution. Therefore, we have

Hence, the generator for can be given by the following algorithm:

(1)Generate ~.(2)Use (27) and obtain an outcome of by .

By using the quantile function (27), we can examine the Bowley skewness [14] and Moors kurtosis [15] for HLGW as follows:

Table 1 illustrates the values of skewness and kurtosis for the HLGW model for some values of , , and . It can be noted that the skewness and kurtosis are free of parameter and they are decreasing functions of the parameters and .

Table 1: Skewness and kurtosis of HLGW for different values of , , and .

0.5
0.5
0.6530
2.4642
0.3366
1.4540
0.1303
1.2492

1.0
0.4634
1.6286
0.1808
1.2395
0.0197
1.2074

1.5
0.3913
1.4463
0.1264
1.1952
−0.0196
1.2045

2.0
0.3541
1.3700
0.0989
1.1779
−0.0396
1.2054

2.5
0.3314
1.3286
0.0824
1.1691
−0.0517
1.2067

1.5
0.5
0.6530
2.4642
0.3366
1.4540
0.1303
1.2492

1.0
0.4634
1.6286
0.1808
1.2395
0.0197
1.2074

1.5
0.3913
1.4463
0.1264
1.1952
−0.0196
1.2045

2.0
0.3541
1.3700
0.0989
1.1779
−0.0396
1.2054

2.5
0.3314
1.3286
0.0824
1.1691
−0.0517
1.2067

2.0
0.5
0.6530
2.4642
0.3366
1.4540
0.1303
1.2492

1.0
0.4634
1.6286
0.1808
1.2395
0.0197
1.2074

1.5
0.3913
1.446
0.1264
1.1952
−0.0196
1.2045

2.0
0.3541
1.3700
0.0989
1.1779
−0.0396
1.2054

2.5
0.3314
1.3286
0.0824
1.1691
−0.0517
1.2067

3.0
0.5
0.6530
2.4642
0.3366
1.4540
0.1303
1.2492

1.0
0.4634
1.6286
0.1808
1.2395
0.0197
1.2074

1.5
0.3913
1.4463
0.1264
1.1956
−0.0196
1.2045

2.0
0.3541
1.3700
0.0989
1.1779
−0.0396
1.2054

2.5
0.3314
1.3286
0.0824
1.1691
−0.0517
1.2067

4. Estimation

The log-likelihood function is expressed as

Taking the first partial derivatives of with respect to , , and and letting them equal zero, we obtain a nonlinear system of equations.

The above equations cannot be solved analytically, and statistical software can be used to solve them numerically via iterative methods and get the maximum likelihood estimate (MLE) of , , and .

4.1. Asymptotic Distribution

In order to have approximate confidence intervals (CIs) of the involved parameters, we require the estimated values of the elements of variance-covariance matrix of the MLEs. The variance-covariance matrix is estimated by the observed information matrix , where

where , are the second partial derivatives of (29) with respect to , , and . They are the entries of Fisher’s information matrix analogous to , , and , respectively, which are given in Appendix. The diagonal of matrix in (31) gives the variances of the MLEs of , and , respectively.

Approximation by a standard normal (SN) distribution of the distribution of , where , results in an approximate % confidence interval for as

where is the upper th percentile of SN distribution. We can use the likelihood ratio (LR) test to compare the fit of the HLGW distribution with its submodels for a given data set. For example, to test , the LR statistic is , where , and are the unrestricted estimates and are the restricted estimates. The LR test rejects the null hypothesis if , where denotes the upper point of the -distribution with 1 degree of freedom.

5. Simulation Study

The MLEs can be checked out by a simulation study. The following steps can be followed:

(1)By using (6), 5,000 samples of size are achieved. The variates of the HLGW distribution are developed using

 for ~.(2)The MLEs are computed for the samples, say for .(3)The mean square errors (MSEs) are calculated for every parameter.

The above steps were repeated for , and with and . Table 2 shows the bias and MSEs of , and . It can be deduced through the table that MSEs for individual parameters fall to zero when sample size increases.

Table 2: Bias and MSE for the HLGW parameters.

Bias ()
MSE ()
Bias ()
MSE ()
Bias ()
MSE ()

1.5
5
0.5
20
0.2103
0.0392
1.1935
3.2415
0.3454
0.9104

40
0.1732
0.0314
1.1431
2.9381
0.2931
0.5830

60
0.1526
0.0283
1.1063
2.4897
0.2480
0.3148

80
0.1473
0.0231
0.9733
2.1738
0.1812
0.1908

100
0.1332
0.0207
0.8910
1.9318
0.1317
0.0813

250
0.1171
0.0113
0.5337
1.3877
0.0811
0.0031

500
0.1010
0.0061
0.2701
1.0814
0.0213
0.0008

1.5
10
1
20
0.3118
0.0433
2.6754
15.9312
0.6125
10.2311

40
0.2918
0.0395
2.3487
12.3471
0.2451
7.6401

60
0.2554
0.0365
2.1174
9.2877
0.1615
4.3881

80
0.2375
0.0335
2.0968
6.2514
0.0941
1.9722

100
0.2114
0.0317
2.0532
4.9934
0.0532
1.5531

250
0.1783
0.0307
0.0065
1.9951
0.0113
0.4899

500
0.1265
0.0299
0.0013
0.2164
0.0095
0.0989

6. Applications

This section deals with the applications of the HLGW model to two lifetime data sets, that is, the data of 213 observed values of intervals between failures of air conditioning system of Boeing 720 jet airplanes and the data of waiting times (min) of 100 bank customers. Estimates of the parameters of HLGW distribution (standard errors in parentheses) and Cramer-von Mises , Anderson Darling , and K-S statistics are presented for the data sets. In general, the smaller the values of , , and K-S statistics, the better the fit. We compare the proposed model with other models for the same data sets to check the potentiality and flexibility of new model.

6.1. Air Conditioning Systems Failure Data

The first application concerns 213 observed values of intervals between failures of air conditioning system of Boeing 720 jet airplanes firstly analysed by Proschan [16]. We have compared the performance of the HLGW distribution with its submodels as well as with some other well-known models given below.

The Weibull Poisson (WP) distribution by Lu and Shi [17] with the PDF:

for , , and .

The complementary Weibull geometric (CWG) distribution by Tojeiro et al. [18] with the PDF:

for , and .

The Weibull (W) distribution with PDF:

for , , and .

The estimated values of the parameters with standard errors (SE) are found using the method of maximum likelihood estimation. Table 3 gives the estimated values of the parameters along with their standard errors and the test statistics have the smallest values of , , and K-S statistic for the data set under HLGW distribution as compared to the other models. Based on these values, it is concluded that HLGW distribution is the best model as compared to the other models to fit this data set. This conclusion can also be made by the CDF plots for empirical and fitted HLGW distributions using data of 213 values of intervals between failures of air conditioning system in Figure 3(a) for the data. In Figure 3(b), the TTT plot is shown for the data set. The failure rate shape for this data set is decreasing, as its TTT plot is convex.

Table 3: Estimates of models for the air conditioning systems failure data.

Distributions
Estimates

Statistics

K-S
value

HLGW
0.3470
1.3814
0.0296
2349.674
0.0339
0.2587
0.0394
0.8957

(0.0885)
(0.2462)
(0.0127)

HLW

0.7804
0.0442
2360.630
0.2099
1.2851
0.0577
0.4770

(0.0426)
(0.0098)

HLE

0.0144
2383.325
0.3446
2.0606
0.1407
0.0004

(0.0009)

HLNH
0.5596

0.0482
2352.930
0.0850
0.5534
0.0445
0.7926

(0.0515)

(0.0108)

Weibull
0.0158
0.9226

2355.171
0.1373
0.8552
0.0514
0.6270

(0.0038)
(0.0459)

CWG
0.2890
0.0079
4.1050
2362.329
0.1519
1.0253
0.0502
0.6563

(0.0238)
(0.0032)
(2.8847)

WP
0.4015
0.4967
−8.0978
2351.850
0.0440
0.3305
0.0412
0.8624

(0.0584)
(0.1838)
(2.5766)

(a) TTT plot

(b) CDF Plot

Figure 3: TTT plot and CDF plot for the air conditioning systems failure data.

The LR test statistics of hypotheses : HLW versus : HLGW, : HLE versus : HLGW, and : HLNH versus : HLGW are ( value = ), ( value = ), and ( value = ), respectively. We conclude that there is a significant difference between HLW and HLGW distributions, HLE and HLGW, and also between HLNH and HLGW distributions at the level.

6.2. Waiting Time Data

The data encountered in the second application involves data of waiting times (min) of 100 bank customers used by Bidram and Nadarajah [4]. We have compared the performance of the HLGW distribution with its submodels as well as with some other well-known distributions, such as the Weibull and the complementary Weibull geometric distribution. We also use the LR test to compare the HLGW distribution and its submodels.

The estimated values of the parameters, the standard errors, and the goodness-of-fit test statistics , , and K-S statistic for the data of waiting time are given in Table 4. These tables illustrate that HLGW model shows a good fit for this data set as compared to the other distributions. The CDF plot for empirical and HLGW distributions using the data set of waiting time in Figure 4(a) also confirms the fitness of HLGW distribution. Figure 4(b) expresses the TTT plot for this data set. Since the TTT plot of the data set is concave, the data set has increasing hazard rate shape. The LR test statistics of the hypotheses : HLW versus : HLGW, : HLE versus : HLGW, and : HLNH versus : HLGW are 5.501 ( value = ), 11.85 ( value = ), and 9.271 ( value = ). The HLGW distribution is significantly better than HLW, HLE, and HLNH distributions. There is no difference between HLGW and Weibull distribution based on the LR test; however the goodness-of-fit statistics , , and K-S statistic clearly show that HLGW distribution is better than Weibull distribution for the data.

Table 4: Estimates of models for the waiting time data.

Distributions
Estimates

Statistics

K-S
value

HLGW
0.3610
2.1425
0.0773
634.226
0.0176
0.1277
0.0368
0.9992

(0.1268)
(0.5257)
(0.0340)

HLW

1.2336
0.0772
639.730
0.0934
0.5740
0.0627
0.8267

(0.0978)
(0.0216)

HLE

0.1446
646.076
0.0583
0.3624
0.1213
0.1053

(0.0119)

HLNH
1.4044

0.0846
643.497
0.1032
0.6355
0.0877
0.4257

(0.3411)

(0.0316)

Weibull
0.0305
1.4582

637.461
0.0629
0.3961
0.0577
0.8938

(0.0095)
(0.1089)

CWG
0.5066
0.0126
2.2254
645.742
0.0745
0.4653
0.0485
0.9726

(0.0779)
(0.0097)
(1.7142)

(a) TTT plot

(b) CDF Plot

Figure 4: TTT plot and CDF plot for the waiting time data.

7. Conclusion

We have proposed a three-parameter lifetime generalized distribution, referred to as the half-logistic generalized Weibull (HLGW) distribution. The HLGW distribution has three other distributions like the half-logistic exponential, the half-logistic Weibull, and the half-logistic Nadarajah-Haghighi distributions, as its submodels. The new model exhibits a variety of shapes of the failure rate function, that is, increasing, increasing and then constant, decreasing, bathtub, unimodal, and decreasing-increasing-decreasing (DID) shapes. Various statistical properties of the HLGW distribution are derived and studied in detail. We estimated the parameters involved in the model by using the method of MLEs. Two real data sets are used to illustrate the flexibility, potentiality, and usefulness of HLGW distribution. It is concluded that HLGW model delivers better fitting than the other lifetime models and we hope that HLGW distribution may attract wider range of practical applications and this research may serve as a reference and benefit future research in the subject field of study.

Appendix

Fisher’s Information Matrix

The elements of Fisher’s information matrix analogous to , , and :

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

C. D. Lai, “Constructions and applications of lifetime distributions,” Applied Stochastic Models in Business and Industry, vol. 29, no. 2, pp. 127–140, 2013.

A. W. Marshall and I. Olkin, “A new method for adding a parameter to a family of distributions with application to the exponential and Weibull families,” Biometrika, vol. 84, no. 3, pp. 641–652, 1997.

B. Al-Zahrani and H. Sagor, “The Poisson-Lomax distribution,” Revista Colombiana de Estadistica, vol. 37, no. 1, pp. 225–245, 2014.

H. Bidram and S. Nadarajah, “A new lifetime model with decreasing, increasing, bathtub-shaped, and upside-down bathtub-shaped hazard rate function,” Statistics: A Journal of Theoretical and Applied Statistics, 2015.

C. Kus, “A new lifetime distribution,” Computational Statistics & Data Analysis, vol. 51, no. 9, pp. 4497–4509, 2007.

M. R. Gurvich, A. T. DiBenedetto, and S. V. Ranade, “A new statistical distribution for characterizing the random strength of brittle materials,” Journal of Materials Science, vol. 32, no. 10, pp. 2559–2564, 1997.

S. Nadarajah and S. Kotz, “On some recent modifications of Weibull distribution,” IEEE Transactions on Reliability, vol. 54, no. 4, pp. 561-562, 2005.

C. D. Lai, M. Xie, and D. N. P. Murthy, “Reply to "On some recent modifications of Weibull distribution",” IEEE Transactions on Reliability, vol. 54, no. 4, p. 563, 2005.

C. D. Lai, M. Xie, and D. N. P. Murthy, “A modified Weibull distribution,” IEEE Transactions on Reliability, vol. 52, no. 1, pp. 33–37, 2003.

M. Xie, Y. Tang, and T. N. Goh, “A modified Weibull extension with bathtub-shaped failure rate function,” Reliability Engineering & System Safety, vol. 76, no. 3, pp. 279–285, 2002.

T. Dimitrakopoulou, K. Adamidis, and S. Loukas, “A lifetime distribution with an upside-down bathtub-shaped hazard function,” IEEE Transactions on Reliability, vol. 56, no. 2, pp. 308–311, 2007.

G. M. Cordeiro, M. Alizadeh, and P. . Marinho, “The type I half-logistic family of distributions,” Journal of Statistical Computation and Simulation, vol. 86, no. 4, pp. 707–728, 2016.

S. D. Krishnarani, “On a power transformation of half-logistic distribution,” Journal of Probability and Statistics, vol. 2016, Article ID 2084236, 2016.

J. Kenney and E. S. Keeping, Mathematics of Statistics, Van Nostrand, Princeton, NJ, USA, 1962, Van Nostrand., Mathematics of Statistics.

J. J. Moors, “A Quantile Alternative for Kurtosis,” The American Statistician, vol. 37, no. 1, p. 25, 1988.

F. Proschan, “Theoretical Explanation of Observed Decreasing Failure Rate,” Technometrics, vol. 5, no. 3, pp. 375–383, 1963.

W. Lu and D. Shi, “A new compounding life distribution: the Weibull-Poisson distribution,” Journal of Applied Statistics, vol. 39, no. 1, pp. 21–38, 2012.

C. Tojeiro, F. Louzada, M. Roman, and P. Borges, “The complementary Weibull geometric distribution,” Journal of Statistical Computation and Simulation, vol. 84, no. 6, pp. 1345–1362, 2014.