Tutorials › Competing Risks and Multi-State Models › A Worked Example With Two Constant Cause-Specific Hazards

Cause-Specific and Subdistribution Hazards · Tutorial 9 of 699

A Worked Example With Two Constant Cause-Specific Hazards

Learn to add two constant cause-specific hazards, interpret the resulting overall rate, and distinguish event shares from finite-time probabilities.

Intermediate 10 min read

What You'll Learn

  • Add two cause-specific hazards to obtain the overall first-event hazard for mutually exclusive causes.
  • Interpret each hazard as a conditional rate among people still event-free.
  • Check units and convert rates consistently before combining them.
  • Use constant hazards to distinguish a cause’s share of first events from its finite-time probability.
  • Recognize when a rate calculation is not justified by the assumptions.

From Cause-Specific Rates to an Overall Rate

When a first event can occur for one of two mutually exclusive reasons, each cause has its own cause-specific hazard. The hazards describe rates among people who are still event-free, not probabilities among everyone who began follow-up. This distinction matters when interpreting a reported rate or using it to calculate how quickly first events occur overall.

This tutorial works through a population model in which the two cause-specific hazards are constant over time. The main calculation is simple: add the cause-specific hazards to obtain the overall first-event hazard. The interpretation requires more care, especially when translating a rate into statements about the fraction of people who experience each cause by a particular time.

Definition: Suppose the possible first-event causes are \(A\) and \(B\), and a person is event-free at time \(t\). The cause-specific hazard \(\lambda_A(t)\) is the instantaneous rate of a first event from cause \(A\), conditional on being event-free just before \(t\); \(\lambda_B(t)\) is defined similarly. The overall first-event hazard is the instantaneous rate of any first event among people still event-free. For mutually exclusive first-event causes, \(\lambda(t)=\lambda_A(t)+\lambda_B(t)\).

A hazard has units of inverse time, such as events per person-year. It is not itself a probability. Informally, over a sufficiently short interval of length \(\Delta t\), a cause-specific hazard \(\lambda_A(t)\) corresponds to a conditional probability of a cause-\(A\) event of approximately \(\lambda_A(t)\Delta t\), provided that interval is short enough for the hazard to be nearly constant. The exact conditional probability over a longer interval depends on the hazards throughout that interval.

Assumptions for the Worked Calculations

The examples use a population in which every person starts event-free, and the first event must be either cause \(A\) or cause \(B\). The causes are mutually exclusive: a person can have at most one of them as their first event. Assume continuous event times, so simultaneous first events have probability zero, and suppose both cause-specific hazards are constant during the follow-up period being considered. A constant hazard means the cause-specific rate among those still event-free does not change with time; it does not mean the number of event-free people stays fixed.

These are stipulated population rates for the examples, not estimates from a sample. Therefore no standard errors, confidence intervals, or tests are needed. In an application, the constant-hazard assumption would need to be supported by the study design and data, rather than inferred just because a single rate has been reported.

Why the Two Rates Add

Over a very short interval, a first event can be from \(A\) or from \(B\), but not both. The probability of any first event is therefore the sum of the two disjoint cause-specific event probabilities. Dividing by the interval length and taking the small-interval limit gives the overall hazard.

$$ \lambda(t) = \lim_{\Delta t\to 0} \frac{P(\text{any first event in }[t,t+\Delta t)\mid\text{event-free at }t)}{\Delta t} = \lambda_A(t)+\lambda_B(t). $$

When the hazards are constant, write \(\lambda_A\) and \(\lambda_B\) for their values. Then \(\lambda=\lambda_A+\lambda_B\) is constant as well. This describes the instantaneous rate of any first event among the people who remain event-free at each time. As events accumulate, the event-free group can become smaller even though its hazard remains constant.

Worked Example: Adding Two Constant Annual Hazards

Suppose a hypothetical population has constant cause-specific hazards of \(\lambda_A=0.03\) per year and \(\lambda_B=0.02\) per year. Find the overall first-event hazard and interpret all three rates.

Calculation. The causes are mutually exclusive first-event types, so their hazards add:

$$ \lambda=0.03+0.02=0.05\text{ per year}. $$

Interpretation. Among people who are event-free at a given time, the instantaneous rate of a cause-\(A\) first event is \(0.03\) per year, and the corresponding rate for cause \(B\) is \(0.02\) per year. The instantaneous rate of any first event is \(0.05\) per year among that same event-free group.

The ratio \(0.03/0.05=0.60\) is the cause-\(A\) share of the total instantaneous hazard. Under these constant hazards, it is also the share of first events from cause \(A\) among first events occurring over any specified follow-up interval. It does not mean that 60% of all people initially at risk have cause \(A\) by a finite follow-up time.

What the Constant-Hazard Model Implies Over Time

The finite-time implication can be derived directly. Let \(S(t)\) be the probability that a person who starts event-free has had no first event by time \(t\). Over a short interval, the conditional probability of any first event is approximately \(\lambda\Delta t\), so the probability of remaining event-free is approximately \(1-\lambda\Delta t\). Thus \(S(t+\Delta t)\) is approximately \(S(t)(1-\lambda\Delta t)\). Taking the limit gives the differential equation \(S'(t)=-\lambda S(t)\), with \(S(0)=1\), whose solution is \(S(t)=\exp(-\lambda t)\).

For cause \(A\), the probability of a cause-\(A\) first event in a short interval after time \(t\) is approximately \(S(t)\lambda_A\Delta t\): the person must still be event-free, and then the cause-\(A\) hazard applies. Integrating this rate over follow-up gives the cumulative incidence for cause \(A\). With constant hazards, the resulting expression is:

$$ F_A(\tau) = \int_0^\tau S(t)\lambda_A\,dt = \frac{\lambda_A}{\lambda_A+\lambda_B} \left(1-\exp\left[-(\lambda_A+\lambda_B)\tau\right]\right). $$

The total probability of any first event by time \(\tau\) is \(1-\exp[-(\lambda_A+\lambda_B)\tau]\). Dividing \(F_A(\tau)\) by that total gives \(\lambda_A/(\lambda_A+\lambda_B)\) for \(\tau>0\). Therefore, because the hazards are constant here, 60% of first events by any specified follow-up time are from cause \(A\). The probability that an initially event-free person has cause \(A\) by that time is \(0.6(1-\exp(-0.05\tau))\), not 60%.

Worked Example: A Cause’s Share Is Not Its Finite-Time Probability

Use the annual hazards from the first example to calculate the probability of a cause-\(A\) event by 10 years and compare it with the total probability of any first event by 10 years. Assume the stated hazards remain constant for all 10 years.

The overall hazard is \(0.03+0.02=0.05\) per year. The probability of any first event by 10 years is:

$$ 1-\exp[-(0.05)(10)] = 1-\exp(-0.5) \approx 1-0.6065 = 0.3935. $$

The cause-\(A\) cumulative incidence is:

$$ F_A(10) = \frac{0.03}{0.05}\left(1-\exp(-0.5)\right) = 0.60(0.3935) \approx 0.2361. $$

Thus, about 23.61% of people initially event-free have a cause-\(A\) first event by 10 years under this model. The share of first events that are from cause \(A\) is \(0.2361/0.3935\approx0.6000\), or 60.00%. These are different quantities: the first is a probability among everyone initially event-free, while the second is a proportion among those who experience any first event by year 10.

For comparison, the cause-\(B\) cumulative incidence is \(0.40(0.3935)\approx0.1574\). The two cause-specific probabilities sum to \(0.2361+0.1574=0.3935\), the total probability of any first event by 10 years.

Worked Example: Converting Units Before Adding Hazards

In another hypothetical population, the constant cause-specific hazards are \(0.004\) per month for cause \(A\) and \(0.006\) per month for cause \(B\). Report the overall hazard per year and interpret the cause-specific rates. Assume the monthly rates apply throughout the year.

First add the rates in their given common unit:

$$ \lambda=0.004+0.006=0.010\text{ per month}. $$

There are 12 months in a year, so multiply each monthly rate by 12:

$$ \lambda_A=12(0.004)=0.048\text{ per year},\qquad \lambda_B=12(0.006)=0.072\text{ per year}, $$ $$ \lambda=12(0.010)=0.120\text{ per year}. $$

Among people still event-free, the cause-\(A\) instantaneous rate is \(0.048\) per year and the cause-\(B\) instantaneous rate is \(0.072\) per year; their sum, \(0.120\) per year, is the overall first-event rate. Cause \(A\)'s share of the total hazard is \(0.048/0.120=0.40\), and cause \(B\)'s share is \(0.60\). A rate per year is not an exact one-year event probability: converting it to a probability requires accounting for survival event-free during the year.

Checks, Failure Modes, and Alternatives

A short set of checks catches many errors in this calculation. Confirm that both cause-specific rates use the same time unit before adding them. Confirm that they refer to mutually exclusive first-event causes and to the event-free risk set. Finally, check that the total hazard is no smaller than either component and equals their sum. If either cause-specific hazard varies over time, the hazards still add at each time, but a single constant overall rate cannot represent the whole follow-up.

  • Adding probabilities instead of rates. Cause-specific hazards are rates. Exact probabilities over a finite interval are not generally found by simply adding hazards or multiplying a hazard by the entire interval length.
  • Using the wrong denominator. Each cause-specific hazard is conditional on being event-free just before the time in question. It is not a rate among all people initially enrolled regardless of whether they have already had an event.
  • Mixing units. A rate per month cannot be added directly to a rate per year. Convert both to the same unit first.
  • Overextending the cause share. A 60% share of first events does not say that 60% of the initial population experiences cause \(A\) by a specified finite time. The finite-time probability also depends on the total event rate and the length of follow-up.

The overall hazard is useful when the question concerns the instantaneous rate of any first event. Keeping cause-specific hazards separate is useful when the event type matters. Neither rate alone is a finite-time event probability. The constant-hazard formulas above provide a model-based translation between rates and probabilities; when hazards vary with time, finite-time probabilities require accounting for that time variation rather than inserting one average rate without justification.

Key takeaway: For two mutually exclusive first-event causes, add the cause-specific hazards to obtain the overall hazard at each time. Interpret each rate among people who remain event-free. With constant hazards, the cause-\(A\) share among first events over any specified follow-up interval is \(\lambda_A/(\lambda_A+\lambda_B)\), whereas the probability of cause \(A\) by that time is that share multiplied by the probability of any first event.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions.

  1. Constant hazards are \(0.07\) per year for cause \(A\) and \(0.03\) per year for cause \(B\). What is the overall hazard, and what is cause \(A\)'s share of the total hazard?
  2. For \(\lambda_A=0.02\) per year and \(\lambda_B=0.03\) per year, explain why saying “40% of the initial population has cause \(A\) by year 5” is not justified by the cause-\(A\) share alone.
  3. One hazard is \(0.005\) per month and another is \(0.020\) per year. What must be done before calculating their sum?
  4. If the cause-specific hazards change over time, does the overall hazard still equal their sum? What conclusion about a single constant overall rate follows?

Solutions. (1) The overall hazard is \(0.07+0.03=0.10\) per year; cause \(A\)'s share is \(0.07/0.10=0.70\), or 70%. (2) The 40% is cause \(A\)'s share among first events under the constant-hazard model, not the probability among everyone initially event-free. The latter is \(0.40(1-\exp[-(0.05)(5)])\approx0.0885\), or 8.85%. (3) Convert one rate so both use the same time unit; \(0.005\) per month is \(0.060\) per year, so the sum is \(0.080\) per year. (4) Yes. At each time \(t\), \(\lambda(t)=\lambda_A(t)+\lambda_B(t)\). If the component hazards vary, their sum may also vary, so one constant overall rate need not describe the follow-up.

Clinical Biostats