11 research outputs found
Efficient Computation of Location Depth Contours by Methods of Computational Geometry
The concept of location depth was introduced as a way to extend the univariate notion of ranking to a bivariate configuration of data points. It has been used successfully for robust estimation, hypothesis testing, and graphical display. The depth contours form a collection of nested polygons, and the center of the deepest contour is called the Tukey median. The only available implemented algorithms for the depth contours and the Tukey median are slow, which limits their usefulness. In this paper we describe an optimal algorithm which computes all bivariate depth contours in O(n 2) time and space, using topological sweep of the dual arrangement of lines. Once these contours are known, the location depth of any point can be computed in O(log 2 n) time with no additional preprocessing or in O(log n) time after O(n 2) preprocessing. We provide fast implementations of these algorithms to allow their use in everyday statistical practice
Fast computation of Tukey trimmed regions and median in dimension
Given data in , a Tukey -trimmed region is the set of
all points that have at least Tukey depth w.r.t. the data. As they are
visual, affine equivariant and robust, Tukey regions are useful tools in
nonparametric multivariate analysis. While these regions are easily defined and
interpreted, their practical use in applications has been impeded so far by the
lack of efficient computational procedures in dimension . We construct
two novel algorithms to compute a Tukey -trimmed region, a na\"{i}ve
one and a more sophisticated one that is much faster than known algorithms.
Further, a strict bound on the number of facets of a Tukey region is derived.
In a large simulation study the novel fast algorithm is compared with the
na\"{i}ve one, which is slower and by construction exact, yielding in every
case the same correct results. Finally, the approach is extended to an
algorithm that calculates the innermost Tukey region and its barycenter, the
Tukey median
Efficient Computation of Location Depth Contours by Methods of Computational Geometry
The concept of location depth was introduced as a way to extend the univariate notion of ranking to a bivariate configuration of data points. It has been used successfully for robust estimation, hypothesis testing, and graphical display. The depth contours form a collection of nested polygons, and the center of the deepest contour is called the Tukey median. The only available implemented algorithms for the depth contours and the Tukey median are slow, which limits their usefulness. In this paper we describe an optimal algorithm which computes all bivariate depth contours in O(n²) time and space, using topological sweep of the dual arrangement of lines. Once these contours are known, the location depth of any point can be computed in O(log² n) time with no additional preprocessing or in O(log n) time after O(n²) preprocessing. We provide fast implementations of these algorithms to allow their use in everyday statistical practice