Tutorials › Multiple Imputation by Chained Equations (MICE) › Subgroup-Specific Delta Adjustments

Sensitivity of Conclusions to Imputation Assumptions · Tutorial 428 of 497

Subgroup-Specific Delta Adjustments

Specify and evaluate subgroup-specific departures from a MAR reference while keeping subgroup estimates, contrasts, and uncertainty interpretable.

Advanced 12 min read

What You'll Learn

  • Define subgroup-specific delta parameters relative to a MAR reference.
  • Apply additive shifts to continuous missing outcomes and log-odds shifts to binary outcomes.
  • Calculate how subgroup deltas change subgroup means and contrasts.
  • Check subgroup definitions, missingness counts, model support, and shifted distributions.
  • Distinguish subgroup-specific sensitivity analyses from a single shared delta scenario.
  • Interpret multiple-imputation summaries separately at each fixed delta setting.

Why Allow Different Departures by Subgroup?

A single delta adjustment assumes that, after conditioning on the variables in an imputation model, the departure of missing outcomes from their MAR reference is the same for everyone to whom that adjustment applies. That assumption may be implausible when a scientifically defined subgroup has a distinct reason for missing outcomes. For example, participants with greater baseline symptom burden might have missing follow-up outcomes that are systematically worse than their MAR predictions, while the corresponding departure in a lower-burden subgroup could be smaller or even point in the other direction.

A subgroup-specific adjustment makes those departures explicit. It does not identify which departure is true from the observed data alone. Instead, it defines a set of plausible assumptions, applies each to the missing outcomes in the relevant subgroup, and examines how the subgroup estimates and contrasts respond. As in the earlier material on delta adjustments, the MAR imputation model is the reference; the delta parameters describe departures from that reference for missing outcomes.

Definition: Let \(G\) denote a prespecified subgroup with levels \(g=1,\ldots,K\), and let \(\delta_g\) be the sensitivity parameter for level \(g\). A subgroup-specific delta adjustment modifies the imputation distribution for an outcome that is missing in subgroup \(g\), relative to its MAR reference distribution. The observed outcomes are not shifted. The scale and meaning of \(\delta_g\) depend on the outcome and the chosen adjustment model.

The subgroup should be defined in a way that is meaningful for the analysis and available when the adjustment is assigned. A baseline characteristic such as a prespecified age category can often serve this purpose. If subgroup membership itself is missing, it is not generally valid to assign a fixed subgroup adjustment using an unknown group label. The uncertainty in membership must instead be handled in a coherent imputation strategy. Subgroups should not be selected after inspecting which definition produces the most favorable result.

Specification and Assumptions

For a continuous outcome on its original scale, an additive adjustment is \(Y^{\mathrm{mis}}_{\delta}=Y^{\mathrm{mis}}_{\mathrm{MAR}}+\delta_g\) for a missing outcome in subgroup \(g\). Here \(\delta_g\) is measured in the outcome's units. For a binary outcome, a log-odds adjustment can be written as \(\operatorname{logit}(p_{\delta,g})=\operatorname{logit}(p_{\mathrm{MAR}})+\delta_g\), where \(p_{\mathrm{MAR}}\) is the reference probability for a missing outcome given the observed information. Thus, a delta of zero recovers the MAR reference, and a negative log-odds delta lowers the probability of the event.

These adjustments require more than correct arithmetic. The MAR reference model must be a defensible model for the observed data, including relevant predictors and the subgroup structure. The chosen scale must make sense for the outcome, and the resulting distributions must be plausible and respect the outcome's support. The subgroup definition and the rationale for its delta should be specified before interpreting the sensitivity results. A delta does not correct a poor MAR model, and a different delta in each subgroup does not by itself guarantee that the assumed departures are realistic.

For each fixed vector \((\delta_1,\ldots,\delta_K)\), apply the appropriate adjustment to missing outcomes in each subgroup in every imputed data set, then calculate the analysis estimate in each completed data set. Re-estimate subgroup estimates and their contrasts under that scenario. If \(Q_m\) is an estimate from imputation \(m\), and \(U_m\) its estimated variance, the usual multiple-imputation summaries are \(\bar Q=M^{-1}\sum_{m=1}^M Q_m\), \(\bar U=M^{-1}\sum_{m=1}^M U_m\), \(B=(M-1)^{-1}\sum_{m=1}^M(Q_m-\bar Q)^2\), and \(T=\bar U+(1+M^{-1})B\). The sensitivity parameter is held fixed within a scenario; uncertainty across imputations reflects imputation and analysis uncertainty under that scenario, not uncertainty about which delta is correct.

How Subgroup Adjustments Change Estimates

A useful justification is available for an additive adjustment when the target is a subgroup mean. Let subgroup \(g\) contain \(N_g\) individuals, of whom \(n_{\mathrm{mis},g}\) have missing outcomes. If every missing outcome is shifted by the same \(\delta_g\), the sum of the completed outcomes increases by \(n_{\mathrm{mis},g}\delta_g\), while the observed outcomes remain unchanged. Dividing by the subgroup size gives

$$ \bar Y_{\delta,g} = \bar Y_{\mathrm{MAR},g} + \frac{n_{\mathrm{mis},g}}{N_g}\delta_g. $$

The formula shows why the effect depends both on the delta and on the fraction missing in the subgroup. For two subgroups \(A\) and \(B\), the adjusted contrast \(B-A\) changes by \(\{n_{\mathrm{mis},B}/N_B\}\delta_B-\{n_{\mathrm{mis},A}/N_A\}\delta_A\). This is a point-estimate identity for the stated additive setup. With stochastic imputations, analysis-model estimates, nonlinear estimands, or other adjustment scales, calculate the estimand in each completed data set rather than assuming this simple identity applies.

Worked Example: Binary Outcome With Different Log-Odds Deltas

Consider a hypothetical follow-up study with two prespecified baseline symptom subgroups, A and B. In subgroup A, 6 of 12 observed outcomes are events and 8 of 20 outcomes are missing. The MAR reference event probability for those missing outcomes is 0.30. In subgroup B, 12 of 15 observed outcomes are events and 10 of 25 outcomes are missing; their MAR reference probability is 0.60. Suppose the sensitivity scenario uses log-odds deltas \(\delta_A=-0.50\) and \(\delta_B=0.40\). Treat the reference probabilities as given for this illustration and apply the same subgroup delta to each missing outcome.

For subgroup A, the reference odds are \(0.30/(1-0.30)=0.4286\). The adjusted odds are \(0.30/0.70\times\exp(-0.50)=0.2599417\ldots\), so the adjusted event probability is \(0.2599417/(1+0.2599417)\approx0.2063\).

The expected completed event proportion under MAR is \((6+8(0.30))/20=8.4/20=0.4200\). Under the adjusted probability it is \((6+8(0.206312\ldots))/20\approx7.6505/20=0.3825\). The negative delta lowers the expected subgroup A event proportion by about \(0.0375\), or 3.75 percentage points.

For subgroup B, the reference odds are \(0.60/0.40=1.50\). The adjusted odds are \(1.50\exp(0.40)=2.237737\ldots\), giving probability \(2.237737/(1+2.237737)\approx0.6911\). The expected completed proportion is \((12+10(0.691142\ldots))/25\approx18.9114/25=0.7565\), compared with its MAR value \((12+10(0.60))/25=18/25=0.7200\).

The \(B-A\) contrast is \(0.7565-0.3825=0.3739\), using unrounded probabilities for the calculation. Its MAR-reference value is \(0.7200-0.4200=0.3000\); this scenario increases the contrast by about 0.0739, or 7.39 percentage points. These are expected proportions for the stated illustrative calculations, not a claim that fractional events are observed. In actual multiple imputation, binary outcomes are drawn and the estimates and uncertainty are calculated across completed data sets. The contrast describes a difference between subgroups, not a causal effect.

Worked Example: Additive Adjustment for a Continuous Subgroup Mean

Suppose, hypothetically, a subgroup has 10 people: 6 observed outcomes with values 10, 12, 14, 16, 18, and 20 units, and 4 missing outcomes. The MAR imputation model gives conditional mean predictions of 13, 15, 17, and 19 units for the missing outcomes. Consider an additive delta of \(-2\) units for this subgroup. Assume the continuous outcome can validly take all shifted values and that the MAR predictions are appropriate for this illustration.

The observed sum is \(10+12+14+16+18+20=90\) units. The MAR predictions sum to \(13+15+17+19=64\) units, so the MAR completed mean is \((90+64)/10=15.4\) units. Applying the subgroup delta gives adjusted missing values \(11,13,15,17\), whose sum is 56 units. The adjusted subgroup mean is \((90+56)/10=14.6\) units.

The change is \(-0.8\) units, which agrees with the identity \((4/10)(-2)=-0.8\). This checks both the direction and the magnitude: only four of the ten outcomes were shifted, so the subgroup mean changes by less than the full two-unit shift. For real imputations, these predictions would ordinarily be subject to imputation uncertainty; a fixed delta changes the assumed missing-outcome distribution, not the observed values.

Contrasts, Diagnostics, and Alternatives

Separate subgroup estimates are often reported because an overall average can hide different sensitivity patterns. The contrast must be computed from estimates produced under the same delta scenario. A useful sensitivity analysis includes the MAR reference vector, with all deltas set to zero, and a set of scientifically defensible alternatives. For continuous outcomes, the preceding identity helps anticipate how the contrast changes. For binary outcomes, odds shifts alter probabilities nonlinearly, so recompute the subgroup proportions or model-based estimands at each delta setting.

Before interpreting results, check the subgroup-specific missingness counts and fractions, the distribution of observed outcomes and predictors, and whether the MAR model has adequate information within each subgroup. Sparse data can make subgroup-specific reference predictions unstable. Inspect the shifted distributions, especially for bounded outcomes or unusual tail behavior, and verify that group assignments are applied consistently in every imputation. If the subgroup-specific deltas are so large that shifted values become implausible, that is evidence to reconsider the range or the scale, not a reason to overlook the problem.

A shared delta is simpler and may be appropriate when the same departure is scientifically plausible across subgroups. Pattern-specific adjustments instead distinguish missingness patterns; a subgroup adjustment distinguishes defined population groups, and those classifications need not coincide. Subgroup-specific deltas are more flexible but require more assumptions and can yield poorly supported estimates in small subgroups. They are not a substitute for a richer, well-specified MAR model or for reporting uncertainty.

Common pitfall: Applying a subgroup delta to every member of that subgroup, including people whose outcomes were observed, changes the data rather than the missing-outcome assumption. Apply the adjustment only to missing outcomes. Also avoid interpreting a changed subgroup contrast as proof that the subgroups truly differ by that amount: it is a consequence of the chosen sensitivity assumptions.

A second common error is to compare subgroup estimates calculated under different delta settings as if they formed one coherent contrast. For a scenario-based comparison, specify the full delta vector and apply it consistently. Report the subgroup definitions, delta scale and values, the MAR reference, and the resulting estimates and contrasts. This makes clear which conclusions are supported across assumptions and which depend on a particular departure.

Worked Example: A Two-Subgroup Continuous Contrast

In a hypothetical study, subgroup A has 20 people: 12 observed outcomes with mean 50 units and 8 missing outcomes with MAR predicted mean 48 units. Subgroup B has 25 people: 15 observed outcomes with mean 60 units and 10 missing outcomes with MAR predicted mean 58 units. Suppose the additive deltas are \(\delta_A=-3\) units and \(\delta_B=2\) units. Assume both subgroups are defined from observed baseline information and the shifted outcome values are plausible.

For A, the observed sum is \(12(50)=600\), and the MAR missing-outcome sum is \(8(48)=384\). The MAR subgroup mean is \((600+384)/20=49.2\) units. After shifting the missing outcomes, their expected sum is \(8(48-3)=360\), giving adjusted mean \((600+360)/20=48.0\) units.

For B, the observed sum is \(15(60)=900\), and the MAR missing-outcome sum is \(10(58)=580\). The MAR subgroup mean is \((900+580)/25=59.2\) units. Under the delta, the missing-outcome sum is \(10(58+2)=600\), giving adjusted mean \((900+600)/25=60.0\) units.

The MAR \(B-A\) contrast is \(59.2-49.2=10.0\) units. Under the subgroup-specific scenario it is \(60.0-48.0=12.0\) units, a 2.0-unit increase. The contrast-change identity confirms this: \((10/25)(2)-(8/20)(-3)=0.8+1.2=2.0\) units. The overall mean also changes, but by a different amount: the MAR total is \(984+1480=2464\) units across 45 people, or \(54.76\) units; the adjusted total is \(960+1500=2460\) units, or \(54.67\) units. Subgroup contrasts and overall averages answer different questions, so one should not be substituted for the other.

Check Your Understanding

Answer each question, then compare with the solutions.

  1. A subgroup has 30 people, 9 missing continuous outcomes, and an additive delta of 4 units. By how much does its completed mean change from the MAR mean?
  2. For a binary outcome, a subgroup's MAR probability is 0.20 and its log-odds delta is \(\log(2)\). Find the adjusted probability.
  3. Why should observed outcomes not receive a subgroup-specific delta adjustment?
  4. Two subgroup-specific deltas produce a larger contrast than the MAR reference. Does this establish that the subgroups truly differ by the larger amount? Explain.

Solutions. (1) The mean changes by the missing fraction times the delta: \((9/30)(4)=1.2\) units. (2) The MAR odds are \(0.20/0.80=0.25\). Multiplying by \(\exp\{\log(2)\}=2\) gives adjusted odds \(0.50\), and the probability is \(0.50/(1+0.50)=0.3333\). (3) The delta represents an assumption about outcomes that were not observed, relative to their MAR reference; shifting observed values would alter known data rather than describe missing outcomes. (4) No. The larger contrast is the estimate under the specified sensitivity assumptions. Observed data alone do not establish that those deltas, or the resulting contrast, are true.

Clinical Biostats