Tutorials › Multiple Imputation by Chained Equations (MICE) › Reporting Delta Assumptions and Tipping Results

Sensitivity of Conclusions to Imputation Assumptions · Tutorial 461 of 497

Reporting Delta Assumptions and Tipping Results

Specify the assumptions and analysis decisions that make a delta-adjusted sensitivity analysis reproducible and its reported tipping point interpretable.

Advanced 12 min read

What You'll Learn

  • State what quantity delta shifts, its scale, and which missing values receive the adjustment.
  • Describe the reference scenario and justify the explored delta range.
  • Report how estimates and uncertainty were pooled at each delta value.
  • Define a tipping point using an explicit estimand, interval, and decision rule.
  • Distinguish sensitivity-analysis results from conclusions that depend on unreported assumptions.

Why Reporting Choices Matter

A tipping-point number is not interpretable on its own. A reported value such as \(\delta=-1.8\) could mean a shift of 1.8 outcome units, a log-odds shift, or something else entirely. Its direction also depends on which group and which missing values receive the adjustment. To reproduce or assess the result, a reader needs the assumptions that give delta its meaning, the reference scenario, the range examined, the analysis used to pool estimates, and the exact rule defining a “tip.”

This tutorial focuses on documenting those choices. It does not treat a delta as an empirically identified feature of the data: delta encodes a sensitivity assumption about missing outcomes that cannot be checked from the observed outcomes alone. The earlier material on Assessing the Plausibility of Extreme Delta Scenarios concerns the substantive plausibility of scenarios; reporting should make the assumptions and the basis for the chosen range visible rather than imply that the data selected a uniquely correct delta.

Definition: A delta assumption specifies how the imputation distribution for a targeted set of missing values is shifted from a stated reference imputation model. A reference scenario is the specified baseline, commonly the unadjusted imputation model with \(\delta=0\). A tipping-point definition is a prespecified rule for identifying when the reported conclusion changes as delta varies.

What a Reproducible Report Specifies

State the missing variable, the missingness pattern or group to which the adjustment applies, and the scale and direction of delta. For a continuous outcome, an additive shift might replace a reference imputation \(Y_{\mathrm{mis}}^{(0)}\) with \(Y_{\mathrm{mis}}^{(\delta)}=Y_{\mathrm{mis}}^{(0)}+\delta\); delta then has the outcome’s units. For a binary outcome, a logit shift can be written as \(\operatorname{logit}\{P(Y=1\mid X,\mathrm{mis})\}=\operatorname{logit}\{P(Y=1\mid X,\mathrm{ref})\}+\delta\). Here delta is dimensionless: the corresponding odds are multiplied by \(e^\delta\). These are different models, so a report must not describe both merely as “a shift by delta.”

Name the reference scenario and explain what \(\delta=0\) means in this analysis. Zero means no additional delta adjustment relative to the stated reference imputation model; it does not, by itself, establish that the model is true or that missingness is ignorable. Report the explored range, its increment or search method, and the substantive reasoning for its bounds. If ranges differ across groups or outcomes, specify each one separately. A range should include enough detail for a reader to see whether the tipping point lies inside the examined region or whether the analysis only establishes that it was not found there.

At each delta, describe how the imputed data were analyzed and how the resulting estimates and variances were combined. With \(M\) imputations, let \(\hat Q_j(\delta)\) be the estimate and \(U_j(\delta)\) its estimated variance in imputation \(j\). The pooled estimate is \(\bar Q(\delta)\), the mean of the \(M\) estimates. Rubin’s combining rules use the average within-imputation variance \(\bar U(\delta)\), the between-imputation variance \(B(\delta)\), and the total variance \(T(\delta)\):

$$ \bar Q(\delta)=\frac{1}{M}\sum_{j=1}^{M}\hat Q_j(\delta),\qquad \bar U(\delta)=\frac{1}{M}\sum_{j=1}^{M}U_j(\delta),\qquad B(\delta)=\frac{1}{M-1}\sum_{j=1}^{M}\{\hat Q_j(\delta)-\bar Q(\delta)\}^2, \qquad T(\delta)=\bar U(\delta)+\left(1+\frac{1}{M}\right)B(\delta). $$

Report the number of imputations and the pooling and interval procedure, including the degrees-of-freedom method when it affects the interval. State whether each delta scenario was refitted and pooled independently, or whether any quantities were held fixed. In general, estimates and variances should be recomputed under each scenario; carrying a reference-scenario standard error across a delta grid is an extra approximation, not an automatic property of delta adjustment.

Worked Examples of Reporting Decisions

Worked Example: Describing an Additive Delta and Its Direction

Suppose a hypothetical trial estimates a treatment-minus-control difference in a continuous symptom score, measured in points. Only missing outcomes in the treatment group receive an additive delta shift; control-group imputations remain at their reference model. The fraction of treatment-group outcomes that is missing is \(q_T=0.30\). At the reference scenario, the estimated contrast is 2.0 points.

If the reference imputed mean is shifted by \(\delta\) points for each missing treatment outcome, the all-participant treatment mean changes by \(q_T\delta\), assuming the observed treatment outcomes and control mean are unchanged. Therefore, in this simplified setup, the contrast changes by the same amount:

$$ \widehat{\Delta}(\delta)=2.0+0.30\delta\ \text{points}. $$

For example, at \(\delta=-2.0\) points, the adjusted contrast is \(2.0+0.30(-2.0)=1.4\) points. A negative delta makes the imputed treatment outcomes lower than their reference imputations, and the estimated treatment-minus-control contrast falls. This relationship follows because \(0.30\) of the treatment-group outcomes receive the shift: the mean changes by \(0.30\times(-2.0)=-0.60\) points. In an actual analysis, model refitting and uncertainty calculations may change other quantities too; this algebra describes the specified simple mean contrast, not a universal formula for every estimand.

A reproducible report would say, for instance: “We added \(\delta\) points to imputed symptom scores in the treatment group only; \(\delta=0\) denotes the reference imputation model. We examined \(-3\) to \(2\) points in 0.1-point increments, with negative values representing worse imputed treatment outcomes.” The units, target, reference, direction, and range are all explicit.

Worked Example: Reporting Pooling at a Binary-Outcome Scenario

Consider a hypothetical binary outcome analysis in which delta is a logit shift for missing outcomes in one prespecified group. At a particular delta, suppose five completed-data analyses give log-odds-ratio estimates \(0.30, 0.26, 0.34, 0.22,\) and \(0.28\), each with estimated within-imputation variance \(0.040\). The log-odds-ratio scale is used for pooling.

The pooled estimate is \(\bar Q=(0.30+0.26+0.34+0.22+0.28)/5=1.40/5=0.280\). Since all five within-imputation variances are \(0.040\), \(\bar U=0.040\). The deviations from \(0.280\) are \(0.020,-0.020,0.060,-0.060,0.000\). Thus:

$$ B=\frac{0.020^2+(-0.020)^2+0.060^2+(-0.060)^2+0^2}{5-1} =\frac{0.0080}{4}=0.0020, \qquad T=0.040+1.2(0.0020)=0.0424, \qquad SE=\sqrt{0.0424}=0.2059. $$

For a transparent large-sample illustration using a 1.96 multiplier, the 95% interval is \(0.280\pm1.96(0.2059)=[-0.124,0.684]\) on the log-odds-ratio scale. Exponentiating gives an odds-ratio estimate of \(e^{0.280}=1.323\) and interval \([e^{-0.124},e^{0.684}]=[0.883,1.982]\). The interval includes an odds ratio of 1. A report should identify the interval method actually used; with finite degrees of freedom, a \(t\) multiplier may be appropriate instead of 1.96. Five imputations are used here only to make the arithmetic short, not as a general recommendation.

To make the scenario reproducible, the report must also state that delta shifts the log odds, which missing group receives it, and what zero means. For example, \(\delta=0.405\) multiplies the odds for the targeted missing outcomes by \(e^{0.405}\approx1.50\) relative to the reference imputation model. Reporting only “delta was 0.405” would leave the scale and its practical interpretation unclear.

Worked Example: Defining a Statistical Tipping Point

Suppose the estimand is a treatment-minus-control mean difference in points, with a two-sided 95% confidence interval. The prespecified statistical conclusion is that the difference is positive when the entire interval is above zero. At each delta, the analysis is rerun and pooled. In a hypothetical search, the lower interval limit is 0.08 points at \(\delta=0.18\) and \(-0.03\) points at \(\delta=0.21\).

Under the stated rule, the positive conclusion holds at 0.18 but no longer holds at 0.21 because the interval then includes zero. If a finer search estimates the first crossing at \(\delta=0.20\), the report should define this as the smallest delta in the examined direction at which the lower 95% limit is at or below zero. It should also say whether the value is a grid result, an interpolated estimate, or a root-finding result. The number is tied to this confidence-interval rule; another rule, such as a practical-effect threshold, can yield a different tipping point.

This example does not imply that the null hypothesis becomes true at \(\delta=0.20\). It identifies the sensitivity assumption at which the specified statistical conclusion changes. A complete report gives the endpoint, contrast direction, interval level, rule, delta direction and scale, and uncertainty or numerical precision for the estimated crossing when relevant.

Justifying the Mean-Shift Relationship

The arithmetic in the continuous example can be derived directly. Suppose \(n_T\) treatment outcomes contribute to the target mean, \(m_T\) are missing, and every missing outcome is increased by \(\delta\), while observed outcomes are unchanged. The sum of the completed treatment outcomes increases by \(m_T\delta\). Dividing by \(n_T\), the treatment mean increases by \((m_T/n_T)\delta=q_T\delta\). If the control mean is not shifted, subtracting it leaves the contrast change \(q_T\delta\). This derivation depends on a simple difference in means and a common additive shift; for regression coefficients, transformed estimands, or other patterns, the effect need not be this direct.

The underlying analysis also requires that the imputation and analysis models be specified for the intended estimand, that the delta adjustment be applied as described, and that uncertainty be pooled appropriately at each scenario. The sensitivity analysis does not make an untestable missing-outcome assumption testable. Instead, it reveals how conclusions vary across the stated alternatives. Diagnostics should therefore check that the generated values obey the outcome’s support and that the adjustment has the intended direction and magnitude. For example, an additive shift should not silently create impossible scores, and a logit shift should be interpreted on the odds scale rather than as a fixed probability increment.

Common Reporting Failures and Alternatives

A frequent mistake is to report a tipping point without defining the rule: “The tipping point was \(-1.8\)” does not tell the reader whether the estimate crossed zero, a confidence interval included zero, or a practical threshold was crossed. Another is to report only a delta range without specifying its units or the target missing values. Both errors prevent interpretation even if the numerical calculation itself is correct.

Useful diagnostics include checking the reference result against the primary analysis, inspecting estimates and interval limits across the full explored range, and confirming that the reported crossing is consistent with neighboring scenarios. Report if a tipping point was not found within the range; do not imply that it cannot occur outside it. Also disclose grid resolution or numerical tolerance, the number of imputations, and any Monte Carlo uncertainty that could materially affect the crossing. A coarse grid can locate a transition only to its resolution; a more precise-looking decimal does not overcome that limitation.

The main alternative is to report a set of prespecified scenarios rather than a single crossing. Scenario results can be easier to interpret when stakeholders can judge plausible assumptions directly, while a tipping-point summary compactly describes where a stated decision rule changes. These approaches answer different questions and can be reported together. Neither replaces the reference analysis, and neither establishes which untestable delta assumption is correct. The earlier material on Separating Statistical and Practical Tipping Points is relevant when the decision rule concerns both statistical evidence and a meaningful effect size.

Reporting checklist: Identify the estimand and missing values targeted; define delta’s scale, direction, and reference; give the range and rationale; state the imputation and pooling procedure, number of imputations, and interval method; define the tipping rule; and report the crossing, numerical resolution, and any important limitations.

Check Your Understanding

Answer each question using the reporting principles above. Brief solutions follow the questions.

  1. A continuous outcome is shifted by \(\delta\) only among missing outcomes in one group. What information is needed to interpret a reported value of \(\delta=-1.5\)?
  2. For a binary outcome, delta is added to the logit of the event probability. By what factor do the odds change when \(\delta=\log(2)\)?
  3. An analysis reports a tipping point at \(\delta=0.4\) but gives no decision rule. What essential definition is missing?
  4. A search from \(-2\) to \(2\) finds no tipping point. What may the report conclude, and what should it avoid claiming?

Solutions. (1) State the outcome units, which missing values or group are shifted, the direction of the adjustment, and the reference scenario; also give the explored range and analysis context. (2) The odds are multiplied by \(e^{\log(2)}=2\), so they double relative to the reference model. (3) The report must define what changes at the tip, including the estimand, interval or threshold, and crossing rule. (4) It may state that no tipping point was found within the examined range under the specified rule. It should not claim that no tipping point exists outside that range or under other assumptions.

Clinical Biostats