Tutorials › Pattern-Mixture and Selection Models › Diagnosing Implausible Delta-Adjusted Predictions

Delta Adjustment and Its Assumptions · Tutorial 335 of 499

Diagnosing Implausible Delta-Adjusted Predictions

You will learn how to audit delta-adjusted predictions for logical and clinical plausibility without mistaking a diagnostic warning for evidence that identifies the correct delta.

Advanced 13 min read

What You'll Learn

  • Define the support and conditional reference distribution needed to assess adjusted predictions.
  • Calculate how much adjusted predictive probability falls outside an outcome’s possible range.
  • Compare reference and adjusted predictions across clinically distinct subgroups.
  • Diagnose implausible changes in level and rate across repeated measurements.
  • Distinguish a useful warning from proof that a delta value or missingness assumption is wrong.
  • Choose graphical and numerical checks that complement, rather than replace, clinical judgment.

Why Inspect Adjusted Predictions?

A delta adjustment expresses a sensitivity assumption about outcomes that are missing. Even when a chosen delta is defensible, its consequences may not be. Predictions can fall outside the outcome’s possible range, change sharply across a subgroup boundary, or imply a longitudinal course that is difficult to reconcile with the clinical setting. These are reasons to investigate the model and the assumptions—not automatic proof that a particular delta is invalid.

The earlier material on checking the reference model provides the starting point: assume the reference prediction model has been examined for its intended use. The task here is different. We inspect what happens after applying the adjustment, including how that result behaves across outcome values, subgroups, and time. The observed data do not reveal the missing outcomes, so these checks assess implications and coherence rather than determine the true delta.

Definition: Let \(Y\) be an outcome, \(H\) the observed covariates or history used for prediction, and \(R=0\) indicate that \(Y\) is missing. Write \(F_{\text{ref}}(y\mid H)\) for the reference predictive distribution. Under an additive location delta \(\delta(H)\), the delta-adjusted predictive distribution is defined by \(Y_{\text{adj}}=Y_{\text{ref}}+\delta(H)\), conditional on \(H\). Thus a positive delta shifts predictions upward, in the units of \(Y\). The support of an outcome is the set of values it can logically take.

This definition describes a distributional shift, not merely an instruction to add delta to a reported average. A check based only on adjusted means can miss implausible tails or individual trajectories. Conversely, a randomly drawn imputed value outside a plausible central range is not automatically impossible: it matters whether it lies outside the outcome’s actual support, and how much predictive probability is assigned there.

Start With Support and Scale

Before plotting, state what values the outcome can take and what its units mean. A score may be bounded by a test scale, a count cannot be negative, and a proportion lies between 0 and 1. Some limits are logical; others are clinical or operational and require justification. Keep these categories distinct. A value can be rare or clinically concerning without being mathematically impossible.

For a continuous outcome with support \([L,U]\), an additive delta shifts the reference distribution. The adjusted probability below or above the support is therefore

$$ P(Y_{\text{adj}}<L\mid H)=F_{\text{ref}}(L-\delta(H)\mid H),\qquad P(Y_{\text{adj}}>U\mid H)=1-F_{\text{ref}}(U-\delta(H)\mid H). $$

These expressions make the diagnostic concrete: evaluate the reference distribution at the shifted boundaries. For discrete outcomes, use the corresponding probabilities at values outside the allowed set, taking account of whether the boundary values themselves are allowed. If a predictive model already respects the outcome’s support, a location shift may not. A large out-of-support probability is a warning about combining that reference distribution with that additive shift.

Worked Example: Out-of-Range Predictions for a Bounded Score

Suppose a hypothetical symptom score is measured from 0 to 100. For a particular history \(H\), the reference predictive distribution is \(N(82,8^2)\), and the proposed additive delta is \(+18\) score points. Treat the normal distribution as a working predictive model for this calculation, while recognizing that it itself has unbounded support.

The adjusted distribution has mean \(82+18=100\) and standard deviation 8. The probability of exceeding the score maximum is

$$ P(Y_{\text{adj}}>100\mid H) =P(Y_{\text{ref}}>100-18\mid H) =P\left(Z>\frac{82-82}{8}\right)=P(Z>0)=0.500. $$

Thus the adjusted model assigns 50.0% probability to scores above 100. That is not a handful of extreme draws; it is a major incompatibility between this normal reference distribution, the additive shift, and the score’s support. The calculation does not show that the delta is clinically wrong. It shows that this specification cannot be interpreted literally as a predictive distribution on the 0-to-100 score scale without further modeling decisions.

One should not silently truncate predictions at 100 and call the result the same additive-delta model. Truncation changes the distribution and generally changes its mean. A bounded-outcome model or a different adjustment may be more coherent, but the new specification must be stated and checked on its own terms.

Look for Abrupt Subgroup Changes

When the adjustment depends on treatment arm or a clinically distinct subgroup, compare the reference and adjusted predictions side by side. Useful summaries include subgroup means, selected quantiles, predicted probabilities above a clinical threshold, and contrasts between subgroups. Plotting these summaries often makes a discontinuity easier to see than reading a table of delta parameters.

A discontinuity is not necessarily a flaw. If the delta is deliberately different between groups, the adjusted predictions may also differ. The diagnostic question is whether the size and direction of that change have a credible interpretation given the scientific rationale for the subgroup-specific assumption. Check that the change is not driven by a coding error, a sparse subgroup, or an unintended interaction.

Worked Example: A Sudden Difference Across Subgroups

In a hypothetical study, the reference predicted mean outcome is 70 points in subgroup A and 68 points in subgroup B. The proposed deltas are \(-2\) points for A and \(+15\) points for B. Assume both means and deltas refer to the same outcome scale and comparable prediction histories.

The adjusted subgroup means are \(70+(-2)=68\) for A and \(68+15=83\) for B. The reference contrast, B minus A, is \(68-70=-2\) points. The adjusted contrast is \(83-68=15\) points. The change in the contrast is therefore

$$ 15-(-2)=17\text{ points}. $$

The adjustment reverses the ordering and creates a 17-point change in the subgroup contrast. This is a diagnostic flag to review the clinical justification and the implementation of the subgroup-specific deltas. It does not establish that the groups’ missing outcomes should be equal, nor does it establish that their adjusted means must be similar. If the shift was intended to represent a real difference in missing outcomes, the resulting contrast should be described and defended in those terms.

A plot of the two reference means and two adjusted means would show both the reversal and its size. A table should also report the number of individuals and amount of missing information in each subgroup: a sharp estimate in a small subgroup may be particularly sensitive to the reference predictions.

Inspect Trajectories, Not Only Time-Specific Means

For repeated outcomes, a collection of individually plausible time-specific predictions can still form an implausible path. Plot adjusted trajectories against time for representative histories, and compare their increments with the reference trajectories. Examine whether shifts create sudden jumps, reversals, or sustained rates of change that conflict with the outcome’s context. Do not assume a smooth trajectory is always correct; some outcomes change abruptly. The relevant benchmark is clinical and design knowledge, not smoothness for its own sake.

For equally spaced times, a simple numerical check is to calculate adjacent increments. If adjusted means at times \(t_1,\ldots,t_k\) are \(m_1,\ldots,m_k\), the increment from \(t_j\) to \(t_{j+1}\) is \(m_{j+1}-m_j\). Large or alternating increments can flag a pattern for review. With unequal time intervals, compare rates of change by dividing each increment by the elapsed time.

Worked Example: A Change in the Implied Longitudinal Course

Consider hypothetical reference means for an outcome measured every four weeks: 50, 55, 60, and 64 points. Suppose the time-specific deltas are 0, 2, 5, and 12 points. The adjusted means are 50, 57, 65, and 76 points.

The reference increments are \(55-50=5\), \(60-55=5\), and \(64-60=4\) points per four weeks. The adjusted increments are \(57-50=7\), \(65-57=8\), and \(76-65=11\) points per four weeks. The adjusted increase in the last interval is \(11-4=7\) points larger than the reference increase for that interval.

A second-difference check summarizes changes in increments. For the adjusted means, the second differences are \(8-7=1\) and \(11-8=3\) points per four weeks between successive increments. These values show an accelerating adjusted increase in this example. Whether that acceleration is implausible depends on what the score measures and the time scale; the arithmetic identifies the feature to discuss, not a universal cutoff.

A trajectory plot should show the reference and adjusted paths together, with the time intervals and units labeled. If the path is inconsistent with the intended meaning of the delta, inspect the time-specific adjustments and the prediction histories used at each visit.

A Practical Diagnostic Sequence

A disciplined review combines simple numerical summaries with plots. Use the same prediction records for reference and adjusted summaries so that differences reflect the adjustment rather than a changed target population.

1
Specify the scale and support.
Record the outcome units, logical limits, and any separately justified clinical limits. State which limits are hard constraints and which are warning thresholds.
2
Compare predictive distributions.
For relevant histories or groups, compare reference and adjusted means, quantiles, and probabilities near or beyond boundaries. Report the fraction of predictions outside support as well as the magnitude of any excursions.
3
Plot subgroup and time patterns.
Show reference and adjusted summaries on the same axes. Label groups, times, outcome units, and uncertainty where available; do not let a plot’s scale hide a consequential change.
4
Trace warnings to their assumptions.
Check the delta’s sign, units, subgroup or time indexing, and whether the pattern was intended. Ask whether the reference predictions and the adjustment together tell a coherent scientific story.

Numerical checks should be chosen to match the concern. For a bounded outcome, report out-of-support probability. For a subgroup adjustment, report contrasts before and after adjustment. For longitudinal predictions, report adjacent increments and, when useful, second differences. No single number substitutes for inspecting the distribution and the context.

What These Checks Can and Cannot Establish

An implausible adjusted prediction can arise from several sources: an overly strong delta, an unsuitable reference distribution, a support mismatch, an indexing or coding mistake, or a real feature of the scientific scenario that needs clearer explanation. The diagnostic alone does not distinguish these causes. In particular, observed outcomes cannot verify the assumed distribution of outcomes among those missing. A plot that looks reasonable is not evidence that the unobserved outcomes follow the adjusted distribution.

The main alternatives answer different questions. Reference-model diagnostics assess the reference predictions before adjustment. Delta sensitivity analysis compares predictions under alternative assumptions about missing outcomes. Distributional approaches, such as a bounded-outcome model or an adjustment of a different distributional feature, may better respect the scale in some settings. None removes the need to state the assumption. Mechanically clipping values, smoothing subgroup changes, or deleting uncomfortable predictions can conceal rather than solve a mismatch.

Treat diagnostics as a reasoned review, not a significance test. There is generally no universal cutoff for an acceptable subgroup jump or trajectory curvature. Document the feature, its size and units, the groups or histories affected, and the scientific rationale for retaining or revising the specification. If a different adjustment is considered, recompute and inspect its implications rather than assuming it fixes the problem.

Key takeaway: A delta-adjusted prediction is a consequence of both the reference model and a sensitivity assumption. Check its support, subgroup contrasts, and longitudinal behavior; use warnings to clarify or revise assumptions, not to claim that the missing outcomes have been discovered.

Check Your Understanding

Answer each question, then compare with the solutions below.

  1. A continuous outcome is restricted to 0–40. For one history, the reference prediction is \(N(30,5^2)\) and the additive delta is \(+10\). What proportion of adjusted predictions exceeds 40, and what does this indicate?
  2. Two subgroups have reference means 12 and 16 units, with deltas \(+3\) and \(-2\), respectively. Calculate the adjusted means and the change in the second-minus-first contrast.
  3. Why is a prediction above a clinical warning threshold not necessarily impossible, and why is a value outside a hard logical bound different?
  4. Adjusted repeated-outcome means increase by 4, 4, and 13 units over three equal intervals. What numerical feature should be flagged, and what can it not establish by itself?

Solutions. (1) The adjusted mean is \(30+10=40\), with standard deviation 5. Thus \(P(Y_{\text{adj}}>40)=P(Y_{\text{ref}}>30)=0.500\). Half the predictive distribution exceeds the scale maximum, indicating a serious support mismatch for this model and shift, not proving which component is wrong. (2) The adjusted means are \(12+3=15\) and \(16-2=14\). The original contrast is \(16-12=4\); the adjusted contrast is \(14-15=-1\), a change of \(-1-4=-5\) units. (3) A clinical warning threshold may identify a rare or concerning but possible outcome; a hard bound rules values out by definition. (4) The final increment is 9 units larger than the previous increment, so the path has a sharp acceleration. This flags a trajectory to review but does not, alone, show that the delta or the trajectory is scientifically wrong.

Clinical Biostats