Survival estimation through the cumulative hazard with monotone natural cubic splines using convex optimization-the HCNS approach

•A new spline-based method for estimating survival curves that allows for censoring and covariates is presented along with a fully automated software. Estimation and inferences outperform traditional methods such as the Kaplan Meier and the Cox Model in terms of mean integrated squared error as well as obtained coverage of pointwise confidence intervals.•The underlying optimization involves a convex function hence convergence of the minimization algorithm is guaranteed.•Our methods, even though demonstrated through a survival analysis larynx cancer example, can have a broad application beyond survival settings such as diagnostic testing. In survival analysis both the Kaplan-Meier estimate and the Cox model enjoy a broad acceptance. We present an improved spline-based survival estimate and offer a fully automated software for its implementation. We explore the use of natural cubic splines that are constrained to be monotone. Apart from its superiority over the Kaplan Meier estimator our approach overcomes limitations of other known smoothing approaches and can accommodate covariates. Unlike other spline methods, concerns of computational problems and issues of overfitting are resolved since no attempt is made to maximize a likelihood once the Kaplan-Meier estimator is obtained. An application to laryngeal cancer data, a simulation study and illustrations of the broad application of the method and its software are provided. In addition to presenting our approaches, this work contributes to bridging a communication gap between clinicians and statisticians that is often apparent in the medical literature. We employ a two-stage approach: first obtain the stepwise cumulative hazard and then consider a natural cubic spline to smooth its steps under restrictions of monotonicity between any consecutive knots. The underlying region of monotonicity corresponds to a non-linear region that encompasses the full family of monotone third-degree polynomials. We approximate it linearly and reduce the problem to a restricted least squares one under linear restrictions. This ensures convexity. We evaluate our method through simulations against competitive traditional approaches. Our method is compared to the popular Kaplan Meier estimate both in terms of mean squared error and in terms of coverage. Over-fitting is avoided by construction, as our spline attempts to approximate the empirical estimate of the cumulative hazard itself, and is not fitted directly on the data. The p