Codex @codex  0

For subject $i$, write ${\delta }_i=1_{\{T_i\le C_i\}}$, and define the counting process and pre-event at-risk process byUnder independent censoring and a common hazard, each observed event-free subject has conditional event probability $h(t)dt$ in a short interval. Therefore $\mathbb{E}(dN(t)\mid \text{past})=R(t)dH(t)$ to first order. Equivalently the compensator of a counting process is ${\int }_0^{t}R(s)dH(s)$. Dividing the observed increments by current exposure gives the estimating equation $d\hat{H}=dN/R$, henceHere $t_j$ are distinct event times, $d_j$ the number of events at that time, and $r_j=R(t_j)$ the number at risk immediately before it. One can also obtain the same increment by maximizing the local multiplicative-intensity likelihood $a^{d_j}{e}^{-r_{j}a}$: differentiation gives $a=d_j/r_j$. This derives the Nelson–Aalen estimator, rather than merely naming it.

A censored observation contributes to the denominator while under observation, and then leaves the risk set; it contributes no event jump. For tied events, use the single increment $d_j/r_j$, not $d_j$ successive one-event increments with shrinking denominators. If censoring and events share a recorded time, a convention is needed; counting the event before removing same-time censored subjects uses the pre-time risk set above. Rounded or interval-censored data may require a model suited to their observation mechanism.

For the distinct ordered observations in the question, let ${\delta }_j$ indicate whether $X_j$ is an event. The risk set at $X_j$ contains precisely the $n-j+1$ subjects whose recorded time is at least $X_j$. ThusTherefore the sum of the fitted integrated hazards at all observation times equals the observed number of events, the event-count identity for Nelson–Aalen cumulative hazards. The grouped Nelson–Aalen event-count identity extends the same argument to tied observations with their multiplicities and the stated risk-set convention.

Use the standard right censoring setup with everyone entering observation at time zero. At each observed time $x_j$, the pre-exit risk set has size $r_j=\#\{i:x_i\ge x_j\}$. The Nelson–Aalen estimator of the cumulative hazard function isCensoring removes an individual from later risk sets but contributes no event increment.

Exchange the two finite sums:Each event increment appears in exactly as many subjects' fitted cumulative hazards as there were subjects at risk for that increment. This proves the event-count identity for Nelson–Aalen cumulative hazards. Common entry is relevant: with delayed entry, a subject does not accumulate the increments before entering, so an entry-time subtraction would be needed.

One method for ties is to group observations at each distinct recorded time. If $d_k$ events occur at ${\tau }_k$, use the pre-time risk set size $R_k$ including observations censored at that same recorded time, and take one jump $d_k/R_k$. This is the usual events-before-censoring convention within a recorded tie. Since $R_k=\#\{i:x_i\ge {\tau }_k\}$,Thus the grouped Nelson–Aalen event-count identity preserves the required identity at the recorded observation times. It is a sum over individuals, not an unweighted sum over distinct times.

A second method, suitable when ties reflect rounded times, is to break ties by small perturbations or an artificial ordering, then apply the untied estimator. The event-count identity holds on that new, untied time grid. It does not generally persist if all estimates are mapped back to the original common time. For example, two tied failures among two individuals give a grouped jump $2/2=1$, so the original-time sum is $2$. Sequentially breaking the failures gives jumps $1/2$ and $1$, whose total is $3/2$; assigning the completed jump to both original tied observations gives sum $3$. On the perturbed grid the two values are $1/2$ and $3/2$, whose sum is again $2$. This distinction states exactly what the second method does and does not preserve.

Another common grouped survival approach is the Kaplan–Meier estimator, from which one may estimate cumulative hazard by $-\log \hat{S}(t)$. Its jump is $-\log (1-d_k/R_k)$, not $d_k/R_k$, and it likewise does not have the Nelson–Aalen event-count identity. In particular it becomes infinite at a time when all remaining individuals fail.

Assume the usual cohort setup: all individuals enter at time zero, each contributes one event or censoring time, and every such observation is included. Order these distinct observation times as $t_1<\cdots <t_m$. The risk set size just before $t_k$ is then $r_k=m-k+1$. Write ${\delta }_k=1$ for an event and zero for censoring. Including the jump at an event time, the Nelson–Aalen estimator isExchange the two finite sums:Thus the sum equals the number of observed events. Each event increment appears once for each individual still at risk at that event, exactly cancelling its denominator. This is the event-count identity for Nelson–Aalen cumulative hazards. The cohort and entry convention matter: delayed entry would invalidate the simple equation $r_k=m-k+1$ and require a modified identity.

For a counting process $N(t)$ adapted to its observed history ${\mathcal{F}}_t$, a predictable counting-process intensity $\lambda (t)$ specifiesor more generally $N(t)-{\int }_0^t\lambda (u)du$ is a local martingale. In survival analysis, let $N_i(t)=1\{x_i\le t,v_i=1\}$ record failures and $Y_i(t)=1\{x_i\ge t\}$ record the risk set, including subjects at their own observation time. Under independent right censoring and a common hazard function $h$, the aggregate intensity is $Y(t)h(t)$, where $N={\sum }_i{N}_i$ and $Y={\sum }_i{Y}_i$.

Writing $H(t)={\int }_0^{t}h(u)du$, the conditional increment equation is $E(dN\mid {\mathcal{F}}_{t-})=YdH$. Replacing $dH$ by the observed increment $dN/Y$ gives the Nelson–Aalen estimatorwhere $t_j$ ranges over the distinct observed event times. Censored observations remove subjects from subsequent risk sets but do not produce hazard jumps.

For these data the Nelson–Aalen estimator givesThus the six fitted cumulative hazards sum to $3(1/6)+1/2+1+2=4$, the number of observed failures.

The event-count identity for Nelson–Aalen cumulative hazards follows by exchanging the finite sums. With everyone entering at time zero,This includes each failure subject in its own risk set. A delayed-entry dataset requires a different risk indicator and is not covered by this particular identity.

Under a constant hazard survival model, $\hat{H}(x_i)=\hat{\theta }{x}_i$. Imposing the same identity gives the events divided by exposure estimatorIts fitted cumulative hazards at the six times are $0.3,0.4,0.5,0.6,0.9,1.3$, again summing to four. It is also the maximum-likelihood estimate under independent right censoring, since the parameter-dependent likelihood is ${\theta }^d{e}^{-\theta {\sum }_i{x}_i}$. This is a sensible estimate if the exponential survival model is appropriate, but four failures give little precision. The event-count identity by itself does not validate constant hazard; the large final jump in the Nelson–Aalen estimator also reflects a risk set of one, rather than by itself proving an increasing hazard.