Multiple step-stress accelerated life test: the tampered failure rate model

Battacharyya and Soejoeti (1989) proposed the tampered failure rate model for step-stress accelerated life testing. In this note, their model is generalized from the simple (2-step) step-stress setting to the multiple (k-step, k > 2) setting. For the parametric setting where the life distribution under constant stress is Weibull, maximum likelihood estimation is investigated and the situation where the different stress levels are equispaced is looked at.     

Generalized exponential distributions

The three-parameter gamma and three-parameter Weibull distributions are commonly used for analysing any lifetime data or skewed data. Both distributions have several desirable properties, and nice physical interpretations. Because of the scale and shape parameters, both have quite a bit of flexibility for analysing different types of lifetime data. They have increasing as well as decreasing hazard rate depending on the shape parameter. Unfortunately both distributions also have certain drawbacks. This paper considers a three-parameter distribution which is a particular case of the exponentiated Weibull distribution originally proposed by Mudholkar, Srivastava & Freimer (1995) when the location parameter is not present. The study examines different properties of this model and observes that this family has some interesting features which are quite similar to those of the gamma family and the Weibull family, and certain distinct properties also. It appears this model can be used as an alternative to the gamma model or the Weibull model in many situations. One dataset is provided where the three-parameter generalized exponential distribution fits better than the three-parameter Weibull distribution or the three-parameter gamma distribution.     

The generalized inverse Weibull distribution

The inverse Weibull distribution has the ability to model failure rates which are quite common in reliability and biological studies. A three-parameter generalized inverse Weibull distribution with decreasing and unimodal failure rate is introduced and studied. We provide a comprehensive treatment of the mathematical properties of the new distribution including expressions for the moment generating function and the rth generalized moment. The mixture model of two generalized inverse Weibull distributions is investigated. The identifiability property of the mixture model is demonstrated. For the first time, we propose a location-scale regression model based on the log-generalized inverse Weibull distribution for modeling lifetime data. In addition, we develop some diagnostic tools for sensitivity analysis. Two applications of real data are given to illustrate the potentiality of the proposed regression model.     

Reliability statistical analysis about products of Birnbaum–Saunders fatigue life distribution for life test and progressive stress accelerated life test

Birnbaum–Saunders fatigue life distribution is an important failure model in the probability physical methods. It is more suitable for describing the life rules of fatigue failure products than common life distributions such as Weibull distribution and lognormal distribution. Besides, it is mainly applied to analytical research about fatigue failure and degradation failure of electronic product performance. The characteristic properties such as numerical characteristics and image features of density function and failure rate function are studied for generalized BS fatigue life distribution GBS(α, β, m) in this paper. Then the point estimates and approximate interval estimates of parameters are proposed for generalized BS fatigue life distribution GBS(α, β, m), and the precision of estimates are investigated by Monte Carlo simulations. Finally, when the scale parameter satisfies inverse power law model, the failure distribution model is given for the products of two-parameter BS fatigue life distribution BS(α, β) under progressive stress accelerated life test according to the time conversion idea of famous Nelson assumption, and then the points estimates of parameters are given.     

Estimation of Weibull distribution parameters for irregular interval group failure data with unknown failure times

SUMMARY This paper presents three methods for estimating Weibull distribution parameters for the case of irregular interval group failure data with unknown failure times. The methods are based on the concepts of the piecewise linear distribution function (PLDF), an average interval failure rate (AIFR) and sequential updating of the distribution function (SUDF), and use an analytical approach similar to that of Ackoff and Sasieni for regular interval group data. Results from a large number of simulated case problems generated with specified values of Weibull distribution parameters have been presented, which clearly indicate that the SUDF method produces near-perfect parameter estimates for all types of failure pattern. The performances of the PLDF and AIFR methods have been evaluated by goodness-of-fit testing and statistical confidence limits on the shape parameter. It has been found that, while the PLDF method produces acceptable parameter estimates, the AIFR method may fail for low and high shape parameter values that represent the cases of random and wear-out types of failure. A real-life application of the proposed methods is also presented, by analyzing failures of hydrogen make-up compressor valves in a petroleum refinery.     

The Weibull–Pareto Composite Family with Applications to the Analysis of Unimodal Failure Rate Data

The Weibull distribution is composited with Pareto model to obtain a flexible, reliable long-tailed parametric distribution for modeling unimodal failure rate data. The hazard function of the composite family accommodates decreasing and unimodal failure rates, which are separated by the boundary line of the space of shape parameter, gamma, when it equals to a known constant. The least square and maximum likelihood parameter estimation techniques are discussed. The advantages of using the proposed family are demonstrated and compared by illustrating well-known examples: guinea pigs survival time data, head and neck cancer data, and nasopharynx cancer survival data.     

Inverse Weibull power series distributions: properties and applications

Inverse Weibull (IW) distribution is one of the widely used probability distributions for nonnegative data modelling, specifically, for describing degradation phenomena of mechanical components. In this paper, by compounding IW and power series distributions we introduce a new lifetime distribution. The compounding procedure follows the same set-up carried out by Adamidis and Loukas [A lifetime distribution with decreasing failure rate. Stat Probab Lett. 1998;39:35–42]. We provide mathematical properties of this new distribution such as moments, estimation by maximum likelihood with censored data, inference for a large sample and the EM algorithm to determine the maximum likelihood estimates of the parameters. Furthermore, we characterize the proposed distributions using a simple relationship between two truncated moments and maximum entropy principle under suitable constraints. Finally, to show the flexibility of this type of distributions, we demonstrate applications of two real data sets.     

Approximate Bayesian analysis of doubly censored samples from mixture of two Weibull distributions

The purpose of the paper is to estimate the parameters of the two-component mixture of Weibull distribution under doubly censored samples using Bayesian approach. The choice of Weibull distribution is made due to its (i) capability to model failure time data from engineering, medical and biological sciences (ii) added advantages over the well-known lifetime distributions such as exponential, Raleigh, lognormal and gamma distribution in terms of flexibility, increasing and decreasing hazard rate and closed-form distribution function and hazard rate. The proposed two-component mixture of Weibull distribution is even more flexible than its conventional form. However, the estimation of the parameters from the proposed mixture is more complex. Further, we have assumed couple of loss functions under non informative prior for the Bayesian analysis of the parameters from the mixture model. As the resultant Bayes estimators and associated posterior risks cannot be derived in the closed form, we have used the importance sampling and Lindley’s approximation to obtain the approximate estimates for the parameters of the mixture model. The comparison between the performances of approximation techniques has been made on the basis of simulation study and real-life data analysis. The importance sampling is found to be better than Lindley’s approximation as it gives better estimation for shape and mixing parameters of the mixture model and computations under this technique are much easier/shorter than those under Lindley’s approximation.     

The role of weighted distributions in stochastic modeling

C. R. Rao pointed out that “The role of statistical methodology is to extract the relevant information from a given sample to answer specific questions about the parent population” and raised the question “What population does a sample represent”? Wrong specification can lead to invalid inference giving rise to a third kind of error. Rao introduced the concept of weighted distributions as a method of adjustment applicable to many situations.

In this paper, we study the relationship between the weighted distributions and the parent distributions in the context of reliability and life testing. These relationships depend on the nature of the weight function and give rise to interesting connections between the different ageing criteria of the two distributions. As special cases, the length biased distribution, the equilibrium distribution of the backward and forward recurrence times and the residual life distribution, which frequently arise in practice, are studied and their relationships with the original distribution are examined. Their survival functions, failure rates and mean residual life functions are compared and some characterization results are established.     

Likelihood-based confidence interval for the ratio of scale parameters of two independent Weibull distributions

The Weibull distribution is widely used in lifetime data analysis. For example, in studies on the time to the occurrence of tumors in human populations or in laboratory animals, the time of occurrence of tumors is generally assumed to be distributed as a Weibull distribution. Moreover, in engineering, the voltage levels at which failure occurred in electrical cable insulation has been shown to be distributed as a Weibull distribution. When comparing two independent Weibull distributions, it is often assumed that only the scale parameter is altered. In this paper, we propose a simple and accurate procedure to obtain inference concerning the ratio of the two scale parameters of two independent distributions. The performance of the proposed method is assessed through Monte Carlo simulation studies. The numerical results show that the proposed method is extremely accurate even for very small samples. The method is applied to a set of real-life data.     

A GLM approach to step-stress accelerated life testing with interval censoring

In this paper, we present a statistical inference procedure for the step-stress accelerated life testing (SSALT) model with Weibull failure time distribution and interval censoring via the formulation of generalized linear model (GLM). The likelihood function of an interval censored SSALT is in general too complicated to obtain analytical results. However, by transforming the failure time to an exponential distribution and using a binomial random variable for failure counts occurred in inspection intervals, a GLM formulation with a complementary log-log link function can be constructed. The estimations of the regression coefficients used for the Weibull scale parameter are obtained through the iterative weighted least square (IWLS) method, and the shape parameter is updated by a direct maximum likelihood (ML) estimation. The confidence intervals for these parameters are estimated through bootstrapping. The application of the proposed GLM approach is demonstrated by an industrial example.     

The Kumaraswamy modified Weibull distribution: theory and applications

A five-parameter extension of the Weibull distribution capable of modelling a bathtub-shaped hazard rate function is introduced and studied. The beauty and importance of the new distribution lies in its ability to model both monotone and non-monotone failure rates that are quite common in lifetime problems and reliability. The proposed distribution has a number of well-known lifetime distributions as special sub-models, such as the Weibull, extreme value, exponentiated Weibull, generalized Rayleigh and modified Weibull (MW) distributions, among others. We obtain quantile and generating functions, mean deviations, Bonferroni and Lorenz curves and reliability. We provide explicit expressions for the density function of the order statistics and their moments. For the first time, we define the log-Kumaraswamy MW regression model to analyse censored data. The method of maximum likelihood is used for estimating the model parameters and the observed information matrix is determined. Two applications illustrate the potentiality of the proposed distribution.     

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

Lifetimes of modern mechanic or electronic units usually exhibit bathtub-shaped failure rates. An appropriate probability distribution to model such data is the modified Weibull distribution proposed by Lai et al. [15]. This distribution has both the two-parameter Weibull and type-1 extreme value distribution as special cases. It is able to model lifetime data with monotonic and bathtub-shaped failure rates, and thus attracts some interest among researchers because of this property. In this paper, the procedure of obtaining the maximum likelihood estimates (MLEs) of the parameters for progressively type-2 censored and complete samples are studied. Existence and uniqueness of the MLEs are proved.     

Weibull accelerated life tests when there are competing causes of failure

Accelerated life testing of a product under more severe than normal conditions is commonly used to reduce test time and costs. Data collected at such accelerated conditions are used to obtain estimates of the parameters of a stress translation function. This function is then used to make inference about the product's life under normal operating conditions. We consider the problem of accelerated life tests when the product of interest is a p component series system. Each of the components is assumed to have an independent Weibull time to failure distribution with different shape parameters and different scale parameters which are increasing functions stress. A general model i s used for the scale parameter includes the standard engineering models as special This model also has an appealing biological interpretation     

A new compounding life distribution: the Weibull–Poisson distribution

In this paper, a new compounding distribution, named the Weibull–Poisson distribution is introduced. The shape of failure rate function of the new compounding distribution is flexible, it can be decreasing, increasing, upside-down bathtub-shaped or unimodal. A comprehensive mathematical treatment of the proposed distribution and expressions of its density, cumulative distribution function, survival function, failure rate function, the kth raw moment and quantiles are provided. Maximum likelihood method using EM algorithm is developed for parameter estimation. Asymptotic properties of the maximum likelihood estimates are discussed, and intensive simulation studies are conducted for evaluating the performance of parameter estimation. The use of the proposed distribution is illustrated with examples.     

Optimum Design for Type-I Step-stress Accelerated Life Tests of Two-parameter Weibull Distributions

In this article, we focus on the general k-step step-stress accelerated life tests with Type-I censoring for two-parameter Weibull distributions based on the tampered failure rate (TFR) model. We get the optimum design for the tests under the criterion of the minimization of the asymptotic variance of the maximum likelihood estimate of the pth percentile of the lifetime under the normal operating conditions. Optimum test plans for the simple step-stress accelerated life tests under Type-I censoring are developed for the Weibull distribution and the exponential distribution in particular. Finally, an example is provided to illustrate the proposed design and a sensitivity analysis is conducted to investigate the robustness of the design.     

Bayes prediction of number of failures in weibull samples

This work is concerned with the Bayesian prediction problem of the number of components which will fail in a future time interval, when the failure times are Weibull distributed. Both the 1-sample and the 2-sample prediction problems are dealed with, and some choices of the prior densities on the distribution parameters are discussed which are relatively easy to work with and allow different degrees of knowledge on the failure mechanism to be incorporated in the predictive procedure. Useful relations between the predictive distribution on the number of future failures and the predictive distribution on the future failure times are derived. Numerical examples are also given.     

Optimal simple step stress accelerated life test design for reliability prediction

A step stress accelerated life testing model is presented to obtain the optimal hold time at which the stress level is changed. The experimental test is designed to minimize the asymptotic variance of reliability estimate at time ζ$\zeta$. A Weibull distribution is assumed for the failure time at any constant stress level. The scale parameter of the Weibull failure time distribution at constant stress levels is assumed to be a log-linear function of the stress level. The maximum likelihood function is given for the step stress accelerated life testing model with Type I censoring, from which the asymptotic variance and the Fisher information matrix are obtained. An optimal test plan with the minimum asymptotic variance of reliability estimate at time ζ$\zeta$ is determined.     

Analysis of hybrid censored competing risks data

In this paper, we consider the analysis of hybrid censored competing risks data, based on Cox's latent failure time model assumptions. It is assumed that lifetime distributions of latent causes of failure follow Weibull distribution with the same shape parameter, but different scale parameters. Maximum likelihood estimators (MLEs) of the unknown parameters can be obtained by solving a one-dimensional optimization problem, and we propose a fixed-point type algorithm to solve this optimization problem. Approximate MLEs have been proposed based on Taylor series expansion, and they have explicit expressions. Bayesian inference of the unknown parameters are obtained based on the assumption that the shape parameter has a log-concave prior density function, and for the given shape parameter, the scale parameters have Beta–Gamma priors. We propose to use Markov Chain Monte Carlo samples to compute Bayes estimates and also to construct highest posterior density credible intervals. Monte Carlo simulations are performed to investigate the performances of the different estimators, and two data sets have been analysed for illustrative purposes.     

A comparison of parametric and semi-parametric survival models with artificial neural networks

Survival models are used to examine data in the event of an occurrence. These are discussed in various types including parametric, non-parametric and semi-parametric models. Parametric models require a clear distribution of survival time, and semi-parametric models assume proportional hazards. Among these models, the non-parametric model of artificial neural network has the fewest assumptions and can be often replaced by other models. Given the importance of distribution Weibull survival models in this study of simulation shape parameter of the Weibull distribution have been assumed as 1, 2 and 3, and also the average rate at levels of 0%–75% have been censored. The values predicted by the neural network forecasting model with parametric survival and Cox regression models were compared. This comparison considering levels of complexity due to the hazard model using the ROC curve and the corresponding tests have been carried out.