Tutorials › Multi-Arm Multi-Stage (MAMS) Designs › Limits of the Jackknife for Selected-Arm Estimates

Treatment-Effect Estimation After Selection · Tutorial 279 of 499

Limits of the Jackknife for Selected-Arm Estimates

Learn how selection boundaries disrupt ordinary jackknife bias estimation, how to diagnose that instability, and what its results can and cannot tell you.

Advanced 11 min read

What You'll Learn

  • Define the ordinary leave-one-out jackknife bias estimate and corrected estimate.
  • Explain how arm switching at a selection boundary disrupts jackknife approximations.
  • Distinguish an unstable estimate from a jackknife replicate that is not defined.
  • Calculate jackknife quantities in examples with arm selection and a fixed selected arm.
  • Diagnose instability using selection changes, boundary margins, and jackknife replicates.
  • Compare the ordinary jackknife with validation-based and selection-aware estimation.

Why the Ordinary Jackknife Can Mislead After Selection

In a multi-arm, multi-stage trial, the reported treatment estimate can depend on a sequence of decisions: an arm may be retained or dropped at an interim analysis, and one or more arms may be selected for reporting. The ordinary leave-one-out jackknife estimates bias by repeatedly removing one observation and recalculating the statistic. That procedure is most informative when removing one observation makes only a small, smooth change to the statistic. A selection rule can violate that condition: a small data change may switch the selected arm or change whether an arm is reported at all.

The issue is not that a jackknife calculation is always impossible after selection. It is that its usual bias interpretation can be unreliable when leave-one-out samples cross selection boundaries. Even when every replicate produces a number, some replicates may estimate different arms’ effects. Their average then need not describe bias for one fixed treatment effect.

Definition: Let \(X_1,\ldots,X_n\) be the observational units used by a statistic \(T\), and let \(\hat\theta=T(X_1,\ldots,X_n)\). Define \(\hat\theta_{(-i)}\) as the statistic recalculated after deleting unit \(i\), including rerunning any selection and dropping rules that form part of \(T\). The ordinary jackknife bias estimate and bias-corrected estimate are
$$ \widehat{\operatorname{Bias}}_J=(n-1)\left(\bar{\theta}_{(-)}-\hat\theta\right), \qquad \hat\theta_J=\hat\theta-\widehat{\operatorname{Bias}}_J =n\hat\theta-(n-1)\bar{\theta}_{(-)}, \qquad \bar{\theta}_{(-)}=\frac{1}{n}\sum_{i=1}^{n}\hat\theta_{(-i)}. $$
These formulas require every \(\hat\theta_{(-i)}\) to be defined on a compatible scale and to target the same estimand. They do not, by themselves, make a nonsmooth selected-arm statistic suitable for jackknife bias correction.

For a MAMS application, the relevant statistic must be specified before computing these quantities. If it is the treatment-effect estimate for the arm selected by a prespecified rule, each deletion replicate should rerun that rule. Holding the observed selection fixed instead answers a different, conditional question: how the estimate changes when an observation is removed, given the selected arm. Neither calculation repairs an unstable selection boundary; they simply describe different procedures.

Why Smoothness Matters

A short derivation shows why the max of two estimates is problematic near a tie. Suppose two arms have independent interim estimates \(A\) and \(B\), each distributed as \(N(\mu,\sigma^2/m)\), where \(m\) is the number of observations contributing to each estimate. The reported estimate for the arm with the larger interim estimate is \(\max(A,B)\). The identity

$$ \max(A,B)=\frac{A+B+|A-B|}{2} $$

gives the expected selected estimate directly. Since \(A-B\sim N(0,2\sigma^2/m)\), and a mean-zero normal variable with variance \(v\) has expected absolute value \(\sqrt{2v/\pi}\),

$$ E\{\max(A,B)\} =\mu+\frac{1}{2}E|A-B| =\mu+\frac{\sigma}{\sqrt{\pi m}}. $$

Thus, when the true arm effects are equal, selecting the larger noisy estimate creates a positive bias of order \(m^{-1/2}\). For many smooth statistics, the ordinary jackknife is motivated by small, gradual changes in the statistic under one-observation perturbations; it is not a general correction for a boundary effect of this kind. At a tie, the function \(\max(A,B)\) has a kink: an arbitrarily small change can change which arm determines the reported estimate. A drop rule has a related boundary, where a small change can change the set of arms still under study.

This calculation describes a simple selection rule, not every MAMS design. Multi-stage designs can have several boundaries, and their estimates may use information from different stages. The earlier material on bootstrap failure at selection boundaries considers the same broad difficulty for a different resampling method; here the point is specifically that leave-one-out perturbations can also cross a boundary and make an ordinary jackknife bias estimate unstable.

Worked Example: A Leave-One-Out Deletion Switches the Selected Arm

Consider a hypothetical interim comparison with four independent, equally weighted effect contributions for each of two active arms. The contributions, in treatment-effect units, are A: \(6,6,6,2\) and B: \(5,5,5,5\). The prespecified rule selects the arm with the larger mean, breaking an exact tie in favor of A. Calculate the ordinary jackknife bias estimate for the selected mean, rerunning the selection after every deletion.

The full-sample means are \(\bar X_A=20/4=5\) and \(\bar X_B=20/4=5\). The tie rule selects A, so \(\hat\theta=5\). There are \(n=8\) contributions.

Deleting any one of A’s three contributions equal to 6 leaves A’s mean at \(14/3=4.667\), so B is selected and the reported mean is 5. Deleting A’s contribution equal to 2 leaves A’s mean at \(18/3=6\), so A is selected and the reported mean is 6. Deleting any one of B’s four contributions leaves B’s mean at \(15/3=5\); the tie rule selects A, and the reported mean remains 5.

Thus seven leave-one-out reported means are 5 and one is 6:

$$ \bar{\theta}_{(-)}=\frac{7(5)+6}{8}=\frac{41}{8}=5.125. $$

The jackknife bias estimate is \((8-1)(5.125-5)=0.875\) treatment-effect units, giving \(\hat\theta_J=5-0.875=4.125\). The full estimate selected A, but three deletions select B. The numeric correction is therefore based on replicates that do not all estimate the same arm’s effect. Its value should not be read as a dependable estimate of the bias for A. The example also illustrates how one unusual replicate can materially affect the jackknife average.

When a Jackknife Replicate Is Not Defined

A separate problem occurs when a deletion changes a retain/drop decision and the selected-arm estimator is defined only for arms that remain in the trial. In that case, a replicate may not be a noisy version of the original estimate: it may have no estimate at all under the stated reporting rule. Assigning an arbitrary number to that missing replicate does not solve the problem; it changes the statistic being jackknifed.

Worked Example: A Deletion Crosses a Drop Boundary

Suppose a hypothetical arm is retained when its interim mean is at least 4.05 treatment-effect units. Its four interim contributions are \(3.5, 4.1, 4.4,\) and \(4.4\), so its full-sample mean is \(16.4/4=4.1\). The arm is retained, and its reported estimate is 4.1.

The four leave-one-out means are

$$ \frac{12.9}{3}=4.3,\qquad \frac{12.3}{3}=4.1,\qquad \frac{12.0}{3}=4.0,\qquad \frac{12.0}{3}=4.0. $$

Deleting either contribution equal to 4.4 gives a mean of 4.0, below the 4.05 retention boundary. Those two deletion replicates drop the arm. If the statistic is defined only as the reported estimate for a retained arm, these replicates are undefined. The jackknife average \(\bar\theta_{(-)}\) and bias estimate therefore cannot be computed for this statistic without an additional, scientifically meaningful definition for dropped-arm outcomes. Coding each dropped replicate as zero, for example, would impose a new reporting rule rather than estimate bias under the original one.

When Deletions Are Informative—and When They Are Not

The jackknife is not automatically useless after a selection decision. If arm selection is made using interim data and a separate, independent validation sample is used to estimate the chosen arm’s effect, then deleting observations only from that validation sample leaves the selection fixed. The resulting jackknife describes variation in the validation estimator conditional on the selected arm. It does not estimate the full procedure’s selection bias, because the selection data and decision were not perturbed.

Worked Example: Jackknifing a Validation Mean With Selection Frozen

In a hypothetical trial, interim data select arm A. An independent validation sample then supplies four effect contributions for A: \(1,2,3,\) and \(4\) points. Assume these observations are independent and that the estimate of A’s effect is their arithmetic mean. Selection is based only on the separate interim data and is held fixed for this calculation.

The validation estimate is \(\hat\theta=(1+2+3+4)/4=2.5\) points. Deleting each observation in turn gives leave-one-out means \(3,\;8/3,\;7/3,\) and \(2\) points, whose average is

$$ \bar{\theta}_{(-)} =\frac{3+8/3+7/3+2}{4} =2.5\ \text{points}. $$

The jackknife bias estimate is \((4-1)(2.5-2.5)=0\) points, and the jackknife-corrected estimate is also 2.5 points. This is expected for the arithmetic mean: its leave-one-out average equals the full-sample mean. It says nothing about bias from choosing A using interim data, because the validation observations are independent of that choice and the jackknife did not redo selection.

Diagnostics, Alternatives, and Limitations

Before interpreting an ordinary jackknife correction, check what the deletion replicates actually do. Record the selected arm or retained set for each replicate, not only the resulting estimate. If many deletions switch arms, or if an estimate changes sharply when the selected arm changes, the average may mix distinct effects. Also record each replicate’s distance from relevant selection boundaries. A full-sample result just beyond a boundary, combined with leave-one-out results on both sides, is direct evidence that the result is locally sensitive to individual observations.

Jackknife pseudo-values can expose extreme contributions:

$$ P_i=n\hat\theta-(n-1)\hat\theta_{(-i)}, \qquad \frac{1}{n}\sum_{i=1}^{n}P_i=\hat\theta_J. $$

Large or negative pseudo-values deserve investigation, especially when they coincide with changes in the selected arm. They do not establish that the jackknife is invalid in every setting, but they help reveal how strongly the correction depends on particular deletions. If deletion can make the estimator undefined, that is a more basic failure: the stated statistic does not supply all the replicates required by the ordinary formula.

The main alternatives answer different questions. A separate validation sample can estimate an effect conditional on the arm chosen by independent interim data. Sample splitting likewise separates selection information from estimation information, at the cost of using less data for each task. Selection-aware estimation instead models or conditions on the selection rule. These approaches require their own design and modeling assumptions; they are not interchangeable corrections to a jackknife number. The ordinary jackknife may still be useful as a local sensitivity check when selection is stable and all replicates target the same estimand, but a small correction does not prove that selection bias is absent.

Key takeaway: Ordinary leave-one-out jackknife bias estimation relies on meaningful, comparable deletion replicates. When deleting one observation changes which arm is selected, crosses a drop boundary, or changes the estimand, the resulting average may be unstable, misleading, or undefined. Diagnose the selection paths before interpreting the correction.

Check Your Understanding

Answer each question, then compare your reasoning with the brief solutions.

  1. For \(n=10\), a statistic is \(\hat\theta=3.2\) and the average leave-one-out value is 3.3. What are the jackknife bias estimate and corrected estimate?
  2. Why can the maximum of two arm estimates have a larger selection bias near equal true effects than a smooth-statistic jackknife approximation is designed to address?
  3. A deletion causes the selected arm to switch, but the reported numeric estimate changes very little. Does that alone establish that the jackknife estimates bias for the originally selected arm?
  4. What should be concluded if two leave-one-out replicates drop an arm and the estimator is defined only for retained arms?
  5. When can jackknifing a validation-sample mean be informative after interim selection, and what does it not measure?

Solutions. (1) The bias estimate is \(9(3.3-3.2)=0.9\), and the corrected estimate is \(3.2-0.9=2.3\), in the statistic’s units. (2) Near a tie, small changes can switch the selected arm; for two independent normal estimates with equal means, the expected maximum exceeds the common mean by \(\sigma/\sqrt{\pi m}\). (3) No. The replicates may refer to different arm effects even if their numeric values are close. (4) The ordinary jackknife average is undefined for that estimator unless a defensible value is specified for the dropped-arm outcome; an arbitrary code changes the statistic. (5) It can describe the validation estimator conditional on the interim-selected arm when validation data are independent and selection is held fixed. It does not measure bias caused by the selection process.

Clinical Biostats