Neighborliness of randomly projected simplices in high dimensions

  1. David L. Donoho* and
  2. Jared Tanner
  1. Department of Statistics, Sequoia Hall, Serra Mall 370, Stanford University, Stanford, CA 94305-4065
  1. Contributed by David L. Donoho, March 30, 2005

Abstract

Let A be a d × n matrix and T = Tn -1 be the standard simplex in Rn. Suppose that d and n are both large and comparable: d ≈ δn, δ ∈ (0, 1). We count the faces of the projected simplex AT when the projector A is chosen uniformly at random from the Grassmann manifold of d-dimensional orthoprojectors of Rn. We derive ρN(δ) > 0 with the property that, for any ρ < ρN(δ), with overwhelming probability for large d, the number of k-dimensional faces of P = AT is exactly the same as for T, for 0 ≤ k ≤ ρd. This implies that P is Formula-neighborly, and its skeleton Formula is combinatorially equivalent to Formula. We also study a weaker notion of neighborliness where the numbers of k-dimensional faces fk(P)fk(T)(1 - ε). Vershik and Sporyshev previously showed existence of a threshold ρVS(δ) > 0 at which phase transition occurs in k/d. We compute and display ρVS and compare with ρN. Corollaries are as follows. (1) The convex hull of n Gaussian samples in Rd, with n large and proportional to d, has the same k-skeleton as the (n - 1) simplex, for k < ρN (d/n)d(1 + oP(1)). (2) There is a “phase transition” in the ability of linear programming to find the sparsest nonnegative solution to systems of underdetermined linear equations. For most systems having a solution with fewer than ρVS(d/n)d(1 + o(1)) nonzeros, linear programming will find that solution.

Footnotes

  • * To whom correspondence should be addressed. E-mail: donoho{at}stat.stanford.edu.

  • Author contributions: D.L.D. designed research; D.L.D. and J.T. performed research; J.T. analyzed data, developed high-accuracy computational tools, and performed high-accuracy calculations; and D.L.D. and J.T. wrote the paper.

« Previous | Next Article »Table of Contents