Tutorials › Statistical Theory I - Probability and Distribution Theory › Reconstructing Survival from Cumulative Hazard

Reliability and Extreme-Value Laws · Tutorial 671 of 994

Reconstructing Survival from Cumulative Hazard

Learn to move precisely between a continuous lifetime’s hazard, cumulative hazard, and survival function, including the conditions needed for each direction.

Intermediate 9 min read

What You'll Learn

  • Define cumulative hazard as the time integral of a nonnegative hazard rate.
  • Derive survival as the exponential of negative cumulative hazard.
  • Recover cumulative hazard from survival and hazard from its derivative.
  • Convert specified and piecewise hazard functions into survival functions.
  • Check when a hazard function defines a proper finite-lifetime distribution.
  • Distinguish cumulative hazard from a probability of failure.

From an Instantaneous Rate to Survival

For a lifetime \(T\), the survival function \(S(t)=P(T>t)\) describes the probability that the lifetime exceeds \(t\). The hazard \(h(t)\) describes the instantaneous failure rate at time \(t\), among units that have survived up to that time. These quantities answer different questions: survival is a probability, while hazard is a rate conditional on continued survival. Integrating the hazard over time accumulates that rate, and the accumulated quantity determines survival exactly in the continuous-lifetime setting developed here.

Assume \(T\) is nonnegative, has a density \(f\), and has no atom at zero. Then \(S(0)=1\), and on times for which \(S(t)>0\), the hazard is \(h(t)=f(t)/S(t)\). We will assume the hazard is measurable and integrable on every bounded interval where it is used. These conditions allow the integral and the logarithm in the derivation below to be treated rigorously.

Definition: For a continuous lifetime with hazard \(h\), its cumulative hazard up to time \(t\) is $$ H(t)=\int_0^t h(u)\,du, $$ whenever this integral is defined. The variable \(u\) is the integration time; \(t\) is the time at which the accumulated hazard is evaluated.

Cumulative hazard starts at zero and is nondecreasing because \(h(u)\geq 0\). It does not have to stay below one: it is not a probability. The key relation is that survival decays exponentially in this accumulated quantity.

Deriving Survival from Cumulative Hazard

Theorem: Suppose \(T\) has a density on \([0,\infty)\), has no atom at zero, and \(S(t)>0\) on the time interval under consideration. Define \(h(t)=f(t)/S(t)\) there, and suppose \(h\) is integrable on each bounded subinterval. Then, for every \(t\) in that interval, $$ S(t)=\exp\{-H(t)\} \qquad\text{where}\qquad H(t)=\int_0^t h(u)\,du. $$ Equivalently, \(H(t)=-\log S(t)\).

Proof. Since \(T\) has density \(f\), its survival function is absolutely continuous on bounded time intervals, with derivative \(S'(u)=-f(u)\) for almost every \(u\). On any interval \([0,t]\) on which \(S(t)>0\), monotonicity of \(S\) implies \(S(u)\geq S(t)>0\). Thus \(\log S(u)\) is well defined throughout that interval. The chain rule for absolutely continuous functions gives $$ \frac{d}{du}\log S(u)=\frac{S'(u)}{S(u)} =-\frac{f(u)}{S(u)} =-h(u) $$ for almost every \(u\in[0,t]\). Integrating both sides from \(0\) to \(t\), and using \(S(0)=1\), yields $$ \log S(t)-\log S(0)=-\int_0^t h(u)\,du, \qquad \log S(t)=-H(t). $$ Exponentiating gives \(S(t)=\exp\{-H(t)\}\), and rearranging gives \(H(t)=-\log S(t)\). \(\square\)

The proof shows why the initial condition matters: the integration constant is \(\log S(0)=0\). If a lifetime has positive probability of being zero, this version of the identity cannot be used from time zero without accounting separately for that atom.

Worked Example: A Hazard That Increases Linearly

Suppose the hazard is \(h(t)=0.2+0.1t\) for \(t\geq 0\), with time measured in years. Find the cumulative hazard and survival function, then calculate survival at five years.

Accumulate the hazard. Integrating each term gives $$ H(t)=\int_0^t (0.2+0.1u)\,du =0.2t+0.05t^2. $$

Recover survival. The identity gives $$ S(t)=\exp\{-(0.2t+0.05t^2)\}. $$ At \(t=5\), \(H(5)=0.2(5)+0.05(25)=1+1.25=2.25\), so $$ S(5)=e^{-2.25}\approx 0.1054. $$

This is a survival probability, so it lies between zero and one. The cumulative hazard \(2.25\) is not itself a probability. As a check, differentiating the displayed survival function gives \(S'(t)=-(0.2+0.1t)S(t)\), consistent with the specified hazard.

Recovering Hazard from Survival

The relation can also be read in the reverse direction. If a survival function is positive and differentiable, its logarithmic rate of decrease determines the hazard. Under absolute continuity, differentiability is needed only almost everywhere.

Theorem: Suppose \(S(0)=1\), \(S(t)>0\) on an interval, and \(S\) is absolutely continuous there. Define $$ H(t)=-\log S(t). $$ Then \(H\) is absolutely continuous and $$ H'(t)=-\frac{S'(t)}{S(t)} $$ for almost every \(t\). If the lifetime has density \(f(t)=-S'(t)\) almost everywhere, this derivative equals its hazard \(h(t)=f(t)/S(t)\).

Proof. On any bounded interval where \(S\) is positive, it is bounded away from zero by its value at the right endpoint, since survival is nonincreasing. The logarithm is continuously differentiable on the range of \(S\) there, so the chain rule for absolutely continuous functions applies. It gives \(H'(t)=-S'(t)/S(t)\) almost everywhere. When \(S'(t)=-f(t)\) almost everywhere, substitution gives \(H'(t)=f(t)/S(t)=h(t)\). Since \(H(0)=-\log 1=0\), integrating its derivative recovers \(H(t)=\int_0^t h(u)\,du\). \(\square\)

Conversely, a nonnegative, measurable hazard that is integrable on bounded intervals can be used to construct a survival function by setting \(H(t)=\int_0^t h(u)\,du\) and \(S(t)=e^{-H(t)}\). This function is nonincreasing, starts at one, and is absolutely continuous on bounded intervals. Almost everywhere, \(S'(t)=-h(t)S(t)\). For this construction to describe a lifetime that is finite with probability one, it must also satisfy \(S(t)\to 0\) as \(t\to\infty\), equivalently \(H(t)\to\infty\). If \(H\) instead approaches a finite limit, the constructed survival function approaches a positive limit; it does not describe a finite lifetime with probability one.

Worked Example: A Piecewise Hazard

Suppose a device has hazard \(h(t)=0.4\) for \(0\leq t\leq 2\), and \(h(t)=0.4+0.3(t-2)\) for \(t>2\). Find its cumulative hazard and survival at \(t=2\) and \(t=4\).

For \(0\leq t\leq 2\), direct integration gives \(H(t)=0.4t\). In particular, \(H(2)=0.8\) and \(S(2)=e^{-0.8}\approx 0.4493\).

For \(t>2\), split the integral at the point where the hazard changes: $$ H(t)=\int_0^2 0.4\,du+\int_2^t [0.4+0.3(u-2)]\,du =0.8+0.4(t-2)+0.15(t-2)^2. $$ Consequently \(H(4)=0.8+0.4(2)+0.15(2^2)=0.8+0.8+0.6=2.2\), and \(S(4)=e^{-2.2}\approx 0.1108\).

The two expressions for \(H\) agree at \(t=2\), so the cumulative hazard is continuous there. The survival function is continuous there as well, and its value from either side is \(e^{-0.8}\). A change in hazard does not create a jump in survival in this continuous-lifetime model.

Worked Example: Recovering Both Quantities from Survival

Let \(S(t)=(1+t/3)^{-2}\) for \(t\geq 0\). Find the cumulative hazard and hazard, and verify that this is a proper survival function for a finite lifetime.

First, \(S(0)=1\), and \(S(t)>0\). Taking the negative logarithm gives $$ H(t)=-\log S(t)=2\log(1+t/3). $$ Differentiating, \(H'(t)=2/(3+t)\), so the hazard is \(h(t)=2/(3+t)\).

To check the reverse calculation directly, integrate the hazard: $$ \int_0^t \frac{2}{3+u}\,du =2\bigl[\log(3+u)\bigr]_0^t =2\log\left(\frac{3+t}{3}\right) =2\log(1+t/3). $$ Thus \(\exp\{-H(t)\}=(1+t/3)^{-2}\), as specified. Finally, \(S(t)\to 0\) as \(t\to\infty\), so the survival function does not leave positive probability at an infinite lifetime.

What Cumulative Hazard Does—and Does Not—Measure

The value \(H(t)\) measures accumulated hazard, not the probability of failure by time \(t\). The probability of surviving to \(t\) is \(e^{-H(t)}\); therefore the probability of failing by \(t\) is \(1-e^{-H(t)}\). This distinction follows from the survival identity, rather than from treating hazard as a probability.

A common mistake is to add hazards over time and call the result a failure probability. For example, if \(H(t)=0.8\), the failure probability by time \(t\) is \(1-e^{-0.8}\), approximately \(0.5507\), not \(0.8\). The correction is to accumulate the hazard first and then apply the exponential relation.

The ordinary integral formula in this tutorial applies to continuous lifetimes with a density. A lifetime with atoms can have jumps in its survival function; an ordinary density-based hazard integral alone does not account for those jumps. Within the continuous setting, however, the inverse relations are exact: integrate hazard to obtain cumulative hazard, exponentiate its negative to obtain survival, or take the negative logarithm of survival to recover cumulative hazard.

Key takeaway: For a continuous lifetime with \(S(0)=1\) and positive survival, $$ H(t)=\int_0^t h(u)\,du=-\log S(t), \qquad S(t)=e^{-H(t)}. $$ Hazard is the instantaneous conditional rate; cumulative hazard is its integral; survival is the exponential transformation of that integral.

Check Your Understanding

Answer each question using the continuous-lifetime assumptions stated above.

  1. If \(h(t)=0.3\) for \(0\leq t\leq 4\), what are \(H(4)\) and \(S(4)\)?
  2. Given \(S(t)=e^{-0.5t^2}\) for \(t\geq 0\), find \(H(t)\) and \(h(t)\).
  3. Why must survival be positive on the interval when using \(H(t)=-\log S(t)\)?
  4. Suppose \(H(t)\) approaches \(1.7\) as \(t\to\infty\). What does the corresponding limiting survival value imply about a finite lifetime?
  5. Explain why a cumulative hazard of \(0.8\) does not mean an 80% probability of failure.