Tutorials › Pattern-Mixture and Selection Models › Worked Example: Longitudinal Delta Adjustment

Delta Adjustment and Its Assumptions · Tutorial 344 of 499

Worked Example: Longitudinal Delta Adjustment

Work through sequential delta-adjusted predictions and imputations across visits, including how prior adjusted values affect later predictions.

Advanced 11 min read

What You'll Learn

  • Define visit-specific delta adjustments for post-dropout outcomes in a longitudinal study.
  • Distinguish reference-model predictions from delta-adjusted conditional predictions.
  • Calculate how an adjustment at one visit propagates to later visits.
  • Compare adjusted trajectories for different dropout patterns.
  • Check simulated imputations for clinical plausibility and trajectory coherence.
  • Recognize what the delta adjustment assumes and what it does not establish.

Why Longitudinal Delta Adjustment Needs a Trajectory

For a single missing outcome, a delta adjustment shifts a reference-model prediction by a specified amount. In a longitudinal study, however, a missing outcome at one visit can also become part of the history used to predict later outcomes. Applying a shift at each visit therefore does not generally produce a trajectory that is just the reference trajectory plus the same fixed amount. Earlier shifts can be carried forward through the sequential model.

This tutorial works through that propagation for monotone dropout: once a participant drops out, all later outcomes are missing. The calculations concern hypothetical continuous symptom scores measured at baseline and at weeks 4, 8, and 12. Higher scores indicate worse symptoms. The numerical trajectories are conditional means, useful for understanding the adjustment; actual multiple imputations also include random variation.

Definition: Let \(Y_t\) denote the outcome at visit \(t\), and let \(H_{t-1}\) denote the participant's outcome and covariate history available before that visit. A sequential reference model specifies the conditional distribution of \(Y_t\) given \(H_{t-1}\) and treatment group. A longitudinal delta adjustment replaces its conditional mean \(m_t(H_{t-1})\) for a missing outcome with \(m_t(H_{t-1})+\delta_t\), where \(\delta_t\) is a specified shift in outcome units. The adjusted missing value is then used in the history for predicting later missing visits.

Because higher symptom scores are worse in the examples, a positive \(\delta_t\) represents a sensitivity assumption that missing outcomes are worse than their reference-model predictions. The sign and size of a delta are substantive choices, not quantities identified by the observed outcomes alone. The earlier material on calibrating delta values and anchoring deltas to clinical thresholds motivates how such choices can be made.

Assumptions Behind the Sequential Calculation

The arithmetic below relies on a clearly specified reference model and a clearly specified sensitivity assumption. In practice, both are essential: an adjustment cannot repair a poorly fitting reference model, and a delta value does not by itself describe every feature of a missing-outcome distribution.

  • Reference model adequacy. The conditional models should represent the observed longitudinal data reasonably well, including important predictors and dependence on prior outcomes. The earlier material on checking the reference model and diagnosing implausible predictions addresses this prerequisite.
  • Sequential conditioning. At each visit, the model conditions on the history available through the preceding visit. For a participant who has dropped out, this history includes earlier imputed outcomes as well as observed outcomes.
  • Specified delta scale and pattern. Each \(\delta_t\) is expressed in the outcome's units and applies to the missing outcomes to which it was assigned. It does not alter outcomes already observed.
  • Distributional choice. A mean shift alone does not specify the full imputation distribution. A common implementation retains the reference model's residual distribution and variance while shifting its location. That is an additional assumption; another distributional change would need to be defined explicitly.
  • Clinical feasibility. Adjusted predictions and draws should be compatible with the outcome's possible range and with plausible within-person change. For a bounded scale, an unconstrained normal model can generate impossible values even when its mean is plausible.

The examples use specified conditional means to make the propagation transparent. They do not claim that a particular delta is correct, nor that a mean trajectory alone captures the distribution of outcomes. Sensitivity analysis compares conclusions over defensible assumptions; it does not identify the unobserved outcomes.

Worked Example: Propagating a Constant Shift Across Three Missing Visits

Suppose a hypothetical symptom scale ranges from 0 to 30, with higher values indicating worse symptoms. A participant has dropped out before week 4, so weeks 4, 8, and 12 are all missing. The reference-model conditional means are

$$ m_4=20,\qquad m_8=4.8+0.60Y_4,\qquad m_{12}=1+0.70Y_8+0.15Y_4. $$

Use a constant delta of \(+2\) points at each missing visit. For this calculation, take the reference-only trajectory to be the recursively calculated conditional means; then compare it with the recursively adjusted conditional means.

At week 4, the reference mean is 20 points, so the adjusted mean is \(20+2=22\) points. At week 8, the reference-only mean is

$$ m_8=4.8+0.60(20)=4.8+12=16.8. $$

With delta adjustment, the week-8 prediction uses the adjusted week-4 value of 22:

$$ m_8^{\delta}=4.8+0.60(22)+2=4.8+13.2+2=20.0. $$

Thus the week-8 difference is \(20.0-16.8=3.2\) points, not 2 points: the week-4 shift contributes \(0.60(2)=1.2\) points to the next conditional mean, in addition to the new 2-point shift.

At week 12, the reference-only mean is

$$ m_{12}=1+0.70(16.8)+0.15(20) =1+11.76+3 =15.76. $$

The adjusted mean uses both earlier adjusted values and the new week-12 shift:

$$ m_{12}^{\delta}=1+0.70(20.0)+0.15(22)+2 =1+14+3.3+2 =20.3. $$

The week-12 difference is therefore \(20.3-15.76=4.54\) points. Equivalently, if \(d_t\) is the adjusted-minus-reference difference, the recurrence gives

$$ d_4=2,\qquad d_8=0.60(2)+2=3.2,\qquad d_{12}=0.70(3.2)+0.15(2)+2=2.24+0.30+2=4.54. $$

All three adjusted means lie within the 0–30 scale. The growing difference is coherent with the specified conditional model: earlier worse symptom predictions contribute to later predictions, while each new missing visit receives its own 2-point shift.

Dropout Timing Determines Which Visits Are Adjusted

The delta is applied only to missing outcomes. A participant who remains observed through week 8 has no week-8 delta adjustment, even if another participant with the same treatment assignment dropped out earlier. The later conditional prediction can nevertheless differ between these participants because their observed histories differ. If delta values vary by dropout pattern, that pattern must be specified along with the visit.

Worked Example: Comparing Two Dropout Patterns

Consider another hypothetical study on the same 0–30 symptom scale. For a participant observed through week 4 with \(Y_4=18\), use the reference conditional means

$$ m_8=4+0.70Y_4,\qquad m_{12}=1+0.65Y_8+0.10Y_4. $$

Suppose the sensitivity schedule assigns \(\delta_8=1.5\) points and \(\delta_{12}=2\) points to missing visits. First consider dropout after week 4, so both later outcomes are missing. The reference-only week-8 mean is \(4+0.70(18)=16.6\) points. The adjusted week-8 mean is

$$ m_8^\delta=4+0.70(18)+1.5=18.1. $$

Using that adjusted week-8 mean as history, the reference-only week-12 mean is

$$ m_{12}=1+0.65(16.6)+0.10(18) =1+10.79+1.8 =13.59. $$

The adjusted week-12 mean is

$$ m_{12}^\delta=1+0.65(18.1)+0.10(18)+2 =1+11.765+1.8+2 =16.565. $$

The difference at week 12 is \(16.565-13.59=2.975\) points. It is also \(0.65(1.5)+2=0.975+2=2.975\): the week-8 shift propagates through the coefficient on \(Y_8\), and the week-12 shift is added.

Now consider a different participant observed through week 8, with \(Y_4=18\) and \(Y_8=15\), who drops out before week 12. The observed week-8 outcome is not shifted. The reference week-12 mean is \(1+0.65(15)+0.10(18)=12.55\) points, and the adjusted mean is \(12.55+2=14.55\) points. This participant receives no week-8 delta because that outcome is observed. The comparison illustrates why dropout timing and observed history both matter when interpreting adjusted trajectories.

From Adjusted Means to Imputed Values

In multiple imputation, the adjusted conditional mean is generally not substituted as the missing outcome. Instead, a random draw is made from a specified adjusted predictive distribution. For example, if the reference conditional distribution is normal with standard deviation \(\sigma_t\), and the chosen assumption shifts its mean but keeps its standard deviation, then a missing value is drawn from \(N(m_t(H_{t-1})+\delta_t,\sigma_t^2)\). That draw becomes part of the history used at the next visit. Repeating this sequentially preserves both model-based uncertainty and the dependence induced by the longitudinal model.

Worked Example: Generating One Sequential Imputation

Use the conditional models from the first example and suppose each visit's reference residual standard deviation is 1 point. Keep the \(+2\)-point shift at each missing visit, and use the following hypothetical standard-normal residual draws in sequence: \(z_4=-0.5\), \(z_8=0.4\), and \(z_{12}=-0.25\). The scale remains 0–30.

The week-4 adjusted mean is 22, so the first imputed draw is \(Y_4^\delta=22+1(-0.5)=21.5\) points. For week 8, use this draw in the conditional mean:

$$ m_8^\delta=4.8+0.60(21.5)+2=4.8+12.9+2=19.7, \qquad Y_8^\delta=19.7+1(0.4)=20.1. $$

At week 12, condition on both imputed values:

$$ m_{12}^\delta=1+0.70(20.1)+0.15(21.5)+2 =1+14.07+3.225+2 =20.295, \qquad Y_{12}^\delta=20.295+1(-0.25)=20.045. $$

The imputed values, 21.5, 20.1, and 20.045 points, are all within the scale's bounds. They are one possible draw, not the adjusted conditional means. In repeated imputations the residual draws vary, and subsequent predictions change because each draw contributes to later history.

Diagnostics, Alternatives, and Limitations

A trajectory should be checked as a sequence, not only one visit at a time. Useful checks include plotting observed and imputed outcome profiles by visit and dropout pattern, comparing adjusted means and quantiles with clinical thresholds, and checking the frequency of values outside the scale's bounds. For a clinical scale, also inspect changes between adjacent visits: plausible individual-level values can still form an implausible trajectory if they imply abrupt, sustained deterioration or improvement unsupported by the scenario.

Look especially for these warning signs:

  • Adjusted means, imputed draws, or clinically important quantiles exceed the outcome's possible range.
  • The same delta is applied to observed visits, or a visit is shifted despite being observed for that participant.
  • Later predictions are calculated from reference-only earlier values when the intended sequential procedure requires adjusted imputed history.
  • A delta schedule gives a pattern that conflicts with the stated clinical direction or produces implausible change across visits.
  • Predictions depend heavily on poorly supported histories, such as extreme imputed values in a sparse part of the data.

Longitudinal delta adjustment differs from applying an independent fixed shift to every missing value: sequential conditioning lets earlier adjusted values affect later conditional means. It also differs from leaving the reference model unchanged, which corresponds to no delta shift under the specified location-shift setup. A different sensitivity approach could alter variance, distributional shape, or the evolution of change over time; those are different assumptions and should not be silently combined with a mean shift.

The method is limited by the reference model, the chosen delta schedule, and the assumed distribution for imputation. If dropout patterns are sparse, a detailed pattern-specific schedule may be weakly supported and unstable. If outcomes are bounded, shifting a normal mean does not ensure valid draws; a model respecting the bounds or a carefully justified constrained procedure may be needed. Finally, a coherent adjusted trajectory is not evidence that it is the true trajectory. It shows the implications of a stated sensitivity assumption.

Key takeaway: For each missing visit, shift the reference conditional mean by the specified delta, generate or calculate the adjusted outcome using the chosen distributional assumption, and carry that adjusted history forward to subsequent visits. Then check both the numerical propagation and the clinical plausibility of the resulting trajectories.

Check Your Understanding

Answer the questions, then compare your work with the solutions.

  1. A participant has an observed week-4 score of 10. The reference week-8 mean is \(3+0.5Y_4\), and the week-8 delta is \(+1.2\). What is the adjusted week-8 mean?
  2. In the first worked example, what is the adjusted-minus-reference difference at week 12, and why is it larger than the week-12 delta?
  3. A participant's week-8 outcome is observed, and week 12 is missing. Should the week-8 outcome receive a delta adjustment? State what is shifted at week 12.
  4. Give one reason an adjusted mean can be clinically plausible while some imputed draws are not.
  5. If the reference model's conditional residual standard deviation is 2 points and the assumed location shift is \(+1\) point, what distributional assumption is made if the adjusted predictive distribution is normal with unchanged variance?

Solutions. (1) The reference mean is \(3+0.5(10)=8\), so the adjusted mean is \(8+1.2=9.2\) points. (2) The difference is 4.54 points. It includes the current 2-point shift plus propagated earlier differences: \(0.70(3.2)+0.15(2)+2=4.54\). (3) No. An observed outcome is retained as observed; the week-12 reference conditional mean, based on that observed history, receives the specified week-12 delta. (4) A distribution has variation around its mean, so some draws can cross a bound even when the mean is inside the scale's range. (5) The assumed distribution is normal with its mean increased by 1 point and variance unchanged at \(2^2=4\) points squared; the residual spread is retained while the location is shifted.

Clinical Biostats