Iterated residuals and time‐varying covariate effects in Cox regression

Summary. The Cox proportional hazards model, which is widely used for the analysis of treatment and prognostic effects with censored survival data, makes the assumption that the hazard ratio is constant over time. Nonparametric estimators have been developed for an extended model in which the hazard ratio is allowed to change over time. Estimators based on residuals are appealing as they are easy to use and relate in a simple way to the more restricted Cox model estimator. After fitting a Cox model and calculating the residuals, one can obtain a crude estimate of the time‐varying coefficients by adding a smooth of the residuals to the initial (constant) estimate. Treating the crude estimate as the fit, one can re‐estimate the residuals. Iteration leads to consistent estimation of the nonparametric time‐varying coefficients. This approach leads to clear guidelines for residual analysis in applications. The results are illustrated by an analysis of the Medical Research Council's myeloma trials, and by simulation.

[1]  D. Speed,et al.  Melphalan in the Treatment of Myelomatosis , 1964, British medical journal.

[2]  James M. Ortega,et al.  Iterative solution of nonlinear equations in several variables , 2014, Computer science and applied mathematics.

[3]  David R. Cox,et al.  Regression models and life tables (with discussion , 1972 .

[4]  L. Schumaker.  Spline Functions: Basic Theory , 1981 .

[5]  J. Cameron,et al.  Lupus Nephritis , 1959, Clinics in rheumatic diseases.

[6]  S. Geman,et al.  Nonparametric Maximum Likelihood Estimation by the Method of Sieves , 1982 .

[7]  David A. Schoenfeld,et al.  Partial residuals for the proportional hazards regression model , 1982 .

[8]  R Kay,et al.  Goodness of fit methods for the proportional hazards regression model: a review. , 1984, Revue d'epidemiologie et de sante publique.

[9]  S. Pocock,et al.  Regression Models and Non‐Proportional Hazards in the Analysis of Breast Cancer Survival , 1984 .

[10]  I. Maclennan,et al.  Objective evaluation of the role of vincristine in induction and maintenance therapy for myelomatosis. Medical Research Council Working Party on Leukaemia in Adults. , 1985, British Journal of Cancer.

[11]  J. Cuzick,et al.  Results of the mrc myelomatosis trials for patients entered since 1980 , 1988, Hematological oncology.

[12]  F. O’Sullivan.  Nonparametric Estimation of Relative Risk Using Splines and Cross-Validation , 1988 .

[13]  A. Karr,et al.  Nonparametric Survival Analysis with Time-Dependent Covariate Effects: A Penalized Partial Likelihood Approach , 1990 .

[14]  Jack Cuzick,et al.  Long‐term prognostic value of serum β2 microglobulin in myelomatosis , 1990 .

[15]  A. Pettitt,et al.  Investigating Time Dependence in Cox's Proportional Hazards Model , 1990 .

[16]  D. Gamerman.  Dynamic Bayesian Models for Survival Data , 1991 .

[17]  P. Sen,et al.  Time-dependent coefficients in a Cox-type regression model , 1991 .

[18]  I. Maclennan,et al.  Combined chemotherapy with ABCM versus melphalan for treatment of myelomatosis , 1992, The Lancet.

[19]  Robert Gray,et al.  Flexible Methods for Analyzing Survival Data Using Splines, with Applications to Breast Cancer Prognosis , 1992 .

[20]  S. Maclennan,et al.  Screening blood donations for HCV , 1992, The Lancet.

[21]  Niels Keiding,et al.  Statistical Models Based on Counting Processes , 1993 .

[22]  Peter J Green,et al.  Lead discussion on 'Varying coefficient models' by Hastie & Tibshirani , 1993 .

[23]  K R Hess,et al.  Assessing time-by-covariate interactions in proportional hazards regression models using cubic spline functions. , 1994, Statistics in medicine.

[24]  P. Grambsch,et al.  Proportional hazards tests and diagnostics based on weighted residuals , 1994 .

[25]  D. J. Donnell,et al.  Analysis of Additive Dependencies and Concurvities Using Smallest Additive Principal Components , 1994 .

[26]  H C van Houwelingen,et al.  Time-dependent effects of fixed covariates in Cox regression. , 1995, Biometrics.

[27]  James F. Watkins,et al.  Analysing Survival Data from Clinical Trials and Observational Studies. , 1995 .

[28]  Michal Abrahamowicz,et al.  Time-Dependent Hazard Ratio: Modeling and Hypothesis Testing with Application in Lupus Nephritis , 1996 .

[29]  R. Strawderman,et al.  ON CONSISTENCY IN PARAMETER SPACES OF EXPANDING DIMENSION: AN APPLICATION OF THE INVERSE FUNCTION THEOREM , 1996 .

[30]  P. Sasieni,et al.  Evaluation of long-term survival: use of diagnostics and robust estimators with Cox's proportional hazards model. , 1996, Statistics in medicine.

[31]  E. Mammen.  The Bootstrap and Edgeworth Expansion , 1997 .

[32]  D J Sargent,et al.  A Flexible Approach to Time-varying Coefficients in the Cox Regression Setting , 1997, Lifetime data analysis.

[33]  P. Sasieni,et al.  Miscellanea. A note on scaled Schoenfeld residuals for the proportional hazards model , 2001 .

[34]  P. Sasieni,et al.  A note on scaled Schoenfeld residuals for the proportional hazards model , 2001 .