A new modified Weibull distribution
Introduction
The Weibull distribution [28] has been used in many different fields with many applications, see for example [18]. The hazard function of the Weibull distribution can only be increasing, decreasing or constant. Thus it cannot be used to model lifetime data with a bathtub shaped hazard function, such as human mortality and machine life cycles. For many years, researchers have been developing various extensions and modified forms of the Weibull distribution, with number of parameters ranging from 2 to 5. The two-parameter flexible Weibull extension of Bebbington et al. [5] has a hazard function that can be increasing, decreasing or bathtub shaped. Zhang and Xie [31] studied the characteristics and application of the truncated Weibull distribution which has a bathtub shaped hazard function. A three-parameter model, called exponentiated Weibull distribution, was introduced by Mudholkar and Srivastave [17]. Another three-parameter model is by Marshall and Olkin [15] and called extended Weibull distribution. Xie et al. [30] proposed a three-parameter modified Weibull extension with a bathtub shaped hazard function. The modified Weibull (MW) distribution of Lai et al. [13] multiplies the Weibull cumulative hazard function by , which was later generalized to exponentiated form by Carrasco et al. [6]Recent studies of the modified Weibull include [11], [26], [27].
Among the four-parameter distributions, the additive Weibull distribution (AddW) of Xie and Lai [29] with cumulative distribution function (CDF) has a bathtub-shaped hazard function consisting of two Weibull hazards, one increasing () and one decreasing (). The modified Weibull distribution of Sarhan and Zaindin (SZMW) [21] can be derived from the additive Weibull distribution by setting . A four-parameter beta Weibull distribution was proposed by Famoye et al. [10]. Cordeiro et al. [8] introduced another four-parameter called the Kumaraswamy Weibull distribution.
Five-parameter modified Weibull distributions include Phani's modified Weibull [20], the beta modified Weibull (BMW) introduced by Silva et al. [24] and further studied by Nadarajah et al. [19]. The latest examples include the beta generalized Weibull distribution by Singla et al. [25], exponentiated generalized linear exponential distribution by Sarhan et al. [22] and the generalized Gomprtz distribution by El-Gohary et al. [9].
We propose a new lifetime distribution based on the Weibull and the modified Weibull (MW) distributions by combining them in a serial system. The hazard function of the new distribution is the sum of a Weibull hazard function and a modified Weibull hazard function. Section 2 gives definition, motivation and usefulness of this model and lists its sub-models. Section 3 considers properties of the new distribution such as hazard, moments and order statistics. Section 4 discusses estimation of the parameters. Two real data sets are analyzed in Section 5 and the results are compared with existing distributions. Section 6 concludes the paper.
Section snippets
Definition
We define a new modified Weibull distribution (NMW) by the following CDF:where , , , and are non-negative, with and being shape parameters and and being scale parameters and acceleration parameter.
The probability density function (PDF) isIt can be rewritten aswhere SW, hW, SMW and hMW are survival and hazard functions of the Weibull and modified Weibull
The hazard function
The hazard function can have many different shapes, including bathtub, as shown in Fig. 2. We can deduce from (6) that it is increasing if , , decreasing if , and and bathtub shaped otherwise.
It is desirable for a bathtub shaped hazard function to have a long useful life period [12], with relatively constant failure rate in the middle. A few distributions have this property, so does the NMW as shown in Fig. 3.
The moments
It is customary to derive the moments when a new distribution is
Parameter estimation
Given a random sample x1, …, xn from the NMW with parameters , the usual method of estimation is by maximum likelihood [7]. Other possible approaches include Bayesian estimation using Lindley approximation [14] or MCMC [26], [27].
The log-likelihood function is given bySetting the first partial derivatives of with respect to , , , and to zero, the likelihood equations are
Applications
In this section we provide results of fitting the NMW to two well-known data sets and compare its goodness-of-fit with other modified Weibull distributions using Kolmogorov–Smirnov (K–S) statistic, as well as Akaike information criterion (AIC) [2] and Bayesian information criterion (BIC) [23] values.
Sub-model of the NMW with
To simplify the statistical inference, it is always a good idea to reduce the number of parameters of any distribution and investigate how that affects the ability of the reduced model to fit the data. In this section we reduce the number of parameters from five to four, by setting . We test the reduced model against the original model . For each data set, Table 6 shows ML estimates of the four parameter NMW, the log-likelihood value under H0, likelihood ratio statistic (LRT)
Conclusions
A new distribution, based on Weibull and modified Weibull distributions, has been proposed and its properties studied. The idea is to combine two components in a serial system, so that the hazard function is either increasing or more importantly, bathtub shaped. Using a modified Weibull component, the distribution has flexibility to model the second peak in a distribution. We have shown that the new modified Weibull distribution fits certain well-known data sets better than existing
Acknowledgments
We would like to thank the referees for their comments and suggestions which improved the presentation of the paper. The first author wishes to thank the Saudi Arabia Culture Bureau in the UK and the Taif University for their financial support.
References (31)
- et al.
A flexible Weibull extension
Reliability Engineering & System Safety
(2007) - et al.
The Kumaraswamy Weibull distribution with application to failure data
Journal of the Franklin Institute
(2010) - et al.
The beta generalized Weibull distributionproperties and applications
Reliability Engineering & System Safety
(2012) - et al.
Modified Weibull model: a Bayes study using MCMC approach based on progressive censoring data
Reliability Engineering & System Safety
(2012) - et al.
A modified Weibull extension with bathtub-shaped failure rate function
Reliability Engineering System Safety
(2002) - et al.
On the upper truncated Weibull distribution and its reliability implications
Reliability Engineering System Safety
(2011) How to identify bathtub hazard rate
IEEE Transactions on Reliability
(1987)A new look at the statistical model identification
IEEE Transactions on Automatic Control
(1974)- et al.
A first course in order statistics
(1992) Analysis for the linear failure-rate life-testing distribution
Technometrics
(1974)