Continuous threshold models with two-way interactions in survival analysis

Proportional hazards model with the biomarker–treatment interaction plays an important role in the survival analysis of the subset treatment effect. A threshold parameter for a continuous biomarker variable defines the subset of patients who can benefit or lose from a certain new treatment. In this article, we focus on a continuous threshold effect using the rectified linear unit and propose a gradient descent method to obtain the maximum likelihood estimation of the regression coefficients and the threshold parameter simultaneously. Under certain regularity conditions, we prove the consistency, asymptotic normality and provide a robust estimate of the covariance matrix when the model is misspecified. To illustrate the finite sample properties of the proposed methods, we simulate data to evaluate the empirical biases, the standard errors and the coverage probabilities for both the correctly specified models and misspecified models. The proposed continuous threshold model is applied to a prostate cancer data with serum prostatic acid phosphatase as a biomarker.