A semiparametric extended hazard regression model with time-dependent covariates

We introduce a general class of semiparametric hazard regression models, called extended hazard (EH) models, that are designed to accommodate various survival schemes with time-dependent covariates. The EH model contains both the Cox model and the accelerated failure time (AFT) model as its subclasses so that we can use this nested structure to perform model selection between the Cox model and the AFT model. A class of estimating equations using counting process and martingale techniques is developed to estimate the regression parameters of the proposed model. The performance of the estimating procedure and the impact of model misspecification are assessed through simulation studies. Two data examples, Stanford heart transplant data and Mediterranean fruit flies, egg-laying data, are used to demonstrate the usefulness of the EH model.     

Bayesian analysis of generalized odds-rate hazards models for survival data

In the analysis of censored survival data Cox proportional hazards model (1972) is extremely popular among the practitioners. However, in many real-life situations the proportionality of the hazard ratios does not seem to be an appropriate assumption. To overcome such a problem, we consider a class of nonproportional hazards models known as generalized odds-rate class of regression models. The class is general enough to include several commonly used models, such as proportional hazards model, proportional odds model, and accelerated life time model. The theoretical and computational properties of these models have been re-examined. The propriety of the posterior has been established under some mild conditions. A simulation study is conducted and a detailed analysis of the data from a prostate cancer study is presented to further illustrate the proposed methodology.     

Extended Hazard Regression Model for Reliability and Survival Analysis

We propose an extended hazard regression model which allows the spread parameter to be dependent on covariates. This allows a broad class of models which includes the most common hazard models, such as the proportional hazards model, the accelerated failure time model and a proportional hazards/accelerated failure time hybrid model with constant spread parameter. Simulations based on sub-classes of this model suggest that maximum likelihood performs well even when only small or moderate-size data sets are available and the censoring pattern is heavy. The methodology provides a broad framework for analysis of reliability and survival data. Two numerical examples illustrate the results.     

Modelling accelerated life test data by using a Bayesian approach

Summary. Because of the high reliability of many modern products, accelerated life tests are becoming widely used to obtain timely information about their time-to-failure distributions. We propose a general class of accelerated life testing models which are motivated by the actual failure process of units from a limited failure population with a positive probability of not failing during the technological lifetime. We demonstrate a Bayesian approach to this problem, using a new class of models with non-monotone hazard rates, the hazard model with potential scope for use far beyond accelerated life testing. Our methods are illustrated with the modelling and analysis of a data set on lifetimes of printed circuit boards under humidity accelerated life testing.     

M-estimation of the accelerated failure time model under a convex discrepancy function

Based on the Kaplan–Meier weight functions, we introduce a class of M-estimators of regression parameters for the accelerated failure time (AFT) model with right censored data. The proposed M-estimator is root-n consistent and asymptotically normal under appropriate assumptions. For robustness analysis, we also derive the corresponding influence functions. Appropriate criteria are developed for tests of hypotheses concerning regression parameters. The results are applied to several particular cases. We evaluate the finite-sample performance of the proposed methods through extensive simulation studies.     

Maximum likelihood estimation of the parameters of a multiple step-stress model from the Birnbaum-Saunders distribution under time-constraint: A comparative study

The cumulative exposure model (CEM) is a commonly used statistical model utilized to analyze data from a step-stress accelerated life testing which is a special class of accelerated life testing (ALT). In practice, researchers conduct ALT to: (1) determine the effects of extreme levels of stress factors (e.g., temperature) on the life distribution, and (2) to gain information on the parameters of the life distribution more rapidly than under normal operating (or environmental) conditions. In literature, researchers assume that the CEM is from well-known distributions, such as the Weibull family. This study, on the other hand, considers a p-step-stress model with q stress factors from the two-parameter Birnbaum-Saunders distribution when there is a time constraint on the duration of the experiment. In this comparison paper, we consider different frameworks to numerically compute the point estimation for the unknown parameters of the CEM using the maximum likelihood theory. Each framework implements at least one optimization method; therefore, numerical examples and extensive Monte Carlo simulations are considered to compare and numerically examine the performance of the considered estimation frameworks.     

Testing Goodness of Fit of Parametric AFT and PH Models with Residuals

In this article, we consider the problems of testing the goodness of fit of the parametric accelerated failure time model and the Cox proportional hazards model. We consider omnibus test statistics based on residuals. The statistical distributions of Kolmogorov, Cramer-von Mises–Smirnov, and Anderson–Darling statistics are all investigated by means of Monte Carlo simulations. Type-I, Type-II, and independent random censoring situations are all considered in this study. A Monte Carlo power study has also been carried out for these tests to distinguish between various baseline models, which reveals that the Anderson–Darling test performs better than the others.     

Generalized M-estimation for the accelerated failure time model

The accelerated failure time (AFT) model is an important regression tool to study the association between failure time and covariates. In this paper, we propose a robust weighted generalized M (GM) estimation for the AFT model with right-censored data by appropriately using the Kaplan–Meier weights in the GM–type objective function to estimate the regression coefficients and scale parameter simultaneously. This estimation method is computationally simple and can be implemented with existing software. Asymptotic properties including the root-n consistency and asymptotic normality are established for the resulting estimator under suitable conditions. We further show that the method can be readily extended to handle a class of nonlinear AFT models. Simulation results demonstrate satisfactory finite sample performance of the proposed estimator. The practical utility of the method is illustrated by a real data example.     

Linear Signed-Rank Tests for Paired Survival Data Subject to a Common Censoring Time

In this paper we consider rank-based tests for paired survival data, in which pair members are subject to the same right censoring time. Linear signed-rank tests have already been developed for the two-treatment problem in which pair members receive the opposite treatments. Assuming a bivariate accelerated failure time model, we extend this class of linear signed-rank tests to the case of multiple covariates, making this methodology applicable to more complicated experimental designs. These tests can be reformulated as weighted sums of contigency table measures, giving an alternative method of computation and intuitive view of how these tests work. A simulation study of their small-sample performance relative to other tests demonstrates that the linear signed-rank tests have greater power in cases of moderately to highly correlated data.     

Tests of independence for censored bivariate failure time data

Bivariate failure time data is widely used in survival analysis, for example, in twins study. This article presents a class of chi2-type tests for independence between pairs of failure times after adjusting for covariates. A bivariate accelerated failure time model is proposed for the joint distribution of bivariate failure times while leaving the dependence structures for related failure times completely unspecified. Theoretical properties of the proposed tests are derived and variance estimates of the test statistics are obtained using a resampling technique. Simulation studies show that the proposed tests are appropriate for practical use. Two examples including the study of infection in catheters for patients on dialysis and the diabetic retinopathy study are also given to illustrate the methodology.     

Bayes Analysis and Comparison of Accelerated Weibull and Accelerated Birnbaum–Saunders Models

Several models are proposed in the literature for modeling fatigue data resulting from materials subject to cyclic stress and strain. Accelerated Weibull and accelerated Birnbaum–Saunders distributions are most commonly used models. Whereas the accelerated Weibull model is easier compared to accelerated Birnbaum–Saunders, it fails to represent the situation equally well. The present article focuses on Bayes analysis of the two models and provides a comparison based on some important Bayesian tools. Model compatibility study using predictive simulation ideas is preceded by the said comparison. Throughout, the posterior simulations are carried out by Markov chain Monte Carlo procedure.     

On the estimators and tests for the semiparametric hazards regression model

In the accelerated hazards regression model with censored data, estimation of the covariance matrices of the regression parameters is difficult, since it involves the unknown baseline hazard function and its derivative. This paper provides simple but reliable procedures that yield asymptotically normal estimators whose covariance matrices can be easily estimated. A class of weight functions are introduced to result in the estimators whose asymptotic covariance matrices do not involve the derivative of the unknown hazard function. Based on the estimators obtained from different weight functions, some goodness-of-fit tests are constructed to check the adequacy of the accelerated hazards regression model. Numerical simulations show that the estimators and tests perform well. The procedures are illustrated in the real world example of leukemia cancer. For the leukemia cancer data, the issue of interest is a comparison of two groups of patients that had two different kinds of bone marrow transplants. It is found that the difference of the two groups are well described by a time-scale change in hazard functions, i.e., the accelerated hazards model.     

Exact inference for a simple step-stress model with competing risks for failure from exponential distribution under Type-II censoring

In reliability analysis, accelerated life-testing allows for gradual increment of stress levels on test units during an experiment. In a special class of accelerated life tests known as step-stress tests, the stress levels increase discretely at pre-fixed time points, and this allows the experimenter to obtain information on the parameters of the lifetime distributions more quickly than under normal operating conditions. Moreover, when a test unit fails, there are often more than one fatal cause for the failure, such as mechanical or electrical. In this article, we consider the simple step-stress model under Type-II censoring when the lifetime distributions of the different risk factors are independently exponentially distributed. Under this setup, we derive the maximum likelihood estimators (MLEs) of the unknown mean parameters of the different causes under the assumption of a cumulative exposure model. The exact distributions of the MLEs of the parameters are then derived through the use of conditional moment generating functions. Using these exact distributions as well as the asymptotic distributions and the parametric bootstrap method, we discuss the construction of confidence intervals for the parameters and assess their performance through Monte Carlo simulations. Finally, we illustrate the methods of inference discussed here with an example.     

Exponential progressive step-stress life-testing with link function based on Box–Cox transformation

In order to quickly extract information on the life of a product, accelerated life-tests are usually employed. In this article, we discuss a k-stage step-stress accelerated life-test with M-stress variables when the underlying data are progressively Type-I group censored. The life-testing model assumed is an exponential distribution with a link function that relates the failure rate and the stress variables in a linear way under the Box–Cox transformation, and a cumulative exposure model for modelling the effect of stress changes. The classical maximum likelihood method as well as a fully Bayesian method based on the Markov chain Monte Carlo (MCMC) technique is developed for inference on all the parameters of this model. Numerical examples are presented to illustrate all the methods of inference developed here, and a comparison of the ML and Bayesian methods is also carried out.     

Accelerated Rates Regression Models for Recurrent Failure Time Data

In this article, we formulate a semiparametric model for counting processes in which the effect of covariates is to transform the time scale for a baseline rate function. We assume an arbitrary dependence structure for the counting process and propose a class of estimating equations for the regression parameters. Asymptotic results for these estimators are derived. In addition, goodness of fit methods for assessing the adequacy of the accelerated rates model are proposed. The finite-sample behavior of the proposed methods is examined in simulation studies, and data from a chronic granulomatous disease study are used to illustrate the methodology.     

Diagnostic tools in generalized Weibull linear regression models

We propose some statistical tools for diagnosing the class of generalized Weibull linear regression models [A.A. Prudente and G.M. Cordeiro, Generalized Weibull linear models, Comm. Statist. Theory Methods 39 (2010), pp. 3739–3755]. This class of models is an alternative means of analysing positive, continuous and skewed data and, due to its statistical properties, is very competitive with gamma regression models. First, we show that the Weibull model induces ma-ximum likelihood estimators asymptotically more efficient than the gamma model. Standardized residuals are defined, and their statistical properties are examined empirically. Some measures are derived based on the case-deletion model, including the generalized Cook's distance and measures for identifying influential observations on partial F-tests. The results of a simulation study conducted to assess behaviour of the global influence approach are also presented. Further, we perform a local influence analysis under the case-weights, response and explanatory variables perturbation schemes. The Weibull, gamma and other Weibull-type regression models are fitted into three data sets to illustrate the proposed diagnostic tools. Statistical analyses indicate that the Weibull model fitted into these data yields better fits than other common alternative models.     

On a variance stabilizing model and its application to genomic data

In this paper, we propose a model based on a class of symmetric distributions, which avoids the transformation of data, stabilizes the variance of the observations, and provides robust estimation of parameters and high flexibility for modeling different types of data. Probabilistic and statistical aspects of this new model are developed throughout the article, which include mathematical properties, estimation of parameters and inference. The obtained results are illustrated by means of real genomic data.     

Bayesian parametric accelerated failure time spatial model and its application to prostate cancer

Prostate cancer (PrCA) is the most common cancer diagnosed in American men and the second leading cause of death from malignancies. There are large geographical variation and racial disparities existing in the survival rate of PrCA. Much work on the spatial survival model is based on the proportional hazards (PH) model, but few focused on the accelerated failure time (AFT) model. In this paper, we investigate the PrCA data of Louisiana from the Surveillance, Epidemiology, and End Results program and the violation of the PH assumption suggests that the spatial survival model based on the AFT model is more appropriate for this data set. To account for the possible extra-variation, we consider spatially referenced independent or dependent spatial structures. The deviance information criterion is used to select a best-fitting model within the Bayesian frame work. The results from our study indicate that age, race, stage, and geographical distribution are significant in evaluating PrCA survival.     

Models Comparison for Step-Stress Accelerated Life Testing

In this article, four basic models for step-stress accelerated life testing are introduced and compared: cumulative exposure model (CEM), linear cumulative exposure model (LCEM), tampered random variable model (TRVM), and tampered failure rate model (TFRM). Limitations of the four models are also introduced for better use of the models.     

Smooth Semi‐nonparametric Analysis for Mixture Cure Models and Its Application to Breast Cancer

Mixture cure models are widely used when a proportion of patients are cured. The proportional hazards mixture cure model and the accelerated failure time mixture cure model are the most popular models in practice. Usually the expectation–maximisation (EM) algorithm is applied to both models for parameter estimation. Bootstrap methods are used for variance estimation. In this paper we propose a smooth semi‐nonparametric (SNP) approach in which maximum likelihood is applied directly to mixture cure models for parameter estimation. The variance can be estimated by the inverse of the second derivative of the SNP likelihood. A comprehensive simulation study indicates good performance of the proposed method. We investigate stage effects in breast cancer by applying the proposed method to breast cancer data from the South Carolina Cancer Registry.