Tutorials › Joint Modeling of Longitudinal and Time-to-Event Data › Omitting Survival to the Landmark From the Update

Dynamic Prediction · Tutorial 541 of 697

Omitting Survival to the Landmark From the Update

Work through patient-specific predictions that update latent risk using both the biomarker history and survival to the landmark.

Advanced 10 min read

What You'll Learn

  • Write the posterior distribution of patient-specific random effects conditional on biomarker history and survival to the landmark.
  • Derive the conditional event-risk prediction over a future horizon.
  • Calculate how omitting landmark survival changes a patient’s posterior risk and predicted event probability.
  • Distinguish conditioning on survival to the landmark from predicting survival from baseline.
  • Recognize assumptions and diagnostics relevant to the survival contribution in the update.

Why Survival to the Landmark Is Part of the Update

A dynamic prediction at landmark time \(s\) concerns a patient known to be event-free at \(s\), given the biomarker measurements available up to \(s\). The event-free status is information about the patient, just as the biomarker history is. If patients with higher latent risk are more likely to have experienced the event before \(s\), the patients who remain event-free are a selected group. A patient-specific prediction that updates latent risk using only the biomarkers ignores that selection.

The consequence is not merely a technical detail in the likelihood. Omitting survival to the landmark can leave too much posterior weight on latent risk profiles that are less compatible with the patient having remained event-free. The resulting future risk can be too high, sometimes substantially so. The direction and size of the error depend on how survival to \(s\) differs across latent profiles.

Definition: Let \(b_i\) denote patient \(i\)'s latent random effects, \(\mathcal{Y}_i(s)\) the biomarker history observed through landmark \(s\), and \(T_i\) the event time. A dynamic prediction at \(s\) conditions on both \(\mathcal{Y}_i(s)\) and \(T_i>s\). Under a joint model, the posterior distribution of \(b_i\) therefore includes the likelihood of the patient remaining event-free through \(s\).

The Correct Update and Prediction

Suppose the joint model specifies a prior density \(p(b_i)\), a longitudinal likelihood \(L_Y\{\mathcal{Y}_i(s)\mid b_i\}\), and a hazard \(h_i(t\mid b_i)\). Define the event-free survival function conditional on \(b_i\) by

$$ S_i(t\mid b_i)=\exp\left\{-\int_0^t h_i(v\mid b_i)\,dv\right\}. $$

Assume that, conditional on \(b_i\) and the model covariates, the longitudinal and event processes are modeled jointly in the usual way: the observed marker history contributes \(L_Y\), and the event-free status through \(s\) contributes \(S_i(s\mid b_i)\). Then Bayes' rule gives

$$ p\{b_i\mid \mathcal{Y}_i(s),T_i>s\} = \frac{ L_Y\{\mathcal{Y}_i(s)\mid b_i\}\,S_i(s\mid b_i)\,p(b_i) }{ \int L_Y\{\mathcal{Y}_i(s)\mid b\}\,S_i(s\mid b)\,p(b)\,db }. $$

The survival factor is a likelihood contribution, not a second future-risk adjustment. Once this posterior has been formed, the probability of an event by \(s+w\), conditional on being event-free at \(s\) and on \(b_i\), is

$$ P(s<T_i\le s+w\mid T_i>s,b_i) = 1-\frac{S_i(s+w\mid b_i)}{S_i(s\mid b_i)}. $$

Averaging this conditional risk over the updated posterior gives the patient-specific prediction:

$$ \widehat{R}_i(s,w) = \int \left\{1-\frac{S_i(s+w\mid b)}{S_i(s\mid b)}\right\} p\{b\mid \mathcal{Y}_i(s),T_i>s\}\,db. $$

The denominator \(S_i(s\mid b)\) in the conditional risk matters: the target is risk over the future interval among those already known to be event-free at \(s\), not event probability from time zero. An equivalent expression makes the role of the landmark update especially clear:

$$ \widehat{R}_i(s,w) = \frac{ \int L_Y\{\mathcal{Y}_i(s)\mid b\} \{S_i(s\mid b)-S_i(s+w\mid b)\}p(b)\,db }{ \int L_Y\{\mathcal{Y}_i(s)\mid b\}S_i(s\mid b)p(b)\,db }. $$

Why the Survival Factor Belongs in the Posterior

The update follows directly from conditioning. The joint likelihood contribution for the observed marker history and event-free status, conditional on \(b_i\), is \(L_Y\{\mathcal{Y}_i(s)\mid b_i\}S_i(s\mid b_i)\). Multiplying by the prior density and normalizing yields the posterior above. If the survival factor is left out, the calculation instead produces \(p\{b_i\mid\mathcal{Y}_i(s)\}\): a posterior informed by markers but not by the fact that the patient survived event-free to \(s\).

When latent profiles with greater event hazard also have smaller \(S_i(s\mid b)\), the correct update downweights those profiles relative to the marker-only update. This is a selection effect, not a claim that the patient’s survival caused a change in their latent risk. The prediction conditions on the information observed about the patient.

Worked Example: A Marker Update Is Not Enough

Consider a hypothetical patient with two possible latent risk profiles: high-risk \(H\), with prior probability \(0.40\), and low-risk \(L\), with prior probability \(0.60\). The patient’s biomarker history is twice as likely under \(H\) as under \(L\), so the marker likelihoods are \(2\) and \(1\). By the landmark, the probability of remaining event-free is \(0.50\) under \(H\) and \(0.90\) under \(L\). The one-year event risks after the landmark are \(0.40\) under \(H\) and \(0.10\) under \(L\).

A marker-only update assigns unnormalized weights \(0.40(2)=0.80\) to \(H\) and \(0.60(1)=0.60\) to \(L\). Its high-risk posterior probability is \(0.80/(0.80+0.60)=0.5714\).

The patient is also known to have remained event-free to the landmark. Including that information gives weights \(0.40(2)(0.50)=0.40\) and \(0.60(1)(0.90)=0.54\). Thus the correct posterior probability of \(H\) is \(0.40/(0.40+0.54)=0.4255\).

The correct one-year predicted risk is \(0.4255(0.40)+0.5745(0.10)=0.2277\), or 22.77%. Using the marker-only posterior instead gives \(0.5714(0.40)+0.4286(0.10)=0.2714\), or 27.14%. In this hypothetical case, omitting survival through the landmark overstates the predicted one-year risk by 4.37 percentage points.

Worked Example: Exponential Hazards Over Two Time Intervals

Suppose a hypothetical patient has prior probabilities \(0.50\) for each of two latent profiles. The biomarker history has likelihood \(3\) under the high-risk profile and \(1\) under the low-risk profile. The hazards are constant: \(0.20\) per year for the high-risk profile and \(0.05\) per year for the low-risk profile. The patient is event-free at \(s=2\) years, and the prediction horizon is the following year.

Under a constant hazard \(\lambda\), \(S(t)=e^{-\lambda t}\). Survival to year 2 is \(e^{-0.40}=0.6703\) for the high-risk profile and \(e^{-0.10}=0.9048\) for the low-risk profile. Including the marker and survival contributions, the unnormalized weights are \(0.50(3)(0.6703)=1.0055\) and \(0.50(1)(0.9048)=0.4524\). The posterior probability of the high-risk profile is therefore \(1.0055/(1.0055+0.4524)=0.6897\).

Over the next year, the conditional risks are \(1-e^{-0.20}=0.1813\) and \(1-e^{-0.05}=0.04877\). Averaging over the updated posterior gives \(0.6897(0.1813)+0.3103(0.04877)=0.1402\), or 14.02%.

If survival to year 2 is omitted, the marker-only high-risk probability is \(0.50(3)/\{0.50(3)+0.50(1)\}=0.7500\). The resulting risk is \(0.7500(0.1813)+0.2500(0.04877)=0.1482\), or 14.82%. The difference is smaller than in the previous example, but it arises from the same omitted likelihood contribution.

When the Omission Can Matter Most

The omission is most consequential when survival through \(s\) is quite different across the latent profiles that remain plausible after observing the marker history. A strong marker signal does not guarantee that survival information is negligible. In the next example, the marker initially favors the high-risk profile, but survival to the landmark changes that assessment sharply.

Worked Example: Survival Information Can Reverse the Update

For a hypothetical patient, let the prior probabilities of high and low risk be \(0.20\) and \(0.80\). The biomarker history has likelihood ratio \(8\) in favor of high risk. The probability of being event-free at the landmark is \(0.10\) under high risk and \(0.90\) under low risk. The future-horizon risks are \(0.60\) and \(0.05\), respectively.

Using markers alone, the weights are \(0.20(8)=1.60\) and \(0.80(1)=0.80\), so the marker-only posterior probability of high risk is \(1.60/(1.60+0.80)=0.6667\).

After including landmark survival, the weights are \(0.20(8)(0.10)=0.16\) and \(0.80(1)(0.90)=0.72\). The correct posterior probability of high risk is \(0.16/(0.16+0.72)=0.1818\). The correct predicted future risk is \(0.1818(0.60)+0.8182(0.05)=0.1500\), or 15.00%. The marker-only prediction would be \(0.6667(0.60)+0.3333(0.05)=0.4167\), or 41.67%. These figures illustrate how strongly survival to the landmark can alter the update when latent profiles have very different survival probabilities.

Assumptions, Diagnostics, and Alternatives

The calculation relies on the fitted joint model being an appropriate description of the relevant data. In particular, the hazard \(h_i(t\mid b_i)\), the distribution \(p(b_i)\), and the longitudinal likelihood must be specified coherently. The survival contribution \(S_i(s\mid b_i)\) assumes the patient’s event-free status through \(s\) is known and correctly represented by the event model. If a patient was censored before \(s\), that does not establish \(T_i>s\); the observed censoring information must be handled according to the censoring model rather than treated as known event-free survival through \(s\).

A practical diagnostic is to compare the posterior distribution of the latent effects before and after adding the survival contribution. If survival probabilities differ meaningfully across latent profiles, but the two posteriors are identical, the update may have omitted the event-free likelihood. Also check that the landmark survival is included exactly once: including it in the posterior and then treating survival from baseline as the future conditional risk would not target the stated prediction. Numerical integration methods such as quadrature or Monte Carlo integration approximate the required integrals; they do not correct a missing survival factor.

The principal alternative is a marker-only posterior update, which can be appropriate for a different target that conditions on biomarker history but not on event-free survival to \(s\). It is not the patient-specific prediction for someone known to be event-free at \(s\). A population-level marginal survival probability is different again: it averages over a population distribution and does not, by itself, incorporate this patient’s biomarker history and landmark status. The distinction between marginal survival and conditional risk is useful here, but the key point is that the landmark information must be included in the patient-specific update.

Common mistake: “The marker history already tells us the patient’s risk profile, so we can update using the marker likelihood alone.” This reasoning ignores that the patient also remained event-free to \(s\). The correction is to multiply the marker likelihood by \(S_i(s\mid b_i)\) when forming the posterior, then average the conditional future-horizon risk over that updated posterior.

Check Your Understanding

Answer each question, then compare with the solutions below.

  1. In the posterior update, what information does \(S_i(s\mid b_i)\) contribute?
  2. A patient has two equally likely latent profiles before observing data. The marker likelihoods are \(2\) and \(1\), and survival probabilities to the landmark are \(0.60\) and \(0.90\). What is the posterior probability of the first profile when both sources of information are used?
  3. Why is \(1-S_i(s+w\mid b_i)\) not generally the conditional risk from landmark \(s\) to \(s+w\)?
  4. Under what circumstance can omitting survival to the landmark have little effect on the posterior update?

Solutions. (1) It is the likelihood of the patient remaining event-free through \(s\), conditional on the latent profile. (2) The weights are \(0.50(2)(0.60)=0.60\) and \(0.50(1)(0.90)=0.45\), so the posterior probability is \(0.60/(0.60+0.45)=0.5714\). (3) It is event probability from time zero through \(s+w\), not risk conditional on already being event-free at \(s\). The conditional risk is \(1-S_i(s+w\mid b_i)/S_i(s\mid b_i)\). (4) The omission may have little effect when survival to \(s\) is approximately the same across the latent profiles supported by the marker history, so the survival factor changes their relative posterior weights very little.

Clinical Biostats