Tutorials › Joint Modeling of Longitudinal and Time-to-Event Data › A Causal Diagram for Biomarker, Dropout, and Event

Why Joint Modeling · Tutorial 18 of 697

A Causal Diagram for Biomarker, Dropout, and Event

Learn to represent the assumed pathways among underlying health, repeated biomarker measurements, dropout, and terminal events in a time-ordered directed acyclic graph.

Beginner 10 min read

What You'll Learn

  • Define time-indexed nodes for latent health, biomarkers, dropout, event occurrence, and biomarker observation.
  • Choose arrows that reflect temporal order and the scientific assumptions being made.
  • Distinguish a biomarker that measures latent health from one that may affect later health or event risk.
  • Use a graph to explain how dropout and terminal events can be related without causing one another.
  • Check a proposed diagram for cycles, ambiguous timing, and assumptions that the data alone cannot establish.

Why Draw the Data-Generating Story?

A repeated biomarker is not observed in isolation. It reflects a person’s underlying health at the time of measurement, while that health may also help determine who remains in follow-up and who experiences a terminal event. A directed acyclic graph (DAG) makes these proposed relationships explicit. It helps separate a biomarker’s role as a measurement of health from a possible role as a cause of later health or event risk.

The diagram is a statement of assumptions, not a finding produced automatically by the data. In particular, an arrow does not show that an effect has been estimated, and the absence of an arrow is also an assumption. The diagram is useful because it forces the analyst to specify which variables precede others and which pathways could connect them.

Definition: A directed acyclic graph is a set of nodes representing variables and directed arrows representing proposed direct dependencies. “Acyclic” means that following arrows forward cannot lead back to the same node. In a time-indexed graph, an arrow \(A \to B\) states that \(A\) is a direct parent of \(B\) in the proposed data-generating structure; it does not, by itself, establish a causal effect from observational data.

Define Nodes and Order Them in Time

Let \(i\) index individuals and let \(k\) index scheduled biomarker visits at times \(t_k\). Use the following nodes:

  • \(H_{ik}\): latent health at visit time \(t_k\), where larger values are defined here to mean worse health.
  • \(Y_{ik}\): the biomarker value at \(t_k\), which may be observed only for some individuals.
  • \(D_{i,k+1}\): an indicator that the individual drops out between \(t_k\) and \(t_{k+1}\), before a terminal event is known to have occurred.
  • \(E_{i,k+1}\): an indicator that the terminal event occurs in the interval after \(t_k\) and by \(t_{k+1}\).
  • \(R_{i,k+1}\): an indicator that a biomarker is observed at \(t_{k+1}\).

The timing matters. \(H_{ik}\) precedes the next-interval indicators \(D_{i,k+1}\) and \(E_{i,k+1}\). If dropout or the terminal event has occurred, the next scheduled biomarker is not observed, so the observation indicator \(R_{i,k+1}\) depends on those prior follow-up outcomes. A simple version of the graph has these arrows:

$$ H_{ik} \to Y_{ik}, \qquad H_{ik} \to H_{i,k+1}, \qquad H_{ik} \to D_{i,k+1}, \qquad H_{ik} \to E_{i,k+1}, \qquad D_{i,k+1} \to R_{i,k+1}, \qquad E_{i,k+1} \to R_{i,k+1}. $$

The arrow \(H_{ik} \to Y_{ik}\) represents the biomarker as a measurement influenced by current latent health. The arrow \(H_{ik} \to H_{i,k+1}\) represents persistence in health over time. The arrows from current health to subsequent dropout and event represent the possibility that worsening health makes either outcome more likely. The final arrows represent observation availability, not biological effects: dropout and a terminal event prevent a later scheduled measurement from being recorded.

This is a deliberately compact graph. It omits baseline factors and other causes that may matter in a particular application. It also does not assert that dropout causes the terminal event. A shared cause, such as latent health, can make dropout and event occurrence dependent even when neither causes the other. The graph must be expanded if the scientific setting calls for additional nodes or arrows.

Key distinction: \(D_{i,k+1}\) describes loss from observation; \(E_{i,k+1}\) describes the terminal event itself. A dropout is not automatically a terminal event, and an unobserved event after dropout is not evidence that no event occurred. Keep the underlying event process distinct from whether its status is observed.

Constructing a Coherent Diagram

A reliable construction begins with time rather than with a search for every conceivable association. First list the variables and their time points. Then add arrows only when a direct dependence is part of the proposed scientific story. Finally, add observation indicators so that the graph distinguishes a value that exists from a value that is recorded.

1
Place nodes at their relevant times.
Separate current health and the current biomarker from next-interval dropout, event occurrence, and measurement availability.
2
State the substantive pathways.
For example, latent health may influence both the biomarker and the probability of dropout or a terminal event.
3
Represent observation separately.
Use an observation indicator to show why a later biomarker may be missing. Do not treat a missing value as if it were a measured biomarker value.
4
Check ordering and scope.
Confirm that arrows point forward in time and that the diagram includes the causes needed for the question being asked. Record important omitted causes as limitations.

The biomarker’s position requires particular care. If \(Y_{ik}\) is only a measurement of latent health, the association between \(Y_{ik}\) and a later event can arise because both depend on health. A direct arrow from \(Y_{ik}\) to \(E_{i,k+1}\) would instead assert an additional pathway: changing the biomarker itself, while holding its causes fixed, could affect the event. That claim needs scientific justification. If the biomarker could affect subsequent health, a possible pathway is \(Y_{ik} \to H_{i,k+1} \to E_{i,k+1}\). Whether that arrow belongs in the graph is a substantive assumption, not a decision made from the biomarker’s predictive strength alone.

Worked Example: A Marker of Worsening Health

Hypothetical setting. A cohort has scheduled blood tests every three months. Larger values of a particular marker are thought to indicate worse underlying health. Patients in worse health may be more likely to stop attending visits and more likely to experience a terminal event.

At month 3, represent underlying health by \(H_{i1}\) and the blood test by \(Y_{i1}\). Represent dropout and event occurrence during the following interval by \(D_{i2}\) and \(E_{i2}\). The core assumptions give \(H_{i1} \to Y_{i1}\), \(H_{i1} \to D_{i2}\), and \(H_{i1} \to E_{i2}\). If either dropout or the event prevents a month-6 blood test, include arrows from \(D_{i2}\) and \(E_{i2}\) to \(R_{i2}\).

The graph does not require an arrow \(Y_{i1} \to E_{i2}\). Under the stated measurement-only interpretation, the marker and event may be associated because both are influenced by \(H_{i1}\). Adding a direct arrow would change the scientific story by proposing a marker-to-event pathway beyond the shared influence of health. The assumptions are coherent with the timing because the month-3 health and marker precede the next-interval outcomes.

What the Graph Implies About Dependence

A DAG can show why two outcomes may be associated without a direct arrow between them. In the compact graph, \(H_{ik}\) is a common cause of \(D_{i,k+1}\) and \(E_{i,k+1}\). If the graph includes no other shared causes and no direct dependence between those outcomes, it represents them as conditionally independent given \(H_{ik}\). In symbols, the proposed model has \(P(D_{i,k+1}, E_{i,k+1}\mid H_{ik}) = P(D_{i,k+1}\mid H_{ik})P(E_{i,k+1}\mid H_{ik})\). This conditional independence is a model assumption, not a general fact about dropout and events.

The following calculation shows how marginal dependence can still arise. For simplicity, let current health have two categories: \(H=1\) means worse health and \(H=0\) means better health. Suppose the graph’s conditional independence assumption holds for dropout and a later event given \(H\).

Worked Example: Shared Health Can Link Dropout and Event

Hypothetical probabilities. Suppose \(P(H=1)=0.30\), \(P(D=1\mid H=1)=0.40\), and \(P(D=1\mid H=0)=0.10\). Suppose the later event probabilities are \(P(E=1\mid H=1)=0.20\) and \(P(E=1\mid H=0)=0.05\). Assume \(D\) and \(E\) are conditionally independent given \(H\).

The marginal dropout probability is

$$ P(D=1)=0.30(0.40)+0.70(0.10)=0.12+0.07=0.19. $$

The marginal event probability is

$$ P(E=1)=0.30(0.20)+0.70(0.05)=0.06+0.035=0.095. $$

Using conditional independence within each health category, the probability of both outcomes is

$$ P(D=1,E=1)=0.30(0.40)(0.20)+0.70(0.10)(0.05) =0.024+0.0035=0.0275. $$

If dropout and event were marginally independent, their joint probability would be \(0.19(0.095)=0.01805\), not \(0.0275\). In this hypothetical model, their joint occurrence is more common than that marginal-independence calculation suggests because worse health raises the probability of each. The calculation supports a shared-cause interpretation; it does not show that dropout causes the event.

The probability calculation is an instance of the graph’s factorization. With \(H\) as the only shared parent and conditional independence given \(H\), sum the product of the conditional probabilities over the possible health states:

$$ P(D,E) =\sum_{h \in \{0,1\}}P(H=h)P(D\mid H=h)P(E\mid H=h). $$

This follows by conditioning on \(H\), applying the assumed conditional independence within each health state, and adding the two mutually exclusive cases \(H=0\) and \(H=1\). If additional common causes or direct links are scientifically plausible, this factorization must be changed to match the expanded graph.

Worked Example: Calculating Who Reaches the Next Visit

Continue with the same hypothetical probabilities. Suppose an individual reaches the next scheduled visit only if neither dropout nor the terminal event has occurred, and assume conditional independence of \(D\) and \(E\) given health. Among those with worse health, the probability of reaching the visit is

$$ P(R=1\mid H=1)=(1-0.40)(1-0.20)=0.60(0.80)=0.48. $$

Among those with better health, it is

$$ P(R=1\mid H=0)=(1-0.10)(1-0.05)=0.90(0.95)=0.855. $$

Averaging over the health categories gives

$$ P(R=1)=0.30(0.48)+0.70(0.855) =0.144+0.5985=0.7425. $$

Thus, under these hypothetical assumptions, 74.25% of the cohort reaches the next measurement occasion. The probability differs by health category, so the set with an observed next biomarker is not a random subset with respect to \(H\). This is a consequence of the specified graph and probabilities, not a claim that the numerical values apply to an actual study.

Diagnostics, Alternatives, and Limits

A graph is a disciplined description of assumptions, but it cannot verify unobserved causes by itself. Use the following checks before relying on one:

  • Check temporal order. A later event cannot cause an earlier biomarker. If a proposed cycle appears, separate the variables by time, such as \(H_{ik}\) and \(H_{i,k+1}\), and reconsider the ordering.
  • Distinguish event occurrence from event observation. A terminal event can stop later measurements; dropout can stop event ascertainment. The recorded event indicator may therefore need a separate observation or ascertainment node.
  • Ask what each arrow means. An arrow from latent health to a measured biomarker is different from an arrow asserting that changing the biomarker would change future health.
  • Look for implausible omissions. If treatment, baseline severity, or another factor plausibly affects several nodes, omitting it may make the displayed conditional independencies misleading.
  • Do not equate fit with truth. Observed data can be inconsistent with some implications of a proposed graph, but agreement with observed patterns cannot establish that all unmeasured causes or arrows have been captured.

A DAG and a joint statistical model serve different purposes. The DAG clarifies the assumed structure and time ordering; it does not specify probability distributions, estimate a longitudinal trajectory, or provide an event-risk estimate. A joint model gives a statistical framework for the longitudinal and event processes, with parameters and distributional assumptions. Conversely, analyzing each process separately can be simpler, but may fail to represent dependencies suggested by a shared-health graph. The graph helps explain why dependence may matter; choosing and fitting a statistical model requires additional assumptions.

A common mistake is to read every association as a direct causal arrow. For example, “the biomarker predicts the event, therefore the biomarker causes the event” confuses predictive information with a causal pathway. The correction is to ask whether the biomarker is a measurement of latent health, a cause of subsequent health, or both, and to draw only the arrows justified by the scientific question. A second mistake is to collapse dropout and the terminal event into one generic “missingness” node. They have different meanings and may require different observation mechanisms.

Finally, a graph is not a substitute for careful data definitions. The timing of visits, the definition of dropout, and whether an event occurring after dropout is ascertained all affect the appropriate nodes and arrows. Write these definitions down before interpreting the diagram.

Check Your Understanding

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

  1. If latent health affects both the current biomarker and next-interval event risk, which arrows belong in the graph even if there is no assumed biomarker effect on the event?
  2. In a simplified model, \(P(H=1)=0.20\), \(P(D=1\mid H=1)=0.50\), and \(P(D=1\mid H=0)=0.10\). Calculate \(P(D=1)\).
  3. Why should dropout and the terminal event be represented as distinct nodes?
  4. A proposed diagram contains \(H_{k+1}\to H_k\) and \(H_k\to H_{k+1}\). What is the problem, and how could the health process be represented without a cycle?

Solutions. (1) Include \(H_k\to Y_k\) and \(H_k\to E_{k+1}\). Do not add \(Y_k\to E_{k+1}\) unless a direct biomarker pathway is part of the assumed scientific story. (2) Average over the two health categories: \(P(D=1)=0.20(0.50)+0.80(0.10)=0.10+0.08=0.18\). (3) Dropout means follow-up is lost; the terminal event is an outcome. Dropout may prevent observing a later event without being the event itself. (4) The two arrows create a directed cycle. Represent temporal persistence as \(H_k\to H_{k+1}\), with no arrow back from the later health state to the earlier one.

Clinical Biostats