Definition
A counting process N(t) is a right-continuous process that records the cumulative number of events of a point process in [0,t] (nondecreasing, integer-valued, unit jumps). A multivariate counting process N=(N1,…,Nk) collects several such processes. Its intensity process is λi(t)=limh↓0h1E[Ni(t+h)−Ni(t)∣Ft], the instantaneous event rate given the observed history Ft. Counting-process theory (Aalen 1978) makes this the unifying language for survival, event-history, and duration data, with censoring handled naturally through the filtration.
Key Ideas
- Multiplicative intensity model. Aalen's model sets λi(t)=αi(t)Yi(t): an unknown nonnegative rate function αi (the hazard) times an observable predictable at-risk process Yi (e.g., the number of individuals still at risk). It is nonparametric when αi is unrestricted apart from regularity, and it nests censored survival data, finite-state continuous-time Markov chains, and birth-and-death processes.
- Compensator and martingale. The compensator is Λi(t)=∫0tλi(s)ds, and Mi(t)=Ni(t)−Λi(t) is a square-integrable martingale with predictable variation ⟨Mi⟩(t)=∫0tλi(s)ds. "Observed minus expected, given the past" is mean-zero — the property that drives all the inference.
- Nelson-Aalen estimator. The cumulative hazard Bi(t)=∫0tαi(s)ds is estimated by B^i(t)=∫0tYi(s)Ji(s)dNi(s) — add 1/Yi at each observed event while anyone is at risk. It generalizes the empirical cumulative hazard and is the additive-hazard companion of the Kaplan-Meier product-limit (survival) estimator.
- Inference by stochastic integration. Estimators and test statistics are written as stochastic integrals of predictable processes against the martingale Mi; their variances are optional/predictable-variation estimators, and the martingale central limit theorem yields consistency and weak convergence to Gaussian processes — with no i.i.d. or fixed-censoring assumption.
- Two-sample and rank tests. Weighted differences of estimated increments give two-sample tests that generalize the log-rank / Mantel-Haenszel and other rank tests to arbitrary censoring and time-dependent risk sets.
How It Works
- Model the event counts Ni and specify the predictable at-risk process Yi from the observation/censoring scheme.
- Assume the multiplicative intensity λi=αiYi; form the martingale Mi=Ni−∫αiYi.
- Estimate the cumulative hazard by the Nelson-Aalen stochastic integral B^i(t)=∫0t(Ji/Yi)dNi.
- Estimate its variance from the predictable variation and invoke the martingale CLT for confidence bands.
- For group comparisons, integrate a predictable weight against the difference of increment estimates to obtain a rank-type test statistic with a martingale-based variance.
Why It Matters
- Unifies survival analysis. Kaplan-Meier estimation, the log-rank test, and Cox's proportional-hazards partial likelihood all become special cases or close relatives of the counting-process construction — one theory instead of many ad hoc derivations.
- Handles arbitrary censoring. Because censoring enters only through the predictable Yi and the filtration, right-censoring, staggered entry, and time-dependent risk sets are accommodated without special-casing.
- Rigorous, general asymptotics. The martingale CLT replaces i.i.d.-based limit theory, giving valid inference for dependent, time-evolving event data — the probabilistic backbone of duration and event-history modeling.
Open Questions
- Mathematical overhead. The square-integrable-martingale / stochastic-integral machinery is heavy relative to the applied payoff, and downstream methods usually hide it.
- Beyond multiplicative intensity. Additive-hazard, frailty (random-effects), and interval-censored extensions require enlarging the model, and identifiability can become delicate.
- Small-sample behavior. The asymptotics are large-sample; with few events the Gaussian approximation to Nelson-Aalen bands and rank tests can be poor.
Related