Hanqing Li Hanqing Li 李汉青


Survival analysis begins with an event time, but its distinctive mathematics comes from observing that time through censoring, truncation, and evolving risk sets. This page fixes the distributional and observational objects used throughout the note. Measure-theoretic probability and likelihood theory are developed in the notes on probability theory and statistical inference.

Event-time distributions

Definition. Let \(T:\Omega\to[0,\infty]\) be an event time. Its distribution and survival functions are \[ F(t)=P(T\le t),\qquad S(t)=P(T>t)=1-F(t),\qquad t\ge0. \] The right endpoint of observable survival is \(\tau_T=\sup\{t:S(t)>0\}\). If \(F\) is absolutely continuous with density \(f\), its hazard rate is \[ \lambda(t)=\lim_{h\downarrow0} \frac{P(t\le T<t+h\mid T\ge t)}{h} =\frac{f(t)}{S(t-)} \] whenever the limit exists. The cumulative hazard is the Lebesgue-Stieltjes integral \[ \Lambda(t)=\int_{(0,t]}\frac{\mathrm{d}F(u)}{S(u-)}. \]

The measure definition of \(\Lambda\) covers discrete and mixed event-time laws. The familiar integral of a hazard rate represents only its absolutely continuous part.

Theorem (product-integral representation). For every \(t<\tau_T\), \[ S(t)=\prod_{(0,t]}\{1-\mathrm{d}\Lambda(u)\} =\exp\{-\Lambda^c(t)\}\prod_{0<u\le t}\{1-\Delta\Lambda(u)\}, \] where \(\Lambda^c\) is the continuous part and \(\Delta\Lambda(u)=P(T=u\mid T\ge u)\). In particular, if \(F\) is absolutely continuous, \[ S(t)=e^{-\Lambda(t)},\qquad \Lambda(t)=\int_0^t\lambda(u)\,\mathrm{d}u,\qquad f(t)=\lambda(t)S(t). \]

Proof. At a jump time \(u\), \[ \frac{S(u)}{S(u-)} =1-\frac{F(u)-F(u-)}{S(u-)} =1-\Delta\Lambda(u). \] On an interval containing no jumps, \(\mathrm{d}F=-\mathrm{d}S\), hence \(\mathrm{d}\Lambda=-\mathrm{d}S/S\) and \(S(b)/S(a)=\exp\{-[\Lambda^c(b)-\Lambda^c(a)]\}\). Multiplying the continuous changes and all jump ratios gives the product integral. The absolutely continuous identities follow by differentiation.

Example (Weibull law). If \[ S(t)=\exp\{-(t/\eta)^\alpha\},\qquad \alpha,\eta>0, \] then \(\Lambda(t)=(t/\eta)^\alpha\) and \(\lambda(t)=\alpha t^{\alpha-1}/\eta^\alpha\). Thus \(\alpha<1\), \(\alpha=1\), and \(\alpha>1\) give decreasing, constant, and increasing hazards. A survival distribution can therefore be easier to interpret through its local hazard than through its density.

Observation schemes

Definition. The principal incomplete-data schemes are:

  • right censoring: observe \(X=T\wedge C\) and \(\Delta=1\{T\le C\}\);
  • left censoring: observe whether \(T\le C\), possibly together with \(T\) when \(T>C\);
  • interval censoring: observe \(L<T\le R\), where \(R\) may equal \(+\infty\);
  • left truncation: include a subject only when \(T\ge A\), and then observe the subject from entry time \(A\);
  • administrative censoring: \(C\) is fixed by the study calendar rather than generated by the subject.

Censoring makes an included observation incomplete. Truncation changes which units enter the sample.

Definition (independent censoring). With baseline covariates \(Z\), elementary conditional independent censoring is \[ T\perp C\mid Z. \] In a filtration \((\mathcal F_t)\), the operational condition is that the event counting process has the same \(\mathcal F_t\)-intensity after the censoring history is included. This latter condition permits some censoring mechanisms that are not marginally independent of \(T\), provided all common predictors are observed.

Proposition (nonidentifiability without a censoring condition). The law of \((X,\Delta)\) does not in general identify the marginal law of \(T\).

Proof. Suppose \(X\equiv1\) and \(\Delta\equiv0\). This observation law is compatible with every distribution supported on \((1,\infty)\): set \(C\equiv1\) and choose any such \(T\). Hence distinct event-time laws produce exactly the same observed data. More generally, dependence between \(T\) and \(C\) can alter the unobserved tail while preserving the subdistributions \(P(X\le t,\Delta=j)\).

Likelihood under right censoring

Theorem (factorization under conditional independence). Suppose \(T\perp C\mid Z\), both conditional laws are absolutely continuous, and \(O=(X,\Delta,Z)\). Conditional on \(Z=z\), \[ p_O(x,\delta\mid z) =\{f_T(x\mid z)S_C(x\mid z)\}^{\delta} \{f_C(x\mid z)S_T(x\mid z)\}^{1-\delta}. \] Consequently, if the event and censoring laws have variation-independent parameters, the event-time likelihood is \[ L_T =\prod_{i=1}^n f_T(x_i\mid z_i)^{\delta_i}S_T(x_i\mid z_i)^{1-\delta_i} =\prod_{i=1}^n \lambda_T(x_i\mid z_i)^{\delta_i} e^{-\Lambda_T(x_i\mid z_i)} \] in the absolutely continuous case.

Proof. On \(\{\Delta=1\}\), an event occurs at \(x\) while censoring exceeds \(x\), giving \(f_T(x\mid z)S_C(x\mid z)\). On \(\{\Delta=0\}\), censoring occurs at \(x\) while the event exceeds \(x\), giving \(f_C(x\mid z)S_T(x\mid z)\). Multiplication over subjects separates the factors involving the event law from those involving the censoring law. Finally \(f_T=\lambda_TS_T\) and \(S_T=e^{-\Lambda_T}\).

Example (exponential likelihood). If \(T_i\sim\operatorname{Exp}(\lambda)\), the event-time log-likelihood is \[ \ell(\lambda)=D\log\lambda-\lambda\sum_{i=1}^nX_i, \qquad D=\sum_{i=1}^n\Delta_i. \] When \(D>0\), \(\widehat\lambda=D/\sum_iX_i\): the numerator is the number of events and the denominator is total exposure time. This event-over-exposure form reappears in multiplicative intensity models.

Delayed entry and conditional sampling

Proposition. Suppose a subject enters at \(A\) and is sampled only if \(T\ge A\). Conditional on \(A=a\) and inclusion, the event-time density is \[ t\longmapsto \frac{f_T(t)}{S_T(a)}1\{t\ge a\}. \] With right censoring at \(C\), the event-law contribution is therefore \[ \frac{f_T(x)^\Delta S_T(x)^{1-\Delta}}{S_T(a)}, \qquad x\ge a. \]

Proof. This is ordinary conditional probability: \[ P(T\in\mathrm{d}t\mid T\ge a) =\frac{1\{t\ge a\}P(T\in\mathrm{d}t)}{S_T(a)}. \] The censored contribution follows by applying the right-censoring factorization to this conditional law.

In process notation the observable risk indicator is \[ Y_i(t)=1\{A_i\le t\le X_i\}. \] Every later risk-set estimator remains valid with this \(Y_i\), provided entry is conditionally independent in the corresponding filtration. Ignoring \(A_i\) creates length-biased risk sets: individuals who fail before entry could never have appeared.

Support and estimability

Remark. Even under independent censoring, \(S_T(t)\) is nonparametrically estimable only where subjects can remain uncensored and at risk. A convenient endpoint is \[ \tau=\sup\{t:P(X\ge t)>0\}. \] Statements of consistency and Gaussian limits below are therefore made on \([0,\tau_0]\) with \(\tau_0<\tau\), or under an explicit positive-risk condition. Extrapolating beyond the last informative risk set is a model assumption, not a consequence of censoring theory.