Tutorials › Competing Risks and Multi-State Models › Why Transition Hazard Ratios Are Not Overall Risk Ratios

Regression for Multi-State Data · Tutorial 386 of 699

Why Transition Hazard Ratios Are Not Overall Risk Ratios

Learn to distinguish conditional transition hazard ratios from horizon-specific probabilities of experiencing or occupying a state.

Advanced 10 min read

What You'll Learn

  • Define a transition hazard ratio as a conditional comparison for one transition.
  • Derive cumulative incidence under constant competing transition hazards.
  • Calculate state-occupation probabilities in an illness-death model.
  • Compare a transition hazard ratio with risk ratios and risk differences at a fixed time.
  • Identify assumptions and diagnostics that matter when interpreting transition regression.

Why a Transition Effect Is Not an Overall Risk Effect

A transition hazard ratio describes a relative rate for a specified move between states, among individuals currently eligible to make that move. A cumulative incidence or state-occupation probability instead describes how much probability has accumulated along all relevant paths by a chosen time. Those are different quantities: other outgoing transitions, later transitions, and time spent in each state all affect the overall probability.

This distinction matters in interpreting transition regression. A coefficient can be a clear, transition-specific association without being a risk ratio for a clinical endpoint. To compare overall probabilities, one must specify the endpoint and time horizon, and account for the full transition structure. The calculations below use simple continuous-time models so that the difference is visible directly.

Definition: For a transition from state \(r\) to state \(s\), the transition hazard \(\lambda_{rs}(t\mid z)\) is the instantaneous rate of moving from \(r\) to \(s\) at time \(t\), conditional on being in state \(r\) just before \(t\) and on covariate value \(z\). A transition hazard ratio compares these rates for two covariate values at the same time. A cumulative incidence is the probability of having made a specified transition by a stated time, accounting for competing transitions. A state-occupation probability is the probability of being in a specified state at that time.

For a binary covariate \(z\), a proportional transition-hazards model for one transition can be written as \(\lambda_{rs}(t\mid z)=\lambda_{rs,0}(t)\exp(\beta_{rs}z)\). Its hazard ratio comparing \(z=1\) with \(z=0\) is \(\exp(\beta_{rs})\). This comparison is conditional on occupying the origin state; it does not compare the marginal probabilities of reaching or occupying the destination state by time \(t\).

A Competing-Destinations Calculation

Start with a person in state 0, from which two absorbing destinations, state 1 and state 2, are possible. Assume constant transition hazards \(\lambda_{01}\) and \(\lambda_{02}\), no return transitions, and no censoring in the probability calculation. The time to leave state 0 then has total rate \(\lambda_{01}+\lambda_{02}\). The cumulative incidence of transition \(0\to1\) by time \(t\) is the probability of still being in state 0 at time \(u\), multiplied by the transition rate to state 1, integrated over possible transition times \(u\).

$$ F_1(t) =\int_0^t \exp\{-(\lambda_{01}+\lambda_{02})u\}\lambda_{01}\,du =\frac{\lambda_{01}}{\lambda_{01}+\lambda_{02}} \left[1-\exp\{-(\lambda_{01}+\lambda_{02})t\}\right]. $$

The rate \(\lambda_{02}\) enters this probability even though it is not the hazard for transition \(0\to1\). Increasing the hazard of the competing transition can reduce the chance of reaching state 1 first. As a result, even a proportional change in \(\lambda_{01}\) does not generally produce the same proportional change in \(F_1(t)\).

Worked Example: A Twofold Transition Hazard, but a Smaller Risk Ratio

Suppose the constant transition hazards per year are \(\lambda_{01}=0.40\) and \(\lambda_{02}=0.60\) for \(z=0\). For \(z=1\), let \(\lambda_{01}=0.80\) and \(\lambda_{02}=0.60\). These values correspond to a \(0\to1\) hazard ratio of \(0.80/0.40=2\), with the competing transition unchanged. Calculate the cumulative incidence of state 1 by one year in each group.

For \(z=0\), the total rate out of state 0 is \(0.40+0.60=1.00\) per year. Thus

$$ F_1(1\mid z=0) =\frac{0.40}{1.00}\left(1-e^{-1.00}\right) =0.40(0.6321) \approx 0.2528. $$

For \(z=1\), the total rate is \(0.80+0.60=1.40\) per year, so

$$ F_1(1\mid z=1) =\frac{0.80}{1.40}\left(1-e^{-1.40}\right) \approx 0.5714(0.7534) \approx 0.4305. $$

The one-year cumulative-incidence risk ratio is \(0.4305/0.2528\approx1.703\), not 2. The corresponding risk difference is \(0.4305-0.2528=0.1777\), or about 17.77 percentage points. The transition hazard ratio compares instantaneous rates among those still in state 0; the cumulative-incidence ratio compares probabilities of reaching state 1 by one year, with state 2 competing.

State Occupation in an Illness-Death Model

The distinction is even more pronounced when the destination state is not absorbing. Consider three states: 0 is healthy, 1 is ill, and 2 is dead. The allowed transitions are \(0\to1\), \(0\to2\), and \(1\to2\), with constant rates \(a\), \(c\), and \(d\), respectively. Assume everyone starts in state 0, the rates are constant, and the transition intensities depend only on current state and the fixed covariate. This is a time-homogeneous Markov illness-death model. Death is absorbing.

To occupy state 1 at time \(t\), a person must enter state 1 at some time \(u\leq t\), then avoid the \(1\to2\) transition until \(t\). The probability of remaining in state 0 until \(u\) is \(e^{-(a+c)u}\), the rate of entering state 1 at \(u\) is \(a\), and the probability of remaining in state 1 from \(u\) to \(t\) is \(e^{-d(t-u)}\). Integrating these contributions gives the state-occupation probability.

$$ P_1(t) =\int_0^t e^{-(a+c)u}a e^{-d(t-u)}\,du =\frac{a\left(e^{-dt}-e^{-(a+c)t}\right)}{a+c-d}, \qquad a+c\ne d. $$

The expression includes both entry into illness and subsequent exit from illness. It therefore is not simply the cumulative incidence of \(0\to1\). If \(a+c=d\), the integral instead gives \(P_1(t)=at e^{-dt}\), the continuous limiting case. In either form, changing \(a\) changes the time of entry and the amount of time available for the \(1\to2\) transition.

Worked Example: A Twofold Entry Hazard and Illness Occupation

Use the same two covariate groups as above, and suppose the \(1\to2\) death hazard is \(d=0.50\) per year in both groups. For \(z=0\), \(a+c=0.40+0.60=1.00\); for \(z=1\), \(a+c=0.80+0.60=1.40\). Calculate the probability of occupying state 1 at one year.

For \(z=0\), substitute \(a=0.40\), \(a+c=1.00\), \(d=0.50\), and \(t=1\):

$$ P_1(1\mid z=0) =\frac{0.40(e^{-0.50}-e^{-1.00})}{1.00-0.50} =0.80(0.6065-0.3679) \approx 0.1909. $$

For \(z=1\), substitute \(a=0.80\), \(a+c=1.40\), \(d=0.50\), and \(t=1\):

$$ P_1(1\mid z=1) =\frac{0.80(e^{-0.50}-e^{-1.40})}{1.40-0.50} =\frac{0.80}{0.90}(0.6065-0.2466) \approx 0.3199. $$

The one-year state-occupation ratio is approximately \(0.3199/0.1909=1.676\), not 2. Unlike the previous example, the state-occupation ratio here is greater than the cumulative-incidence ratio because the transition \(1\to2\) allows people to leave illness before the one-year assessment. The size and even direction of the difference depend on the transition structure and the time horizon.

The stated values in the preceding example require care: applying the formula gives \(0.80/0.90\) multiplied by \(0.3599\), which is approximately \(0.3199\). For the baseline group, the same formula gives \(0.1909\). These figures illustrate that the illness-occupation ratio is about 1.676 under the specified rates. They do not imply a general ordering between a transition hazard ratio and an occupation ratio.

Worked Example: Overall Probability of Death by One Year

In the same illness-death model, state 0 occupation is \(P_0(t)=e^{-(a+c)t}\). Since the three states exhaust all possibilities, the probability of being dead at \(t\) is \(P_2(t)=1-P_0(t)-P_1(t)\). Use the one-year values above to compare death probabilities.

For \(z=0\), \(P_0(1)=e^{-1.00}\approx0.3679\) and \(P_1(1)\approx0.1909\). Therefore

$$ P_2(1\mid z=0) =1-0.3679-0.1909 \approx0.4412. $$

For \(z=1\), \(P_0(1)=e^{-1.40}\approx0.2466\) and \(P_1(1)\approx0.3199\). Therefore

$$ P_2(1\mid z=1) =1-0.2466-0.3199 \approx0.4335. $$

The one-year death-risk ratio is about \(0.4335/0.4412=0.983\), while the \(0\to1\) transition hazard ratio remains 2. In this illustration, the two death probabilities are quite similar, despite the higher entry rate into illness. This is possible because death can occur directly from state 0 or after occupying state 1, and changing the entry rate redistributes time across states. A coefficient for \(0\to1\) alone is not a death-risk ratio.

What the Calculations Do, and Do Not, Establish

The equations demonstrate a general point rather than a universal numerical relationship. The overall probability for a chosen endpoint depends on all transition rates along paths that can lead to, or away from, that endpoint. It also depends on the time horizon. A hazard ratio is an instantaneous, conditional comparison; a risk ratio is a ratio of probabilities accumulated by a specified time. A state-occupation ratio concerns presence in a state at that time, not necessarily whether that state was ever visited.

The assumptions used above are deliberately restrictive: constant rates, a fixed baseline covariate, the stated allowed transitions, and a Markov model in the illness-death example. Under these assumptions, the formulas follow directly from the transition process. In applications, transition hazards may vary over time, covariates may change, and history beyond the current state may matter. The transition-specific hazard interpretation still needs a well-defined risk set and model, while overall probabilities then require calculations appropriate to the fitted transition structure.

A useful way to avoid conflation is to identify the estimand before interpreting a coefficient:

  • Transition effect: the hazard ratio for a named transition, conditional on occupying its origin state.
  • Cumulative incidence effect: a comparison of probabilities of having experienced a specified event by a named time, accounting for competing transitions.
  • State-occupation effect: a comparison of probabilities of being in a specified state at a named time, accounting for entry and exit paths.

These are not interchangeable summaries. If the scientific question concerns one-year incidence or occupation, the fitted transition model must be translated into that probability and the resulting contrast calculated directly. The transition hazard ratio alone does not supply it.

Diagnostics, Alternatives, and Common Pitfalls

Before interpreting transition regression, check that the transition definitions and risk sets match the scientific question. For each modeled transition, confirm who is eligible to make it and how censoring and state history are represented. A proportional-hazards specification also requires the covariate effect on that transition’s hazard to be adequately represented as proportional over time; a single hazard ratio can be misleading if the effect varies materially over time. In a Markov model, check whether conditioning on current state and included covariates is a defensible description of transition rates. Dependence on time since entry or earlier history may require a richer model.

A common mistake is to read \(\exp(\beta_{01})=2\) as “twice the probability of illness” or “twice the one-year risk.” The reasoning is wrong because the coefficient compares instantaneous \(0\to1\) rates among people currently in state 0, not probabilities at a horizon. The correction is to specify the probability of interest, include the competing and subsequent transitions, and calculate that probability under each covariate value before forming a risk ratio or difference.

Transition-specific regression is useful when the scientific question concerns how a covariate relates to a particular move, and it allows different effects for different transitions. An alternative is to report model-based cumulative incidence or state-occupation probabilities at clinically meaningful times, with absolute differences or ratios as appropriate. These probability summaries answer different questions and depend on the whole transition model. Neither approach makes the other unnecessary; the appropriate summary is determined by the estimand, not by the form of a regression coefficient.

Finally, a mismatch between a transition hazard ratio and an overall risk ratio is not evidence that the regression is faulty. It is often the expected consequence of competing paths and transitions out of the destination state. Model diagnostics address whether the transition model is credible; they do not turn a conditional rate contrast into a marginal probability contrast.

Key takeaway: A transition hazard ratio is conditional on occupying the transition’s origin state. Cumulative incidence and state-occupation probabilities depend on the complete set of relevant transitions and the time horizon, so they must be calculated and compared separately.

Check Your Understanding

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

  1. In the competing-destinations example, why does the cumulative incidence of state 1 depend on \(\lambda_{02}\), even though \(\lambda_{02}\) is not the hazard for entering state 1?
  2. For the illness-death model, distinguish the probability of ever entering state 1 by time \(t\) from the probability of occupying state 1 at time \(t\).
  3. If a transition coefficient is \(\beta_{01}=\log(2)\), what is the transition hazard ratio, and what additional information is needed to calculate a one-year risk ratio?
  4. Why can a covariate double the \(0\to1\) hazard while having little association with one-year probability of death?

Solutions. (1) Transition \(0\to2\) removes people from state 0 before they can make transition \(0\to1\), so it changes the probability of remaining eligible for that transition. (2) Ever entering state 1 counts anyone who made \(0\to1\) by time \(t\), including those who later left state 1; occupying state 1 at time \(t\) counts only those still there then. (3) The hazard ratio is \(\exp\{\log(2)\}=2\). A one-year risk ratio also requires the other relevant transition rates and the model-based one-year probabilities for both covariate values. (4) Overall death probability combines direct death from state 0 and death after entering state 1; changes in the timing and route through states can offset or dilute the effect of one transition.

Clinical Biostats