Confidence bands for cumulative incidence curves under the additive risk model

Fitting additive risk models using auxiliary information

There has been a growing interest in incorporating auxiliary summary information from external studies into the analysis of internal individual-level data. In this paper, we propose an adaptive estimation procedure for an additive risk model to integrate auxiliary subgroup survival information via a penalized method of moments technique. Our approach can accommodate information from heterogeneous data. Parameters to quantify the magnitude of potential incomparability between internal data and external auxiliary information are introduced in our framework while nonzero components of these parameters suggest a violation of the homogeneity assumption. We further develop an efficient computational algorithm to solve the numerical optimization problem by profiling out the nuisance parameters. In an asymptotic sense, our method can be as efficient as if all the incomparable auxiliary information is accurately acknowledged and has been automatically excluded from consideration. The asymptotic normality of the proposed estimator of the regression coefficients is established, with an explicit formula for the asymptotic variance-covariance matrix that can be consistently estimated from the data. Simulation studies show that the proposed method yields a substantial gain in statistical efficiency over the conventional method using the internal data only, and reduces estimation biases when the given auxiliary survival information is incomparable. We illustrate the proposed method with a lung cancer survival study.

Doubly robust estimation of the hazard difference for competing risks data

We consider the conditional treatment effect for competing risks data in observational studies. We derive the efficient score for the treatment effect using modern semiparametric theory, as well as two doubly robust scores with respect to (1) the assumed propensity score for treatment and the censoring model, and (2) the outcome models for the competing risks. An important property regarding the estimators is rate double robustness, in addition to the classical model double robustness. Rate double robustness enables the use of machine learning and nonparametric methods in order to estimate the nuisance parameters, while preserving the root-n $$ n $$ asymptotic normality of the estimated treatment effect for inferential purposes. We study the performance of the estimators using simulation. The estimators are applied to the data from a cohort of Japanese men in Hawaii followed since 1960s in order to study the effect of mid-life drinking behavior on late life cognitive outcomes. The approaches developed in this article are implemented in the R package "HazardDiff".

Competing risks predictions with different time scales under the additive risk model

In the analysis for competing risks data, regression modeling of the cause-specific hazard functions has been usually conducted using the same time scale for all event types. However, when the true time scale is different for each event type, it would be appropriate to specify regression models for the cause-specific hazards on different time scales for different event types. Often, the proportional hazards model has been used for regression modeling of the cause-specific hazard functions. However, the proportionality assumption may not be appropriate in practice. In this article, we consider the additive risk model as an alternative to the proportional hazards model. We propose predictions of the cumulative incidence functions under the cause-specific additive risk models employing different time scales for different event types. We establish the consistency and asymptotic normality of the predicted cumulative incidence functions under the cause-specific additive risk models specified on different time scales using empirical processes and derive consistent variance estimators of the predicted cumulative incidence functions. Through simulation studies, we show that the proposed prediction methods perform well. We illustrate the methods using stage III breast cancer data obtained from the Surveillance, Epidemiology, and End Results (SEER) program of the National Cancer Institute.

EM algorithm for the additive risk mixture cure model with interval-censored data

Interval-censored failure time data arise in a number of fields and many authors have recently paid more attention to their analysis. However, regression analysis of interval-censored data under the additive risk model can be challenging in maximizing the complex likelihood, especially when there exists a non-ignorable cure fraction in the population. For the problem, we develop a sieve maximum likelihood estimation approach based on Bernstein polynomials. To relieve the computational burden, an expectation-maximization algorithm by exploiting a Poisson data augmentation is proposed. Under some mild conditions, the asymptotic properties of the proposed estimator are established. The finite sample performance of the proposed method is evaluated by extensive simulations, and is further illustrated through a real data set from the smoking cessation study.

On a Shape-Invariant Hazard Regression Model with application to an HIV Prevention Study of Mother-to-Child Transmission

In survival analysis, Cox model is widely used for most clinical trial data. Alternatives include the additive hazard model, the accelerated failure time (AFT) model and a more general transformation model. All these models assume that the effects for all covariates are on the same scale. However, it is possible that for different covariates, the effects are on different scales. In this paper, we propose a shape-invariant hazard regression model that allows us to estimate the multiplicative treatment effect with adjustment of covariates that have non-multiplicative effects. We propose moment-based inference procedures for the regression parameters. We also discuss the risk prediction and the goodness of fit test for our proposed model. Numerical studies show good finite sample performance of our proposed estimator. We applied our method to the HIVNET 012 study, a milestone trial of single-dose nevirapine in prevention of mother-to-child transmission of HIV. From the HIVNET 012 data analysis, single-dose nevirapine treatment is shown to improve 18-month infant survival significantly with appropriate adjustment of the maternal CD4 counts and the virus load.

Analysis of Generalized Semiparametric Regression Models for Cumulative Incidence Functions with Missing Covariates

The cumulative incidence function quantifies the probability of failure over time due to a specific cause for competing risks data. The generalized semiparametric regression models for the cumulative incidence functions with missing covariates are investigated. The effects of some covariates are modeled as non-parametric functions of time while others are modeled as parametric functions of time. Different link functions can be selected to add flexibility in modeling the cumulative incidence functions. The estimation procedures based on the direct binomial regression and the inverse probability weighting of complete cases are developed. This approach modifies the full data weighted least squares equations by weighting the contributions of observed members through the inverses of estimated sampling probabilities which depend on the censoring status and the event types among other subject characteristics. The asymptotic properties of the proposed estimators are established. The finite-sample performances of the proposed estimators and their relative efficiencies under different two-phase sampling designs are examined in simulations. The methods are applied to analyze data from the RV144 vaccine efficacy trial to investigate the associations of immune response biomarkers with the cumulative incidence of HIV-1 infection.

Smoothed Rank Regression for the Accelerated Failure Time Competing Risks Model with Missing Cause of Failure

This paper examines the accelerated failure time competing risks model with missing cause of failure using the monotone class rank-based estimating equations approach. We handle the non-smoothness of the rank-based estimating equations using a kernel smoothed estimation method, and estimate the unknown selection probability and the conditional expectation by non-parametric techniques. Under this setup, we propose three methods for estimating the unknown regression parameters based on 1) inverse probability weighting, 2) estimating equations imputation and 3) augmented inverse probability weighting. We also obtain the associated asymptotic theories of the proposed estimators and investigate the estimators' small sample behaviour in a simulation study. A direct plug-in method is suggested for estimating the asymptotic variances of the proposed estimators. A real data application based on a HIV vaccine efficacy trial study is considered.

Instrumental variable with competing risk model

In this paper, we discuss causal inference on the efficacy of a treatment or medication on a time-to-event outcome with competing risks. Although the treatment group can be randomized, there can be confoundings between the compliance and the outcome. Unmeasured confoundings may exist even after adjustment for measured covariates. Instrumental variable methods are commonly used to yield consistent estimations of causal parameters in the presence of unmeasured confoundings. On the basis of a semiparametric additive hazard model for the subdistribution hazard, we propose an instrumental variable estimator to yield consistent estimation of efficacy in the presence of unmeasured confoundings for competing risk settings. We derived the asymptotic properties for the proposed estimator. The estimator is shown to be well performed under finite sample size according to simulation results. We applied our method to a real transplant data example and showed that the unmeasured confoundings lead to significant bias in the estimation of the effect (about 50% attenuated). Copyright © 2017 John Wiley & Sons, Ltd.

On the choice of time scales in competing risks predictions

In the standard analysis of competing risks data, proportional hazards models are fit to the cause-specific hazard functions for all causes on the same time scale. These regression analyses are the foundation for predictions of cause-specific cumulative incidence functions based on combining the estimated cause-specific hazard functions. However, in predictions arising from disease registries, where only subjects with disease enter the database, disease-related mortality may be more naturally modeled on the time since diagnosis time scale while death from other causes may be more naturally modeled on the age time scale. The single time scale methodology may be biased if an incorrect time scale is employed for one of the causes and an alternative methodology is not available. We propose inferences for the cumulative incidence function in which regression models for the cause-specific hazard functions may be specified on different time scales. Using the disease registry data, the analysis of other cause mortality on the age scale requires left truncating the event time at the age of disease diagnosis, complicating the analysis. In addition, standard Martingale theory is not applicable when combining regression models on different time scales. We establish that the covariate conditional predictions are consistent and asymptotically normal using empirical process techniques and propose consistent variance estimators for constructing confidence intervals. Simulation studies show that the proposed two time scales method performs well, outperforming the single time-scale predictions when the time scale is misspecified. The methods are illustrated with stage III colon cancer data obtained from the Surveillance, Epidemiology, and End Results program of National Cancer Institute.

Simultaneous confidence bands for Cox regression from semiparametric random censorship

Cox regression is combined with semiparametric random censorship models to construct simultaneous confidence bands (SCBs) for subject-specific survival curves. Simulation results are presented to compare the performance of the proposed SCBs with the SCBs that are based only on standard Cox. The new SCBs provide correct empirical coverage and are more informative. The proposed SCBs are illustrated with two real examples. An extension to handle missing censoring indicators is also outlined.

A Proportional Hazards Regression Model for the Sub-distribution with Covariates Adjusted Censoring Weight for Competing Risks Data

With competing risks data, one often needs to assess the treatment and covariate effects on the cumulative incidence function. Fine and Gray proposed a proportional hazards regression model for the subdistribution of a competing risk with the assumption that the censoring distribution and the covariates are independent. Covariate-dependent censoring sometimes occurs in medical studies. In this paper, we study the proportional hazards regression model for the subdistribution of a competing risk with proper adjustments for covariate-dependent censoring. We consider a covariate-adjusted weight function by fitting the Cox model for the censoring distribution and using the predictive probability for each individual. Our simulation study shows that the covariate-adjusted weight estimator is basically unbiased when the censoring time depends on the covariates, and the covariate-adjusted weight approach works well for the variance estimator as well. We illustrate our methods with bone marrow transplant data from the Center for International Blood and Marrow Transplant Research (CIBMTR). Here cancer relapse and death in complete remission are two competing risks.

Semiparametric methods for center effect measures based on the ratio of survival functions

The survival function is often of chief interest in epidemiologic studies of time to an event. We develop methods for evaluating center-specific survival outcomes through a ratio of survival functions. The proposed method assumes a center-stratified additive hazards model, which provides a convenient framework for our purposes. Under the proposed methods, the center effects measure is cast as the ratio of subject-specific survival functions under two scenarios: the scenario in which the subject is treated at center [Formula: see text]; and that wherein the subject is treated at a hypothetical center with survival function equal to the population average. The proposed measure reduces to the ratio of baseline survival functions, but is invariant to the choice of baseline covariate level. We derive the asymptotic properties of the proposed estimators, and assess finite-sample characteristics through simulation. The proposed methods are applied to national kidney transplant data.

Population-based absolute risk estimation with survey data

Absolute risk is the probability that a cause-specific event occurs in a given time interval in the presence of competing events. We present methods to estimate population-based absolute risk from a complex survey cohort that can accommodate multiple exposure-specific competing risks. The hazard function for each event type consists of an individualized relative risk multiplied by a baseline hazard function, which is modeled nonparametrically or parametrically with a piecewise exponential model. An influence method is used to derive a Taylor-linearized variance estimate for the absolute risk estimates. We introduce novel measures of the cause-specific influences that can guide modeling choices for the competing event components of the model. To illustrate our methodology, we build and validate cause-specific absolute risk models for cardiovascular and cancer deaths using data from the National Health and Nutrition Examination Survey. Our applications demonstrate the usefulness of survey-based risk prediction models for predicting health outcomes and quantifying the potential impact of disease prevention programs at the population level.

Modeling cumulative incidence function for competing risks data

A frequent occurrence in medical research is that a patient is subject to different causes of failure, where each cause is known as a competing risk. The cumulative incidence curve is a proper summary curve, showing the cumulative failure rates over time due to a particular cause. A common question in medical research is to assess the covariate effects on a cumulative incidence function. The standard approach is to construct regression models for all cause-specific hazard rate functions and then model a covariate-adjusted cumulative incidence curve as a function of all cause-specific hazards for a given set of covariates. New methods have been proposed in recent years, emphasizing direct assessment of covariate effects on cumulative incidence function. Fine and Gray proposed modeling the effects of covariates on a subdistribution hazard function. A different approach is to directly model a covariate-adjusted cumulative incidence function, including a pseudovalue approach by Andersen and Klein and a direct binomial regression by Scheike, Zhang and Gerds. In this paper, we review the standard and new regression methods for modeling a cumulative incidence function, and give the sources of computer packages/programs that implement these regression models. A real bone marrow transplant data set is analyzed to illustrate various regression methods.

Proportional hazards model for competing risks data with missing cause of failure

We consider the semiparametric proportional hazards model for the cause-specific hazard function in analysis of competing risks data with missing cause of failure. The inverse probability weighted equation and augmented inverse probability weighted equation are proposed for estimating the regression parameters in the model, and their theoretical properties are established for inference. Simulation studies demonstrate that the augmented inverse probability weighted estimator is doubly robust and the proposed method is appropriate for practical use. The simulations also compare the proposed estimators with the multiple imputation estimator of Lu and Tsiatis (2001). The application of the proposed method is illustrated using data from a bone marrow transplant study.

A proportional hazards regression model for the subdistribution with right-censored and left-truncated competing risks data

With competing risks failure time data, one often needs to assess the covariate effects on the cumulative incidence probabilities. Fine and Gray proposed a proportional hazards regression model to directly model the subdistribution of a competing risk. They developed the estimating procedure for right-censored competing risks data, based on the inverse probability of censoring weighting. Right-censored and left-truncated competing risks data sometimes occur in biomedical researches. In this paper, we study the proportional hazards regression model for the subdistribution of a competing risk with right-censored and left-truncated data. We adopt a new weighting technique to estimate the parameters in this model. We have derived the large sample properties of the proposed estimators. To illustrate the application of the new method, we analyze the failure time data for children with acute leukemia. In this example, the failure times for children who had bone marrow transplants were left truncated.

Additive mixed effect model for clustered failure time data

We propose an additive mixed effect model to analyze clustered failure time data. The proposed model assumes an additive structure and includes a random effect as an additional component. Our model imitates the commonly used mixed effect models in repeated measurement analysis but under the context of hazards regression; our model can also be considered as a parallel development of the gamma-frailty model in additive model structures. We develop estimating equations for parameter estimation and propose a way of assessing the distribution of the latent random effect in the presence of large clusters. We establish the asymptotic properties of the proposed estimator. The small sample performance of our method is demonstrated via a large number of simulation studies. Finally, we apply the proposed model to analyze data from a diabetic study and a treatment trial for congestive heart failure.

SAS macros for estimation of direct adjusted cumulative incidence curves under proportional subdistribution hazards models

The cumulative incidence function is commonly reported in studies with competing risks. The aim of this paper is to compute the treatment-specific cumulative incidence functions, adjusting for potentially imbalanced prognostic factors among treatment groups. The underlying regression model considered in this study is the proportional hazards model for a subdistribution function [1]. We propose estimating the direct adjusted cumulative incidences for each treatment using the pooled samples as the reference population. We develop two SAS macros for estimating the direct adjusted cumulative incidence function for each treatment based on two regression models. One model assumes the constant subdistribution hazard ratios between the treatments and the alternative model allows each treatment to have its own baseline subdistribution hazard function. The macros compute the standard errors for the direct adjusted cumulative incidence estimates, as well as the standard errors for the differences of adjusted cumulative incidence functions between any two treatments. Based on the macros' output, one can assess treatment effects at predetermined time points. A real bone marrow transplant data example illustrates the practical utility of the SAS macros.

Additive transformation models for clustered failure time data

We propose a class of additive transformation risk models for clustered failure time data. Our models are motivated by the usual additive risk model for independent failure times incorporating a frailty with mean one and constant variability which is a natural generalization of the additive risk model from univariate failure time to multivariate failure time. An estimating equation approach based on the marginal hazards function is proposed. Under the assumption that cluster sizes are completely random, we show the resulting estimators of the regression coefficients are consistent and asymptotically normal. We also provide goodness-of-fit test statistics for choosing the transformation. Simulation studies and real data analysis are conducted to examine the finite-sample performance of our estimators.

Assessing cumulative incidence functions under the semiparametric additive risk model

In analyzing competing risks data, a quantity of considerable interest is the cumulative incidence function. Often, the effect of covariates on the cumulative incidence function is modeled via the proportional hazards model for the cause-specific hazard function. As the proportionality assumption may be too restrictive in practice, we consider an alternative more flexible semiparametric additive hazards model of (Biometrika 1994; 81:501-514) for the cause-specific hazard. This model specifies the effect of covariates on the cause-specific hazard to be additive as well as allows the effect of some covariates to be fixed and that of others to be time varying. We present an approach for constructing confidence intervals as well as confidence bands for the cause-specific cumulative incidence function of subjects with given values of the covariates. Furthermore, we also present an approach for constructing confidence intervals and confidence bands for comparing two cumulative incidence functions given values of the covariates. The finite sample property of the proposed estimators is investigated through simulations. We conclude our paper with an analysis of the well-known malignant melanoma data using our method.

Flexible competing risks regression modeling and goodness-of-fit

In this paper we consider different approaches for estimation and assessment of covariate effects for the cumulative incidence curve in the competing risks model. The classic approach is to model all cause-specific hazards and then estimate the cumulative incidence curve based on these cause-specific hazards. Another recent approach is to directly model the cumulative incidence by a proportional model (Fine and Gray, J Am Stat Assoc 94:496-509, 1999), and then obtain direct estimates of how covariates influences the cumulative incidence curve. We consider a simple and flexible class of regression models that is easy to fit and contains the Fine-Gray model as a special case. One advantage of this approach is that our regression modeling allows for non-proportional hazards. This leads to a new simple goodness-of-fit procedure for the proportional subdistribution hazards assumption that is very easy to use. The test is constructive in the sense that it shows exactly where non-proportionality is present. We illustrate our methods to a bone marrow transplant data from the Center for International Blood and Marrow Transplant Research (CIBMTR). Through this data example we demonstrate the use of the flexible regression models to analyze competing risks data when non-proportionality is present in the data.

Inference for outcome probabilities in multi-state models

In bone marrow transplantation studies, patients are followed over time and a number of events may be observed. These include both ultimate events like death and relapse and transient events like graft versus host disease and graft recovery. Such studies, therefore, lend themselves for using an analytic approach based on multi-state models. We will give a review of such methods with emphasis on regression models for both transition intensities and transition- and state occupation probabilities. Both semi-parametric models, like the Cox regression model, and parametric models based on piecewise constant intensities will be discussed.

Design of the Dialysis Access Consortium (DAC) Aggrenox Prevention Of Access Stenosis Trial

BACKGROUND

Surgically created arteriovenous (AV) grafts are the most common type of hemodialysis vascular access in the United States, but fail frequently due to the development of venous stenosis. The Dialysis Access Consortium (DAC) Aggrenox Prevention of Access Stenosis Trial tests the hypothesis that Aggrenox (containing dipyridamole and aspirin) can prevent stenosis and prolong survival of arteriovenous grafts.

METHODS

This is a multicenter, randomized, double-blind, placebo-controlled trial that will enroll 1056 subjects over four years with one-half year follow-up. Subjects undergoing placement of a new AV graft for hemodialysis are randomized to treatment with Aggrenox or placebo immediately following access surgery. The primary outcome is primary unassisted patency defined as the time from access placement until thrombosis or an access procedure carried out to maintain or restore patency. The major secondary outcome is cumulative access patency. Monthly access flow monitoring is incorporated in the study design to enhance detection of a hemodynamically significant access stenosis before it leads to thrombosis.

RESULTS

This paper describes the key issues in trial design, broadly including: 1) ethical issues surrounding the study of a clinical procedure that, although common, is no longer the clinical intervention of choice; 2) acceptable risk (bleeding) from the primary intervention; 3) inclusion of subjects already receiving a portion of the study intervention; 4) inclusion of subjects with incident rather than prevalent qualifying clinical conditions; 5) timing of the study intervention to balance safety and efficacy concerns; and 6) the selection of primary and secondary study endpoints.

CONCLUSIONS

This is the first, large, multicenter trial evaluating a pharmacologic approach to prevent AV graft stenosis and failure, an important and costly problem in this patient population. Numerous design issues were addressed in implementing the trial and these will form a roadmap for future trials in this area.

Modelling competing risks in cancer studies

Competing risks arise commonly in the analysis of cancer studies. Most common are the competing risks of relapse and death in remission. These two risks are the primary reason that patients fail treatment. In most medical papers the effects of covariates on the three outcomes (relapse, death in remission and treatment failure) are model by distinct proportional hazards regression models. Since the hazards of relapse and death in remission must add to that of treatment failure, we argue that this model leads to internal inconsistencies. We argue that additive models for either the hazard rates or the cumulative incidence functions are more natural and that these models properly partition the effect of a covariate on treatment failure into its component parts. We illustrate the use and interpretation of additive models for the hazard rate or for the cumulative incidence function using data from a study of the efficacy of two preparative regimes for hematopoietic stem cell transplantation.

Extensions and applications of the Cox-Aalen survival model

Cox's regression model is the standard regression tool for survival analysis in most applications. Often, however, the model only provides a rough summary of the effect of some covariates. Therefore, if the aim is to give a detailed description of covariate effects and to consequently calculate predicted probabilities, more flexible models are needed. In another article, Scheike and Zhang (2002, Scandinavian Journal of Statistics 29, 75-88), we suggested a flexible extension of Cox's regression model, which aimed at extending the Cox model only for those covariates where additional flexibility are needed. One important advantage of the suggested approach is that even though covariates are allowed a nonparametric effect, the hassle and difficulty of finding smoothing parameters are not needed. We show how the extended model also leads to simple formulae for predicted probabilities and their standard errors, for example, in the competing risk framework.

Two simulation methods for constructing confidence bands under the additive risk model

With right-censored failure time data, we propose procedures to construct simultaneous confidence bands for the survival curve of a given subject under the additive risk model. The distribution of the subject-specific cumulative hazard function can be approximated with a zero-mean Gaussian process, which can be generated through simulation methods. We propose two different simulation schemes to obtain the distribution of the supremum of the cumulative hazard process over the entire time range. The two simulation approaches are asymptotically equivalent, whereas the numerical forms are different. We construct two types of confidence bands, namely the equal precision band and the Hall-Wellner type band, through choosing suitable weight functions. Monte Carlo simulation studies show that both of the proposed confidence bands are appropriate for finite sample sizes. We illustrate the new proposal with a real example.

Omnibus tests for comparison of competing risks with adjustment for covariate effects

This article develops omnibus tests for comparing cause-specific hazard rates and cumulative incidence functions at specified covariate levels. Confidence bands for the difference and the ratio of two conditional cumulative incidence functions are also constructed. The omnibus test is formulated in terms of a test process given by a weighted difference of estimates of cumulative cause-specific hazard rates under Cox proportional hazards models. A simulation procedure is devised for sampling from the null distribution of the test process, leading to graphical and numerical technques for detecting significant differences in the risks. The approach is applied to a cohort study of type-specific HIV infection rates.