Tutorials › Pattern-Mixture and Selection Models › Additive Delta Shifts for Continuous Outcomes

Delta Adjustment and Its Assumptions · Tutorial 310 of 499

Additive Delta Shifts for Continuous Outcomes

Learn how to shift the conditional mean of missing continuous outcomes by a delta while retaining the MAR reference residual distribution, and how to interpret the resulting sensitivity analysis.

Advanced 11 min read

What You'll Learn

  • Define an additive delta shift relative to a MAR reference distribution.
  • Derive the shifted conditional distribution and identify what remains unchanged.
  • Apply the adjustment to conditional probabilities and imputed values.
  • Calculate how a delta changes completed-data means and treatment contrasts.
  • Check the assumptions and diagnose implausible or incorrectly applied shifts.
  • Distinguish a sensitivity parameter from a quantity estimated from observed outcomes.

Why Shift a Reference Distribution?

For a continuous outcome with missing values, an analysis needs a model for outcomes that were not observed. A common starting point is a MAR reference distribution: conditional on the recorded covariates and observed history, missing outcomes are assigned the distribution specified by the reference model. This is a useful benchmark, but it does not establish that the benchmark is correct for missing outcomes. An additive delta adjustment makes a chosen departure from that reference explicit.

This tutorial focuses on one specific construction: shift the conditional mean of a continuous missing outcome by a fixed amount, while retaining the reference model’s residual distribution. The adjustment applies to missing outcomes only. The observed outcome distribution is not altered.

Definition: Let \(Y\) be a continuous outcome, \(X\) the recorded covariates and history used in the reference model, and \(R\) an observation indicator, with \(R=1\) for an observed outcome and \(R=0\) for a missing outcome. Let \(F_0(y\mid X=x)\) be the MAR reference conditional distribution for a missing outcome, with conditional mean \(m_0(x)\). For a specified additive shift \(\delta\), define the shifted missing-outcome distribution by \(Y_\delta=Y_0+\delta\), where \(Y_0\mid X=x\) has distribution \(F_0(\cdot\mid x)\). Thus the shifted conditional mean is \(m_0(x)+\delta\). The observed outcomes remain as recorded.

Here \(\delta\) is expressed in the original outcome units. For example, if \(Y\) is a symptom score measured in points, \(\delta=2\) shifts the conditional mean for missing outcomes upward by 2 points. This definition uses the convention that a positive \(\delta\) adds to the reference mean; the direction is part of the specification and should be stated when reporting results.

The reference distribution is a model-based benchmark, not a directly observed distribution of missing outcomes. In particular, the data do not generally identify \(\delta\). It is chosen to examine how conclusions change under departures from the MAR reference assumption.

What the Additive Shift Does

The defining operation is a translation of the reference outcome. If \(Y_0\mid X=x\) has reference density \(f_0(\cdot\mid x)\), then adding \(\delta\) gives the shifted density \(f_\delta(y\mid x)=f_0(y-\delta\mid x)\). The same translation gives the shifted cumulative distribution function. These relations hold whether or not the reference residuals are normally distributed.

$$ F_\delta(y\mid X=x)=P(Y_0+\delta\leq y\mid X=x) =F_0(y-\delta\mid X=x), \qquad f_\delta(y\mid x)=f_0(y-\delta\mid x). $$

If the reference model is written as \(Y_0=m_0(X)+\varepsilon\), the additive adjustment can be represented as \(Y_\delta=m_0(X)+\delta+\varepsilon\). The same residual \(\varepsilon\) is retained. In particular, whenever the relevant moments exist,

$$ E(Y_\delta\mid X=x)=m_0(x)+\delta, \qquad \operatorname{Var}(Y_\delta\mid X=x) =\operatorname{Var}(Y_0\mid X=x). $$

The variance result follows because adding a constant does not change deviations from the mean: \(Y_\delta-E(Y_\delta\mid x)=Y_0-E(Y_0\mid x)\). Consequently, the conditional variance and the shape of the residual distribution are retained; only location changes. If the reference model has covariate-dependent residual distributions, the adjustment retains each such conditional distribution rather than imposing a single common variance.

For this construction to have the stated interpretation, the reference distribution must be defined for the covariate and history values at which outcomes are missing, and the same additive shift must be meaningful over the range to which it is applied. The reference model itself must be adequate for its intended conditional prediction. The adjustment does not repair a misspecified reference model; it adds a specified location departure to it.

Worked Example: Shifting a Normal Reference Distribution

Suppose a hypothetical continuous score has reference distribution \(Y_0\mid X=x\sim N(50,100)\) for a particular covariate profile. The second parameter is the variance, so the reference standard deviation is \(10\) score points. Consider \(\delta=3\) points.

The shifted distribution has mean \(50+3=53\) points and variance \(100\), so \(Y_\delta\mid X=x\sim N(53,100)\). To compare upper-tail probabilities, under the reference distribution,

$$ P(Y_0>60\mid x) =1-\Phi\left(\frac{60-50}{10}\right) =1-\Phi(1) \approx 0.1587. $$

After the shift,

$$ P(Y_\delta>60\mid x) =1-\Phi\left(\frac{60-53}{10}\right) =1-\Phi(0.7) \approx 0.2420. $$

Under the specified shifted model, the probability of exceeding 60 is about 0.242 rather than 0.159 for this covariate profile. The increase occurs because the mean moved upward by 3 points; the standard deviation stayed at 10 points. This calculation assumes the normal reference model is appropriate for the profile and that the additive shift is the intended departure. It does not show that the missing outcomes actually follow the shifted distribution.

Applying the Shift to Imputations and Summaries

In a completed-data analysis, the shift is applied to the missing outcomes generated from the reference distribution. If a reference-model draw for a missing record is \(y_0^{(j)}\), the corresponding shifted draw is \(y_\delta^{(j)}=y_0^{(j)}+\delta\). The draw’s residual variation is retained; replacing all draws by a conditional mean would instead remove that variation.

For a sample with \(n\) people, let \(n_1\) outcomes be observed and \(n_0\) be missing, so \(n=n_1+n_0\). If the reference predictions for missing outcomes are \(m_0(X_i)\), the reference completed-data mean is the sum of observed outcomes and predicted missing means divided by \(n\). Under a common additive shift, each missing prediction gains \(\delta\), while observed outcomes do not change. Therefore the completed-data mean changes by \(n_0\delta/n\), not by \(\delta\) unless every outcome is missing.

$$ \bar{Y}_\delta = \frac{\sum_{i:R_i=1}Y_i+\sum_{i:R_i=0}\{m_0(X_i)+\delta\}}{n} = \bar{Y}_0+\frac{n_0}{n}\delta. $$

This identity is for the completed-data mean when the missing outcomes are represented by their conditional means. With stochastic imputations, the same shift adds \(n_0\delta/n\) to each completed-data mean for a fixed set of reference draws. The imputation uncertainty from those draws remains present.

Worked Example: Completed-Data Mean Under a Shift

In a hypothetical sample of 20 patients, 14 scores are observed and sum to 140 points. The six missing scores have reference conditional means that sum to 54 points. Take \(\delta=2\) points for each missing outcome.

The reference completed-data mean is

$$ \bar{Y}_0=\frac{140+54}{20}=\frac{194}{20}=9.70\text{ points}. $$

The shifted missing means sum to \(54+6(2)=66\) points. Thus the shifted completed-data mean is

$$ \bar{Y}_\delta=\frac{140+66}{20} =\frac{206}{20}=10.30\text{ points}. $$

Equivalently, the mean increases by \(6(2)/20=0.60\) points, giving \(9.70+0.60=10.30\) points. The observed sum remains 140 points. This calculation assumes that the same 2-point shift is suitable for all six missing outcomes and that their reference predictions are appropriate for their covariate profiles.

Consequences for Comparisons Between Groups

A common shift does not necessarily cancel from a contrast between groups. Suppose group \(g\) has missing fraction \(q_g\), and its reference completed-data mean is \(\bar{Y}_{g,0}\). When every missing outcome in both groups receives the same \(\delta\), the group mean becomes \(\bar{Y}_{g,0}+q_g\delta\). For the contrast of group A minus group B, the adjusted contrast is therefore the reference contrast plus \((q_A-q_B)\delta\). If missing fractions differ, a common shift changes the contrast.

Worked Example: A Common Shift and a Treatment Contrast

Consider hypothetical trial groups with continuous outcomes measured in points. In group A, 70 of 100 outcomes are observed, with mean 10 points; the 30 missing outcomes have reference mean 8 points. In group B, 75 of 100 outcomes are observed, with mean 9 points; the 25 missing outcomes have reference mean 8.5 points. First calculate the reference completed-data means.

$$ \bar{Y}_{A,0} =\frac{70(10)+30(8)}{100} =\frac{940}{100}=9.40, \qquad \bar{Y}_{B,0} =\frac{75(9)+25(8.5)}{100} =\frac{887.5}{100}=8.875. $$

The reference contrast A minus B is \(9.40-8.875=0.525\) points. Now apply the same \(\delta=2\)-point shift to missing outcomes in both groups. Group A’s mean increases by \((30/100)(2)=0.60\) points, to 10.00. Group B’s mean increases by \((25/100)(2)=0.50\) points, to 9.375. The adjusted contrast is

$$ 10.00-9.375=0.625\text{ points}. $$

The contrast increased by \(0.625-0.525=0.100\) points, which also equals \((0.30-0.25)(2)=0.100\) points. This result follows from the different missing fractions; the common shift does not add equally to the two group means. The calculations describe the point estimates under the stated reference predictions and shift, not the uncertainty interval for the contrast.

Assumptions, Diagnostics, and Alternatives

An additive delta analysis has both a reference-model component and a sensitivity component. The reference model supplies the conditional distributions for missing outcomes; the analyst supplies \(\delta\) to represent a departure from those distributions. A useful analysis states both clearly rather than presenting the shifted result as if it were identified by observed data.

  • Check the reference predictions. Examine whether covariate profiles among missing cases are represented in the data used to fit the reference model. Predictions far outside the observed covariate range may rely heavily on extrapolation.
  • Check the residual assumption. Compare the reference model’s residual spread and shape with observed residuals where outcomes are available. The additive shift retains those features, so a poor residual model remains poor after shifting.
  • Check scale plausibility. A location shift preserves the distribution’s shape and spread. Consider whether a translated distribution remains plausible for the outcome’s possible range and for all relevant covariate profiles.
  • Check the analysis over specified values. Report results across a justified range of \(\delta\) values, including the reference value \(\delta=0\). A single chosen shift can hide how sensitive the conclusion is to the unobserved-outcome assumption.

The additive construction is appropriate when a constant change in the conditional mean, in outcome units, is a meaningful sensitivity departure and retaining the reference residual distribution is defensible. It may be unsuitable when the effect of a departure should vary with the covariates or when a constant translation produces implausible values. In such cases, the model’s assumptions need to be changed deliberately rather than ignored.

The main comparison is with the unshifted MAR reference analysis, which corresponds to \(\delta=0\). Reference-based imputation instead specifies missing outcomes through a reference arm or reference trajectory, as covered in the earlier material on reference-based imputation; it is not generally equivalent to adding a constant to a MAR conditional mean. More broadly, selection models parameterize missingness conditional on outcomes, whereas this pattern-mixture sensitivity construction directly specifies the missing-outcome distribution. These approaches encode different assumptions and should not be treated as interchangeable.

Common Mistake: Shifting More Than the Missing Outcomes

A frequent error is to add \(\delta\) to every outcome, including those observed, or to add the full \(\delta\) to the overall completed-data mean. The first action changes the observed data and no longer represents a departure confined to missing outcomes. The second ignores the proportion missing. For a common shift, the overall mean changes by \(q\delta\), where \(q=n_0/n\), as shown by the completed-data identity above. A second error is to shift the mean and also inflate the residual variance without specifying a different model: that does not retain the reference residual distribution.

Key takeaway: An additive delta adjustment translates the conditional distribution for missing continuous outcomes by \(\delta\), retaining its residual distribution. It leaves observed outcomes unchanged. Its impact on an overall mean or group contrast depends on how many outcomes are missing in each relevant group.

Check Your Understanding

Use the definitions and formulas above. Solutions follow each question.

  1. A reference conditional distribution has mean 24 units and variance 16 units squared. Under an additive shift of \(-3\) units, what are the shifted mean, variance, and standard deviation?
  2. In a sample of 50 people, 10 outcomes are missing. Under a common shift of \(1.5\) units, by how much does the completed-data mean change?
  3. Two groups have missing fractions 0.20 and 0.35. Both receive the same shift \(\delta=2\) units. By how much does the contrast of the first group minus the second group change?
  4. A reference-model draw for a missing outcome is 7.4 units, and the chosen shift is \(-0.6\) units. What is the shifted draw? Should an observed outcome of 8.1 units also be changed?

Solutions. (1) The mean is \(24-3=21\) units. Translation leaves the variance at 16 units squared, so the standard deviation remains \(\sqrt{16}=4\) units. (2) The missing fraction is \(10/50=0.20\), so the mean changes by \(0.20(1.5)=0.30\) units. (3) The contrast changes by \((0.20-0.35)(2)=-0.30\) units; the first-minus-second contrast decreases by 0.30 units. (4) The shifted draw is \(7.4-0.6=6.8\) units. The observed outcome stays at 8.1 units because the adjustment applies only to missing outcomes.

Clinical Biostats