English

Approximating the Simplicial Depth

Computational Geometry 2015-12-29 v2

Abstract

Let PP be a set of nn points in dd-dimensions. The simplicial depth, σP(q)\sigma_P(q) of a point qq is the number of dd-simplices with vertices in PP that contain qq in their convex hulls. The simplicial depth is a notion of data depth with many applications in robust statistics and computational geometry. Computing the simplicial depth of a point is known to be a challenging problem. The trivial solution requires O(nd+1)O(n^{d+1}) time whereas it is generally believed that one cannot do better than O(nd1)O(n^{d-1}). In this paper, we consider approximation algorithms for computing the simplicial depth of a point. For d=2d=2, we present a new data structure that can approximate the simplicial depth in polylogarithmic time, using polylogarithmic query time. In 3D, we can approximate the simplicial depth of a given point in near-linear time, which is clearly optimal up to polylogarithmic factors. For higher dimensions, we consider two approximation algorithms with different worst-case scenarios. By combining these approaches, we compute a (1+ε)(1+\varepsilon)-approximation of the simplicial depth in time O~(nd/2+1)\tilde{O}(n^{d/2 + 1}) ignoring polylogarithmic factor. All of these algorithms are Monte Carlo algorithms. Furthermore, we present a simple strategy to compute the simplicial depth exactly in O(ndlogn)O(n^d \log n) time, which provides the first improvement over the trivial O(nd+1)O(n^{d+1}) time algorithm for d>4d>4. Finally, we show that computing the simplicial depth exactly is #P-complete and W[1]-hard if the dimension is part of the input.

Keywords

Cite

@article{arxiv.1512.04856,
  title  = {Approximating the Simplicial Depth},
  author = {Peyman Afshani and Donald R. Sheehy and Yannik Stein},
  journal= {arXiv preprint arXiv:1512.04856},
  year   = {2015}
}

Comments

25 pages, 4 figures

R2 v1 2026-06-22T12:10:26.918Z