Past exam of the mathematics course of the University of Cambridge / 2006 / iii / Paper 44 / 5

Use the usual semiparametric Cox proportional-hazards model with observed covariates $z_i$:The covariates are implicit in the question's regression fit, even though the displayed data notation only lists exit times and event indicators. Assume entry at a common origin, time-constant covariates, distinct observed times and independent censoring. Write $R_i=\{j:x_j\ge x_i\}$ and $w_j(\beta )=e^{{\beta }^T{z}_j}$. Conditional on a single event at $x_i$ in its risk set, the chance that individual $i$ supplies it isThe baseline hazard cancels. Consequently estimate $\beta$ by maximizing Cox partial likelihoodFor example, the score and information of Cox partial likelihood follow by differentiation. Define $S_i^{(0)}={\sum }_{j\in R_i}{w}_j$, $S_i^{(1)}={\sum }_{j\in R_i}{w}_j{z}_j$ and $S_i^{(2)}={\sum }_{j\in R_i}{w}_j{z}_j{z}_j^T$. ThenWhen a finite identifiable maximum exists, iterate $\beta \leftarrow \beta +I(\beta {)}^{-1}U(\beta )$ until convergence. The observed Fisher information also gives the usual asymptotic coefficient uncertainty. Censored subjects enter the relevant risk sets, even though they supply no numerator event factor.

After fitting $\hat{\beta }$, estimate the cumulative hazard of the reference subject by the Breslow estimatorIts jump at an event is the observed event count divided by the total fitted hazard multiplier at risk. The estimate is an integrated baseline hazard, with hazard represented by a measure placing these jumps at observed event times. If a smooth ordinary hazard function is wanted, smooth these increments with a chosen kernel and bandwidth; differentiating the step function as an ordinary function would incorrectly give zero away from the jumps. The fitted individual cumulative hazard is ${\hat{H}}_i(t)=e^{{\hat{\beta }}^T{z}_i}{\hat{H}}_0(t)$ over that individual's follow-up. The invariances below concern this Cox partial likelihood fit and this Breslow estimator; they need not hold for a fully parametric proportional-hazards fit or an arbitrarily smoothed hazard estimator.

Since the observed times are ordered, the final risk set is $R_m=\{m\}$. Its Cox partial likelihood factor is thereforefor either value of the last event indicator. Every earlier risk set and earlier event factor is unaffected by changing $v_m$. Thus the entire partial likelihood as a function of $\beta$ is unchanged, not merely its value at one fitted coefficient. Its maximizing set, its score function and its observed Fisher information are identical. In particular, any consistently chosen finite fitted $\hat{\beta }$ is unchanged. The lone terminal event conveys no relative-hazard information because there is no other individual with whom to compare it.

Use the right-continuous Breslow estimator, so the increase from the value at $x_{m-1}$ to the value at $x_m$ means the interval $(x_{m-1},x_m]$. There is no event between these consecutive exit times, and at $x_m$ the risk set contains only subject $m$. Henceand the fitted individual cumulative hazard increases byIt is the subject's integrated hazard, rather than the reference baseline increment, that equals the event indicator.

For the final unheaded request, a fitted martingale residual at a subject's exit isThe revised observation leaves the last subject at risk until a later event time, still after every other subject has exited. Every earlier risk set is exactly the same as before, and the new final event again contributes the constant factor one. Thus part a gives unchanged coefficients, and every earlier Breslow estimator increment is unchanged. No additional events occur during the extended interval because all other subjects have already left observation.

For $i<m$, the event indicator and fitted cumulative hazard at $x_i$ remain unchanged, hence so does ${\hat{M}}_i$. For the terminal subject, let $C$ denote its old fitted integrated hazard at its censored exit. That old residual is $-C$. Part b shows that the new final event adds one to its fitted integrated hazard, giving new residualTherefore all fitted martingale residuals remain unchanged. Their invariance does not assert that the fitted survivor function at later times or the terminal event record is unchanged: the new late baseline hazard jump is real. This is the terminal-observation invariance of Cox martingale residuals, which relies on a unique final subject and unchanged covariate history.