B-Spline in the Cox Regression with Application to Cervical Cancer

Recently, Cox proportional hazard (PH) models have played an important role and become increasingly famous in survival analysis. A crucial assumption of the Cox model is the proportional hazards assumption, that is the covariates do not vary over time. One way to check this assumption is to utilize martingale residuals. Martingale residual is an estimate of the overage of events seen in the data but not covered by the model. These residuals are used to examine the best functional form for a given covariate using an assumed Cox model for the remaining covariates. However, one problem that could be occurred when applying martingale residuals is that they tend to be asymmetric and the line does not fall around zero. Hence, in this paper, the main discussion will focus on the use of smoothing martingale residuals, another type of martingale residuals that give a higher rate of flexibility, by using B-spline and the relation to another smoothing technique, locally weighted scatterplot smoothing (LOWESS). An analysis of variables that probably affect the survival rate of patients with cervical cancer is used for illustration.

[1]  D. Harrington,et al.  Regression Splines in the Cox Model with Application to Covariate Effects in Liver Disease , 1990 .

[2]  Rand R. Wilcox,et al.  The Regression Smoother LOWESS: A Confidence Band That Allows Heteroscedasticity And Has Some Specified Simultaneous Probability Coverage , 2017 .

[3]  W. Cleveland,et al.  Locally Weighted Regression: An Approach to Regression Analysis by Local Fitting , 1988 .

[4]  P. Grambsch,et al.  Martingale-based residuals for survival models , 1990 .

[5]  William G. Jacoby Loess: a nonparametric, graphical tool for depicting relationships between variables , 2000 .

[6]  W. Cleveland Robust Locally Weighted Regression and Smoothing Scatterplots , 1979 .

[7]  C. J. Stone,et al.  Additive Regression and Other Nonparametric Models , 1985 .

[8]  Kaili Wang,et al.  Cubic B-Spline Interpolation and Realization , 2011, ICICA.

[9]  J. Klein,et al.  Survival Analysis: Techniques for Censored and Truncated Data , 1997 .

[10]  Michel Artiles,et al.  Cubic B-Spline Curves with Shape Parameter and Their Applications , 2017 .

[11]  N. Breslow Covariance analysis of censored survival data. , 1974, Biometrics.

[12]  I. J. Schoenberg,et al.  On Pólya frequency functions IV: The fundamental spline functions and their limits , 1966 .

[13]  W. Barlow,et al.  Residuals for relative risk regression , 1988 .

[14]  Yongxin Su,et al.  Optimizing Production Scheduling of Steel Plate Hot Rolling for Economic Load Dispatch under Time-of-Use Electricity Pricing , 2017 .

[15]  Shi-Min Hu,et al.  Approximate merging of B-spline curves via knot adjustment and constrained optimization , 2003, Comput. Aided Des..

[16]  Carl de Boor,et al.  A Practical Guide to Splines , 1978, Applied Mathematical Sciences.