Tutorials › Statistical Theory I - Probability and Distribution Theory › Cumulative-Hazard Criteria for Lifetime Comparisons

Stochastic Orders and Comparisons · Tutorial 951 of 994

Cumulative-Hazard Criteria for Lifetime Comparisons

Learn how differences in cumulative hazards determine hazard-rate order, and calculate the comparison for exponential and Weibull lifetimes.

Advanced 10 min read

What You'll Learn

  • Define cumulative hazards for nonnegative, absolutely continuous lifetime distributions.
  • Relate hazard-rate order to monotonicity of a survival ratio.
  • Derive an equivalent criterion using differences of cumulative hazards.
  • Express the comparison through integrated hazard differences over time intervals.
  • Compare exponential lifetimes and Weibull lifetimes with a common shape parameter.
  • Identify why cumulative-hazard pointwise inequalities alone do not establish hazard-rate order.

Why Compare Cumulative Hazards?

A hazard rate describes the instantaneous failure rate among individuals who have survived up to a given time. Comparing two hazard rates directly can be useful, but the functions may be difficult to inspect or may cross. Cumulative hazards provide another route: they aggregate hazard over time, and their difference reveals whether the relative survival of two lifetimes moves consistently in one direction.

The key distinction is between comparing cumulative hazards at a single time and comparing how their difference changes over time. Hazard-rate order is determined by the latter. The resulting criteria are especially convenient for parametric lifetime models, including exponential and Weibull distributions.

Definitions and the Comparison Direction

Let \(X\) and \(Y\) be nonnegative, absolutely continuous lifetimes with densities \(f_X\) and \(f_Y\), survival functions \(S_X(t)=P(X>t)\) and \(S_Y(t)=P(Y>t)\), and \(S_X(t),S_Y(t)>0\) for every finite \(t\geq 0\). Their hazard rates and cumulative hazards are

$$ h_X(t)=\frac{f_X(t)}{S_X(t)},\qquad H_X(t)=\int_0^t h_X(u)\,du=-\log S_X(t), $$

with analogous definitions for \(Y\). The equality \(H_X(t)=-\log S_X(t)\) follows from absolute continuity and \(S_X(0)=1\). Hazard rates are defined almost everywhere; the integral defining \(H_X\) does not depend on their values on a set of measure zero.

Definition: \(X\) is smaller than \(Y\) in hazard-rate order, written \(X\leq_{\mathrm{hr}}Y\), if \(h_X(t)\geq h_Y(t)\) for almost every \(t\geq0\). Under this convention, \(X\) has at least as large an instantaneous failure rate as \(Y\) at each time, apart from a set of times of measure zero.

The comparison direction matters: a larger hazard for \(X\) corresponds to \(X\leq_{\mathrm{hr}}Y\). If hazards are continuous, the almost-everywhere inequality is equivalent to the inequality holding at every time. When the lifetime distributions have finite endpoints or other support complications, criteria must be applied on the relevant common survival interval and the supports must be considered separately. The results below use the stated positive-survival assumption.

The Cumulative-Hazard and Survival-Ratio Criteria

Define the cumulative-hazard difference \(D(t)=H_X(t)-H_Y(t)\). Since \(H_X(0)=H_Y(0)=0\), we have \(D(0)=0\). The survival ratio in the direction that matches \(X\leq_{\mathrm{hr}}Y\) is

$$ R(t)=\frac{S_Y(t)}{S_X(t)} =\exp\{H_X(t)-H_Y(t)\} =\exp\{D(t)\}. $$

Thus a nondecreasing difference \(D\) gives a nondecreasing survival ratio \(R\). The next result shows that these conditions are not just convenient consequences; they are equivalent to the hazard-rate comparison.

Theorem: Under the assumptions above, the following statements are equivalent:
  • \(X\leq_{\mathrm{hr}}Y\), meaning \(h_X(t)\geq h_Y(t)\) almost everywhere.
  • \(H_X(t)-H_Y(t)\) is nondecreasing in \(t\).
  • \(S_Y(t)/S_X(t)\) is nondecreasing in \(t\).

Proof. The cumulative hazards are locally absolutely continuous and satisfy \(H_X'(t)=h_X(t)\) and \(H_Y'(t)=h_Y(t)\) almost everywhere. Therefore \(D\) is locally absolutely continuous, with \(D'(t)=h_X(t)-h_Y(t)\) almost everywhere. A locally absolutely continuous function is nondecreasing exactly when its derivative is nonnegative almost everywhere. Hence \(h_X\geq h_Y\) almost everywhere if and only if \(D\) is nondecreasing.

For every finite \(t\), \(R(t)=\exp\{D(t)\}\). The exponential function is strictly increasing, so \(D\) is nondecreasing if and only if \(R\) is nondecreasing. This proves all three equivalences. \(\square\)

This result uses the difference of cumulative hazards, not their ratio. Equivalently, on any interval where the hazards are continuous, one can differentiate the log survival ratio:

$$ \frac{d}{dt}\log\!\left(\frac{S_Y(t)}{S_X(t)}\right) =h_X(t)-h_Y(t). $$

The derivative form is often the quickest calculation for smooth models. The cumulative-hazard form remains useful when comparing integrated hazards or when derivatives need only be interpreted almost everywhere.

An Interval-by-Interval Criterion

The same comparison can be expressed without checking whether a function is globally increasing. For any \(0\leq s<t\), the survival ratio changes between \(s\) and \(t\) according to the integrated difference in hazards over that interval.

Interval criterion: Under the theorem's assumptions,
$$ \frac{S_Y(t)/S_X(t)}{S_Y(s)/S_X(s)} =\exp\!\left\{\int_s^t \bigl(h_X(u)-h_Y(u)\bigr)\,du\right\}. $$

In particular, \(X\leq_{\mathrm{hr}}Y\) if and only if this ratio of survival ratios is at least \(1\) for every \(0\leq s<t\).

Proof. By the definition of cumulative hazard, \(H_X(t)-H_X(s)=\int_s^t h_X(u)\,du\), and the same identity holds for \(Y\). Taking the difference gives \(D(t)-D(s)=\int_s^t(h_X(u)-h_Y(u))\,du\). Since \(R(v)=\exp\{D(v)\}\), dividing \(R(t)\) by \(R(s)\) yields the displayed identity. The ratio is at least \(1\) for all \(s<t\) exactly when \(D(t)\geq D(s)\) for all such pairs, which is precisely that \(D\) is nondecreasing. By the theorem, this is equivalent to hazard-rate order. \(\square\)

This form emphasizes that the comparison is cumulative over each time interval. A brief interval on which \(h_X<h_Y\) rules out the order if it makes the cumulative-hazard difference decrease there; positive differences elsewhere cannot repair a decrease that has already occurred.

Worked Comparisons for Lifetime Models

Worked Example: Exponential Lifetimes

Suppose \(X\) is exponential with rate \(2\) and \(Y\) is exponential with rate \(1\). For an exponential lifetime with rate \(\lambda>0\), the density is \(f(t)=\lambda e^{-\lambda t}\), the survival function is \(S(t)=e^{-\lambda t}\), the hazard is \(h(t)=\lambda\), and the cumulative hazard is \(H(t)=\lambda t\).

Thus

$$ h_X(t)=2,\quad H_X(t)=2t,\qquad h_Y(t)=1,\quad H_Y(t)=t. $$

The hazard difference is \(h_X(t)-h_Y(t)=1\), which is positive at every time. The cumulative-hazard difference and survival ratio are

$$ H_X(t)-H_Y(t)=t,\qquad \frac{S_Y(t)}{S_X(t)} =\frac{e^{-t}}{e^{-2t}}=e^t. $$

Both criteria are nondecreasing, so \(X\leq_{\mathrm{hr}}Y\). In fact, over any interval \([s,t]\) the survival ratio is multiplied by \(e^{t-s}\), as the interval criterion predicts.

Worked Example: Weibull Lifetimes With a Common Shape

Use the Weibull shape-scale parameterization: for shape \(k>0\) and scale \(\eta>0\), a Weibull lifetime has survival function \(S(t)=\exp\{-(t/\eta)^k\}\), for \(t\geq0\). Its cumulative hazard and hazard are

$$ H(t)=\left(\frac{t}{\eta}\right)^k,\qquad h(t)=\frac{k t^{k-1}}{\eta^k}\quad(t>0). $$

Let \(X\) have shape \(k=2\) and scale \(\eta_X=2\), and let \(Y\) have shape \(k=2\) and scale \(\eta_Y=3\). Then

$$ H_X(t)=\frac{t^2}{4},\qquad H_Y(t)=\frac{t^2}{9}, \qquad h_X(t)=\frac{t}{2},\qquad h_Y(t)=\frac{2t}{9}. $$

For \(t\geq0\), the hazard difference is \(t/2-2t/9=5t/18\geq0\). Equivalently,

$$ H_X(t)-H_Y(t) =t^2\left(\frac14-\frac19\right) =\frac{5t^2}{36}, \qquad \frac{S_Y(t)}{S_X(t)}=\exp\!\left(\frac{5t^2}{36}\right). $$

The cumulative-hazard difference and survival ratio are nondecreasing, so \(X\leq_{\mathrm{hr}}Y\). The hazards agree at \(t=0\), but that single time does not affect the almost-everywhere comparison. For a common shape \(k\), the same calculation shows that the larger scale has the smaller hazard at positive times.

Worked Example: Weibull Shapes Can Reverse the Comparison

Let \(X\) be exponential with rate \(1\), and let \(Y\) be Weibull with shape \(2\) and scale \(1\). The exponential law is also Weibull with shape \(1\) and scale \(1\), but here the shapes differ. Their hazards and cumulative hazards are

$$ h_X(t)=1,\quad H_X(t)=t,\qquad h_Y(t)=2t,\quad H_Y(t)=t^2. $$

The hazard difference is \(h_X(t)-h_Y(t)=1-2t\): it is positive for \(0\leq t<1/2\) and negative for \(t>1/2\). The cumulative-hazard difference is \(D(t)=t-t^2\), which increases up to \(t=1/2\) and decreases afterward. Correspondingly,

$$ \frac{S_Y(t)}{S_X(t)}=\exp(t-t^2). $$

This ratio first increases and then decreases. Therefore \(X\not\leq_{\mathrm{hr}}Y\). The reverse order also fails: it would require \(h_Y\geq h_X\) almost everywhere, but \(2t<1\) for \(0\leq t<1/2\). Thus these two lifetimes have no hazard-rate order in either direction.

Common Pitfall: Pointwise Cumulative-Hazard Inequalities

A frequent mistake is to infer \(X\leq_{\mathrm{hr}}Y\) from \(H_X(t)\geq H_Y(t)\) for every \(t\). That inequality says only that \(D(t)\geq0\); it does not say that \(D\) is nondecreasing. Hazard-rate order requires \(D\) to be nondecreasing, or equivalently requires its almost-everywhere derivative \(h_X-h_Y\) to be nonnegative.

For a concrete check, define two valid cumulative hazards by

$$ H_Y(t)=t,\qquad H_X(t)=t+t e^{-t},\qquad t\geq0. $$

Both start at zero and are increasing: \(H_Y'(t)=1\), while \(H_X'(t)=1+(1-t)e^{-t}>0\) for all \(t\geq0\). They define proper continuous lifetime distributions through \(S_X=e^{-H_X}\) and \(S_Y=e^{-H_Y}\). Although \(H_X(t)-H_Y(t)=t e^{-t}\geq0\) at every time, its derivative is \((1-t)e^{-t}\), which becomes negative for \(t>1\). Hence the survival ratio \(S_Y/S_X=\exp(t e^{-t})\) eventually decreases, and \(X\leq_{\mathrm{hr}}Y\) does not hold. The cumulative-hazard difference, not its sign at each time, is the decisive quantity.

Key takeaway: For continuous lifetimes with positive survival, \(X\leq_{\mathrm{hr}}Y\) exactly when \(H_X-H_Y\) is nondecreasing, or equivalently when \(S_Y/S_X\) is nondecreasing. Pointwise ordering of the cumulative hazards alone is not enough.

Check Your Understanding

Use the hazard-rate convention and criteria developed above.

  1. For \(X\sim\operatorname{Exp}(3)\) and \(Y\sim\operatorname{Exp}(2)\), calculate \(H_X-H_Y\) and \(S_Y/S_X\), then state the hazard-rate order.
  2. For two Weibull lifetimes with common shape \(k>0\) and scales \(\eta_X<\eta_Y\), express \(H_X-H_Y\) and determine its monotonicity.
  3. Why does \(H_X(t)\geq H_Y(t)\) at every \(t\) not, by itself, imply \(X\leq_{\mathrm{hr}}Y\)? State the criterion that is sufficient.
  4. Suppose \(S_Y(t)/S_X(t)\) decreases between two times \(s<t\). Which hazard-rate order does this rule out, and why?
  5. If \(h_X-h_Y\) is nonnegative almost everywhere, what does the interval criterion say about \(S_Y(t)/S_X(t)\) as \(t\) increases?