Analysis of bootstrap and subsampling in high-dimensional regularized regression (code)

We investigate popular resampling methods for estimating the uncertainty of statistical models, such as subsampling, bootstrap and the jackknife, and their performance in high-dimensional supervised regression tasks. We provide a tight asymptotic description of the biases and variances estimated by these methods in the context of generalized linear models, such as ridge and logistic regression, taking the limit where the number of samples n and dimension d of the covariates grow at a comparable fixed rate α = n/d. Our findings are three-fold: i) resampling methods are fraught with problems in high dimensions and exhibit the double-descent-like behavior typical of these situations; ii) only when α is large enough do they provide consistent and reliable error estimations (we give convergence rates); iii) in the over-parametrized regime α < 1 relevant to modern machine learning practice, their predictions are not consistent, even with optimal regularization. This record provides the code to reproduce the numerical experiments of the related paper "Analysis of bootstrap and subsampling in high-dimensional regularized regression".