Tutorials › Joint Modeling of Longitudinal and Time-to-Event Data › Unobserved Health as a Confounder of Biomarker-Event Associations

Why Joint Modeling · Tutorial 21 of 697

Unobserved Health as a Confounder of Biomarker-Event Associations

Learn how unobserved health can make a repeated biomarker appear associated with event risk, and how that association differs from a causal biomarker effect.

Beginner 11 min read

What You'll Learn

  • Describe the shared-cause pathway from latent disease severity to biomarker history and event risk.
  • Calculate how an elevated biomarker changes the probability of unobserved high severity.
  • Show how repeated measurements can strengthen a biomarker’s prognostic association without establishing causation.
  • Distinguish a crude biomarker-event association from a comparison adjusted for severity.
  • Explain how the severity composition of event-free participants can change a marginal hazard over time.
  • Recognize why adjustment or a joint model cannot automatically eliminate confounding by unobserved health.

Why a Biomarker–Event Association Can Mislead

A repeated biomarker may predict an event because it reflects worsening underlying health, not because changing the biomarker itself would change event risk. If latent disease severity affects both the biomarker trajectory and the event process, it is a common cause of the observed association. The biomarker can be useful for prediction while its association with the event does not identify a causal effect.

This distinction matters in longitudinal studies. Measurements taken repeatedly may track a person’s underlying condition more closely than a single baseline value does. But a more informative measure of severity is still not necessarily a cause of the event. The question is whether the biomarker contributes to event risk beyond the health state it reflects.

Definition: A latent confounder is an unobserved variable that influences both an exposure or predictor and an outcome. Here, latent disease severity \(U_i\) influences participant \(i\)’s repeated biomarker values \(M_{ij}\) and event time \(T_i\). A biomarker–event association can therefore arise along the pathway \(M_{ij} \leftarrow U_i \rightarrow T_i\), even if the biomarker has no causal effect on the event.

A simple data-generating model makes the idea explicit. Suppose \(U_i\) is a person’s underlying disease severity, which is not fully observed. Repeated biomarker values may follow

$$ M_{ij} = m(t_{ij}) + aU_i + \epsilon_{ij}, $$

where \(m(t_{ij})\) describes a common time pattern, \(a\) describes how severity shifts the marker, and \(\epsilon_{ij}\) is measurement or within-person variation. Suppose event risk follows a separate model in which severity affects the event, but the biomarker itself does not:

$$ \lambda_i(t \mid U_i) = \lambda_0(t)\exp(\gamma U_i). $$

In this illustration, a higher biomarker may be associated with a higher event rate because both share \(U_i\) as a cause. The model does not say that a biomarker intervention would change the event rate. It describes an association produced by shared latent health.

A Worked Binary Example of Confounding

Consider a hypothetical population followed for one year. Let \(U=1\) denote high latent disease severity and \(U=0\) low severity. Let \(X=1\) indicate an elevated biomarker at the start of follow-up, and let \(Y=1\) indicate an event during the year. Assume that severity causes both biomarker elevation and event risk, while the biomarker has no effect on the event in this example.

The assumptions and probabilities are:

  • \(P(U=1)=0.40\).
  • \(P(X=1\mid U=1)=0.80\) and \(P(X=1\mid U=0)=0.20\).
  • \(P(Y=1\mid U=1,X)=0.30\) and \(P(Y=1\mid U=0,X)=0.05\), for either value of \(X\).
  • Conditional on \(U\), the biomarker and event are independent in the data-generating model.

Worked Example: An Apparent Event-Risk Difference

Calculate the event risk among participants with elevated and non-elevated biomarkers. Then compare those risks, noting that the model gives the biomarker no causal effect.

Elevated biomarker. Bayes’ rule gives the proportion with high severity among those with \(X=1\):

$$ P(U=1\mid X=1) =\frac{P(X=1\mid U=1)P(U=1)} {P(X=1\mid U=1)P(U=1)+P(X=1\mid U=0)P(U=0)} =\frac{0.80(0.40)}{0.80(0.40)+0.20(0.60)} =\frac{0.32}{0.44} =0.7273. $$

Thus, the one-year event risk for \(X=1\) is the severity-weighted average

$$ P(Y=1\mid X=1) =0.7273(0.30)+0.2727(0.05) =0.2318. $$

Non-elevated biomarker. The corresponding high-severity proportion is

$$ P(U=1\mid X=0) =\frac{0.20(0.40)}{0.20(0.40)+0.80(0.60)} =\frac{0.08}{0.56} =0.1429. $$

Therefore,

$$ P(Y=1\mid X=0) =0.1429(0.30)+0.8571(0.05) =0.0857. $$

The crude one-year risk difference is \(0.2318-0.0857=0.1461\), or about 14.6 percentage points. Yet within either severity group, the event risk is identical for \(X=1\) and \(X=0\): 0.30 for high severity and 0.05 for low severity. In this hypothetical model, the crude association is entirely due to the different severity distributions in the two biomarker groups.

Why Repeated Measurements Can Strengthen the Association

Repeated measurements can provide more information about latent severity. That can make a biomarker history a stronger predictor of the event, even when no biomarker value causes the event. The following calculation uses the same hypothetical probabilities as above. Assume two biomarker readings are conditionally independent given \(U\), with each reading having probability 0.80 of being elevated if \(U=1\), and 0.20 if \(U=0\).

Worked Example: Two Elevated Readings as Evidence of Severity

Find the one-year event risk among participants whose two readings are both elevated. This is a calculation about prediction from a repeated history, not about the causal effect of raising the marker.

The probability of two elevated readings is \(0.80^2=0.64\) for high severity and \(0.20^2=0.04\) for low severity. Thus,

$$ P(U=1\mid X_1=1,X_2=1) =\frac{0.40(0.64)}{0.40(0.64)+0.60(0.04)} =\frac{0.256}{0.280} =0.9143. $$

The event risk for this repeated-measurement pattern is

$$ P(Y=1\mid X_1=1,X_2=1) =0.9143(0.30)+0.0857(0.05) =0.2786. $$

For comparison, the population event risk is \(0.40(0.30)+0.60(0.05)=0.1500\). A history of two elevated readings identifies a group with higher risk because it enriches the group for high latent severity. This improvement in prognostic information does not establish that an intervention that lowers the biomarker would lower event risk.

Adjustment: What It Can and Cannot Show

If severity were observed without error in this example, comparing biomarker groups within severity strata would remove the shared-cause association: within \(U=1\), both biomarker groups have risk 0.30; within \(U=0\), both have risk 0.05. Standardizing both biomarker groups to the same severity distribution also gives equal risks. For example, using the population severity proportions, each standardized risk is

$$ 0.40(0.30)+0.60(0.05)=0.1500. $$

This calculation relies on a simplified model with a measured, correctly classified confounder, overlap in biomarker groups at each severity level, and no other source of bias. In practice, disease severity is often only partly recorded. Adjusting for a noisy proxy may leave residual confounding because people with the same recorded proxy value can still differ in their underlying health.

A baseline severity measure may also fail to represent health as it evolves. If a repeated biomarker reflects changes in latent severity over follow-up, adjusting only for baseline information does not necessarily make the evolving biomarker independent of the evolving health state. The relevant question is not simply whether a model includes a baseline covariate, but whether the shared causes of the biomarker history and event process have been adequately represented.

Severity and the Changing Event-Free Risk Set

In time-to-event data, the mixture of latent severity among participants still at risk can change over time. High-severity participants may experience events earlier, leaving a growing proportion of low-severity participants in the event-free risk set. Consequently, a marginal hazard can change over time even when each severity group has a constant hazard. This is one reason that a crude association at one point in follow-up may not summarize the association at another.

Worked Example: A Marginal Hazard Among Survivors

Suppose half the population has high latent severity and half has low severity at the start of follow-up. Their constant event hazards are 0.40 and 0.10 events per person-year, respectively. Assume the biomarker reflects severity but has no effect on the event. Find the marginal hazard at \(t=2\) years among participants who have remained event-free.

Under a constant hazard \(\lambda\), the probability of remaining event-free through time \(t\) is \(e^{-\lambda t}\). Among the initial population, the high-severity survivor weight at time 2 is therefore

$$ w_H(2) =\frac{0.5e^{-0.40(2)}}{0.5e^{-0.40(2)}+0.5e^{-0.10(2)}} =\frac{e^{-0.8}}{e^{-0.8}+e^{-0.2}} \approx \frac{0.449329}{0.449329+0.818731} =0.3543. $$

The low-severity weight is \(1-0.3543=0.6457\) (using unrounded values). The marginal hazard among those still event-free is the weighted average of the two group-specific hazards:

$$ \lambda_{\mathrm{marginal}}(2) =0.3543437(0.40)+0.6456563(0.10) =0.2063031 \approx 0.2063 $$

events per person-year. At baseline the corresponding average is \(0.5(0.40)+0.5(0.10)=0.2500\) events per person-year. The decline in the marginal hazard reflects changing survivor composition; it does not mean that either severity group’s hazard declined. A biomarker associated with severity can inherit this changing risk-set pattern.

Assumptions, Diagnostics, and Alternatives

The examples isolate shared-cause confounding by stipulating that severity causes both the biomarker and event, that the biomarker has no effect on the event, and that no other process changes the stated probabilities or hazards. Those assumptions make the source of the association transparent; they are not facts that can generally be verified from observed data alone. In an applied analysis, a causal interpretation additionally depends on whether relevant confounders have been measured and handled appropriately.

Useful checks include comparing biomarker histories across observed health-severity groups, assessing whether event associations weaken after adjustment for well-measured severity indicators, and checking whether conclusions change under plausible alternative assumptions about unobserved severity. Strong association with multiple measures of health can support the interpretation that the biomarker tracks underlying illness, but it does not prove that all confounding has been removed. A lack of measured imbalance is not proof that latent health is balanced.

The main alternatives answer different questions. A crude time-to-event analysis describes the observed predictive association. Adjustment for measured severity asks whether that association remains conditional on the included variables, subject to model and measurement assumptions. A joint model can represent dependence between a longitudinal process and an event process through shared latent structure. It may therefore be useful when unobserved health plausibly links the two processes, but fitting such a model does not automatically identify a causal biomarker effect: the result depends on the chosen structure and assumptions about the latent effects.

The limitation is fundamental: if severity is unobserved, the data may not distinguish a causal biomarker effect from an association generated by shared latent health without additional information or assumptions. Prediction can remain valuable, but prediction and intervention are different goals. This confounding pathway is also distinct from a biomarker changing in response to impending events; that is a separate explanation for an association.

Common Pitfall

  • Wrong reasoning: “The biomarker is repeatedly associated with event risk, so changing the biomarker would change event risk.” Repeated association may simply show that the biomarker tracks a persistent or worsening latent health state.
  • Correction: Ask what causes both the biomarker trajectory and the event. If latent severity is a shared cause, the observed association need not be causal. State whether the conclusion is predictive or causal, and identify the assumptions needed for the latter.
Key takeaway: Latent disease severity can make repeated biomarker values informative about event risk while generating a biomarker–event association that is not a causal biomarker effect. More repeated measurements can sharpen prediction of underlying health, but they do not by themselves remove confounding.

Check Your Understanding

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

  1. If high severity increases both biomarker elevation and event risk, what path creates the observed association?
  2. In the first worked example, why is the risk higher for the elevated-biomarker group even though the marker has no effect in the model?
  3. What does a high event risk after two elevated readings establish, and what does it not establish?
  4. Why might a marginal hazard decline over time when each severity group has a constant hazard?

Solutions. (1) The shared-cause path is biomarker \( \leftarrow \) latent severity \( \rightarrow \) event. (2) The elevated-biomarker group contains a larger proportion of high-severity participants, who have higher event risk. (3) It establishes predictive information in this hypothetical model, not a causal effect of changing the biomarker. (4) Higher-risk participants tend to experience events earlier, so the remaining event-free group becomes enriched for lower-severity participants and its weighted average hazard falls.

Clinical Biostats