Tutorials › Competing Risks and Multi-State Models › Validating Direct Transitions And Indirect Paths

Transition-Specific Prediction Validation · Tutorial 405 of 699

Validating Direct Transitions And Indirect Paths

Learn to distinguish direct-edge predictions from indirect first-passage predictions, and validate each against the outcome it actually forecasts.

Advanced 12 min read

What You'll Learn

  • Define direct and indirect routes to a target state from a specified starting state.
  • Distinguish first reaching a state from occupying it at a fixed horizon.
  • Match predicted probabilities to route-specific validation outcomes.
  • Check whether direct and total first-passage predictions are coherent.
  • Compare calibration and Brier scores for different route targets.
  • Recognize how competing transitions, censoring, and endpoint definitions can mislead validation.

Why the Route Matters

A prediction that a patient will be in state B at one year does not say how the patient will get there. The patient might move directly from state A to B, or move from A to an intermediate state C and then to B. Those are different outcomes, and a model can predict them with different accuracy. Validation must therefore start by identifying the exact event represented by each predicted probability.

A further distinction is essential: first reaching B by a horizon is not always the same as being in B at that horizon. Someone can enter B and then leave it before the horizon. The earlier material on transition probabilities versus state-occupation probabilities distinguishes these endpoint quantities. Here, route-specific validation uses first passage to B; a state-occupation prediction must instead be checked against occupancy at the specified time.

Definition: Suppose an individual starts in state A at time \(t_0\), and let \(\tau\) be the prediction horizon measured from \(t_0\). Let \(T_1\) be the time of the first transition after \(t_0\), and let \(J_1\) be its destination. Let \(T_B\) be the first time the process enters target state B. Define the direct-route outcome \(D_{AB}(\tau)=1\) if \(T_1\leq t_0+\tau\) and \(J_1=B\), and define the first-passage outcome \(R_{AB}(\tau)=1\) if \(T_B\leq t_0+\tau\). The indirect-route outcome is \(I_{AB}(\tau)=R_{AB}(\tau)-D_{AB}(\tau)\). Thus, a route is indirect when B is first reached by the horizon but the first transition from A was not directly to B.

This definition describes the route from the starting state to the first visit to B. For example, A→C→A→B is classed as an indirect route: the first transition was to C, even though the final transition into B was from A. This convention makes the route outcomes mutually exclusive and exhaustive for whether B was first reached by the horizon: direct, indirect, or not reached. Another edge-specific target—such as whether the edge A→B is ever traversed after any number of earlier transitions—is possible, but it is a different estimand and needs a prediction for that event.

Predictions and Their Matching Outcomes

Let \(q_D(\tau)\) denote the predicted probability of the direct-route outcome and \(q_R(\tau)\) the predicted probability of first reaching B by the horizon. The corresponding indirect-route probability is \(q_I(\tau)=q_R(\tau)-q_D(\tau)\), provided both predictions refer to the same starting state, population, and horizon and use the route definitions above. A coherent set of predictions must satisfy \(0\leq q_D(\tau)\leq q_R(\tau)\leq 1\). If the model predicts the three route categories directly, their probabilities must sum to one after including the probability of not reaching B.

A transition matrix that predicts the probability of occupying B at \(t_0+\tau\) cannot automatically be substituted for \(q_R(\tau)\). It may count people who entered B and then left, and it may omit people who reached B but were elsewhere at the horizon. Unless B is absorbing, or the model separately supplies the first-passage probability, validate occupancy predictions using the observed occupancy outcome—not a route-specific first-passage outcome.

The outcomes for direct and indirect routes are competing categories in a route classification, not independent binary events. A person cannot have both \(D_{AB}=1\) and \(I_{AB}=1\) under these definitions. Likewise, a transition from A to a different state is not censoring for the direct first-transition target: it is a competing outcome. Treating such transitions as ordinary loss to follow-up would change the target probability.

Worked Example: Separating Routes From Occupancy

Consider a hypothetical validation cohort of 100 people starting in A, followed completely through a one-year horizon. By one year, 18 first moved directly from A to B, 16 first reached B after an initial move to another state, and 66 had not reached B. Of the 34 who reached B, four later left B before the one-year assessment.

The observed direct-route proportion is \(18/100=0.18\), and the indirect-route proportion is \(16/100=0.16\). The observed first-passage proportion is therefore \(34/100=0.34\), with \(0.18+0.16=0.34\). But the one-year state-occupation proportion is \(30/100=0.30\), because four people who had reached B were no longer there at the assessment.

Suppose the model predicted \(q_D=0.20\) and \(q_R=0.35\). These are coherent because \(0.20\leq0.35\), and the implied indirect-route prediction is \(q_I=0.35-0.20=0.15\). The observed direct and indirect proportions are 0.18 and 0.16, respectively. The comparison with 0.35 must use first-passage outcomes; comparing it with the occupancy proportion of 0.30 would assess a different target.

Route-Specific Calibration and Scores

For a predicted direct-route probability, assess calibration against the binary outcome \(D_{AB}(\tau)\). For a predicted first-passage probability, use \(R_{AB}(\tau)\). For indirect-route calibration, use \(I_{AB}(\tau)\) and a prediction for that same route. Calibration compares predicted probabilities with observed event frequencies among individuals receiving similar predictions; calibration-in-the-large compares the average prediction with the overall observed event proportion. A good result for one route does not establish good calibration for another.

The Brier score for a binary route outcome \(Y_i\) and its predicted probability \(q_i\) is the mean squared prediction error, \(n^{-1}\sum_{i=1}^{n}(q_i-Y_i)^2\). Smaller values indicate better overall probabilistic accuracy for that outcome and cohort, but they combine aspects of calibration and discrimination. Crucially, the outcome must match the prediction: using first-passage outcomes to score occupancy predictions, or vice versa, gives a number but not a valid assessment of the stated target.

Worked Example: Route-Specific Brier Scores

In a hypothetical validation sample of 10 individuals, consider two groups of five. In group 1, the predicted direct probability is 0.20 and the predicted first-passage probability is 0.40. One person has a direct route, one has an indirect route, and three do not reach B. In group 2, the corresponding predictions are 0.40 and 0.70; two people have a direct route, two an indirect route, and one does not reach B. Predictions are constant within each group.

For the direct-route score, group 1 has one outcome \(D=1\) and four outcomes \(D=0\). Its mean squared error is \((1/5)(1-0.20)^2+(4/5)(0-0.20)^2=(1/5)(0.64)+(4/5)(0.04)=0.16\). Group 2 has two direct outcomes and three non-direct outcomes, giving \((2/5)(1-0.40)^2+(3/5)(0-0.40)^2=(2/5)(0.36)+(3/5)(0.16)=0.24\). The pooled direct-route Brier score is \((5(0.16)+5(0.24))/10=0.20\).

For first passage, group 1 has two people reaching B and three not reaching it. Its score is \((2/5)(1-0.40)^2+(3/5)(0-0.40)^2=0.24\). Group 2 has four reaching B and one not reaching it, so its score is \((4/5)(1-0.70)^2+(1/5)(0-0.70)^2=(4/5)(0.09)+(1/5)(0.49)=0.17\). The pooled first-passage Brier score is \((5(0.24)+5(0.17))/10=0.205\). These two scores concern different outcomes; their small numerical difference does not establish that one route is inherently easier to predict.

Assumptions, Checks, and a Short Justification

Route-specific validation requires predictions and observed outcomes to refer to the same origin state, target state, horizon, and population. The state definitions and event times must be recorded consistently, and the route classification must follow the stated first-passage rule. If the prediction is conditional on covariates or on a landmark history, validation must use the corresponding eligible individuals and prediction time. The observed process need not be Markov for these outcome definitions; however, a model’s own method for generating probabilities may impose assumptions that should be checked separately.

Complete follow-up to the horizon makes the route outcomes directly observable. With right censoring before the horizon, a person censored before a route is observed cannot simply be classified as “did not reach B.” Use an appropriate censoring-adjusted validation method, such as inverse-probability-of-censoring weighting, under its stated assumptions. In particular, the censoring model must support estimation of the relevant route outcome over the horizon; sparse support or highly variable weights can make the assessment unstable. A competing transition is part of the process, not a censoring event.

The key probability identity follows from the outcome definitions. If \(R_{AB}=1\), then either the first transition was directly to B or the route to the first visit to B was indirect; these possibilities are disjoint. If \(R_{AB}=0\), neither route occurred. Therefore \(R_{AB}=D_{AB}+I_{AB}\) for every individual. Taking expectations gives \(P(R_{AB}=1)=P(D_{AB}=1)+P(I_{AB}=1)\), so coherent predictions for these probabilities must obey \(q_R=q_D+q_I\). This is a route-level coherence check, not a substitute for evaluating how well each prediction matches observed outcomes.

1
Fix the target.
Specify the starting state, target state, horizon, and whether the prediction concerns a first transition, a route to first entry, first passage, or occupancy at the horizon.
2
Classify observed routes.
Apply the same rule to each observed trajectory. Keep direct, indirect, and not-reached outcomes distinct; treat competing transitions as outcomes where appropriate.
3
Check compatibility.
Confirm that predicted and observed targets match, and check probability bounds such as \(q_D\leq q_R\). Investigate any violation before calculating performance metrics.
4
Assess each route separately.
Use route-matched calibration summaries and scores. Account for censoring when follow-up does not establish outcomes through the full horizon.

Worked Example: Detecting an Incoherent Route Prediction

A hypothetical model gives a person starting in A a 0.28 probability of the direct route to B by six months, but only a 0.22 probability of reaching B by six months. Under the definitions above, every direct route is also a first passage to B. Hence \(D_{AB}(6\text{ months})\leq R_{AB}(6\text{ months})\), and the corresponding probabilities must satisfy \(q_D\leq q_R\).

Here \(0.28>0.22\), so the implied indirect probability would be \(q_I=0.22-0.28=-0.06\), which is impossible. Do not interpret the two predictions as separately validated estimates or silently truncate the indirect probability to zero. First check for mismatched horizons, different starting populations, inconsistent route definitions, or a prediction-calculation error. If the predictions are intended to describe the same target and cohort, their construction must be reconciled before validation proceeds.

Common Pitfalls and Alternatives

  • Scoring the wrong outcome. A first-passage prediction is not validated by recording only who occupies B at the horizon. Use the first-entry outcome for first-passage predictions and the horizon state for occupancy predictions.
  • Calling every A-to-B edge traversal a direct route. Under the definition here, direct means the first transition from the specified starting state goes to B. If the desired target is any later traversal of A→B, define and predict that event explicitly.
  • Discarding competing transitions. A first transition to C is not ordinary censoring when assessing whether the first transition was A→B; it rules out that direct first-transition outcome.
  • Inferring route performance from total performance. A model can predict total first passage well while overpredicting direct routes and underpredicting indirect routes. Report the route-specific assessments needed for the scientific question.

Route-specific first-passage validation complements, rather than replaces, validation of transition probabilities and state-occupation probabilities. It is particularly useful when the route itself has clinical or operational meaning—for example, whether a patient reaches a complication directly or after an intermediate clinical state. Its limitation is that route categories can be uncommon, making observed frequencies and scores imprecise; a good overall score can also conceal poor calibration in a particular route or subgroup. When follow-up is incomplete, route comparisons additionally depend on appropriate handling of censoring.

Check Your Understanding

Answer using the first-passage route definitions in this tutorial.

  1. Starting in A, an individual moves A→C→B before the horizon. Is this a direct or indirect route to B?
  2. Why might the proportion occupying B at one year differ from the proportion that first reached B by one year?
  3. A model predicts direct-route probability 0.12 and first-passage probability 0.31. What indirect-route probability is implied?
  4. For a direct first-transition target, how should a first transition from A to C be classified: direct event, ordinary censoring, or competing outcome?

Solutions. (1) It is indirect because the first transition from A was to C. (2) Someone may enter B and then leave before the one-year assessment, so first passage does not guarantee occupancy at that time. (3) The implied indirect probability is \(0.31-0.12=0.19\). (4) It is a competing outcome: it prevents the first transition from being A→B, rather than representing loss to follow-up.

Clinical Biostats