Approximating the Simplicial Depth
Abstract
Let be a set of points in -dimensions. The simplicial depth, of a point is the number of -simplices with vertices in that contain 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 time whereas it is generally believed that one cannot do better than . In this paper, we consider approximation algorithms for computing the simplicial depth of a point. For , 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 -approximation of the simplicial depth in time ignoring polylogarithmic factor. All of these algorithms are Monte Carlo algorithms. Furthermore, we present a simple strategy to compute the simplicial depth exactly in time, which provides the first improvement over the trivial time algorithm for . Finally, we show that computing the simplicial depth exactly is #P-complete and W[1]-hard if the dimension is part of the input.
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