Order Parameter Bootstrapping in Model Selection

• Order parameter bootstrapping is a statistical framework that combines Lepski's method, multiple testing, and wild bootstrap calibration to control type I error.
• It rigorously selects the optimal bias-variance trade-off in linear estimators while achieving oracle-like risk bounds in finite samples.
• The method leverages pairwise test statistics and tailored thresholds to adaptively manage heteroscedastic noise without relying on strict parametric assumptions.

The order parameter bootstrapping argument is a statistical framework for model selection within a family of linear estimators, with a focus on controlling type I error in a completely data-driven manner that is robust to unknown and potentially heteroscedastic noise. The method, developed primarily in the context of ordered model selection for smoothing estimators, advances the "smallest accepted" approach, integrating Lepski's method with contemporary multiple testing theory and leveraging wild (multiplier) bootstrap procedures for threshold calibration. The result is a rigorous, finite-sample-valid, heteroscedasticity-robust procedure for selecting the optimal point along a sequence of bias-variance trade-offs, equipped with provable oracle‐type risk bounds and explicit error control (Spokoiny et al., 2015).

1. Framework: Ordered Linear Estimators and Model Selection

The foundational setup involves observed data $Y = \theta + \epsilon$, where $\theta \in \mathbb{R}^n$ is an unknown signal and $\epsilon$ is a zero-mean noise vector with covariance $\Sigma$, which may be heteroscedastic. Candidate estimators are constructed as $\tilde{\theta}_m = S_m Y$ for $m = 1,2,\ldots,M$, with fixed $n \times n$ smoothing matrices $S_m$. These estimators are ordered in terms of variance: $$\operatorname{Var}(\tilde{\theta}_1) \ll \operatorname{Var}(\tilde{\theta}_2) \ll \cdots \ll \operatorname{Var}(\tilde{\theta}_M)$$. Each model $m$ in the sequence offers a specific bias-variance trade-off, with lower-$\theta \in \mathbb{R}^n$0 models generally incurring higher bias and lower variance, and higher-$\theta \in \mathbb{R}^n$1 models the opposite. An acceptance rule formalizes statistical indistinguishability: model $\theta \in \mathbb{R}^n$2 is accepted if it cannot be reliably distinguished from all larger models $\theta \in \mathbb{R}^n$3, based on the observed data.

2. Pairwise Test Statistics and Acceptance Thresholds

For each pair $\theta \in \mathbb{R}^n$4 with $\theta \in \mathbb{R}^n$5, the difference $\theta \in \mathbb{R}^n$6 is considered, with covariance $\theta \in \mathbb{R}^n$7. Under the null hypothesis that model $\theta \in \mathbb{R}^n$8 is true ($\theta \in \mathbb{R}^n$9 in the column space of $\epsilon$0), this difference is zero-mean with covariance $\epsilon$1. The primary diagnostic is the quadratic form

$\epsilon$2

where $\epsilon$3 is a suitable generalized inverse, and under the null, $\epsilon$4 is approximately $\epsilon$5-distributed with $\epsilon$6 degrees of freedom. A model $\epsilon$7 is accepted if $\epsilon$8 for all $\epsilon$9, where $\Sigma$0 are critical values selected to control the family-wise error rate of false rejections via the propagation condition.

3. Wild Bootstrap Calibration for Threshold Selection

A central challenge is that $\Sigma$1 is both unknown and potentially heteroscedastic. The method addresses this through a multiplier (wild) bootstrap to approximate the joint distribution of all test statistics under each model $\Sigma$2:

• A pilot estimator $\Sigma$3 (often $\Sigma$4) defines initial residuals $\Sigma$5.
• For $\Sigma$6 bootstrap iterations, i.i.d. mean-zero multipliers $\Sigma$7 with $\Sigma$8 are generated (e.g., Rademacher or standard normal).
• Bootstrap pseudo-observations are constructed as $\Sigma$9.
• The test statistics $\tilde{\theta}_m = S_m Y$0 are recomputed for all $\tilde{\theta}_m = S_m Y$1 on the bootstrapped data.
• The empirical distribution $\tilde{\theta}_m = S_m Y$2 of $\tilde{\theta}_m = S_m Y$3 is used to set $\tilde{\theta}_m = S_m Y$4, where the levels $\tilde{\theta}_m = S_m Y$5 are chosen so that $\tilde{\theta}_m = S_m Y$6, with $\tilde{\theta}_m = S_m Y$7 the desired overall error level.

This procedure ensures that for any true model $\tilde{\theta}_m = S_m Y$8, the probability of any false rejection among $\tilde{\theta}_m = S_m Y$9 for $m = 1,2,\ldots,M$0 is bounded by $m = 1,2,\ldots,M$1.

4. Selection Rule and Type I Error Control

The final selection rule is

$m = 1,2,\ldots,M$2

Equivalently, the process begins at the lowest index and selects the smallest model $m = 1,2,\ldots,M$3 not rejected by any larger model. By the construction of $m = 1,2,\ldots,M$4 and the bootstrap-based thresholds, it follows that, for any $m = 1,2,\ldots,M$5 in the column space of some $m = 1,2,\ldots,M$6,

$m = 1,2,\ldots,M$7

where $m = 1,2,\ldots,M$8 and $m = 1,2,\ldots,M$9 are constants dependent on the Monte-Carlo and non-Gaussian tail behavior of the test statistics, respectively. This upper bound demonstrates finite-sample family-wise error control for under-smoothing.

5. Finite-Sample Validity and Oracle Inequalities

Under suitable conditions (fixed design, sub-Gaussian noise, moderate dimensionality; see assumptions A1–A3 in the original work), the method guarantees not only type I error control but also risk performance nearly matching the oracle model. Specifically, with probability at least $n \times n$0,

$n \times n$1

where $n \times n$2 is a negligible remainder. The result holds equally for full parameter estimation or for linear functionals. This demonstrates that the method is data-adaptive and achieves a nearly optimal bias-variance balance, subject only to small finite-sample corrections (Spokoiny et al., 2015).

6. Methodological Context and Significance

The order parameter bootstrapping argument builds directly on Lepski's method and the theory of multiple testing, operationalized in a model selection context through the "smallest accepted" approach. Its primary contributions are the rigorous bootstrapped calibration of critical values in the presence of heteroscedastic noise, the explicit finite-sample family-wise error control, and the explicit characterization of the selected index set. The framework is fully data-driven, requiring neither prior variance structure nor parametric assumptions on the noise. These features position it as a robust alternative to classical penalized criteria, particularly in scenarios with non-constant variance and complex dependence structures in the residuals (Spokoiny et al., 2015).