Robust and Efficient Estimation in the Parametric Cox Regression Model under Random Censoring

TL;DR

This paper introduces a robust and efficient estimation method for the parametric Cox regression model using minimum density power divergence, improving robustness against data contamination while maintaining high efficiency.

Contribution

It develops a new robust estimation procedure for the parametric Cox model based on minimum density power divergence, with proven asymptotic properties and practical advantages.

Findings

Robust estimates under data contamination with minimal efficiency loss.

More precise inference than existing likelihood-based methods.

Validated through simulations and real data examples.

Abstract

Cox proportional hazard regression model is a popular tool to analyze the relationship between a censored lifetime variable with other relevant factors. The semi-parametric Cox model is widely used to study different types of data arising from applied disciplines like medical science, biology, reliability studies and many more. A fully parametric version of the Cox regression model, if properly specified, can yield more efficient parameter estimates leading to better insight generation. However, the existing maximum likelihood approach of generating inference under the fully parametric Cox regression model is highly non-robust against data-contamination which restricts its practical usage. In this paper we develop a robust estimation procedure for the parametric Cox regression model based on the minimum density power divergence approach. The proposed minimum density power divergence estimator is seen to produce highly robust estimates under data contamination with only a slight loss in efficiency under pure data. Further, they are always seen to generate more precise inference than the likelihood based estimates under the semi-parametric Cox models or their existing robust versions. We also sketch the derivation of the asymptotic properties of the proposed estimator using the martingale approach and justify their robustness theoretically through the influence function analysis. The practical applicability and usefulness of the proposal are illustrated through simulations and a real data example.