Tutorials › Multiple Imputation by Chained Equations (MICE) › Tipping Points for Effect Direction and Magnitude

Sensitivity of Conclusions to Imputation Assumptions · Tutorial 445 of 497

Tipping Points for Effect Direction and Magnitude

Calculate and interpret separate delta thresholds for an estimated effect reaching zero and for its magnitude crossing a prespecified minimum important effect.

Advanced 9 min read

What You'll Learn

  • Define direction and magnitude tipping points using a sensitivity curve for the estimate.
  • Calculate both thresholds when the curve is affine in the delta parameter.
  • Interpret thresholds in the units and direction of the chosen delta scenario.
  • Identify multiple crossings, tangencies, and cases where a threshold is not reached.
  • Distinguish point-estimate thresholds from statistical-significance and clinical or policy thresholds.

Why Direction and Magnitude Need Separate Thresholds

A sensitivity analysis can ask more than whether a conclusion remains statistically significant. It can ask how much the estimated effect must change before its direction reverses, and how much before its magnitude falls below a prespecified minimum important effect. These questions concern different boundaries: zero for direction, and a nonzero threshold for magnitude. Their tipping points need not coincide.

Let \(Q(\delta)\) denote the estimate of the target effect under a specified delta-adjusted imputation scenario, with \(\delta=0\) denoting the reference scenario. The estimate may be a pooled estimate from the chosen analysis. Keep the analysis, estimand, and delta-adjustment convention fixed while evaluating this sensitivity curve; otherwise, a change in \(Q\) cannot be attributed to the delta scenario alone. The earlier material on pooled sensitivity curves gives the context for interpreting \(Q(\delta)\).

Definition: A direction boundary is a value of \(\delta\) for which \(Q(\delta)=0\). A direction tipping point is a boundary at which the estimate changes sign as \(\delta\) passes through it. A magnitude boundary, for a prespecified minimum important effect \(\kappa>0\), is a value of \(\delta\) for which \(|Q(\delta)|=\kappa\). A magnitude tipping point is a boundary where the estimate's magnitude crosses \(\kappa\). The value \(\kappa\) must be stated in the same units as \(Q\).

A zero boundary and a sign reversal are related but not identical. An estimate may touch zero and return to the same sign; that is a zero boundary without a direction reversal. Similarly, a magnitude boundary may be crossed more than once: an estimate can enter the region \(|Q(\delta)|<\kappa\) and later leave it. Therefore, report which crossing is being described, such as the first boundary reached from \(\delta=0\) in an increasing-\(\delta\) scenario.

Calculating Thresholds When the Curve Is Affine

For some estimands and adjustment schemes, the estimate changes linearly over the sensitivity range. In that case write \(Q(\delta)=q_0+b\delta\), where \(q_0=Q(0)\) and \(b\) is the change in the estimate per unit change in \(\delta\). The units of \(b\) are estimate-units per delta-unit. This affine form is a useful calculation model, not a universal property of delta-adjusted MICE analyses. If the analysis is nonlinear, calculate boundaries from its actual sensitivity curve rather than assuming a constant slope.

$$ Q(\delta)=q_0+b\delta,\qquad \delta_{\text{zero}}=-\frac{q_0}{b}\quad (b\ne 0) $$

The zero-boundary formula follows by setting \(Q(\delta)=0\): \(q_0+b\delta=0\), so \(b\delta=-q_0\), and division by \(b\ne0\) gives \(\delta=-q_0/b\). This is a relevant tipping point only if it lies in the prespecified delta range and in the direction of change being examined. If \(b=0\), the estimate is constant: it has no zero boundary unless \(q_0=0\), in which case it is zero throughout the range.

For magnitude, solve \(|Q(\delta)|=\kappa\). The absolute-value equation is equivalent to the pair of equations \(Q(\delta)=\kappa\) or \(Q(\delta)=-\kappa\). For an affine curve, their candidate solutions are:

$$ \delta_{+\kappa}=\frac{\kappa-q_0}{b}, \qquad \delta_{-\kappa}=\frac{-\kappa-q_0}{b}, \qquad (b\ne0) $$

Not every candidate is a relevant crossing: check whether it lies in the allowed range and whether \(|Q(\delta)|\) actually moves from one side of \(\kappa\) to the other. When \(Q(0)>\kappa\) and the estimate decreases toward zero, the first magnitude threshold is usually the \(+\kappa\) boundary. A later \(-\kappa\) boundary may occur after the estimate has passed through zero and grown in the opposite direction.

Worked Example: A Positive Estimate Decreases Toward Zero

Suppose a hypothetical analysis estimates a treatment contrast in symptom-score points with \(Q(\delta)=4.8-0.6\delta\), over \(0\le\delta\le12\). The positive value indicates a higher score for the first group under the stated contrast convention. Subject-matter experts have prespecified \(\kappa=1.5\) points as the minimum important magnitude. The affine relationship is assumed to describe this range for the example; in an actual analysis, it must be checked against the computed sensitivity curve.

Here \(q_0=4.8\) points and \(b=-0.6\) points per delta-unit. The zero boundary is:

$$ \delta_{\text{zero}}=-\frac{4.8}{-0.6}=8 $$

Substitution confirms \(Q(8)=4.8-0.6(8)=0\) points. The estimate is positive before \(\delta=8\) and negative after it, so this boundary is a direction tipping point.

For the magnitude boundary, solve \(Q(\delta)=1.5\) and \(Q(\delta)=-1.5\). The first gives \(4.8-0.6\delta=1.5\), so \(\delta=3.3/0.6=5.5\). The second gives \(4.8-0.6\delta=-1.5\), so \(\delta=6.3/0.6=10.5\). Thus \(|Q(\delta)|\) first falls below 1.5 points at \(\delta=5.5\), and returns to 1.5 points at \(\delta=10.5\). For example, \(Q(6)=1.2\) points, whose magnitude is below the minimum important effect; \(Q(11)=-1.8\) points, whose magnitude is above it again. The magnitude tipping point at 5.5 and the direction tipping point at 8 answer different questions.

Worked Example: A Negative Estimate Moves Toward Zero

Consider a hypothetical contrast in blood-pressure change, in millimetres of mercury, with \(Q(\delta)=-2.4+0.3\delta\) for \(0\le\delta\le12\). Suppose the minimum important magnitude is \(\kappa=1.0\) mmHg. Here the adjustment increases the estimate as \(\delta\) increases; whether that direction is plausible is a separate subject-matter judgment, not something the arithmetic establishes.

Set the estimate to zero: \(-2.4+0.3\delta=0\), giving \(\delta=2.4/0.3=8\). The estimate is negative before 8 and positive after 8, so its direction reverses there.

For magnitude, the negative-side boundary solves \(-2.4+0.3\delta=-1.0\), giving \(0.3\delta=1.4\) and \(\delta=4.6667\). The positive-side boundary solves \(-2.4+0.3\delta=1.0\), giving \(0.3\delta=3.4\) and \(\delta=11.3333\). Consequently, the estimate's magnitude first becomes less than 1.0 mmHg at about \(\delta=4.667\), then rises above 1.0 mmHg after about \(\delta=11.33\). At \(\delta=6\), \(Q(6)=-0.6\) mmHg, which is within the minimum-magnitude band; at \(\delta=9\), \(Q(9)=0.3\) mmHg, also within that band. These checks agree with the two calculated boundaries and the zero crossing at 8.

Nonlinear Curves, Uniqueness, and First Crossings

When \(Q(\delta)\) is not affine, a boundary is found by solving \(Q(\delta)=0\) or \(|Q(\delta)|=\kappa\) on the stated range. A root is a candidate boundary; the curve on either side tells whether it is a sign reversal or a magnitude crossing. Strict monotonicity is sufficient for an equation such as \(Q(\delta)=0\) to have at most one solution on an interval, but it is not necessary: a non-monotone curve can still have exactly one root. Conversely, weak monotonicity alone does not guarantee a unique solution, because the curve may sit on the boundary over an interval. Uniqueness should be determined from the actual curve or its roots, including any boundary plateau.

Worked Example: A Non-Monotone Curve Has Two Magnitude Boundaries

Suppose a hypothetical estimated effect, in score points, follows \(Q(\delta)=0.3+(\delta-2)^2\) on \(0\le\delta\le4\), and the minimum important magnitude is \(\kappa=1.0\) point. The curve is positive throughout, so it has no zero boundary and no direction reversal. To find magnitude boundaries, solve \(|Q(\delta)|=1\). Since \(Q(\delta)>0\), this reduces to \(0.3+(\delta-2)^2=1\).

$$ (\delta-2)^2=0.7,\qquad \delta=2\pm\sqrt{0.7} $$

Since \(\sqrt{0.7}\approx0.8367\), the two boundaries are \(\delta\approx1.163\) and \(\delta\approx2.837\). At \(\delta=2\), \(Q(2)=0.3\) points, below the 1.0-point minimum; at either endpoint \(\delta=0\) or 4, \(Q(\delta)=4.3\) points, above it. The curve therefore enters the below-minimum region at the first boundary and leaves it at the second. This example has two magnitude crossings despite its single minimum. It also illustrates why one should specify the first crossing when that is the intended summary.

For a direction boundary, a curve can also touch zero without changing sign. For instance, \(Q(\delta)=(\delta-1)^2\) has exactly one zero on \(0\le\delta\le2\), at \(\delta=1\), but is positive on both sides. That is a zero boundary, not a direction reversal. Likewise, if a weakly monotone curve equals \(\kappa\) over a range, every point in that range solves the magnitude-boundary equation; report the plateau rather than claiming a unique tipping point.

Interpretation, Diagnostics, and Alternatives

A direction tipping point describes how far the specified delta scenario must move for the point estimate to reach zero or reverse sign. A magnitude tipping point describes when that point estimate reaches the chosen minimum important magnitude. Neither threshold by itself states whether a confidence interval contains zero, whether a null hypothesis is rejected, or whether the underlying effect is clinically important. Those are distinct questions. In particular, statistical significance and practical importance do not share a boundary: a small effect may be precisely estimated, and an important-sized estimate may be imprecise.

Before interpreting a reported threshold, check the following:

  • Direction and units: Verify the contrast's sign convention, the units of \(Q\) and \(\kappa\), and the meaning and scale of \(\delta\). A threshold without these is difficult to interpret.
  • Range and plausibility: Confirm the solution is inside the delta range being examined and that the corresponding adjusted imputations and implied distributions remain plausible. A mathematical root outside the range is not an observed tipping point for that analysis.
  • Shape and crossings: Check estimates around each candidate boundary. A non-monotone curve can have multiple roots, and a boundary root can be a touch rather than a crossing. A long or unstable plateau can also make a single-number summary misleading.
  • Estimation stability: Check that apparent changes are not dominated by Monte Carlo variability in the imputations or by changes in the analysis specification. A threshold estimated from a noisy curve should not be reported with unjustified precision.

The main alternatives answer related but different questions. A confidence-interval tipping analysis tracks when an interval changes its relation to zero or another boundary, whereas direction and magnitude tipping points here track the point estimate itself. A clinical or policy threshold analysis evaluates a decision boundary defined for the application; a minimum important effect is a prespecified magnitude benchmark and should not be treated as interchangeable with that decision rule. Reporting the sensitivity curve, its relevant boundaries, and a plausible delta range is more informative than reporting a threshold alone.

A common mistake is to say, “The estimate is below the minimum important effect at the direction tipping point.” That need not be true: in the first worked example the estimate reaches the magnitude boundary at \(\delta=5.5\), but reaches zero only at \(\delta=8\). The correction is to solve the two boundary equations separately, inspect all relevant roots in the stated range, and describe each in context. Another mistake is to infer that a single threshold exists because the curve appears broadly increasing or decreasing. Check the actual crossings; monotonicity is not a necessary condition for uniqueness, and weak monotonicity can permit a boundary plateau.

Key takeaway: Set the boundary to zero for direction and to \(+\kappa\) or \(-\kappa\) for magnitude. State the delta direction and range, check each candidate against the sensitivity curve, and distinguish reaching a boundary from crossing it.

Check Your Understanding

Use the stated estimates as sensitivity curves over increasing \(\delta\). Each solution gives a threshold for the point estimate, not a test of statistical significance.

  1. For \(Q(\delta)=3-0.5\delta\), find the zero boundary and the first magnitude boundary for \(\kappa=1\).
  2. For \(Q(\delta)=-1.5+0.25\delta\), find the zero boundary and the first magnitude boundary for \(\kappa=0.5\).
  3. Does a curve that touches zero necessarily have a direction tipping point? Explain using \(Q(\delta)=(\delta-2)^2\).
  4. What should be reported if \(Q(\delta)=\kappa\) throughout an interval rather than at one isolated value?

Solutions. (1) Set \(3-0.5\delta=0\), so \(\delta=6\). Solve \(3-0.5\delta=1\) to obtain \(\delta=4\); this is the first magnitude boundary, before the zero boundary. (2) Set \(-1.5+0.25\delta=0\), so \(\delta=6\). The first magnitude boundary solves \(-1.5+0.25\delta=-0.5\), giving \(\delta=4\); the other boundary, at \(Q=0.5\), is \(\delta=8\). (3) No. The curve has a zero at \(\delta=2\), but is positive on both sides, so it touches zero without reversing direction. (4) Report the boundary plateau as an interval of solutions, not as a unique tipping point; weak monotonicity does not rule out such a plateau.

Clinical Biostats