A new modified Weibull distribution

https://doi.org/10.1016/j.ress.2012.10.018Get rights and content

Abstract

We introduce a new lifetime distribution by considering a serial system with one component following a Weibull distribution and another following a modified Weibull distribution. We study its mathematical properties including moments and order statistics. The estimation of parameters by maximum likelihood is discussed. We demonstrate that the proposed distribution fits two well-known data sets better than other modified Weibull distributions including the latest beta modified Weibull distribution. The model can be simplified by fixing one of the parameters and it still provides a better fit than existing models.

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 αxβ by eλx, which was later generalized to exponentiated form by Carrasco et al. [6]F(x)=(1eαxβeλx)θ,x0.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) F(x)=1eαxθβxγ,x0,has a bathtub-shaped hazard function consisting of two Weibull hazards, one increasing (θ>1) and one decreasing (0<γ<1). The modified Weibull distribution of Sarhan and Zaindin (SZMW) [21] can be derived from the additive Weibull distribution by setting θ=1. 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:F(x)=1eαxθβxγeλx,x0,where α, β, θ, γ and λ are non-negative, with θ and γ being shape parameters and α and β being scale parameters and λ acceleration parameter.

The probability density function (PDF) isf(x)=(αθxθ1+β(γ+λx)xγ1eλx)eαxθβxγeλx,x>0.It can be rewritten asf(x)=[hW(x;α,θ)+hMW(x;β,γ,λ)]SW(x,α,θ)SMW(x;β,γ,λ),where 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 θ, γ1, decreasing if θ, γ<1 and λ=0 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 byL=i=1nln(β(γ+λxi)xiγ1eλxi+αθxiθ1)αi=1nxiθβi=1nxiγeλxi.Setting the first partial derivatives of with respect to α, β, θ, γ and λ to zero, the likelihood equations arei=1nθxiθ1h(xi;α,β,γ,θ,λ)i=1nxiθ=0,i=1n(γ+

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 γ=1

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 γ=1. We test the reduced model H0:γ=1) against the original model Ha:γ1). 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)

  • M. Carrasco et al.

    A generalized modified Weibull distribution for lifetime modeling

    Computational Statistics and Data Analysis

    (2008)
  • R.A. Fisher

    On the mathematical foundation of theoretical statistics

    Philosophical Transactions of the Royal Society A

    (1922)
  • A. El-Gohary et al.

    The generalized Gompertz distribution

    Applied Mathematical Modeling

    (2011)
  • F. Famoye et al.

    The beta-Weibull distribution

    Journal of Statistical Theory and Applications

    (2005)
  • H. Jiang et al.

    On MLEs of the parameters of a modified Weibull distribution for progressively type-2 censored samples

    Journal of Applied Statistical Science

    (2010)
  • Cited by (0)

    View full text