Tutorials › Multiple Imputation by Chained Equations (MICE) › Evaluating Sensitivity Curves for Pooled Estimates

Sensitivity of Conclusions to Imputation Assumptions · Tutorial 441 of 497

Evaluating Sensitivity Curves for Pooled Estimates

Learn to evaluate how a pooled estimate changes across a delta grid, and distinguish meaningful sensitivity from curvature, grid artifacts, and Monte Carlo noise.

Advanced 10 min read

What You'll Learn

  • Define a pooled-estimate sensitivity curve over a one-dimensional delta grid
  • Calculate secant slopes and finite-difference curvature from plotted estimates
  • Interpret curve direction and steepness in the units of the estimand
  • Set a practical-stability tolerance and assess it at evaluated grid points
  • Diagnose sparse grids, Monte Carlo noise, and inconsistent analyses
  • Distinguish point-estimate sensitivity curves from inferential summaries

Why Evaluate a Sensitivity Curve?

A delta adjustment represents a specified departure from the imputation model’s reference assumptions for missing values. When the adjustment varies over a range of plausible values, the resulting pooled estimates show how sensitive the substantive estimate is to that choice. A sensitivity curve makes this pattern visible: its direction indicates how the estimate changes as delta changes, while its steepness and curvature describe how that change behaves across the grid.

The curve is a descriptive summary of analyses under a sequence of assumptions. It does not establish which delta is correct, and it does not by itself summarize uncertainty in the estimates. Its usefulness depends on making the delta’s meaning clear, holding the rest of the analysis fixed, and plotting enough values to reveal the pattern without implying more precision than the grid supports.

Definition: Let \(d_1,\ldots,d_K\) be a one-dimensional grid of delta values. At each \(d_k\), apply the specified delta adjustment, analyze the resulting imputed data, and obtain the pooled point estimate \(\bar{Q}(d_k)\) for the same estimand. The plotted pairs \((d_k,\bar{Q}(d_k))\) form a sensitivity curve for the pooled estimate. The horizontal axis is delta, and the vertical axis is the pooled estimate.

Here, “pooled estimate” means the point estimate obtained by combining the analysis estimates from the imputations at that particular delta value. Pooling is performed separately at each grid value; estimates from different delta values are not pooled together. The curve’s vertical coordinate must refer to the same estimand, analysis model, and units at every point.

Conditions for a Meaningful Curve

A curve is interpretable only if its points differ in the intended way. The delta parameter should have a stated interpretation and scale, the grid should cover a range justified for the application, and all other choices should be held fixed unless the purpose is explicitly to vary them too. For example, changing both delta and the analysis model across points would make the curve unable to isolate sensitivity to delta.

The imputation and analysis procedures must also be applied consistently at each grid value. In practice, this means checking that the delta adjustment was implemented as intended, the chained equations have converged sufficiently for the analysis, and there are enough imputations for the pooled point estimates to be stable. These checks do not prove that the missing-data assumptions are correct. They help ensure that the plotted changes reflect the chosen delta scenarios rather than implementation errors or unstable computation.

The curve is a finite set of evaluated points. Connecting adjacent points with line segments is a visual aid, not proof that the estimate changes linearly between them. A smooth curve fitted through sparse points can suggest unsupported behavior, particularly near the ends of the grid. If the pattern between grid values matters, evaluate additional delta values there.

Direction, Slope, and Curvature

For two evaluated values \(d_a\) and \(d_b\), the secant slope is the change in pooled estimate divided by the change in delta. It describes the average rate of change over that interval, not necessarily the instantaneous rate at any point. Its units are the estimate’s units per unit of delta. If delta and the estimate share units, the slope is often dimensionless.

$$ \text{Secant slope from }d_a\text{ to }d_b = \frac{\bar{Q}(d_b)-\bar{Q}(d_a)}{d_b-d_a} $$

On an equally spaced grid with spacing \(h\), adjacent secant slopes can be compared to assess curvature. An increase in successive slopes indicates that the curve is bending upward over those grid intervals; a decrease indicates downward bending. The finite second difference at an interior grid point is

$$ \Delta^2 \bar{Q}(d_k) = \bar{Q}(d_{k+1})-2\bar{Q}(d_k)+\bar{Q}(d_{k-1}). $$

Dividing this second difference by \(h^2\) gives a finite-difference approximation to a second derivative when the curve is sufficiently smooth. Without that smoothness assumption, the second difference is still a useful description of how neighboring secant slopes change, but it should not be presented as an exact derivative.

Worked Example: Calculating Slopes and Curvature

Suppose a hypothetical study estimates a treatment-control contrast in systolic blood pressure, measured in mmHg. Delta is also measured in mmHg, and all analysis choices are fixed across the grid. The pooled estimates are:

Delta, \(d\) (mmHg)Pooled estimate, \(\bar{Q}(d)\) (mmHg)
-44.80
-25.40
05.90
26.40
46.90

The grid spacing is \(h=2\) mmHg. The four adjacent secant slopes are:

From \(-4\) to \(-2\): \((5.40-4.80)/2=0.60/2=0.300\). From \(-2\) to \(0\): \((5.90-5.40)/2=0.50/2=0.250\). From \(0\) to \(2\): \((6.40-5.90)/2=0.50/2=0.250\). From \(2\) to \(4\): \((6.90-6.40)/2=0.50/2=0.250\).

Thus, over the evaluated range, the estimate increases as delta increases. The slope is largest on the first interval and is the same on the remaining three. At the interior point \(d=-2\), the second difference is \(5.90-2(5.40)+4.80=5.90-10.80+4.80=-0.10\) mmHg. Dividing by \(h^2=4\) gives \(-0.0250\) per mmHg as a finite-difference curvature approximation there. At \(d=0\), the second difference is \(6.40-2(5.90)+5.40=0\); at \(d=2\), it is \(6.90-2(6.40)+5.90=0\). The data therefore show a small change in slope near the lower end and constant adjacent slopes thereafter, not evidence of a strongly curved relationship.

Practical Stability Is a Substantive Judgment

A curve can change without changing enough to matter for a particular application. To describe practical stability, specify a tolerance in the units of the estimand and a reference estimate, such as the estimate at \(d=0\). The tolerance should be justified by clinical, policy, or scientific considerations before interpreting the plotted pattern. It is not a universal statistical cutoff.

For example, one may call estimates practically stable at evaluated grid points when their absolute difference from the reference estimate does not exceed a prespecified tolerance. This statement concerns the point estimates and the chosen tolerance; it is not a claim that the estimates are statistically indistinguishable. The curve alone cannot answer that inferential question.

Worked Example: Applying a Practical-Stability Tolerance

Use the hypothetical estimates in the preceding table, and suppose subject-matter experts specify that differences of at most \(0.50\) mmHg from the estimate at \(d=0\) are practically negligible. The reference estimate is \(\bar{Q}(0)=5.90\) mmHg. The absolute differences at the five grid points are:

At \(d=-4\), \(|4.80-5.90|=1.10\) mmHg. At \(d=-2\), \(|5.40-5.90|=0.50\) mmHg. At \(d=0\), the difference is \(0\) mmHg. At \(d=2\), \(|6.40-5.90|=0.50\) mmHg. At \(d=4\), \(|6.90-5.90|=1.00\) mmHg.

Under the stated rule, the estimates at \(d=-2,0,\) and \(2\) meet the tolerance, while those at \(d=-4\) and \(4\) do not. This supports describing the estimates as practically stable at those three evaluated values relative to the \(d=0\) estimate. It does not establish stability at every value between \(-2\) and \(2\): that would require either a justified interpolation assumption or additional grid evaluations.

Reading the Plot and Checking Its Construction

Plot delta on the horizontal axis and the pooled estimate on the vertical axis. Label both axes with units, identify the reference value if one is used, and show the evaluated points. A line joining the points can help the reader follow the sequence, but the points should remain visible so the grid is not mistaken for a continuous, exactly known function. If comparing prespecified scenarios, use clear labels and ensure each curve represents the same estimand and analysis.

A curve that rises or falls steadily indicates directional sensitivity over the evaluated range, but the direction need not be assumed in advance. The impact of a delta adjustment on an estimand depends on which values are missing, how the adjustment is applied, and how the estimand is defined. Curvature can signal that the effect of changing delta varies across the range. It may also arise from nonlinear analysis or from Monte Carlo variation, so visual shape should be checked rather than interpreted automatically as a substantive feature.

  • Confirm that the plotted delta values match the values used in the imputation procedure.
  • Check that every point uses the same analysis population, estimand, model, and units.
  • Inspect the pattern of point estimates as well as the connecting line; do not infer behavior beyond the evaluated range.
  • If small irregularities are important to the interpretation, assess whether they persist with more imputations or repeated runs.
  • Document the delta range and its rationale so that the curve’s scope is clear.

Worked Example: Distinguishing a Trend from Monte Carlo Irregularity

In a hypothetical analysis, the team repeats the imputation and analysis at three delta values to check whether a small bend is reproducible. The pooled estimates from two runs are:

Delta (units)Run 1 estimateRun 2 estimateMean of the two runs
-25.395.415.40
05.885.925.90
26.376.436.40

The mean estimates at \(-2,0,2\) are \(5.40,5.90,6.40\). The slope from \(-2\) to \(0\) is \((5.90-5.40)/2=0.250\), and the slope from \(0\) to \(2\) is \((6.40-5.90)/2=0.250\). The differences between runs at the three delta values are \(0.02,0.04,\) and \(0.06\), respectively. In this small illustration, the estimates rise consistently, while the run-to-run differences are small relative to the \(0.50\) increase across each interval.

Two runs are not enough to quantify Monte Carlo error reliably. The example only shows why a suspicious bend or jagged point should prompt a computational check rather than an immediate substantive explanation. More imputations or repeated runs may be needed to determine whether the pattern is stable.

Common Pitfalls and Limitations

A frequent mistake is to read a steep slope as proof that the missing-data assumptions are implausible. A steep slope shows that the pooled estimate responds strongly to the chosen delta over that interval; plausibility must be judged using the meaning of delta and relevant subject-matter information. Conversely, a nearly flat curve over a narrow range does not show robustness to values outside that range.

Another mistake is to treat a point-estimate curve as a complete uncertainty analysis. It contains no interval information by itself and cannot establish whether a change is distinguishable from sampling or imputation uncertainty. It also cannot identify a uniquely correct delta. A curve complements, rather than replaces, the substantive justification for the sensitivity range and the appropriate inferential summaries.

A table of estimates is a reasonable alternative when only a few delta values are evaluated and exact values matter most. A plot is usually more informative for seeing direction, curvature, and gaps in the grid. Neither format compensates for a poorly chosen range, inconsistent analyses, or unstable imputations. Where several sensitivity parameters vary together, a one-dimensional curve cannot show their joint behavior; the interpretation must remain limited to the single parameter varied.

Key takeaway: Evaluate a sensitivity curve by checking how the pooled point estimate changes across a justified delta grid. Report direction, interval-specific slopes, visible curvature, and any prespecified practical-stability tolerance in context, while keeping conclusions within the range actually evaluated.

Check Your Understanding

Answer each question using the definitions and calculations in this tutorial. Solutions follow each question.

  1. A pooled estimate is \(12.0\) units at delta \(0\) and \(13.2\) units at delta \(3\). What is the secant slope, and how should its units be stated?
    Solution: The slope is \((13.2-12.0)/(3-0)=1.2/3=0.400\) estimate units per delta unit. If both quantities use the same units, the ratio is dimensionless, but reporting “0.400 estimate units per delta unit” makes the scale explicit.
  2. On an equally spaced grid, three consecutive estimates are \(8.0, 8.5,\) and \(9.2\). What is the second difference at the middle point, and what does its sign indicate?
    Solution: The second difference is \(9.2-2(8.5)+8.0=9.2-17.0+8.0=0.20\). It is positive, indicating that the adjacent slopes increase across these intervals, so the finite differences suggest upward curvature locally.
  3. A reference estimate is \(2.40\), and the prespecified practical tolerance is \(0.30\) units. Does an estimate of \(2.05\) meet the rule?
    Solution: The absolute difference is \(|2.05-2.40|=0.35\) units. Since \(0.35>0.30\), it does not meet the stated practical-stability rule.
  4. Why is it not enough to plot estimates from several delta values if the analysis model also changes at each value?
    Solution: The plotted differences could reflect changes in the model as well as changes in delta. To interpret the curve as sensitivity to delta, hold the analysis model and other relevant choices fixed.
  5. A curve looks flat across the evaluated values. What does that establish, and what does it not establish?
    Solution: It describes little change in the point estimates over the evaluated grid, subject to computational stability. It does not establish robustness outside that range, prove the assumptions are correct, or provide an uncertainty analysis.
Clinical Biostats