Code for the paper
Open a Julia console and run:
julia> using Pkg
julia> Pkg.add(url="https://github.com/SPOC-group/BootstrapAsymptotics")julia> using BootstrapAsymptotics
julia> (; overlaps) = state_evolution(Ridge(λ=0.1), PairBootstrap(), PairBootstrap());
julia> Matrix(overlaps.Q)
2×2 Matrix{Float64}:
1.35642 0.799808
0.799808 1.35642Open the Jupyter notebooks in the experiments folder, select the Julia environment defined in the same folder, and run.
Open the Mathematica notebook experiments/large_alpha_rates.nb
n: number of samplesd: dimensionα: sampling ration/dΔ: noise varianceλ: regularization parameterm,Q,V: overlaps