Alphabeta Math
False statementConstruction: AI-adaptedVerification: AI-generatedPipeline-generatedprecheck pass
How statement and proof provenance work

The first chip identifies the source of the statement or construction; the second identifies the source of its local proof or verification.

  • Literature-sourced: the exact statement appears in a cited source; only wording and notation differ.
  • AI-adapted: a semantically identical restatement of literature-sourced material, modulo indexing, notation, and boundary cases adopted by the library.
  • AI-generated: a genuinely novel statement formulated by AI, with no source for the claim itself.

These labels describe origin, not correctness: citations and verification chips remain separate evidence.

The Cmax-only Wainwright lower signed-support threshold can predict failure when Lasso succeeds

Statement

Assume the Axiom of Choice. Wainwright's planned lower-threshold claim is the following: along a sequence of real Gaussian linear models with S={j:βj∗≠0}, s=∣S∣, N=p−s, s,n−s,N→∞, and independent observations yi=xiTβ∗+wi,xi∼Np(0,Σ),wi∼N(0,σ2), where X has rows xiT, y=(yi), design and noise are independent, there is no separate intercept, and Σ has unit diagonal. Suppose constants Cmin⁡>0, Cmax⁡<∞, Dmax⁡<∞, and 0<ϵ≤1 satisfy λmin⁡(ΣSS)≥Cmin⁡,λmax⁡(Σ)≤Cmax⁡,∥(ΣSS)−1∥∞≤Dmax⁡,∥ΣScS(ΣSS)−1∥∞≤1−ϵ, with ∥M∥∞=max⁡i∑j∣Mij∣. Define θℓ=(Cmax⁡−Cmax⁡−1/Cmax⁡)2Cmax⁡(2−ϵ)2. For every fixed ν>0 and every positive non-increasing sequence λn, if n<2(θℓ−ν)slog⁡N+s+1, then the probability that there exists an optimizer with sgn⁡(β^)=sgn⁡(β∗) for ∥y−Xβ∥22/(2n)+λn∥β∥1 tends to zero.

Facts & Assumptions

Given: The Gaussian model, covariance hypotheses, exact loss normalization, success event, and claimed lower implication in the Statement.

[F1]

For an iid standard-Gaussian matrix of dimension d=o(n), its Gram matrix divided by n converges in operator norm to identity (Gaussian sample covariance converges in operator norm when dimension is little-o of sample size).

[F2]

A sum of q independent squared standard normals divided by n converges in probability to one if q/n→1 and n→∞ (A chi-square variable divided by a comparable dimension converges to one).

[F3]

Under AC, a centered Gaussian vector with covariance C⪰0 can be realized as C1/2G for a vector of independent standard normals (Multivariate normal law, including singular covariance).

[F4]

A nondegenerate multivariate normal has its usual Gaussian density; in particular the density with identity covariance factors into independent standard-normal densities (A positive-definite multivariate normal law has its Gaussian density, Independent random elements have product joint law).

[F5]

A standard normal has density (2π)−1/2e−x2/2 (Standard normal and normal laws). Completing the square gives EetU=evt2/2 for U∼N(0,v); applying Markov to etU with t=1/v, and likewise to e−tU, yields P(∣U∣≥1)≤2e−1/(2v) when v>0 (Markov's inequality for random variables). For v=0 the probability is zero.

[F6]

The Lasso objective is F(β)=∥y−Xβ∥22/(2n)+λ∥β∥1 (Fixed-design least-squares, ridge, lasso, group, and fused objectives).

[A1]

AC is assumed for the finite Gaussian constructions and suppliers (The Axiom of Choice).

Refutation

1.1A1F3F4givenconstruct

Fix δ=ϵ=10−6, c=(1−ϵ)/2, a=1/1000, and ν=1/100000. For integers N≥2 tending to infinity set s=N3,p=s+N,n=s+⌈aslog⁡N⌉,S={1,…,s},λn=1. In each observation let the first 2N active coordinates consist of independent pairs (Uj,Vj), 1≤j≤N, each with covariance B=(1−(1−δ)−(1−δ)1), and let the other s−2N active coordinates be independent standard normals, independent of these pairs. Take independent standard normals Zj, independent also of all active variables, and define the inactive coordinates by Tj=c(Uj−Vj)+rZj,r=1−(1−ϵ)2(1−δ/2). Take independent observations and independent standard-normal noise wi; set βi∗=n on S and zero off S. These finite Gaussian constructions define the model in the Statement, with σ2=1 and no intercept.

2.1step 1.1algebra

The eigenvalues of B are δ and 2−δ, and each row of B−1 has absolute sum 1/δ. Thus Cmin⁡=δ and Dmax⁡=1/δ work. Each inactive coordinate regresses on its pair with coefficients (c,−c), so ∥ΣScSΣSS−1∥∞=2c=1−ϵ. Its variance is one by the definition of r, and 0<r<2ϵ+δ/2<3⋅10−6. The full covariance, after reordering, is a direct sum of positive-definite three-dimensional covariance blocks and singleton identity blocks: each triple is an invertible linear image of (Uj,Vj,Zj) because r>0. Every triple has trace three, so every eigenvalue is at most three. Consequently Cmax⁡=3 bounds the maximum eigenvalue of the full covariance, and all its diagonal entries are one.

3.1givenstep 1.1step 2.1algebra

For these constants, 3−8/3=1/33+8/3>112,3(2−ϵ)2<12, so the precise claimed constant satisfies θℓ>1/1728 and 2(θℓ−ν)>2(1/1728−1/100000)>a. Since ⌈aslog⁡N⌉≤aslog⁡N+1, the claimed failure inequality n<2(θℓ−ν)slog⁡N+s+1 holds. Also s,n−s,N→∞ and s/n→0. The constant penalty is positive and non-increasing. The lower claim puts no restriction on ρn=n.

4.1F1F3F5step 2.1step 3.1constructalgebra

Write XS=ZΣSS1/2, with Z an n×s iid standard-Gaussian matrix, and put G=XSTXS and A=G/n. By [F1] and s/n→0, with probability tending to one the minimum eigenvalue of A is at least δ/2 and ∥A−1∥op≤2/δ. On this event define β^S=n1+η−A−11,η=A−1XSTw/n,β^Sc=0. The shrinkage term satisfies ∥A−11∥∞≤2s/δ=o(n). Conditional on XS, η is centered Gaussian with covariance A−1/n, so [F5] and a finite union bound give P(∥η∥∞≥1∣XS)≤2se−nδ/4⟶0 uniformly on the spectral event. Thus every active candidate coordinate is positive with probability tending to one.

4.2A1F2F3F4step 1.1step 2.1step 3.1algebrachoose

Put b=1s and T=bTΣSS−1b=s−2N+2N/δ, so T/s→1. To compute Q=bTG−1b, put v=ΣSS−1/2b and choose a fixed orthogonal matrix with first column v/∥v∥2. Orthogonal invariance and [F4] make the rotated columns of Z, written [h,H], mutually independent standard-Gaussian columns. They have full column rank almost surely: condition successively on earlier columns; for each fixed proper column subspace, projection onto any fixed nonzero normal vector has a continuous normal law, so a new standard-Gaussian column belongs to that subspace with probability zero. For full-rank H put PH=H(HTH)−1HT,D=hT(I−PH)h. Block elimination in the Gram matrix of [h,H] gives exactly Q=T/D. For each fixed full-rank H, the projector I−PH has rank n−s+1, so orthogonal invariance gives the conditional law D∼χn−s+12. This conditional law is the same for every such H; the projector is a measurable function of H, and no random basis is selected. Since (n−s+1)/n→1, [F2] gives D/n→1 in probability, and therefore Q/(s/n)→1 in probability.

5.1F2F3F4step 1.1step 4.1step 4.2algebra

Let P=XSG−1XST. The candidate residual is u=y−XSβ^S=(I−P)w+nXSG−1b,XSTu/n=b. Its two displayed summands are orthogonal. For inactive column j, the construction gives Xj=XSaj+rzj, where aj has entries (c,−c) on its pair and zero elsewhere, and the zj are mutually independent standard-Gaussian vectors independent of XS,w. Because ajTb=0, the dual variables Vj=XjTu/n=rzjTu/n are, conditional on XS,w, independent centered normals of common variance rM,M=∥u∥22/n2=Q+∥(I−P)w∥22/n2. By [F2], ∥w∥22/n→1 in probability. Since the second term is at most ∥w∥22/n2, step 4.2 and s→∞ show P(M≤4s/n)→1. This event depends only on XS,w and is independent of all inactive residual vectors.

6.1F5step 1.1step 2.1step 4.1step 5.1algebra

On the event M≤4s/n, since n≥aslog⁡N and r<3⋅10−6, rM≤4ralog⁡N<3250log⁡N. Applying [F5] and a finite union bound conditionally on XS,w yields the uniform bound P ⁣(max⁡1≤j≤N∣Vj∣≥1 | XS,w)≤2Nexp⁡ ⁣(−12rM)≤2N1−125/3⟶0. Removing the conditioning and adding the probability of M>4s/n proves max⁡j∣Vj∣<1 with probability tending to one. Intersecting with the event in step 4.1 preserves probability tending to one and gives both active sign correctness and inactive dual feasibility.

7.1F6step 1.1step 3.1step 4.1step 5.1step 6.1algebra∎

On this intersection put z=XTu/n. It equals one on S and satisfies ∣zj∣<1 off S. Since the candidate is positive on S and zero off S, the scalar inequality ∣β^i+hi∣−∣β^i∣≥zihi holds for every coordinate and every real hi. For every h∈Rp, expand [F6], with λn=1, to obtain F(β^+h)−F(β^)=∥Xh∥222n−hTXTu/n+∥β^+h∥1−∥β^∥1≥0. Thus this candidate is an optimizer with sgn⁡(β^)=sgn⁡(β∗). Its existence has probability tending to one, whereas step 3.1 meets the exact hypotheses that predict probability tending to zero. This refutes the stated lower implication.

Source notes

This refutes the conclusion itself, independently of the erroneous inverse-block identity in the 2006 Appendix D.6. The construction uses s/(p−s)→∞, which the stated source and planned hypotheses permit; it does not refute a theorem with an additional dimensional restriction. The 2009 journal article's equation (27) defines ρℓ(C)=12min⁡i≠j(Cii+Cjj−2Cij), a minimum pairwise increment-variance parameter, not a least eigenvalue. Its revised threshold interface is distinct and is not substituted for the claim refuted here.

Depends on

Used by

Nothing in the library uses this result yet.

Dependency tree · two levels

46 results within two dependency steps of this one, each drawn at its shortest distance from it. An arrow runs from a result to what uses it, so the chart reads left to right and ends at this result, which carries a heavier outline. Every node is a link to that result. Click elsewhere on the chart to enlarge it.

Sources