REVIEW 2 major objections 1 minor 57 references
Studentized Cheap Bootstrap: Achieving Higher-Order Coverage Accuracy with Low Computation
T0 review · 2 major / 1 minor · reviewed 2026-06-25 · grok-4.3
Pith's one-line read The studentized cheap bootstrap matches conventional higher-order accuracy using only a few Monte Carlo replications.
desk verdict The studentized cheap bootstrap claims to match full studentization's higher-order coverage with just a few inner MC reps by formalizing a t-distribution link where df equals the MC effort, but the Edgeworth expansions that carry the claim are not shown. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The studentized cheap bootstrap based on limiting t-distributions with degrees of freedom tied to Monte Carlo replications, supported by Edgeworth and Cornish-Fisher expansions.
What would settle it
A numerical experiment where the empirical coverage error of the studentized cheap bootstrap with few replications fails to match the higher-order reduction predicted by the t-distribution expansions.
Extended reading notes
Core claim
By formalizing the link between the bootstrap and the t-distribution where the degrees of freedom equal the Monte Carlo computation effort, the studentized cheap bootstrap achieves the same higher-order coverage error reduction as the conventional studentized bootstrap but with substantially thinned resampling effort.
Load-bearing premise
The higher-order coverage accuracy relies on the explicit calculations and geometric analyses of higher-order terms in the Edgeworth and Cornish-Fisher expansions for limiting t-distributions.
Editorial extensions
If this is right
- The higher-order coverage accuracy is attained with minimal Monte Carlo replications in the additional layer.
- The degrees of freedom in the t-distribution correspond to the Monte Carlo effort rather than data size.
- Explicit calculations of higher-order terms in expansions for t-distributions underpin the accuracy.
- Geometric analyses of these terms confirm the coverage properties.
Reading between the lines
- This method could make higher-order bootstrap inference practical for problems where repeated resampling is costly.
- The insight on degrees of freedom might suggest similar formal links in other statistical procedures involving approximations.
- Further work could explore the finite-sample behavior when the number of replications is small but fixed.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes the studentized cheap bootstrap, which achieves the same higher-order coverage accuracy as conventional studentized bootstrap while reducing the inner resampling layer to only a few Monte Carlo replications. It claims a formal (rather than informal) link between studentization and the t-distribution, with the degrees of freedom equal to the Monte Carlo effort rather than sample size, and supports the coverage claim via explicit calculations and geometric analyses of higher-order terms in Edgeworth and Cornish-Fisher expansions tailored to the limiting t-distribution.
Significance. If the derivations are correct, the result would meaningfully advance bootstrap methodology by making higher-order accurate inference computationally tractable at scale. The reframing of the t-distribution degrees of freedom as tied to Monte Carlo effort rather than n is a potentially useful theoretical insight for resampling-based inference.
major comments (2)
- [theoretical development of Edgeworth/Cornish-Fisher expansions] The central higher-order accuracy claim (abstract) rests on the correctness of the Edgeworth and Cornish-Fisher expansions for the studentized cheap bootstrap statistic under a limiting t-distribution whose degrees of freedom equal the number of inner Monte Carlo replications. The manuscript must exhibit the full expansions (including any cross terms arising from reuse of the same inner replicates) so that readers can verify that the claimed O(n^{-3/2}) error reduction is not undermined by omitted dependence.
- [section presenting the geometric analysis] The geometric analysis of higher-order terms (abstract) is invoked to establish equivalence with conventional studentization; without the explicit geometric arguments or the resulting coverage-error bounds, it is impossible to confirm that the cheap version inherits the same accuracy order.
minor comments (1)
- [Introduction] Clarify in the introduction whether the few inner replications are drawn independently of the outer bootstrap or share randomness, as this affects the dependence structure in the expansions.
Simulated Author's Rebuttal
Thank you for the careful reading and constructive comments. We address each major point below and will revise the manuscript to improve transparency.
read point-by-point responses
-
Referee: [theoretical development of Edgeworth/Cornish-Fisher expansions] The central higher-order accuracy claim (abstract) rests on the correctness of the Edgeworth and Cornish-Fisher expansions for the studentized cheap bootstrap statistic under a limiting t-distribution whose degrees of freedom equal the number of inner Monte Carlo replications. The manuscript must exhibit the full expansions (including any cross terms arising from reuse of the same inner replicates) so that readers can verify that the claimed O(n^{-3/2}) error reduction is not undermined by omitted dependence.
Authors: We agree that the full expansions, including all cross terms from shared inner replicates, must be shown explicitly. The current manuscript derives the leading terms under the t-limit but omits the complete dependence structure for brevity. In revision we will add the full Edgeworth and Cornish-Fisher expansions with every cross term derived, confirming that the O(n^{-3/2}) coverage error is preserved. revision: yes
-
Referee: [section presenting the geometric analysis] The geometric analysis of higher-order terms (abstract) is invoked to establish equivalence with conventional studentization; without the explicit geometric arguments or the resulting coverage-error bounds, it is impossible to confirm that the cheap version inherits the same accuracy order.
Authors: The geometric analysis appears in the higher-order terms section and establishes the equivalence via the limiting t-distribution. To address the request for explicit arguments and bounds, we will expand this section with the full geometric derivations and the resulting coverage-error order statements. revision: yes
Circularity Check
No significant circularity; higher-order claims rest on explicit Edgeworth expansions
full rationale
The paper derives its higher-order coverage accuracy from explicit calculations and geometric analyses of Edgeworth and Cornish-Fisher expansions for limiting t-distributions (with df tied to Monte Carlo replications). No load-bearing steps reduce by definition, by fitting a parameter then relabeling it a prediction, or by self-citation chains. The formal link between studentization and t-limits is presented as a derived insight rather than an input assumption. The derivation is self-contained against standard asymptotic theory.
Assumptions & free parameters
assumptions (1)
- domain assumption Applicability of Edgeworth and Cornish-Fisher expansions to the bootstrap distribution under t-distribution limits
Cite this review
Pith. "Pith review of Studentized Cheap Bootstrap: Achieving Higher-Order Coverage Accuracy with Low Computation." pith.science (2026). https://pith.science/paper/QHOSIFBB
@misc{pith2026260625968,
author = {Pith},
title = {Pith review of: Studentized Cheap Bootstrap: Achieving Higher-Order Coverage Accuracy with Low Computation},
year = {2026},
howpublished = {\url{https://pith.science/paper/QHOSIFBB}},
note = {Machine review of arXiv:2606.25968}
}
read the original abstract
The bootstrap is a versatile method for quantifying statistical uncertainty. Among its variants, a popular approach, the studentized bootstrap, provably achieves higher-order coverage error reduction compared to other benchmarks. However, its implementation typically requires an analytical form of the standard error, or otherwise an additional layer of resampling effort which can be computationally expensive. In this paper, we introduce what we call the studentized cheap bootstrap that achieves the same higher-order coverage accuracy as the conventional studentization, but substantially thinning the computational effort in the additional resampling layer to only very few Monte Carlo replications. Intriguingly, while conventional wisdom views "studentization" as an informal link between the bootstrap and t-distribution, we provide a first recognition that this link is in fact formal, notably with a distinct insight that the degree of freedom in the t-distribution corresponds to the Monte Carlo computation effort in the additional resampling layer, rather than the data size as in traditional thinking. Moreover, our desirable higher-order coverage accuracy builds crucially on this insight, as well as explicit calculations and geometric analyses of higher-order terms in the Edgeworth and Cornish-Fisher expansions tailored to limiting t-distributions.
Reference graph
Works this paper leans on
-
[1]
and SINGH, K
ABRAMOVITCH, L. and SINGH, K. (1985). Edgeworth corrected pivotal statistics and the bootstrap.The Annals of Statistics13116–132
1985
-
[2]
BABU, G. J. and SINGH, K. (1983). Inference on means using the bootstrap.The Annals of Statistics11 999–1003
1983
-
[3]
BARNDORFF-NIELSEN, O. E. (1983). On a formula for the distribution of the maximum likelihood estima- tor.Biometrika70343–365
1983
-
[4]
BARNDORFF-NIELSEN, O. E. (1991). Modified signed log likelihood ratio.Biometrika78557–563
1991
-
[5]
BARNDORFF-NIELSEN, O. E. and COX, D. R. (1984). Bartlett adjustments to the likelihood ratio statistic and the distribution of the maximum likelihood estimator.Journal of the Royal Statistical Society. Series B (Methodological)46483–495
1984
-
[6]
BARNDORFF-NIELSEN, O. E. and COX, D. R. (1994).Inference and Asymptotics. Chapman & Hall
1994
-
[7]
BARTLETT, M. S. (1937). Properties of sufficiency and statistical tests.Proceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences160268–282
1937
-
[8]
BERAN, R. (1987). Prepivoting to reduce level error of confidence sets.Biometrika74457–468
1987
Show all 57 references
-
[9]
BHATTACHARYA, R. N. and GHOSH, J. K. (1978). On the validity of the formal Edgeworth expansion.The Annals of Statistics6434–451
1978
-
[10]
BHATTACHARYA, R. N. and RAO, R. R. (2010).Normal Approximation and Asymptotic Expansions. SIAM
2010
-
[11]
J., GÖTZE, F
BICKEL, P. J., GÖTZE, F. andVANZWET, W. R. (1997). Resampling fewer thannobservations: Gains, losses, and remedies for losses.Statistica Sinica71–31
1997
-
[12]
BOOTH, J. G. and HALL, P. (1994). Monte Carlo approximation and the iterated bootstrap.Biometrika81 331–340
1994
-
[13]
R., DAVISON, A
BRAZZALE, A. R., DAVISON, A. C. and REID, N. (2007).Applied Asymptotics: Case Studies in Small- Sample Statistics. Cambridge University Press. 74
2007
-
[14]
CHOWDHURY, A. R. and LAM, H. (2025). Efficient uncertainty quantification of bagging via the cheap bootstrap. InProceedings of the 2025 Winter Simulation Conference (WSC)3442–3453. IEEE
2025
-
[15]
COX, D. R. and REID, N. (1987). Parameter orthogonality and approximate conditional inference.Journal of the Royal Statistical Society. Series B (Methodological)491–39
1987
-
[16]
DANIELS, H. E. (1954). Saddlepoint approximations in statistics.The Annals of Mathematical Statistics25 631–650
1954
-
[17]
DAVISON, A. C. and HINKLEY, D. V. (1997).Bootstrap Methods and Their Application. Cambridge Uni- versity Press
1997
-
[18]
and EFRON, B
DICICCIO, T. and EFRON, B. (1992). More accurate confidence intervals in exponential families. Biometrika79231–245
1992
-
[19]
and ROMANO, J
DICICCIO, T., HALL, P. and ROMANO, J. (1991). Empirical likelihood is Bartlett-correctable.The Annals of Statistics191053–1061
1991
-
[20]
DICICCIO, T. J. and EFRON, B. (1996). Bootstrap confidence intervals.Statistical Science11189–228
1996
-
[21]
EFRON, B. (1981). Nonparametric estimates of standard error: The jackknife, the bootstrap and other meth- ods.Biometrika68589–599
1981
-
[22]
EFRON, B. (1987). Better bootstrap confidence intervals.Journal of the American Statistical Association82 171–185
1987
-
[23]
EFRON, B. (2003). Second thoughts on the bootstrap.Statistical Science18135–140
2003
-
[24]
and NARASIMHAN, B
EFRON, B. and NARASIMHAN, B. (2020). The automatic construction of bootstrap confidence intervals. Journal of Computational and Graphical Statistics29608–619
2020
-
[25]
and TIBSHIRANI, R
EFRON, B. and TIBSHIRANI, R. J. (1994).An Introduction to the Bootstrap. CRC Press
1994
-
[26]
and BRODERICK, T
GIORDANO, R., STEPHENSON, W., LIU, R., JORDAN, M. and BRODERICK, T. (2019). A Swiss army infinitesimal jackknife. InProceedings of the 22nd International Conference on Artificial Intelligence and Statistics1139–1147. PMLR
2019
-
[27]
HALL, P. (1983). Inverting an Edgeworth expansion.The Annals of Statistics11569–576
1983
-
[28]
HALL, P. (1986). On the bootstrap and confidence intervals.The Annals of Statistics141431–1452
1986
-
[29]
HALL, P. (1988). Theoretical comparison of bootstrap confidence intervals.The Annals of Statistics16 927–953
1988
-
[30]
HALL, P. (1988). On symmetric bootstrap confidence intervals.Journal of the Royal Statistical Society. Series B (Methodological)5035–45
1988
-
[31]
(2013).The Bootstrap and Edgeworth Expansion
HALL, P. (2013).The Bootstrap and Edgeworth Expansion. Springer
2013
-
[32]
and MARTIN, M
HALL, P. and MARTIN, M. A. (1988). On bootstrap resampling and iteration.Biometrika75661–671
1988
-
[33]
and LAM, H
HE, S. and LAM, H. (2024). Higher-order coverage errors of batching methods via Edgeworth expansions on t-statistics.The Annals of Statistics521360–1383
2024
-
[34]
and ZHANG, H
HUANG, Z., LAM, H. and ZHANG, H. (2023). Efficient uncertainty quantification and reduction for over- parameterized neural networks.Advances in Neural Information Processing Systems3664428–64467
2023
-
[35]
JENSEN, J. L. (1995).Saddlepoint Approximations. Oxford University Press, Oxford
1995
-
[36]
and JORDAN, M
KLEINER, A., TALWALKAR, A., SARKAR, P. and JORDAN, M. I. (2014). A scalable bootstrap for massive data.Journal of the Royal Statistical Society: Series B (Statistical Methodology)76795–816
2014
-
[37]
LAM, H. (2022). A cheap bootstrap method for fast inference. arXiv:2202.00090
2022
-
[38]
LAM, H. (2022). Cheap bootstrap for input uncertainty quantification. InProceedings of the 2022 Winter Simulation Conference (WSC)2318–2329. IEEE
2022
-
[39]
and LIU, Z
LAM, H. and LIU, Z. (2023). Bootstrap in high dimension with low computation. InProceedings of the 40th International Conference on Machine Learning18419–18453. PMLR
2023
-
[40]
and WANG, Z
LAM, H. and WANG, Z. (2026). Cheap bootstrap for fast uncertainty quantification of stochastic gradient descent.Journal of Machine Learning Research. To appear
2026
-
[41]
LAWLEY, D. N. (1956). A general method for approximating to the distribution of likelihood ratio criteria. Biometrika43295–303
1956
-
[42]
LEE, S. M. S. and YOUNG, G. A. (1995). Asymptotic iterated bootstrap confidence intervals.The Annals of Statistics231301–1330
1995
-
[43]
LOPES, M. E. (2019). Estimating the algorithmic variance of randomized ensembles via the bootstrap.The Annals of Statistics471088–1112
2019
-
[44]
and SHA, F
LU, Z., IE, E. and SHA, F. (2020). Uncertainty estimation with infinitesimal jackknife, its distribution and mean-field approximation. arXiv:2006.07584
2020
-
[45]
and RICE, S
LUGANNANI, R. and RICE, S. (1980). Saddlepoint approximation for the distribution of the sum of inde- pendent random variables.Advances in Applied Probability12475–490
1980
-
[46]
S., MUNCH, A., SØRENSEN, K
OHLENDORFF, J. S., MUNCH, A., SØRENSEN, K. K. and GERDS, T. A. (2025). Cheap subsampling bootstrap confidence intervals for fast and robust inference. arXiv:2501.10289
2025
-
[47]
OWEN, A. B. (1988). Empirical likelihood ratio confidence intervals for a single functional.Biometrika75 237–249. STUDENTIZED CHEAP BOOTSTRAP75
1988
-
[48]
OWEN, A. B. (1990). Empirical likelihood ratio confidence regions.The Annals of Statistics1890–120
1990
-
[49]
N., ROMANO, J
POLITIS, D. N., ROMANO, J. P. and WOLF, M. (1999).Subsampling. Springer
1999
-
[50]
QUMSIYEH, M. B. (1990). Edgeworth expansion in regression models.Journal of Multivariate Analysis35 86–101
1990
-
[51]
REID, N. (2003). Asymptotics and the theory of inference.The Annals of Statistics311695–1731
2003
-
[52]
and SARIA, S
SCHULAM, P. and SARIA, S. (2019). Can you trust this prediction? Auditing pointwise reliability after learning. InProceedings of the 22nd International Conference on Artificial Intelligence and Statistics 1022–1031. PMLR
2019
-
[53]
and SHAO, X
SENGUPTA, S., VOLGUSHEV, S. and SHAO, X. (2016). A subsampled double bootstrap for massive data. Journal of the American Statistical Association1111222–1232
2016
-
[54]
SEVERINI, T. A. (2000).Likelihood Methods in Statistics. Oxford University Press, Oxford
2000
-
[55]
and TU, D
SHAO, J. and TU, D. (2012).The Jackknife and Bootstrap. Springer. [56]VAN DERVAART, A. W. and WELLNER, J. A. (1996).Weak Convergence and Empirical Processes: With Applications to Statistics. Springer
2012
-
[56]
and EFRON, B
WAGER, S., HASTIE, T. and EFRON, B. (2014). Confidence intervals for random forests: The jackknife and the infinitesimal jackknife.Journal of Machine Learning Research151625–1651
2014
-
[57]
WITHERS, C. S. (1984). Asymptotic expansions for distributions and quantiles with power series cumulants. Journal of the Royal Statistical Society. Series B (Methodological)46389–396
1984
Reviewed June 25, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.