Tutorials › Statistical Theory I - Probability and Distribution Theory › Hazard-Rate Comparisons for Exponential and Weibull Laws

Stochastic Orders and Comparisons · Tutorial 966 of 994

Hazard-Rate Comparisons for Exponential and Weibull Laws

Learn to compare exponential and Weibull lifetimes by calculating their hazards, survival ratios, and any hazard-crossing time.

Advanced 10 min read

What You'll Learn

  • Calculate exponential and Weibull hazard functions from their survival functions.
  • Relate the survival ratio to the difference between hazards.
  • Determine the hazard-rate order for exponential laws and equal-shape Weibull laws.
  • Show why Weibull laws with different shapes have crossing hazards.
  • Calculate and interpret the unique crossing time for two Weibull hazards.
  • Distinguish pointwise hazard statements for smooth laws from claims about arbitrary density versions.

Why Compare Hazards?

For a lifetime \(T\), the hazard at time \(t\) describes the instantaneous event rate among individuals who have survived to \(t\). Comparing hazards can therefore reveal how the relative event rates evolve over time. Exponential and Weibull laws make this comparison especially transparent: their hazards are constants or simple power functions.

The direction of the comparison matters. If \(X\) has a higher hazard than \(Y\) at every positive time, then \(X\) is smaller in hazard-rate order: its survival is relatively less favorable as time passes. If the hazards cross, neither lifetime has a hazard advantage that holds at all times. The calculations below make these claims precise for the two distribution families.

Hazards and Survival Ratios

Write \(S_T(t)=P(T>t)\) for the survival function of a nonnegative lifetime \(T\). For the smooth survival functions used here, define its hazard by the negative derivative of the log survival function. This definition specifies a pointwise hazard without ambiguity about which version of a density is chosen.

Definition: If \(S_T(t)>0\) and \(S_T\) is differentiable at \(t>0\), its hazard is \[ h_T(t)=-\frac{d}{dt}\log S_T(t)=-\frac{S'_T(t)}{S_T(t)}. \] For the exponential and Weibull laws considered below, the density is smooth on \(t>0\), and \(h_T(t)=f_T(t)/S_T(t)\) there. For general absolutely continuous distributions, \(f_T/S_T\) agrees with this derivative only almost everywhere unless additional regularity is imposed.

For two positive, differentiable survival functions, the survival ratio \(S_Y(t)/S_X(t)\) connects directly to the hazards. This is the pointwise criterion used throughout the tutorial.

Theorem: Suppose \(S_X\) and \(S_Y\) are positive and differentiable on \((0,\infty)\), and define \(h_X=-S'_X/S_X\) and \(h_Y=-S'_Y/S_Y\). Then \[ \frac{S_Y(t)}{S_X(t)}\text{ is nondecreasing on }(0,\infty) \quad\Longleftrightarrow\quad h_X(t)\ge h_Y(t)\text{ for every }t>0. \] In this setting, \(X\le_{\mathrm{hr}}Y\) means that \(X\) is smaller than \(Y\) in hazard-rate order.

To prove the criterion, differentiate the ratio using the product and quotient rules:

$$ \left(\frac{S_Y}{S_X}\right)' =\frac{S_Y}{S_X}\left(\frac{S'_Y}{S_Y}-\frac{S'_X}{S_X}\right) =\frac{S_Y}{S_X}(h_X-h_Y). $$

The ratio is positive, so its derivative is nonnegative exactly when \(h_X-h_Y\ge0\). If \(h_X\ge h_Y\) everywhere, the derivative is nonnegative everywhere, and the mean value theorem implies that the ratio is nondecreasing. Conversely, a differentiable nondecreasing function has a nonnegative derivative wherever its derivative exists; applying this fact to \(S_Y/S_X\) gives \(h_X\ge h_Y\). \(\square\)

Exponential and Weibull Hazard Functions

Use the rate parameterization \(T\sim\operatorname{Exp}(\lambda)\), with \(\lambda>0\), for the exponential law. Its survival function is \(S_T(t)=e^{-\lambda t}\) for \(t\ge0\). A Weibull law will be written with shape \(k>0\) and scale \(\eta>0\), so its survival function is \(S_T(t)=\exp[-(t/\eta)^k]\) for \(t\ge0\).

Result: For \(t>0\), the hazards are \[ h_{\operatorname{Exp}(\lambda)}(t)=\lambda, \qquad h_{\operatorname{Weibull}(k,\eta)}(t)=\frac{k}{\eta^k}t^{k-1}. \] The exponential law is the Weibull law with \(k=1\) and \(\eta=1/\lambda\).

For the exponential case, differentiating \(\log S_T(t)=-\lambda t\) gives the constant hazard \(\lambda\). For the Weibull case, \(\log S_T(t)=-t^k/\eta^k\), so differentiating and changing sign gives \(k t^{k-1}/\eta^k\). Equivalently, differentiating the Weibull distribution function gives \[ f_T(t)=\frac{k}{\eta^k}t^{k-1}\exp[-(t/\eta)^k], \] and division by its positive survival function gives the same hazard. These expressions are pointwise valid for every \(t>0\), including when the hazard tends to zero or infinity as \(t\) approaches zero.

Parameter Cases for Hazard-Rate Order

For two exponential lifetimes, the hazard comparison is determined entirely by their rates. For two Weibull lifetimes with the same shape, it is determined by the coefficients multiplying the common power of time. If the shapes differ, the powers themselves determine the behavior near zero and at large times.

Theorem: Let \(X\sim\operatorname{Exp}(\lambda_X)\) and \(Y\sim\operatorname{Exp}(\lambda_Y)\), with positive rates. Then \(X\le_{\mathrm{hr}}Y\) if and only if \(\lambda_X\ge\lambda_Y\). For \(X\sim\operatorname{Weibull}(k,\eta_X)\) and \(Y\sim\operatorname{Weibull}(k,\eta_Y)\) with the same positive shape \(k\), \(X\le_{\mathrm{hr}}Y\) if and only if \(\eta_X\le\eta_Y\).

For exponentials, the hazards are the constants \(\lambda_X\) and \(\lambda_Y\); the hazard-rate criterion therefore gives the first equivalence. The survival ratio makes the same comparison explicit:

$$ \frac{S_Y(t)}{S_X(t)} =\exp[-(\lambda_Y-\lambda_X)t] =\exp[(\lambda_X-\lambda_Y)t]. $$

This ratio is nondecreasing exactly when \(\lambda_X\ge\lambda_Y\). For equal-shape Weibull laws, the hazard difference is \[ h_X(t)-h_Y(t)=k t^{k-1}\left(\frac{1}{\eta_X^k}-\frac{1}{\eta_Y^k}\right). \] Since \(k t^{k-1}>0\) for \(t>0\), it is nonnegative everywhere exactly when \(\eta_X\le\eta_Y\). This proves the second equivalence.

Worked Example: Comparing Two Exponential Lifetimes

Let \(X\sim\operatorname{Exp}(0.4)\) and \(Y\sim\operatorname{Exp}(0.25)\), with rates measured per year. Find their hazards and survival ratio, and decide their hazard-rate order.

The hazards are constant: \[ h_X(t)=0.4,\qquad h_Y(t)=0.25. \] Thus \(h_X(t)-h_Y(t)=0.15>0\) for every \(t>0\), so \(X\le_{\mathrm{hr}}Y\).

Their survival functions are \(S_X(t)=e^{-0.4t}\) and \(S_Y(t)=e^{-0.25t}\). Therefore, \[ \frac{S_Y(t)}{S_X(t)} =\frac{e^{-0.25t}}{e^{-0.4t}} =e^{0.15t}. \] Its derivative is \(0.15e^{0.15t}>0\), confirming that the ratio is strictly increasing. The larger rate for \(X\) means its survival decays faster.

Worked Example: Equal-Shape Weibull Laws

Let \(X\sim\operatorname{Weibull}(2,2)\) and \(Y\sim\operatorname{Weibull}(2,3)\), where the second parameter is the scale. Calculate the hazards and ratio, then compare them.

The hazards are \[ h_X(t)=\frac{2}{2^2}t=\frac{t}{2}, \qquad h_Y(t)=\frac{2}{3^2}t=\frac{2t}{9}. \] Their difference is \(t/2-2t/9=5t/18>0\) for \(t>0\), so \(X\le_{\mathrm{hr}}Y\).

The survival functions are \(S_X(t)=e^{-t^2/4}\) and \(S_Y(t)=e^{-t^2/9}\). Hence \[ \frac{S_Y(t)}{S_X(t)} =\exp\left(\frac{t^2}{4}-\frac{t^2}{9}\right) =\exp\left(\frac{5t^2}{36}\right). \] Its derivative is \((5t/18)\exp(5t^2/36)>0\) for \(t>0\), consistent with the hazard comparison. Here the smaller scale for \(X\) produces the higher hazard at every positive time.

Different Weibull Shapes: Hazards Must Cross

Write the Weibull hazard as \(h(t)=A t^{k-1}\), where \(A=k/\eta^k>0\). For two different shapes, the ratio of the hazards is a nonconstant power of \(t\). That ratio tends either to zero or to infinity near zero and has the opposite limiting behavior as \(t\) grows without bound. Consequently, the hazard ordering reverses exactly once.

Result: Suppose \(X\sim\operatorname{Weibull}(k_X,\eta_X)\) and \(Y\sim\operatorname{Weibull}(k_Y,\eta_Y)\), with positive parameters and \(k_X\ne k_Y\). Their hazards cross once, at \[ t_*=\left(\frac{k_Y/\eta_Y^{k_Y}}{k_X/\eta_X^{k_X}}\right)^{1/(k_X-k_Y)}. \] For \(k_X>k_Y\), \(h_X(t)<h_Y(t)\) when \(0<t<t_*\), and \(h_X(t)>h_Y(t)\) when \(t>t_*\). For \(k_X<k_Y\), these inequalities are reversed. Thus neither law is smaller in hazard-rate order.

To derive the crossing time, set the hazards equal and define \(A_X=k_X/\eta_X^{k_X}\), \(A_Y=k_Y/\eta_Y^{k_Y}\): \[ A_X t^{k_X-1}=A_Y t^{k_Y-1} \quad\Longleftrightarrow\quad t^{k_X-k_Y}=\frac{A_Y}{A_X}. \] The right side is positive, and the nonzero power has exactly one positive solution, \(t_*\) as displayed. Moreover, \[ \frac{h_X(t)}{h_Y(t)} =\frac{A_X}{A_Y}t^{k_X-k_Y}. \] If \(k_X>k_Y\), this ratio increases strictly from zero to infinity, so it equals one once and gives the stated direction on either side. If \(k_X<k_Y\), it decreases strictly from infinity to zero, giving the reversed direction. Since each hazard is higher on one side of the crossing and lower on the other, neither can dominate for every \(t>0\). \(\square\)

Worked Example: An Exponential and a Weibull Hazard Cross

Let \(X\sim\operatorname{Weibull}(1,1)\), which is \(\operatorname{Exp}(1)\), and let \(Y\sim\operatorname{Weibull}(2,1)\). Find the crossing time and the behavior of the survival ratio.

The hazards are \(h_X(t)=1\) and \(h_Y(t)=2t\). Setting them equal gives \(1=2t\), so \(t_*=1/2\). For \(0<t<1/2\), \(h_X(t)>h_Y(t)\); for \(t>1/2\), \(h_X(t)<h_Y(t)\). The hazards cross, so neither lifetime is smaller in hazard-rate order.

The survivals are \(S_X(t)=e^{-t}\) and \(S_Y(t)=e^{-t^2}\). Their ratio and its derivative are \[ \frac{S_Y(t)}{S_X(t)}=e^{t-t^2}, \qquad \left(\frac{S_Y}{S_X}\right)'=e^{t-t^2}(1-2t). \] The derivative is positive before \(1/2\), zero at \(1/2\), and negative afterward. The ratio first increases and then decreases, agreeing with the hazard crossing.

Common Mistake and Interpretation

Common mistake: “Both hazards are positive, so one lifetime must be smaller in hazard-rate order.” Positivity only says that each hazard is nonnegative; it does not compare the hazards. In the last example, both hazards are positive for every \(t>0\), yet their difference changes sign at \(1/2\). The correction is to compare \(h_X(t)-h_Y(t)\) over the full time range, or equivalently check whether \(S_Y(t)/S_X(t)\) is nondecreasing over that range.

A second point of care concerns pointwise claims. An arbitrary density can be altered at an individual point without changing its distribution, so an equality between \(f/S\) and a derivative-based hazard need not hold pointwise for every density version. Here the exponential and Weibull survival functions are differentiable on \((0,\infty)\), and their displayed hazards are defined as \(-S'/S\); their smooth density versions agree with these hazards there. This regularity is why the pointwise comparisons and crossing statements are justified.

The parameter patterns are useful when selecting a comparison: exponential rates decide the constant-hazard case, a common Weibull shape reduces the comparison to scales, and different Weibull shapes force a crossing. The hazard-rate order records an all-times comparison; a crossing means that this particular order does not hold in either direction.

Check Your Understanding

Use the hazard functions and survival-ratio criterion to answer each question.

  1. If \(X\sim\operatorname{Exp}(0.3)\) and \(Y\sim\operatorname{Exp}(0.5)\), which lifetime is smaller in hazard-rate order? Verify using the survival ratio.
  2. For Weibull laws with common shape \(k=3\) and scales \(\eta_X=4\), \(\eta_Y=2\), which hazard is larger for every positive time?
  3. Why must two Weibull hazards with distinct shapes cross at exactly one positive time?
  4. For \(X\sim\operatorname{Weibull}(1,2)\) and \(Y\sim\operatorname{Weibull}(2,1)\), calculate the hazard-crossing time.
  5. Why is defining hazard as \(-S'/S\) useful when making a pointwise comparison?