Tutorials › Joint Modeling of Longitudinal and Time-to-Event Data › When Separate Longitudinal and Event Analyses Give Conflicting Results

Why Joint Modeling · Tutorial 25 of 697

When Separate Longitudinal and Event Analyses Give Conflicting Results

Learn why separate summaries of biomarker change and event risk can tell conflicting stories when health-related dropout selectively removes participants from longitudinal follow-up.

Beginner 11 min read

What You'll Learn

  • Distinguish a change in the observed-sample mean from within-person biomarker change.
  • Calculate how health-related dropout changes the composition of a measured sample.
  • Compare baseline and landmark biomarker–event associations in a hypothetical cohort.
  • Identify the assumptions behind complete-case and available-case summaries.
  • Use simple sensitivity calculations to show how unobserved biomarker values affect conclusions.

Why Separate Analyses Can Disagree

A longitudinal analysis asks how a biomarker changes over time. An event analysis asks how event occurrence relates to the biomarker. When both are conducted separately, each may use a different subset of participants and a different time reference. If participants with worsening health are more likely to miss later biomarker visits, the observed measurements need not represent the original cohort. The biomarker average can then move in a direction that does not describe typical within-person change. The biomarker–event association can also differ depending on whether it is calculated at baseline or only among people with a later measurement.

The central issue is selection into the observed sample. A measured value is available only for people who remain under biomarker observation at that time. If continued observation is related to health, the measured participants may differ systematically from those without a measurement. This is not fixed merely by analyzing the longitudinal data and event data in separate models: the two analyses can still represent different selected populations and answer different questions.

Definition: Let \(Y(t)\) denote a participant’s biomarker value at time \(t\), and let \(R(t)=1\) indicate that the value is observed at that time; \(R(t)=0\) indicates it is not observed. The observed-sample mean is \(E\{Y(t)\mid R(t)=1\}\). It is generally not the same as the cohort mean \(E\{Y(t)\}\). When observation depends on health or biomarker history, the missingness is informative for describing the observed trend.

An event can end biomarker measurement because the participant dies or experiences another terminal event. Health-related dropout is different: the participant stops providing biomarker measurements, but event status may still be ascertained afterward. Keeping these situations distinct matters. The example below uses health-related dropout before a landmark visit, while event follow-up continues. The landmark biomarker values for those who dropped out are therefore unobserved, not values recorded after a terminal event.

A Hypothetical Cohort

Consider eight hypothetical participants followed from baseline to a scheduled six-month visit, then for events after that visit. A higher biomarker value indicates worse health. Four participants begin at 16 and four begin at 8. By six months, the marker would be 18 in the first group and 10 in the second if measured. However, three of the four participants in the first group stop biomarker follow-up before the visit because their health is worsening. All four participants in the second group and one participant in the first group have a measured six-month value. Event follow-up continues for everyone.

GroupParticipantsBaseline markerSix-month marker if measuredSix-month values observedEvents after baseline
Higher baseline marker416 for each18 for each18 for 1 participant; 3 missing3
Lower baseline marker48 for each10 for each10 for all 41

Worked Example: The Observed Average Versus the Cohort Trajectory

Calculate the baseline mean, the observed six-month mean, and the six-month mean if all eight scheduled values were available. Then compare the change in the observed mean with within-person change.

Baseline mean. All eight have baseline values, so

$$ \bar{Y}(0)=\frac{4(16)+4(8)}{8} =\frac{96}{8}=12. $$

Observed six-month mean. Only five values are measured: one value of 18 and four values of 10. Thus

$$ \bar{Y}_{\mathrm{observed}}(6) =\frac{18+4(10)}{5} =\frac{58}{5}=11.6. $$

Comparing these two means gives an apparent decrease of \(11.6-12=-0.4\) marker units. Yet every participant’s stipulated six-month value is two units higher than their own baseline value: \(18-16=2\) in the higher group and \(10-8=2\) in the lower group. If all scheduled values were available, the six-month mean would be

$$ \bar{Y}(6)=\frac{4(18)+4(10)}{8} =\frac{112}{8}=14, $$

which is two units above the baseline mean. The apparent decline is caused by a change in who contributes a measured value: the observed sample includes just one participant from the higher-marker group, compared with four at baseline. It is not evidence that the measured participants’ markers declined. In fact, among the five people with paired measurements, each marker rises by two units.

Why Selection Changes a Mean

The difference between an observed-sample mean and a cohort mean can be expressed directly. Let \(Y\) be a biomarker value at a fixed time, and let \(R\) indicate whether it is observed. Since \(E(R)=P(R=1)\) and \(E(YR)=P(R=1)E(Y\mid R=1)\),

$$ \operatorname{Cov}(Y,R) =E(YR)-E(Y)E(R) =P(R=1)\{E(Y\mid R=1)-E(Y)\}. $$

Therefore,

$$ E(Y\mid R=1)-E(Y) =\frac{\operatorname{Cov}(Y,R)}{P(R=1)}. $$

This identity applies when \(P(R=1)>0\). It shows that the observed mean equals the cohort mean when observation is unrelated to the value in the covariance sense. If higher values are less likely to be observed, then \(\operatorname{Cov}(Y,R)<0\), and the observed mean is below the cohort mean. In the example, participants with a higher six-month marker are less likely to have a measurement, so the observed mean, 11.6, is below the all-participant mean of 14.

This is a statement about representation, not a claim that missing values can always be recovered. If the missing values are unknown, the full cohort mean is not identified from the observed mean alone. The value 14 in the example follows from stipulated hypothetical values; an analyst who only saw the five observed values would need additional assumptions or information to estimate the mean for all eight.

1
State the target.
Decide whether the question concerns change among participants with measurements or change in the original cohort. These are different summaries when follow-up is incomplete.
2
Describe who remains observed.
Compare observation rates across health or biomarker histories that are available. A strong difference is a warning that the measured sample may be selected.
3
Keep the event time frame explicit.
State whether the biomarker is measured at baseline or at a later landmark, and which subsequent events are counted.
4
Interpret each result for its actual sample.
Do not treat a complete-case or landmark association as if it automatically described the original cohort.

The Biomarker–Event Association Can Also Change

The same hypothetical cohort illustrates why separate event analyses can appear to conflict. At baseline, the higher-marker group has three events among four participants, while the lower-marker group has one among four. This is a positive baseline association between the marker group and subsequent event risk. But an analysis at the six-month landmark can include only participants who have a measured six-month marker. The three high-marker participants who dropped out are then excluded from the biomarker–event comparison, even though their event follow-up continues.

Worked Example: Baseline Association and Landmark Association

Calculate the event risks by baseline marker group. Then compare event risks after the six-month landmark among participants with a measured six-month marker. All events in this example occur after the landmark, so the comparison uses participants at risk at that time.

Baseline grouping. In the higher-marker group, the event risk is \(3/4=0.75\). In the lower-marker group, it is \(1/4=0.25\). The risk difference is

$$ 0.75-0.25=0.50, $$

and the risk ratio is

$$ \frac{0.75}{0.25}=3. $$

In this small hypothetical cohort, the baseline higher-marker group has an event risk 0.50 higher, or three times the risk, than the lower-marker group over the stated follow-up.

Six-month measured participants. The one observed higher-marker participant has no subsequent event, giving risk \(0/1=0\). The four observed lower-marker participants have one event, giving risk \(1/4=0.25\). The measured-participant risk difference, higher versus lower marker, is \(0-0.25=-0.25\). Thus the landmark complete-case comparison points in the opposite direction from the baseline comparison.

This reversal does not show that a higher marker protects against events. The landmark comparison has removed three high-marker participants whose health-related dropout preceded the visit; their event outcomes are not represented in that biomarker analysis. The calculations also involve very few participants, so they are illustrations of selection, not stable estimates of a population association.

Baseline and landmark analyses answer different questions. The baseline comparison uses the original groups and counts events after baseline. The landmark comparison relates a later measured marker to events after the landmark among a selected subset with measurements. A landmark analysis may be useful when the scientific question is explicitly about prognosis among people reaching that time point with a measurement. It does not automatically estimate the original cohort association. Likewise, an analysis using only participants with complete longitudinal data estimates a selected-sample pattern unless the selection process can be justified as irrelevant to the target.

The earlier material on complete-case comparisons, survivor averages, and health-related attrition motivates these distinctions. Here the practical lesson is to report the sampling set and time origin alongside each result. A numerical biomarker–event association without that information is incomplete.

Assumptions, Diagnostics, and Limits

A complete-case longitudinal summary describes people who have measurements at the required times. To interpret it as representing a broader cohort trajectory requires a defensible reason that selection into complete follow-up does not distort the target summary, or a valid adjustment based on adequate observed information. Merely observing that the available measurements rise or fall does not establish this condition. Health-related dropout makes it particularly doubtful when health history is associated with both the biomarker and continued observation.

Useful diagnostics begin with the observed data rather than a single model coefficient. Tabulate the number measured at each visit; describe how dropout varies by baseline marker, earlier marker values, and recorded health indicators; and compare the baseline characteristics of participants with and without later measurements. For the event analysis, state how many people are included, the time origin, and whether participants who stopped biomarker follow-up remain included in event follow-up. These checks can reveal selection patterns, but they cannot prove that unobserved marker values follow any particular trajectory.

Several alternatives answer distinct questions. A paired analysis uses only participants measured at both times and describes their average within-person change. An available-case mean uses everyone measured at each visit but may compare different participant mixtures over time. A baseline event analysis preserves the original cohort grouping, while a landmark analysis asks about later event risk conditional on reaching the landmark and having the required data. None of these choices is universally correct: the appropriate analysis depends on the target question and the assumptions that can be supported. A joint modeling approach is motivated in part by the need to represent longitudinal and event processes together, but the examples here do not establish a particular joint model or guarantee that one resolves all missing-data problems.

A common mistake is to say, “The average marker fell, so participants’ health improved.” That conclusion confuses a change in the observed-sample mean with within-person change. In the worked cohort, the observed mean falls by 0.4 units, while every paired value rises by two units and the stipulated full-cohort mean rises by two. The correction is to report which participants contribute to each mean, examine how observation relates to health, and avoid assigning the observed-sample trend to the whole cohort without justification.

Key takeaway: Separate longitudinal and event analyses can give conflicting results because health-related dropout changes who contributes biomarker data. The direction of an observed mean trend or association should be interpreted only after stating the time point, analysis sample, and target population.

Check Your Understanding

Answer each question, then compare with the solutions below.

  1. In the hypothetical cohort, why is the observed six-month mean not a valid description of the six-month mean for all eight participants?
  2. Suppose the four missing six-month marker values were instead all 14. What would the full-cohort six-month mean be?
  3. Why does the landmark association not establish a protective effect of a higher marker?
  4. Name one question answered by a paired analysis and one question answered by a baseline event comparison.

Solutions. (1) Only five participants have observed values, and observation is related to health and marker group; the observed group composition differs from the baseline cohort. (2) The four observed low-group values sum to 40, the one observed high-group value is 18, and the four missing values would sum to 56. Thus the full mean would be \((40+18+56)/8=114/8=14.25\) units. (3) The landmark comparison excludes three high-marker participants who dropped out before the visit, has only one observed participant in the higher-marker group, and is based on very small numbers; it is not a causal comparison. (4) A paired analysis describes average within-person change among people measured at both times. A baseline event comparison describes subsequent event risks by baseline marker group in the original cohort.

Clinical Biostats