Past exam of the mathematics course of the University of Cambridge / 2016 / iii / Paper 207 / 5

For a continuous proper event time with cumulative hazard function $H$ and survivor function $F(t)=e^{-H(t)}$, use the assumed inverse of $H$. For $u\ge 0$,Therefore the integrated-hazard transformation is unit exponential:This cumulative hazard probability transformation is the basis for the following survival-model residual diagnostics.

The transformed event time is unit exponential, and the transformed censoring time remains a censoring time because $H$ is increasing. Under independent censoring and a correctly fitted model, the Kaplan–Meier estimator of the survivor function for the Cox–Snell residuals $\hat{H}(x_i)$ should approximate $e^{-u}$.

Thus a logarithmic survivor plot should follow a straight line through the origin with slope minus one:The empirical plot is a step function around that line. Systematic curvature suggests model misspecification; the late part is less stable when the risk set is small. The retained event indicators $v_i$ must be used rather than treating censored residuals as complete event times.

Put $U=H(T)$ and $V=H(C)$, so $Y=\min (U,V)$ and $0\le \mathbb{E}Y\le \mathbb{E}U=1$. Its expectation is not generally one and depends on censoring. Under independent censoring, condition on $V=v$: the exponential distribution givesWriting $\Delta =1_{T\le C}$ for the observed event indicator, it follows that $\mathbb{E}Y=\mathbb{E}\Delta$. The exponential mean imputation under independent censoring therefore gives a modified Cox–Snell residual with known mean:On censoring, the added one is the expected remaining transformed lifetime, by memorylessness of the exponential distribution. This gives a mean-one variable, not necessarily another exponential variable. Equivalently, $M=\Delta -Y=1-Y^{\ast }$ is a martingale residual with mean zero. Independent censoring, or its appropriate conditional version, is essential for these identities; arbitrary informative censoring does not justify the correction.

Fit a model using the relevant covariates, calculate ${\hat{M}}_i=v_i-{\hat{H}}_i(x_i)$, and plot these martingale residuals against a candidate explanatory variable or against an included variable's value. A smooth systematic trend away from zero can indicate an omitted effect or an unsuitable functional form, such as a nonlinear age effect. After fitting an appropriate effect, the residual trend should diminish. The equivalent plot of ${\hat{Y}}_i^{\ast }$ has mean-one reference rather than zero.