skip to main content


Title: Hidden Integrality and Semirandom Robustness of SDP Relaxation for Sub-Gaussian Mixture Model
We consider the problem of estimating the discrete clustering structures under the sub-Gaussian mixture model. Our main results establish a hidden integrality property of a semidefinite programming (SDP) relaxation for this problem: while the optimal solution to the SDP is not integer-valued in general, its estimation error can be upper bounded by that of an idealized integer program. The error of the integer program, and hence that of the SDP, are further shown to decay exponentially in the signal-to-noise ratio. In addition, we show that the SDP relaxation is robust under the semirandom setting in which an adversary can modify the data generated from the mixture model. In particular, we generalize the hidden integrality property to the semirandom model and thereby show that SDP achieves the optimal error bound in this setting. These results together highlight the “global-to-local” mechanism that drives the performance of the SDP relaxation. To the best of our knowledge, our result is the first exponentially decaying error bound for convex relaxations of mixture models. A corollary of our results shows that in certain regimes, the SDP solutions are in fact integral and exact. More generally, our results establish sufficient conditions for the SDP to correctly recover the cluster memberships of [Formula: see text] fraction of the points for any [Formula: see text]. As a special case, we show that under the [Formula: see text]-dimensional stochastic ball model, SDP achieves nontrivial (sometimes exact) recovery when the center separation is as small as [Formula: see text], which improves upon previous exact recovery results that require constant separation.  more » « less
Award ID(s):
2047910 1704828 2233152
NSF-PAR ID:
10327714
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Mathematics of Operations Research
ISSN:
0364-765X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. For generalized Korteweg–De Vries (KdV) models with polynomial nonlinearity, we establish a local smoothing property in [Formula: see text] for [Formula: see text]. Such smoothing effect persists globally, provided that the [Formula: see text] norm does not blow up in finite time. More specifically, we show that a translate of the nonlinear part of the solution gains [Formula: see text] derivatives for [Formula: see text]. Following a new simple method, which is of independent interest, we establish that, for [Formula: see text], [Formula: see text] norm of a solution grows at most by [Formula: see text] if [Formula: see text] norm is a priori controlled. 
    more » « less
  2. null (Ed.)
    Consider a fractional Brownian motion (fBM) [Formula: see text] with Hurst index [Formula: see text]. We construct a probability space supporting both B H and a fully simulatable process [Formula: see text] such that[Formula: see text] with probability one for any user-specified error bound [Formula: see text]. When [Formula: see text], we further enhance our error guarantee to the α-Hölder norm for any [Formula: see text]. This enables us to extend our algorithm to the simulation of fBM-driven stochastic differential equations [Formula: see text]. Under mild regularity conditions on the drift and diffusion coefficients of Y, we construct a probability space supporting both Y and a fully simulatable process [Formula: see text] such that[Formula: see text] with probability one. Our algorithms enjoy the tolerance-enforcement feature, under which the error bounds can be updated sequentially in an efficient way. Thus, the algorithms can be readily combined with other advanced simulation techniques to estimate the expectations of functionals of fBMs efficiently. 
    more » « less
  3. We investigate the problem of recovering jointly [Formula: see text]-rank and [Formula: see text]-bisparse matrices from as few linear measurements as possible, considering arbitrary measurements as well as rank-one measurements. In both cases, we show that [Formula: see text] measurements make the recovery possible in theory, meaning via a nonpractical algorithm. In case of arbitrary measurements, we investigate the possibility of achieving practical recovery via an iterative-hard-thresholding algorithm when [Formula: see text] for some exponent [Formula: see text]. We show that this is feasible for [Formula: see text], and that the proposed analysis cannot cover the case [Formula: see text]. The precise value of the optimal exponent [Formula: see text] is the object of a question, raised but unresolved in this paper, about head projections for the jointly low-rank and bisparse structure. Some related questions are partially answered in passing. For rank-one measurements, we suggest on arcane grounds an iterative-hard-thresholding algorithm modified to exploit the nonstandard restricted isometry property obeyed by this type of measurements. 
    more » « less
  4. Auxetic foams exhibit novel mechanical properties due to their unique microstructure for improved energy-absorption and cavity expansion applications that have fascinated the scientific community since their inception. Given the advancements in material processing and performance of polymer open cell auxetic foams, there is a strong desire to fully understand the nonlinear rate-dependent deformation of these materials. The influence of nonlinear compressibility is introduced here along with relaxation effects to improve model predictions for different stretch rates and finite deformation regimes. The viscoelastic behavior of the material is analyzed by comparing fractional order and integer order calculus models. All results are statistically validated using maximum entropy methods to obtain Bayesian posterior densities for the hyperelastic, auxetic, and viscoelastic parameters. It is shown that fractional order viscoelasticity provides [Formula: see text]–[Formula: see text] improvement in prediction over integer order viscoelastic models when the model is calibrated at higher stretch rates where viscoelasticity is more significant.

     
    more » « less
  5. We consider a [Formula: see text] system of hyperbolic balance laws that is the converted form under inverse Hopf–Cole transformation of a Keller–Segel type chemotaxis model. We study Cauchy problem when Cauchy data connect two different end-states as [Formula: see text]. The background wave is a diffusive contact wave of the reduced system. We establish global existence of solution and study the time asymptotic behavior. In the special case where the cellular population initially approaches its stable equilibrium value as [Formula: see text], we obtain nonlinear stability of the diffusive contact wave under smallness assumption. In the general case where the population initially does not approach to its stable equilibrium value at least at one of the far fields, we use a correction function in the time asymptotic ansatz, and show that the population approaches logistically to its stable equilibrium value. Our result shows two significant differences when comparing to Euler equations with damping. The first one is the existence of a secondary wave in the time asymptotic ansatz. This implies that our solutions converge to the diffusive contact wave slower than those of Euler equations with damping. The second one is that the correction function logistically grows rather than exponentially decays. 
    more » « less