A Cox–Snell residual is the fitted cumulative hazard function evaluated at an observed event or censoring time: . Its reference probability distribution comes from the probability integral transform. If the true event time has a proper continuous survival function , then satisfiesFor example, an exponential distribution with rate has , so rescaling its event time by produces a unit-rate exponential distribution.
With right censoring, and the transformed observation is together with its original event indicator . Under independent censoring, a Kaplan–Meier estimator applied to these transformed observations should approximately follow ; its logarithm should approximately follow a straight line of slope minus one. Censored residuals must remain censored in this comparison.
The exact exponential identity uses the true model, whereas is estimated from data: fitted residuals are approximate diagnostics and need not be independent exact exponentials. Atoms in a discrete event-time probability distribution, or a cure fraction with positive probability of no finite event, invalidate the unmodified proper-continuous reference law. Such models require diagnostics respecting their jumps or surviving mass.
New to topics? Read the docs here!