Sam Babichenko

Reading notes / High-dimensional SGD

The information exponent

Ben Arous et al., 2021; Damian et al., 2023

Context. Online SGD estimating a unit vector from a stream of samples starts from a random point, whose correlation with the truth is about N−1/2N^{-1/2}: on the equator, where per-sample gradient noise dwarfs the population gradient. The ODE method assumes a warm start and scaling limits fix m0>0m_0\gt0, so it was open how many samples SGD needs just to leave the equator, and what sets that power of NN. Ben Arous et al. (2021) answer for online SGD, and Damian et al. (2023) show how to do better. Only the population loss near the equator matters. If ϕ(m)\phi(m) has its first nonzero derivative at m=0m=0 in order kk (the information exponent), the drift at the start is mk−1≈N−(k−1)/2m^{k-1}\approx N^{-(k-1)/2} against order-one noise per step, and escaping costs about Nk−1N^{k-1} samples for k≥3k\ge3, Nlog⁡NN\log N to Nlog⁡2NN\log^2N for k=2k=2, and NN for k=1k=1; once the correlation is macroscopic the descent is an ODE and takes a vanishing fraction of the data. Smoothing the loss over a cap of radius about λ\lambda around the iterate picks up signal from the far side of the cap and reduces the noise faster than the signal, which lowers the cost to dk/2d^{k/2} at λ=d1/4\lambda=d^{1/4}, the correlational-query lower bound.

Setting. Samples are i.i.d. from PθP_\theta, M=αNM=\alpha N of them, with population loss Φ(x)=ϕ(m(x))\Phi(x)=\phi(m(x)), m=x⋅θm=x\cdot\theta. Spherical online SGD takes X~t=Xt−1−δN∇L(Xt−1;Yt)\tilde X_t=X_{t-1}-\tfrac{\delta}{N}\nabla L(X_{t-1};Y^t) and renormalizes, one fresh sample per step. Assumption A: ϕ′<0\phi'\lt0 on (0,1)(0,1). Assumption B: the sample-wise error H=L−ΦH=L-\Phi has gradient noise of order one along θ\theta and order N\sqrt N in norm (4+ι4+\iota moments), as for an isotropic subgaussian vector. For the single-index model y=f(a⋅θ)+εy=f(a\cdot\theta)+\varepsilon with Gaussian aa and square loss, the noise operator is diagonal on Hermite polynomials, so ϕ(m)=2∑j≥0uj(f)2(1−mj)+Eε2,uj(f)=E[f(Z)hj(Z)],\phi(m)=2\sum_{j\ge0}u_j(f)^2(1-m^j)+\mathbb E\varepsilon^2,\qquad u_j(f)=\mathbb E[f(Z)h_j(Z)], and kk is the index of the first nonzero Hermite coefficient with j≥1j\ge1. Linear, ReLU, sigmoid, x3x^3 and canonical-link GLMs have k=1k=1; phase retrieval, online PCA and the two-Gaussian mixture have k=2k=2; the spiked pp-tensor has k=pk=p; x3−3x=He3x^3-3x=\mathrm{He}_3 has k=3k=3. With a misspecified link the sign of uk(f)uk(g)u_k(f)u_k(g) decides whether SGD moves toward θ\theta at all. Damian et al. write dd for NN, α=w⋅w⋆\alpha=w\cdot w^\star for mm, k⋆k^\star for k≥3k\ge3, a known link σ=∑kckHek/k!\sigma=\sum_kc_k\mathrm{He}_k/k!, and the correlation loss 1−y σ(w⋅x)1-y\,\sigma(w\cdot x).

Ben Arous et al.'s theorems. Let αc(N,k)=1, log⁡N, Nk−2\alpha_c(N,k)=1,\ \log N,\ N^{k-2} for k=1,2,≥3k=1,2,\ge3. Recovery (Theorem 1.3): under A and B, with α−1≪δ≪α−1/2\alpha^{-1}\ll\delta\ll\alpha^{-1/2} and α≫1\alpha\gg1, log⁡2N\log^2N, Nk−2log⁡2NN^{k-2}\log^2N respectively, SGD from a uniform start on the upper hemisphere has m(XM)→1m(X_M)\to1 in probability; in samples that is NN, Nlog⁡2NN\log^2N and Nk−1log⁡2NN^{k-1}\log^2N. Refutation (Theorem 1.4): if α≪αc\alpha\ll\alpha_c and δ=O(1)\delta=O(1) (k≤2k\le2) or δ=O(α−1/2)\delta=O(\alpha^{-1/2}) (k≥3k\ge3), then sup⁡t≤M∣m(Xt)∣→0\sup_{t\le M}|m(X_t)|\to0, so the thresholds are sharp up to log⁡2N\log^2N; this is about online SGD in these step ranges only, since information-theoretically order NN samples suffice. Descent (Theorems 3.2 and 1.5): from m(X0)=η>0m(X_0)=\eta\gt0 with αδ2→0\alpha\delta^2\to0, SGD tracks projected gradient descent on Φ\Phi, kk plays no role, and the climb to 1−η1-\eta takes order N/δN/\delta samples, small next to the search.

Damian et al.'s theorems. Smoothing averages over the tangent sphere, (Lλg)(w)=Ez∼μw g((w+λz)/∥w+λz∥)(\mathcal L_\lambda g)(w)=\mathbb E_{z\sim\mu_w}\,g\big((w+\lambda z)/\|w+\lambda z\|\big) with μw\mu_w uniform on {z⊥w}\{z\perp w\}, and the algorithm is online spherical SGD on LλL(⋅;xt,yt)\mathcal L_\lambda L(\cdot;x_t,y_t). Theorem 1: from w0⋅w⋆≳d−1/2w_0\cdot w^\star\gtrsim d^{-1/2}, with λ∈[1,d1/4]\lambda\in[1,d^{1/4}] and step O~(d−k⋆/2λ2k⋆−2)\tilde O(d^{-k^\star/2}\lambda^{2k^\star-2}) for T1=O~(dk⋆−1λ−2k⋆+4)T_1=\tilde O(d^{k^\star-1}\lambda^{-2k^\star+4}) steps, then λ=0\lambda=0 and ηt=O((d+t−T1)−1)\eta_t=O((d+t-T_1)^{-1}), the last iterate has loss O(d/(d+T2))O(d/(d+T_2)) with high probability, so n ≳ dk⋆−1λ−2k⋆+4+dε,so at λ=d1/4:n ≳ dk⋆/2+dεn\ \gtrsim\ d^{k^\star-1}\lambda^{-2k^\star+4}+\frac d\varepsilon,\qquad\text{so at }\lambda=d^{1/4}:\quad n\ \gtrsim\ d^{k^\star/2}+\frac d\varepsilon up to polylogs; λ=O(1)\lambda=O(1) recovers the Ben Arous et al. rate. Theorem 2: any correlational statistical query (CSQ) algorithm making qq queries of tolerance τ\tau that reaches loss below 1/21/2 needs τ≲(log⁡(qd)/d)k⋆/4\tau\lesssim(\log(qd)/d)^{k^\star/4}, which with τ≈n−1/2\tau\approx n^{-1/2} means n≳dk⋆/2n\gtrsim d^{k^\star/2}.

Step 1

The escape time is set by a drift budget and a noise budget. Track only mtm_t. To first order, mt≈mt−1+δN akmt−1k−1−δN∇Ht(Xt−1)⋅θ,m_t\approx m_{t-1}+\tfrac{\delta}{N}\,a_k m_{t-1}^{k-1}-\tfrac{\delta}{N}\nabla H^t(X_{t-1})\cdot\theta, a drift from the Taylor expansion of ϕ′\phi' plus the increment of the directional error martingale. In ODE time s=δt/Ns=\delta t/N the drift m˙=akmk−1\dot m=a_km^{k-1} from m0≈N−1/2m_0\approx N^{-1/2} needs time 11, log⁡N\log N and N(k−2)/2N^{(k-2)/2} for k=1,2,≥3k=1,2,\ge3, and the run has ODE time αδ\alpha\delta. The martingale over all αN\alpha N steps is δα/N\delta\sqrt{\alpha/N} by Doob, and it must stay below the lead N−1/2N^{-1/2}, so αδ2=O(1)\alpha\delta^2=O(1) and the available ODE time is at most α\sqrt\alpha; escape requires α≳log⁡N\sqrt\alpha\gtrsim\log N or N(k−2)/2N^{(k-2)/2}, which is αc(N,k)\alpha_c(N,k) up to logs. During the search the noise is as large as mtm_t, so the trajectory does not concentrate, and the proof works without a law of large numbers.

Step 2

The normalization costs less than half the drift as long as δ≲mk−2\delta\lesssim m^{k-2}. Since ∇L⊥Xt−1\nabla L\perp X_{t-1}, the spherical step is the Step 1 increment divided by rt=∥X~t∥r_t=\|\tilde X_t\|, mt=1rt(mt−1−δN∇Φ⋅θ−δN∇Ht⋅θ),1≤rt≤1+δ2(AN2+LtN),Lt=∥∇Ht∥2N.m_t=\frac1{r_t}\Big(m_{t-1}-\tfrac{\delta}{N}\nabla\Phi\cdot\theta-\tfrac{\delta}{N}\nabla H^t\cdot\theta\Big),\qquad1\le r_t\le1+\delta^2\Big(\tfrac{A}{N^2}+\tfrac{L_t}{N}\Big),\quad L_t=\tfrac{\|\nabla H^t\|^2}{N}. The division shrinks mm by about δ2Lt∣m∣/N\delta^2L_t|m|/N per step, second order but with LtL_t of order one, and it stays below the drift δakmk−1/N\delta a_km^{k-1}/N iff δ≲mk−2\delta\lesssim m^{k-2}. At m≈N−1/2m\approx N^{-1/2} this is the step cap δ≲N−(k−2)/2/log⁡N\delta\lesssim N^{-(k-2)/2}/\log N, a second reason small steps are forced for k≥3k\ge3. The rigorous version truncates LtL_t and uses Markov and Freedman's inequality with a union bound.

Step 3

The martingale stays below half the initial correlation, which gives a one-sided comparison. The proof does not need the noise to average out, only to stay below half the initial correlation, and the drift then compounds on the remaining half. Doob's L2L^2 maximal inequality keeps the directional martingale below m0/10m_0/10 up to tˉ=Nγ2/(D2δ2)\bar t=N\gamma^2/(D^2\delta^2), so with Step 2, for m0≥γ/Nm_0\ge\gamma/\sqrt N, with probability 1−C/D2−o(1)1-C/D^2-o(1), mt ≥ m02+δak8N∑j=0t−1mjk−1for all t≤τη+∧tˉ∧M.m_t\ \ge\ \frac{m_0}{2}+\frac{\delta a_k}{8N}\sum_{j=0}^{t-1}m_j^{k-1}\qquad\text{for all }t\le\tau^+_\eta\wedge\bar t\wedge M. The condition tˉ≥M\bar t\ge M is the noise budget αδ2≲1\alpha\delta^2\lesssim1 of Step 1.

Step 4

Solving the inequality gives the three escape times. It is a discrete version of the drift ODE started at m0/2m_0/2. For k=1k=1 growth is linear, t∗≍N/δt^*\asymp N/\delta; for k=2k=2 discrete Grönwall gives mt≥m02eδa2t/(8N)m_t\ge\frac{m_0}{2}e^{\delta a_2t/(8N)}, so t∗≍Nlog⁡N/δt^*\asymp N\log N/\delta; for k≥3k\ge3 the discrete Bihari-LaSalle inequality compares with the finite-time blow-up of m˙=cmk−1\dot m=cm^{k-1}, mt ≥ m0(1−δak8N(k−2)m0k−2 t)1/(k−2),t∗≍N⋅N(k−2)/2δ.m_t\ \ge\ \frac{m_0}{\big(1-\frac{\delta a_k}{8N}(k-2)m_0^{k-2}\,t\big)^{1/(k-2)}},\qquad t^*\asymp\frac{N\cdot N^{(k-2)/2}}{\delta}. The budget checks t∗≤Mt^*\le M and t∗≤tˉt^*\le\bar t become, with δ≪α−1/2\delta\ll\alpha^{-1/2}, Theorem 1.3's conditions; for k=2k=2 the log⁡2N\log^2N is exactly α≫log⁡N\sqrt\alpha\gg\log N. A uniform start has m0≥γ/Nm_0\ge\gamma/\sqrt N with probability near 1.

Step 5

After weak recovery, a law of large numbers by Grönwall. From m=ηm=\eta the drift is order one per unit ODE time, so noise of total size o(1)o(1) cannot matter. Linearizing SGD and the population iterate, the dropped normalization terms sum to O(αδ2)=o(1)O(\alpha\delta^2)=o(1), Doob controls the martingale, and Lipschitz continuity of the drift with discrete Grönwall gives sup⁡∣mt−mˉt∣≤CεeKT\sup|m_t-\bar m_t|\le C\varepsilon e^{KT} on a window long enough for the population iterate to reach 1−γ/51-\gamma/5, after which Assumption A keeps SGD close. The strong Markov property at τη+\tau^+_\eta connects this to Step 4.

Step 6

The lower bound runs the same comparison upward, excursion by excursion. Upward the normalization only helps, since rt≥1r_t\ge1, so while m>0m\gt0 mt≤mt−1+δN32akmt−1k−1−δN∇Ht⋅θ.m_t\le m_{t-1}+\tfrac{\delta}{N}\tfrac32a_km_{t-1}^{k-1}-\tfrac{\delta}{N}\nabla H^t\cdot\theta. Restart the comparison at every positive excursion: Chebyshev and a union bound keep each noise increment below dN/Nd_N/\sqrt N, so every excursion starts below dNN−1/2d_NN^{-1/2}, and reverse Grönwall or Bihari-LaSalle bounds keep mm below η\eta for order N/δN/\delta, Nlog⁡N/δN\log N/\delta and Nk/2/δN^{k/2}/\delta steps, which exceeds MM under the hypotheses. The step ranges come from the choice of dNd_N: for k=2k=2 exponential growth sees the start only through log⁡dN\log d_N, so dN=N1/4d_N=N^{1/4} covers every δ=O(1)\delta=O(1); for k≥3k\ge3 a larger start shortens the blow-up, which forces αδ2=O(1)\alpha\delta^2=O(1).

Step 7

Online SGD moves the overlap at the rate signal2^2/noise. Damian et al. reread Step 1 as a signal-to-noise ratio. With vtv_t the negative spherical gradient and renormalization, αt+1≈αt+ηt (vt⋅w⋆)−12ηt2∥vt∥2αt,\alpha_{t+1}\approx\alpha_t+\eta_t\,(v_t\cdot w^\star)-\tfrac12\eta_t^2\|v_t\|^2\alpha_t , and with the largest step whose drag stays below the drift the overlap follows α′≈SNR/α\alpha'\approx\mathrm{SNR}/\alpha, SNR=E[v⋅w⋆]2/E∥v∥2\mathrm{SNR}=\mathbb E[v\cdot w^\star]^2/\mathbb E\|v\|^2. Unsmoothed, the signal is of order αk⋆−1\alpha^{k^\star-1} by Hermite orthogonality and the noise of order dd, so α′≈α2k⋆−3/d\alpha'\approx\alpha^{2k^\star-3}/d and escape from d−1/2d^{-1/2} takes d α0−2(k⋆−2)=dk⋆−1d\,\alpha_0^{-2(k^\star-2)}=d^{k^\star-1}. The rest of the proof changes the two halves of this ratio.

Step 8

Below α=λ/d\alpha=\lambda/\sqrt d, smoothing gives the loss a floor of signal. The cap mixes in the random tangent step λz\lambda z, whose overlap with w⋆w^\star is about λ/d\lambda/\sqrt d. Expanding binomially, odd moments of z⋅w⋆z\cdot w^\star vanish and even ones are ≍d−j\asymp d^{-j}: Lλ(αk) ≍ (α1+λ2)k∑j≤k/2(λ2α2d)j.\mathcal L_\lambda(\alpha^k)\ \asymp\ \Big(\frac{\alpha}{\sqrt{1+\lambda^2}}\Big)^k\sum_{j\le k/2}\Big(\frac{\lambda^2}{\alpha^2d}\Big)^j . The sum is dominated by its first term when α≥λ/d\alpha\ge\lambda/\sqrt d and by its last, the far side of the cap, below. Smoothing nearly commutes with differentiation (a spherical Stein identity), so the smoothed drift is λ−1(α/λ)k⋆−1\lambda^{-1}(\alpha/\lambda)^{k^\star-1} above the threshold and, below it, the constant λ−1d−(k⋆−1)/2\lambda^{-1}d^{-(k^\star-1)/2} for odd k⋆k^\star or linear in α\alpha for even k⋆k^\star: near the equator the smoothed landscape looks like information exponent 1 or 2, with a small slope.

Step 9

Smoothing shrinks the gradient noise by λ−2k⋆\lambda^{-2k^\star}, as long as λ≤d1/4\lambda\le d^{1/4}. The per-sample loss depends on all of xx, so Step 8 does not reach the noise. Hermite tensors make smoothing linear, Hek(w⋅x)=⟨Hek(x),w⊗k⟩\mathrm{He}_k(w\cdot x)=\langle\mathrm{He}_k(x),w^{\otimes k}\rangle, and their orthogonality gives Ex(u⋅∇wLλHek(w⋅x))2=k! ∥∇wTk(w)[u]∥F2 ≲ k! k21+λ2 Lλ(αk−1)∣α=1/1+λ2:\mathbb E_x\big(u\cdot\nabla_w\mathcal L_\lambda\mathrm{He}_k(w\cdot x)\big)^2=k!\,\big\|\nabla_wT_k(w)[u]\big\|_F^2\ \lesssim\ \frac{k!\,k^2}{1+\lambda^2}\,\mathcal L_\lambda(\alpha^{k-1})\Big|_{\alpha=1/\sqrt{1+\lambda^2}}: the noise in a direction is the smoothed power at the overlap of the cap center with a typical cap point. That sits in Step 8's first regime exactly when 1/λ≥λ/d1/\lambda\ge\lambda/\sqrt d, and then E∥v∥2≲dλ−2k⋆\mathbb E\|v\|^2\lesssim d\lambda^{-2k^\star}. Above the threshold signal2^2 and noise shrink alike; below it the signal has its floor, a gain of λ2k⋆−2\lambda^{2k^\star-2} in SNR for odd k⋆k^\star. Past λ2≈d\lambda^2\approx\sqrt d the noise falls only like λ−2\lambda^{-2}, no faster than the squared signal, so the gain stops at λ=d1/4\lambda=d^{1/4}.

Step 10

Stage 1 escapes from d−1/2d^{-1/2} to λd−1/2\lambda d^{-1/2}. A one-step progress bound makes Step 7 rigorous, and an induction with the Rosenthal-Burkholder-Pinelis inequality keeps αt≥12(α0+η∑s<tcλ(αs))\alpha_t\ge\frac12(\alpha_0+\eta\sum_{s<t}c_\lambda(\alpha_s)), as in Step 3. The step η=d−k⋆/2λ2k⋆−2/polylog(d)\eta=d^{-k^\star/2}\lambda^{2k^\star-2}/\mathrm{polylog}(d) keeps the drag below half the drift, and for odd k⋆k^\star the constant drift covers the distance in T1≈λd−1/2ηλ−1d−(k⋆−1)/2=dk⋆−1λ−2k⋆+4T_1\approx\frac{\lambda d^{-1/2}}{\eta\lambda^{-1}d^{-(k^\star-1)/2}}=d^{k^\star-1}\lambda^{-2k^\star+4} steps, dk⋆/2d^{k^\star/2} at λ=d1/4\lambda=d^{1/4}; for even k⋆k^\star a Bihari-type comparison gives exponential growth in the same number of steps up to logs.

Step 11

The rest of the trajectory is the unsmoothed problem. Above λd−1/2\lambda d^{-1/2} smoothing no longer changes the SNR, so Ben Arous et al.'s analysis at the same step carries α\alpha to 1−1/k⋆1-1/k^\star, and a contraction for 1−αt1-\alpha_t brings it to d−1/4d^{-1/4} by T1T_1. Then λ=0\lambda=0 and ηt∝1/(d+t)\eta_t\propto1/(d+t): near w⋆w^\star the loss is locally strongly convex, and the Robbins-Monro induction gives 1−αt≲d/(d+t)1-\alpha_t\lesssim d/(d+t), the d/εd/\varepsilon term.

Step 12

The CSQ lower bound comes from many nearly orthogonal directions. There are 12ecϵ2d\frac12e^{c\epsilon^2d} unit vectors with pairwise overlaps at most ϵ\epsilon, and Hermite orthogonality makes their models nearly uncorrelated, ∣E σ(wi⋅x)σ(wj⋅x)∣=∣∑kck2k!(wi⋅wj)k∣≤ϵk⋆.\big|\mathbb E\,\sigma(w_i\cdot x)\sigma(w_j\cdot x)\big|=\Big|\sum_k\frac{c_k^2}{k!}(w_i\cdot w_j)^k\Big|\le\epsilon^{k^\star}. Szörényi's CSQ lemma then demands about s(τ2−ϵk⋆)/2s(\tau^2-\epsilon^{k^\star})/2 queries, and ϵ≍log⁡(qdk⋆/2)/d\epsilon\asymp\sqrt{\log(qd^{k^\star/2})/d} gives the bound. The same ϵk⋆\epsilon^{k^\star} that flattens the landscape at the equator makes eΩ(d)e^{\Omega(d)} directions indistinguishable to correlational queries.

Consequences and limits. For online SGD the difficulty of a rank-one estimation problem is one integer read off the population loss at the equator, for single-index models the degree of the link's lowest nonzero Hermite coefficient. Small changes to the activation move the cost by powers of NN (x3x^3 needs NN samples, x3−3xx^3-3x about N2N^2), while from a warm start every problem is linear, and almost all of the samples used before strong recovery go into the search. The dk⋆−1d^{k^\star-1} barrier is specific to online SGD on the raw loss: smoothing reaches the CSQ rate, and experiments with σ=Hek⋆/k⋆!\sigma=\mathrm{He}_{k^\star}/\sqrt{k^\star!} fit exponents 1.5, 2.1 and 2.5 for k⋆=3,4,5k^\star=3,4,5. For the spiked tensor, α≍Nk−2\alpha\asymp N^{k-2} matches the threshold conjectured for first-order methods without spectral initialization, and as λ→∞\lambda\to\infty the smoothed gradient becomes the partial-trace estimator of tensor PCA, a link the paper sketches without proof. Damian et al. conjecture, without proof, that minibatch SGD gets the same benefit implicitly. Limits: rank one, known norm, one pass, a start at the CLT scale (m0=0m_0=0 can be trapped), log⁡N\log N to log⁡2N\log^2N gaps, and a refutation specific to online SGD in its step range; smoothing needs a known link and isotropic Gaussian inputs, an initial overlap of the right sign (constant probability from a random start), and is optimal only within CSQ, since non-correlational methods (Chen and Meka) learn polynomial links with about dd samples.

Terms

m=x⋅θm=x\cdot\theta
correlation (latitude) with the truth; a random start has m≈N−1/2m\approx N^{-1/2}, the equator
Information exponent kk
order of the first nonzero derivative of ϕ\phi at m=0m=0; drift near the equator is akmk−1a_km^{k-1}
Hermite rule
single-index square loss: ϕ(m)=2∑juj(f)2(1−mj)+const\phi(m)=2\sum_ju_j(f)^2(1-m^j)+\text{const}, so kk is the first j≥1j\ge1 with uj≠0u_j\ne0
αc(N,k)\alpha_c(N,k)
1, log⁡N, Nk−21,\ \log N,\ N^{k-2}; sample complexity αcN\alpha_cN up to log⁡2N\log^2N
Directional error martingale
δN∑∇Hj⋅θ\frac{\delta}{N}\sum\nabla H^{j}\cdot\theta; size δT/N\delta\sqrt T/N, so αδ2=O(1)\alpha\delta^2=O(1) keeps it below N−1/2N^{-1/2}
Step cap δˉN(k)\bar\delta_N(k)
≈N−(k−2)/2/log⁡N\approx N^{-(k-2)/2}/\log N: normalization shrinkage δ2m/N\delta^2m/N must stay below the drift δmk−1/N\delta m^{k-1}/N
Lλ\mathcal L_\lambda
average over (w+λz)/∥w+λz∥(w+\lambda z)/\|w+\lambda z\|, zz uniform on the tangent sphere
λ/d\lambda/\sqrt d
overlap reached by the random part of a cap point; below it smoothing adds signal
E∥v∥2≲dλ−2k⋆\mathbb E\|v\|^2\lesssim d\lambda^{-2k^\star}
smoothed gradient noise, valid for λ≤d1/4\lambda\le d^{1/4}
CSQ bound
τ≲(log⁡(qd)/d)k⋆/4\tau\lesssim(\log(qd)/d)^{k^\star/4}, so n≳dk⋆/2n\gtrsim d^{k^\star/2} with τ≈n−1/2\tau\approx n^{-1/2}

Quiz

  1. Gaussian single-index model, square loss: f(x)=x3f(x)=x^3 has k=1k=1, but f(x)=x3−3xf(x)=x^3-3x has k=3k=3. Why?
  2. What does Theorem 1.4 (the refutation) rule out?
  3. For k≥3k\ge3, where does the dimension enter the escape cost?
  4. Why is the smoothing radius capped at λ=d1/4\lambda=d^{1/4}?

Sources