Tutorials › Competing Risks and Multi-State Models › Changing Target Populations Across Sequential Landmarks

Landmarking and Dynamic Prediction · Tutorial 499 of 699

Changing Target Populations Across Sequential Landmarks

Learn to identify the population behind a landmark prediction, separate changes in cohort composition from changes in conditional predictions, and distinguish landmark-cohort estimates from fixed-population summaries.

Advanced 10 min read

What You'll Learn

  • Define the eligible population at each landmark and the target of its prediction.
  • Explain how survival, absorption, and observation affect landmark-cohort composition.
  • Decompose changes in average predictions into composition and within-stratum components.
  • Standardize landmark predictions to a common reference distribution when appropriate.
  • Distinguish conditional landmark predictions from unconditional state probabilities in a fixed baseline population.
  • Diagnose limited overlap and other reasons landmark comparisons can mislead.

Why the Target Population Can Change

A prediction made at a landmark applies to people who are eligible to be predicted at that time, not automatically to the people who formed the original study population. At a later landmark, some individuals may have entered an absorbing event state, died, been lost to observation, or failed another eligibility rule. Those who remain can differ systematically from those eligible earlier. Consequently, a difference between two average landmark predictions may reflect a changed population as well as changed predictions within comparable people.

This distinction matters whenever sequential predictions are reported side by side. A lower predicted risk at a later landmark does not, by itself, imply that the underlying process has become safer. The later cohort may simply contain a different mix of states, histories, or covariates. The first step in interpreting any such summary is to state precisely who is included and what probability is being averaged.

Define the Landmark Cohort and Its Prediction

Let \(t_l\) denote landmark \(l\), and let \(Y(t)\) be an individual's state at time \(t\). Define \(E_l=1\) when the individual meets the prespecified eligibility criteria at \(t_l\). For example, eligibility might require being alive, not yet in an absorbing event state, and under observation. The definition must be applied consistently at every landmark; otherwise, differences can arise merely because the rules changed.

Definition: The landmark cohort at \(t_l\) is the population with \(E_l=1\). For a future horizon \(u>0\) and state \(k\), an individual landmark prediction is \(q_{lk}(u)=P\{Y(t_l+u)=k\mid E_l=1,H_l\}\), where \(H_l\) is the information available for that individual at \(t_l\). The cohort-average prediction is the average of these individual predictions over a stated target distribution within the eligible population. If that distribution is the observed landmark cohort, the average targets that cohort, which can differ from the cohort at another landmark.

The conditioning information \(H_l\) may include the current state, covariates, and relevant observed history. It must not include information from after \(t_l\). The eligible population and the prediction horizon also need to be explicit: a one-year prediction from \(t_l\) concerns a different calendar endpoint from a one-year prediction from \(t_{l+1}\).

Eligibility is often selective rather than random. For instance, people who remain free of an absorbing event until a later landmark have necessarily had a different event history from those who experienced that event earlier. Conditioning on later eligibility can therefore change the distribution of prognostic factors, even if those factors were balanced at baseline. This is a feature of the target population, not automatically a defect in the prediction calculation.

How Composition Changes an Average

Suppose \(X\) is a categorical characteristic recorded in comparable form at each landmark, and \(p_l(x)\) is the predicted probability of a specified future state for people in stratum \(X=x\) at landmark \(l\). The crude average for that landmark is

$$ \psi_l=\sum_x p_l(x)F_l(x), $$

where \(F_l(x)=P(X=x\mid E_l=1)\) is the stratum proportion among those eligible at \(t_l\). Thus, \(\psi_l\) can change because the stratum-specific predictions \(p_l(x)\) change, because the eligible cohort's composition \(F_l(x)\) changes, or both.

One exact decomposition makes the two components visible. For two landmarks 1 and 2, add and subtract \(\sum_x p_1(x)F_2(x)\):

$$ \psi_2-\psi_1 = \sum_x p_1(x)\{F_2(x)-F_1(x)\} + \sum_x F_2(x)\{p_2(x)-p_1(x)\}. $$

The first term is the change associated with applying the earlier stratum-specific predictions to the later cohort's composition. The second is the change in stratum-specific predictions, averaged over the later composition. This decomposition is an accounting identity, not a causal attribution. Its interpretation requires comparable strata and a common, clearly defined prediction target at both landmarks.

Worked Example: Separating Composition From Within-Stratum Change

Consider a hypothetical study with two eligibility landmarks. At the first landmark, 20% of the eligible cohort is in a high-risk stratum and 80% is in a lower-risk stratum. The predicted one-year probability of a specified event is 0.30 in the high-risk stratum and 0.10 in the lower-risk stratum. At the second landmark, the eligible cohort is 60% high-risk and 40% lower-risk, and the corresponding predictions are 0.32 and 0.12. Assume that strata and the one-year event definition are comparable across landmarks.

At landmark 1, the crude cohort-average prediction is

$$ \psi_1=(0.20)(0.30)+(0.80)(0.10)=0.060+0.080=0.140. $$

At landmark 2, it is

$$ \psi_2=(0.60)(0.32)+(0.40)(0.12)=0.192+0.048=0.240. $$

The crude average rises by \(0.240-0.140=0.100\), or 10 percentage points. Using the decomposition, the composition component is

$$ (0.30)(0.60-0.20)+(0.10)(0.40-0.80) =0.120-0.040=0.080. $$

The within-stratum component, weighted by the second cohort's composition, is

$$ (0.60)(0.32-0.30)+(0.40)(0.12-0.10) =0.012+0.008=0.020. $$

The components sum to \(0.080+0.020=0.100\), the observed change in the crude averages. In this example, most of the increase is associated with the greater share of high-risk individuals in the later cohort. The calculation does not establish why the composition changed or imply that landmark timing caused the risk difference.

Standardizing to a Fixed Reference Distribution

A common way to compare conditional predictions across landmarks is to average them over the same reference distribution \(G(x)\), rather than over each landmark's own cohort distribution. The standardized summary is

$$ \psi_l^{G}=\sum_x p_l(x)G(x), $$

where \(G(x)\) is fixed across the comparisons. This answers a different question from the crude cohort average: it compares the landmark-specific stratum predictions after holding the stratum mix constant. It does not turn a conditional prediction into an unconditional prediction for everyone in the original study population.

Standardization requires that the strata be defined comparably, that predictions be available for the reference strata at each landmark, and that the reference distribution be substantively appropriate. If a stratum has no eligible individuals at one landmark, its prediction there may not be estimable from the data without additional modeling assumptions. Extrapolating into such a stratum should be identified as extrapolation, not presented as a data-supported comparison. When current state or history matters, these features must also be handled in the definition of the reference distribution or in the conditioning set; standardizing only a baseline covariate may leave important composition differences unresolved.

Worked Example: Comparing Predictions With a Common Mix

Use the hypothetical stratum-specific predictions in the preceding example. Choose a fixed reference distribution with 50% high-risk and 50% lower-risk individuals. At landmark 1, the standardized average is

$$ \psi_1^G=(0.50)(0.30)+(0.50)(0.10)=0.150+0.050=0.200. $$

At landmark 2, it is

$$ \psi_2^G=(0.50)(0.32)+(0.50)(0.12)=0.160+0.060=0.220. $$

The standardized difference is \(0.220-0.200=0.020\), or 2 percentage points, rather than the 10-point difference between the crude averages. This comparison holds the high-risk proportion fixed at 50%. It is meaningful only if both landmark-specific predictions can be supported for the two reference strata and if the reference mix is relevant to the question.

Landmark-Conditional Predictions Versus a Fixed Population

A fixed baseline population is a different target. If the question concerns the state occupied at a calendar endpoint among everyone in the baseline population, the summary should retain everyone in its denominator and account for their state at that endpoint, including absorbing states. People who experienced an absorbing event before a later landmark are not simply omitted from this fixed-population question.

By contrast, a landmark-cohort prediction is conditional on eligibility at its landmark. It describes future outcomes among those who meet that landmark's eligibility rules. The two summaries can both be useful, but they answer different questions and their probabilities should not be compared as if they had the same denominator.

Worked Example: Reconstructing a Fixed-Population State Probability

Suppose a hypothetical baseline population has 1,000 people. By landmark 1, 120 have entered an absorbing target-event state and 80 have entered an absorbing competing-event state. The remaining 800 are eligible at landmark 1. Assume that everyone is followed or appropriately accounted for, that the two event states are absorbing, and that the landmark prediction for the eligible group gives probabilities of 0.25 for entering the target-event state by one year after the landmark, 0.10 for entering the competing-event state, and 0.65 for remaining in a transient state at that endpoint.

The predicted target-event count in the eligible group is \(800(0.25)=200\). The competing-event count is \(800(0.10)=80\), and the transient-state count is \(800(0.65)=520\). For the fixed baseline population at this endpoint, the target-event state includes those already absorbed before landmark 1:

$$ P\{Y(t_1+1)=\text{target event}\} = \frac{120+200}{1000} =0.320. $$

Likewise, the competing-event probability is \((80+80)/1000=0.160\), and the transient-state probability is \(520/1000=0.520\). These sum to \(1.000\). The landmark-conditional target-event probability is 0.25, whereas the fixed-baseline target-event probability is 0.32. Neither is a calculation error: the first uses the 800 eligible people as its target, while the second includes the 1,000 baseline individuals and counts earlier target events in the target state.

Assumptions, Diagnostics, and Limitations

Interpretation depends first on clear eligibility and outcome definitions. Eligibility should be set using information available at the landmark, and predictions should use no future information. The target state must be defined consistently, including how competing absorbing states are represented. For cohort-average predictions, the averaging distribution must be stated: the observed eligible cohort, a fixed reference distribution, or another explicit target.

When estimates rely on observed outcomes, censoring and loss to follow-up also matter. If individuals with incomplete follow-up differ in ways related to their future state, simply averaging observed outcomes may not estimate the intended cohort probability. The analysis must use an appropriate method for its observation process and state its assumptions. A model-based prediction does not remove the need to examine whether relevant predictors and histories are represented in the data.

Useful diagnostics include reporting the number eligible at each landmark, the distribution of current states and key covariates, and the number of observations contributing to each stratum-specific estimate. Compare these distributions before interpreting changes in crude averages. Look for strata present at one landmark but absent at another, sparse strata, or marked shifts in state and history composition. These are warnings that standardization may depend heavily on extrapolation or that a simple composition decomposition omits relevant differences.

The main alternatives answer distinct questions. A crude landmark average describes each landmark's own eligible cohort. A standardized average compares predictions using a common reference mix, subject to overlap and modeling support. A fixed-baseline state probability describes occupation in a defined original population at a calendar time and includes outcomes before later landmarks. Reporting one does not substitute for another when the scientific question concerns a different target.

SummaryPopulation being averagedQuestion answered
Crude landmark averageEligible cohort at that landmarkWhat is the average prediction for this landmark's eligible people?
Standardized landmark averageFixed reference distribution, where supportedHow do landmark-specific predictions compare under a common mix?
Fixed-baseline state probabilityOriginal baseline populationWhat proportion of the baseline population occupies each state at the endpoint?

A common mistake is to treat the average at each landmark as if it described the same people, then explain any difference as an improvement or deterioration in prognosis. That reasoning ignores selection into later eligibility. The correction is to name the target population for each average, inspect how its composition changed, and use standardization or a fixed-population summary only when it answers the intended question and its assumptions are defensible.

Key takeaway: Sequential landmark averages generally target sequentially changing eligible cohorts. To compare them fairly, distinguish changes in cohort composition from changes in conditional predictions, and do not confuse a standardized conditional summary with an unconditional state probability for the baseline population.

Check Your Understanding

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

  1. At a later landmark, more eligible individuals are in a high-risk stratum. The stratum-specific predictions are unchanged. Can the crude cohort-average prediction still change? Why?
  2. A later landmark's average predicted event probability is 0.18, compared with 0.12 at an earlier landmark. What additional information is needed before interpreting this as a change in prognosis?
  3. In a hypothetical reference population that is 25% high-risk and 75% lower-risk, stratum-specific predictions are 0.28 and 0.08. Calculate the standardized average.
  4. Why is a fixed-baseline state probability not generally equal to a landmark-conditional probability?

Solutions. (1) Yes. The average weights the stratum-specific predictions by the stratum proportions among those eligible; changing those proportions can change the average even when the predictions within strata are constant. (2) Report eligibility counts and composition at both landmarks, clarify whether horizons and outcomes match, and examine stratum-specific predictions or a common-reference standardized comparison. Also check whether the target population and prediction support are comparable. (3) The standardized average is \((0.25)(0.28)+(0.75)(0.08)=0.070+0.060=0.130\), or 13.0%. (4) The landmark-conditional probability is restricted to people eligible at that landmark. A fixed-baseline probability includes the original population and accounts for states, including absorbing events, reached before the landmark as well as later outcomes.

Clinical Biostats