Why Trace Causal Paths?
A biological system can contain far more processes than a particular QSP decision requires. The model boundary is not a catalogue of everything known about the disease or drug. It is a proposed scope for representing the causal processes needed to connect an intervention to the outcome that informs a decision. Causal path tracing makes that scope explicit: start at the intervention, follow plausible causal links, and ask which processes can transmit its effect to the decision-relevant outcome.
The outcome and decision context must be stated before drawing the boundary. For example, a pathway irrelevant to an efficacy endpoint may matter to a safety endpoint, and a mechanism important for one patient population may be inactive in another. The earlier material on decision questions and model output boundaries establishes why that context matters. Here, the focus is the causal diagram used to trace routes within it.
A causal diagram uses nodes for variables or processes and arrows for hypothesized causal influence. A direct path is a single arrow, such as \(D \rightarrow Y\). A mediated path includes one or more intervening nodes, such as \(D \rightarrow M \rightarrow Y\). Several paths may run in parallel. A node affected by the intervention but with no directed route onward to the specified outcome is a branch that does not transmit an effect to that outcome in the diagram.
The arrows are causal claims, not merely observed associations. A diagram organizes a proposed explanation; it does not establish that the links are true. Its use assumes that the intervention is sufficiently well defined, the relevant context and outcome are specified, the temporal ordering of causes is coherent, and the diagram includes the plausible causal routes needed for the question. If important causes or routes are omitted, a boundary drawn from the diagram may be inadequate.
From a Diagram to a Boundary
Begin at the intervention and trace forward along every plausible directed route that reaches the decision-relevant outcome. The mechanisms on those routes are candidates for representation. Then consider whether other variables are needed to specify how those routes operate—for example, a context variable that changes a mediator’s level or a cause of both a mediator and the outcome. Such a variable is not necessarily a drug mechanism, but omitting it may prevent the model from representing the intervention contrast appropriately.
Do not equate “not on a directed path from the intervention” with “always irrelevant.” A variable may be a common cause of a pathway component and the outcome, may define the population, or may be relevant to another decision outcome. Conversely, a biologically interesting process is not automatically required merely because it changes after treatment. The question is whether it can transmit an effect to the stated outcome under the causal structure and context being modeled.
State the intervention or contrast, the population and context, and the outcome used for the decision.
Separate causal hypotheses from associations and order the nodes consistently with the proposed biology.
Mark every plausible directed route from the intervention to the outcome, including parallel routes and intermediate processes.
Ask what other variables are needed to specify or interpret the routes in the stated population and intervention setting.
Include the required route processes; document branches excluded for this outcome and context rather than treating them as universally unimportant.
Why Path Coefficients Matter
A simple linear structural model illustrates how paths carry an intervention effect. Suppose \(M\) is a mediator and \(Y\) is the outcome, with equations \(M=aD+U\) and \(Y=bM+cD+V\), where \(U\) and \(V\) are fixed for the intervention contrast. Substituting the first equation into the second gives \(Y=abD+bU+cD+V\). Therefore, changing \(D\) by one unit changes the predicted \(Y\) by \(ab+c\): \(ab\) is the contribution along \(D \rightarrow M \rightarrow Y\), and \(c\) is the contribution along the direct arrow \(D \rightarrow Y\). This calculation relies on the stated linear equations and fixed context; it is not a general formula for nonlinear or confounded systems.
More generally, in a linear acyclic structural model, a directed path contributes the product of its link coefficients to the total intervention effect, and the total effect is the sum of contributions across directed paths. This follows by repeatedly substituting each node’s equation into the outcome equation: each substitution generates terms for the routes through that node, and terms reaching the outcome correspond to directed paths. The result is a way to check a specified model’s implications, not evidence that the diagram or coefficient values are correct.
Worked Example: A Mediated Effect With a Context Variable
Consider a hypothetical model of an intervention \(D\), a biological mediator \(M\), a decision-relevant score \(Y\), and a context variable \(U\). Let \(D=0\) denote control and \(D=1\) denote intervention. The proposed structural equations are \(M=8-2D+U\) and \(Y=100+3M-5U\). The diagram has \(D \rightarrow M \rightarrow Y\), \(U \rightarrow M\), and \(U \rightarrow Y\). Calculate the predicted scores and intervention contrast at each context value \(U=0\) and \(U=1\).
At \(U=0\), control gives \(M=8-2(0)+0=8\), so \(Y=100+3(8)-5(0)=124\). Intervention gives \(M=8-2(1)+0=6\), so \(Y=100+3(6)-5(0)=118\). The contrast is \(118-124=-6\) score units.
At \(U=1\), control gives \(M=8-2(0)+1=9\), so \(Y=100+3(9)-5(1)=100+27-5=122\). Intervention gives \(M=8-2(1)+1=7\), so \(Y=100+3(7)-5(1)=100+21-5=116\). The contrast is \(116-122=-6\) score units.
In both contexts the intervention contrast is \((-2)(3)=-6\) score units: the intervention changes \(M\) by \(-2\) units, and each unit of \(M\) changes \(Y\) by \(3\) score units. The direct \(U\)-to-\(Y\) term shifts the scores with context but cancels in the within-context contrast. In this example \(U\) is needed to predict absolute scores, though it does not alter the contrast. The equations and values are hypothetical and do not establish a biological mechanism.
Parallel Routes and Excluded Branches
Path tracing is especially useful when a drug may affect more than one process. If two plausible routes reach the outcome, excluding either route can change the predicted intervention effect. A branch that changes under treatment but does not lead to the selected outcome may be outside the boundary for that question. The conclusion is specific: if the branch becomes relevant to another endpoint, the boundary must be reconsidered.
Worked Example: Two Routes to an Outcome
In a hypothetical linear model, intervention \(D\) affects a target state \(T\) and a second process \(O\), with both influencing an inflammation score \(Y\). Let \(T=10-2D\), \(O=4-D\), and \(Y=50-3T-2O\). Assume these equations correctly describe the specified context and that no other routes affect the intervention contrast. Calculate the contrast and identify the route contributions.
Under control, \(D=0\), so \(T=10\), \(O=4\), and \(Y=50-3(10)-2(4)=12\). Under intervention, \(D=1\), so \(T=8\), \(O=3\), and \(Y=50-3(8)-2(3)=20\). Thus the intervention contrast is \(20-12=8\) inflammation-score units.
Along the target route, \(D\) changes \(T\) by \(-2\), and the coefficient from \(T\) to \(Y\) is \(-3\). Its contribution is \((-2)(-3)=6\) score units. Along the second route, \(D\) changes \(O\) by \(-1\), and the coefficient from \(O\) to \(Y\) is \(-2\). Its contribution is \((-1)(-2)=2\) score units. The contributions sum to \(6+2=8\), matching the direct calculation. A boundary that included only \(T\) would predict a contrast of \(6\), not \(8\), under these assumed equations.
Worked Example: A Treatment-Affected Branch Outside the Outcome Boundary
Suppose a hypothetical intervention changes an immune-cell marker \(B\), while a separate mediator \(M\) changes the selected efficacy outcome \(Y\). The proposed equations are \(B=5+2D\), \(M=10-3D\), and \(Y=40-2M\). The diagram contains \(D \rightarrow B\), \(D \rightarrow M \rightarrow Y\), and no arrow from \(B\) to \(Y\) or to \(M\). Compare the predicted efficacy outcome under control and intervention.
For \(D=0\), \(M=10\), so \(Y=40-2(10)=20\). For \(D=1\), \(M=7\), so \(Y=40-2(7)=26\). The intervention contrast is \(26-20=6\) outcome units. Although \(B\) changes from \(5\) to \(7\), it does not enter the equation for \(Y\); including it would not change this efficacy prediction under the stated diagram.
This exclusion is conditional on the question and causal structure. If \(B\) were itself a safety outcome, or if evidence supported \(B \rightarrow Y\), it would no longer be an irrelevant branch for the corresponding model purpose. A missing arrow cannot be justified merely because the available data do not show an association.
Diagnostics, Alternatives, and Limitations
A path-based boundary is a hypothesis to scrutinize, not a guarantee of adequacy. Warning signs include a claimed intervention effect with no route to the modeled outcome, a route that depends on an unsupported causal link, or a process that changes after treatment but has been excluded without considering whether it affects the chosen endpoint. Also check the direction and timing of each proposed link. An observed association between two nodes does not establish which causes which, and a biomarker can track an outcome without mediating the intervention’s effect.
The boundary should also be checked against the model’s intended use. A single-outcome boundary may be narrower than one needed for a joint efficacy-and-safety decision. A population-specific route may fail in another population. If feedback is important, a simple acyclic diagram at one time point may not represent it; time-indexed states or another explicit representation may be needed. These are reasons to revisit the proposed structure, not reasons to include all biology by default.
A broad mechanistic network can represent many candidate routes, but it increases the number of assumptions and quantities that must be supported. A reduced mechanistic model focuses on routes needed for the decision, while retaining causal interpretation for those routes. An empirical exposure–response fit may summarize a relationship in the data without representing the intervening biological pathways; it can be useful for interpolation within its supported setting, but may not answer a counterfactual that depends on changing a mechanism. These approaches answer different questions, and no diagram alone determines which is adequate.
Path tracing has limits. The diagram may omit a route, misstate a causal direction, or oversimplify interactions and context dependence. A process can influence a path without appearing as a mediator on the path itself, and nonlinear effects can make a single set of path contributions inadequate. Thus, use the diagram to make scope and assumptions inspectable, then compare its implied predictions with biological knowledge and available evidence. A defensible boundary is explicit about what is represented, what is excluded, and for which outcome and context that choice is being made.
Check Your Understanding
Answer each question, then compare with the solutions.
- A diagram contains \(D \rightarrow A \rightarrow Y\) and \(D \rightarrow B\), with no route from \(B\) to \(Y\). Is \(B\) on a directed path from \(D\) to \(Y\)?
- In the linear model \(M=5- D\), \(Y=20+4M+2D\), calculate the predicted \(Y\) at \(D=0\) and \(D=1\), and give the intervention contrast.
- Why does an observed association between a biomarker and an outcome not by itself establish that the biomarker mediates the intervention effect?
- A process has no route to the efficacy outcome but is itself a decision-relevant safety endpoint. Should it automatically be excluded from all models? Explain.
- In a diagram with two directed routes from \(D\) to \(Y\), what is the risk of omitting one route from a model intended to predict the intervention contrast?
Solutions. (1) No. The route ends at \(B\) and does not continue to \(Y\) in the stated diagram. (2) At \(D=0\), \(M=5\) and \(Y=20+4(5)+2(0)=40\). At \(D=1\), \(M=4\) and \(Y=20+4(4)+2(1)=38\). The contrast is \(38-40=-2\) outcome units; the mediated contribution is \((-1)(4)=-4\), and the direct contribution is \(+2\), summing to \(-2\). (3) Association alone does not establish causal direction or rule out common causes; mediation requires a supported causal route through the biomarker. (4) No. Its lack of a path to efficacy does not make it irrelevant to a safety decision, so the model boundary depends on the endpoint and use. (5) The omitted route may contribute to the total intervention effect, leading to a biased or incomplete prediction under the proposed causal structure.