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.     

The Marshall-Olkin logistic-exponential distribution

We introduce a three-parameter extension of the exponential distribution which contains as sub-models the exponential, logistic-exponential and Marshall-Olkin exponential distributions. The new model is very flexible and its associated density function can be decreasing or unimodal. Further, it can produce all of the four major shapes of the hazard rate, that is, increasing, decreasing, bathtub and upside-down bathtub. Given that closed-form expressions are available for the survival and hazard rate functions, the new distribution is quite tractable. It can be used to analyze various types of observations including censored data. Computable representations of the quantile function, ordinary and incomplete moments, generating function and probability density function of order statistics are obtained. The maximum likelihood method is utilized to estimate the model parameters. A simulation study is carried out to assess the performance of the maximum likelihood estimators. Two actual data sets are used to illustrate the applicability of the proposed model.     

Generalized exponential distribution: Existing results and some recent developments

Mudholkar and Srivastava [1993. Exponentiated Weibull family for analyzing bathtub failure data. IEEE Trans. Reliability 42, 299–302] introduced three-parameter exponentiated Weibull distribution. Two-parameter exponentiated exponential or generalized exponential distribution is a particular member of the exponentiated Weibull distribution. Generalized exponential distribution has a right skewed unimodal density function and monotone hazard function similar to the density functions and hazard functions of the gamma and Weibull distributions. It is observed that it can be used quite effectively to analyze lifetime data in place of gamma, Weibull and log-normal distributions. The genesis of this model, several properties, different estimation procedures and their properties, estimation of the stress-strength parameter, closeness of this distribution to some of the well-known distribution functions are discussed in this article.     

Proportional Hazards Model for Multivariate Failure Time Data

Multivariate failure time data also referred to as correlated or clustered failure time data, often arise in survival studies when each study subject may experience multiple events. Statistical analysis of such data needs to account for intracluster dependence. In this article, we consider a bivariate proportional hazards model using vector hazard rate, in which the covariates under study have different effect on two components of the vector hazard rate function. Estimation of the parameters as well as base line hazard function are discussed. Properties of the estimators are investigated. We illustrated the method using two real life data. A simulation study is reported to assess the performance of the estimator.     

Estimating the turning point of the log-logistic hazard function in the presence of long-term survivors with an application for uterine cervical cancer data

The hazard function plays an important role in cancer patient survival studies, as it quantifies the instantaneous risk of death of a patient at any given time. Often in cancer clinical trials, unimodal hazard functions are observed, and it is of interest to detect (estimate) the turning point (mode) of hazard function, as this may be an important measure in patient treatment strategies with cancer. Moreover, when patient cure is a possibility, estimating cure rates at different stages of cancer, in addition to their proportions, may provide a better summary of the effects of stages on survival rates. Therefore, the main objective of this paper is to consider the problem of estimating the mode of hazard function of patients at different stages of cervical cancer in the presence of long-term survivors. To this end, a mixture cure rate model is proposed using the log-logistic distribution. The model is conveniently parameterized through the mode of the hazard function, in which cancer stages can affect both the cured fraction and the mode. In addition, we discuss aspects of model inference through the maximum likelihood estimation method. A Monte Carlo simulation study assesses the coverage probability of asymptotic confidence intervals.     

Estimation de densités unimodales

This paper proposes a new nonparametric unimodal estimator of a unimodal probability density function, in the case where the mode is known. The classical solution to this problem is the maximum-likelihood estimator under monotonicity constraint, considered by Grenander (1956). Our approach is based on a unimodal rearrangement of the kernel estimator of the density. Asymptotic properties of this estimator are studied, and its small-sample behaviour is examined through simulations.     

An extension of the generalized linear failure rate distribution

In this article, we introduce a new extension of the generalized linear failure rate (GLFR) distributions. It includes some well-known lifetime distributions such as extension of generalized exponential and GLFR distributions as special sub-models. In addition, it can have a constant, decreasing, increasing, upside-down bathtub (unimodal), and bathtub-shaped hazard rate function (hrf) depending on its parameters. We provide some of its statistical properties such as moments, quantiles, skewness, kurtosis, hrf, and reversible hrf. The maximum likelihood estimation of the parameters is also discussed. At the end, a real dataset is given to illustrate the usefulness of this new distribution in analyzing lifetime data.     

Exponentiated Weibull regression for time-to-event data

The Weibull, log-logistic and log-normal distributions are extensively used to model time-to-event data. The Weibull family accommodates only monotone hazard rates, whereas the log-logistic and log-normal are widely used to model unimodal hazard functions. The increasing availability of lifetime data with a wide range of characteristics motivate us to develop more flexible models that accommodate both monotone and nonmonotone hazard functions. One such model is the exponentiated Weibull distribution which not only accommodates monotone hazard functions but also allows for unimodal and bathtub shape hazard rates. This distribution has demonstrated considerable potential in univariate analysis of time-to-event data. However, the primary focus of many studies is rather on understanding the relationship between the time to the occurrence of an event and one or more covariates. This leads to a consideration of regression models that can be formulated in different ways in survival analysis. One such strategy involves formulating models for the accelerated failure time family of distributions. The most commonly used distributions serving this purpose are the Weibull, log-logistic and log-normal distributions. In this study, we show that the exponentiated Weibull distribution is closed under the accelerated failure time family. We then formulate a regression model based on the exponentiated Weibull distribution, and develop large sample theory for statistical inference. We also describe a Bayesian approach for inference. Two comparative studies based on real and simulated data sets reveal that the exponentiated Weibull regression can be valuable in adequately describing different types of time-to-event data.     

A Note On Beta Inverse-Weibull Distribution

The inverse Weibull distribution is one of the widely applied distribution for problems in reliability theory. In this article, we introduce a generalization—referred to as the Beta Inverse-Weibull distribution—generated from the logit of a beta random variable. We provide a comprehensive treatment of the mathematical properties of the Beta Inverse-Weibull distribution. The shapes of the corresponding probability density function and the hazard rate function have been obtained and graphical illustrations have been given. The distribution is found to be unimodal. Results for the non central moments are obtained. The relationship between the parameters and the mean, variance, skewness, and kurtosis are provided. The method of maximum likelihood is proposed for estimating the model parameters. We hope that this generalization will attract wider applicability to the problems in reliability theory and mechanical engineering.     

Generalized Linear Failure Rate Distribution

The exponential and Rayleigh are the two most commonly used distributions for analyzing lifetime data. These distributions have several desirable properties and nice physical interpretations. Unfortunately, the exponential distribution only has constant failure rate and the Rayleigh distribution has increasing failure rate. The linear failure rate distribution generalizes both these distributions which may have non increasing hazard function also. This article introduces a new distribution, which generalizes linear failure rate distribution. This distribution generalizes the well-known (1) exponential distribution, (2) linear failure rate distribution, (3) generalized exponential distribution, and (4) generalized Rayleigh distribution. The properties of this distribution are discussed in this article. The maximum likelihood estimates of the unknown parameters are obtained. A real data set is analyzed and it is observed that the present distribution can provide a better fit than some other very well-known distributions.     

The inverted Nadarajah–Haghighi distribution: estimation methods and applications

The inverted (or inverse) distributions are sometimes very useful to explore additional properties of the phenomenons which non-inverted distributions cannot. We introduce a new inverted model called the inverted Nadarajah–Haghighi distribution which exhibits decreasing and unimodal (right-skewed) density while the hazard rate shapes are decreasing and upside-down bathtub. Our main focus is the estimation (from both frequentist and Bayesian points of view) of the unknown parameters along with some mathematical properties of the new model. The Bayes estimators and the associated credible intervals are obtained using Markov Chain Monte Carlo techniques under squared error loss function. The gamma priors are adopted for both scale and shape parameters. The potentiality of the distribution is analysed by means of two real data sets. In fact, it is found to be superior in its ability to sufficiently model the data as compared to the inverted Weibull, inverted Rayleigh, inverted exponential, inverted gamma, inverted Lindley and inverted power Lindley models.     

The compound class of linear failure rate-power series distributions: Model,properties, and applications

In this article, a new class of distributions is introduced, which generalizes the linear failure rate distribution and is obtained by compounding this distribution and power series class of distributions. This new class of distributions is called the linear failure rate-power series distributions and contains some new distributions such as linear failure rate-geometric, linear failure rate-Poisson, linear failure rate-logarithmic, linear failure rate-binomial distributions, and Rayleigh-power series class of distributions. Some former works such as exponential-power series class of distributions, exponential-geometric, exponential-Poisson, and exponential-logarithmic distributions are special cases of the new proposed model. The ability of the linear failure rate-power series class of distributions is in covering five possible hazard rate function, that is, increasing, decreasing, upside-down bathtub (unimodal), bathtub and increasing-decreasing-increasing shaped. Several properties of this class of distributions such as moments, maximum likelihood estimation procedure via an EM-algorithm and inference for a large sample, are discussed in this article. In order to show the flexibility and potentiality, the fitted results of the new class of distributions and some of its submodels are compared using two real datasets.     

Sensitivity analysis of reliability functions of the exponential power series lifetime distribution

ABSTRACT

Hazard rate functions are often used in modeling of lifetime data. The Exponential Power Series (EPS) family has a monotone hazard rate function. In this article, the influence of input factors such as time and parameters on the variability of hazard rate function is assessed by local and global sensitivity analysis. Two different indices based on local and global sensitivity indices are presented. The simulation results for two datasets show that the hazard rate functions of the EPS family are sensitive to input parameters. The results also show that the hazard rate function of the EPS family is more sensitive to the exponential distribution than power series distributions.     

Testing dependence between the failure time and failure modes: An application of enlarged filtration

The model of independent competing risks provides no information for the assessment of competing failure modes if the failure mechanisms underlying these modes are coupled. Models for dependent competing risks in the literature can be distinguished on the basis of the functional behaviour of the conditional probability of failure due to a particular failure mode given that the failure time exceeds a fixed time, as a function of time. There is an interesting link between monotonicity of such conditional probability and dependence between failure time and failure mode, via crude hazard rates. In this paper, we propose tests for testing the dependence between failure time and failure mode using the crude hazards and using the conditional probabilities mentioned above. We establish the equivalence between the two approaches and provide an asymptotically efficient weight function under a sequence of local alternatives. The tests are applied to simulated data and to mortality follow-up data.     

Inferential statistics on the dynamic system model with time-dependent failure-rate

This study focuses on the classical and Bayesian analysis of a k-components load-sharing parallel system in which components have time-dependent failure rates. In the classical set up, the maximum likelihood estimates of the load-share parameters with their standard errors (SEs) are obtained. (1?γ) 100% simultaneous and two bootstrap confidence intervals for the parameters and system reliability and hazard functions have been constructed. Further, on recognizing the fact that life-testing experiments are very time consuming, the parameters involved in the failure time distribution of the system are expected to follow some random variations. Therefore, Bayes estimates along with their posterior SEs of the parameters and system reliability and hazard functions are obtained by assuming gamma and Jeffrey's priors of the unknown parameters. Markov chain Monte Carlo technique such as Gibbs sampler has been used to obtain Bayes estimates and highest posterior density credible intervals.     

A new four-parameter lifetime distribution

Generalizing lifetime distributions is always precious for applied statisticians. In this paper, we introduce a new four-parameter generalization of the exponentiated power Lindley (EPL) distribution, called the exponentiated power Lindley geometric (EPLG) distribution, obtained by compounding EPL and geometric distributions. The new distribution arises in a latent complementary risks scenario, in which the lifetime associated with a particular risk is not observable; rather, we observe only the maximum lifetime value among all risks. The distribution exhibits decreasing, increasing, unimodal and bathtub-shaped hazard rate functions, depending on its parameters. It contains several lifetime distributions as particular cases: EPL, new generalized Lindley, generalized Lindley, power Lindley and Lindley geometric distributions. We derive several properties of the new distribution such as closed-form expressions for the density, cumulative distribution function, survival function, hazard rate function, the rth raw moment, and also the moments of order statistics. Moreover, we discuss maximum likelihood estimation and provide formulas for the elements of the Fisher information matrix. Simulation studies are also provided. Finally, two real data applications are given for showing the flexibility and potentiality of the new distribution.     

The Kumaraswamy generalized gamma distribution with application in survival analysis

We introduce and study the so-called Kumaraswamy generalized gamma distribution that is capable of modeling bathtub-shaped hazard rate functions. The beauty and importance of this distribution lies in its ability to model monotone and non-monotone failure rate functions, which are quite common in lifetime data analysis and reliability. The new distribution has a large number of well-known lifetime special sub-models such as the exponentiated generalized gamma, exponentiated Weibull, exponentiated generalized half-normal, exponentiated gamma, generalized Rayleigh, among others. Some structural properties of the new distribution are studied. We obtain two infinite sum representations for the moments and an expansion for the generating function. We calculate the density function of the order statistics and an expansion for their moments. The method of maximum likelihood and a Bayesian procedure are adopted for estimating the model parameters. The usefulness of the new distribution is illustrated in two real data sets.     

The exponentiated generalized gamma distribution with application to lifetime data

A four-parameter extension of the generalized gamma distribution capable of modelling a bathtub-shaped hazard rate function is defined and studied. The beauty and importance of this distribution lies in its ability to model monotone and non-monotone failure rate functions, which are quite common in lifetime data analysis and reliability. The new distribution has a number of well-known lifetime special sub-models, such as the exponentiated Weibull, exponentiated generalized half-normal, exponentiated gamma and generalized Rayleigh, among others. We derive two infinite sum representations for its moments. We calculate the density of the order statistics and two expansions for their moments. The method of maximum likelihood is used for estimating the model parameters and the observed information matrix is obtained. Finally, a real data set from the medical area is analysed.     

Bayesian nonparametric inference for unimodal skew-symmetric distributions

This paper studies the case where the observations come from a unimodal and skew density function with an unknown mode. The skew-symmetric representation of such a density has a symmetric component which can be written as a scale mixture of uniform densities. A Dirichlet process (DP) prior is assigned to mixing distribution. We also assume prior distributions for the mode and the skewed component. A computational approach is used to obtain the Bayes estimate of the components. An example is given to illustrate the approach.     

Density and hazard rate estimation for right-censored data by using wavelet methods

This paper describes a wavelet method for the estimation of density and hazard rate functions from randomly right-censored data. We adopt a nonparametric approach in assuming that the density and hazard rate have no specific parametric form. The method is based on dividing the time axis into a dyadic number of intervals and then counting the number of events within each interval. The number of events and the survival function of the observations are then separately smoothed over time via linear wavelet smoothers, and then the hazard rate function estimators are obtained by taking the ratio. We prove that the estimators have pointwise and global mean-square consistency, obtain the best possible asymptotic mean integrated squared error convergence rate and are also asymptotically normally distributed. We also describe simulation experiments that show that these estimators are reasonably reliable in practice. The method is illustrated with two real examples. The first uses survival time data for patients with liver metastases from a colorectal primary tumour without other distant metastases. The second is concerned with times of unemployment for women and the wavelet estimate, through its flexibility, provides a new and interesting interpretation.