Exact and approximate computation of the scatter halfspace depth

The scatter halfspace depth (sHD) is an extension of the location halfspace (also called Tukey) depth that is applicable in the nonparametric analysis of scatter. Using sHD, it is possible to define minimax optimal robust scatter estimators for multivariate data. The problem of exact computation of sHD for data of dimension $d \geq 2$ has, however, not been addressed in the literature. We develop an exact algorithm for the computation of sHD in any dimension $d$ and implement it efficiently using C++ for $d \leq 5$, and in R for any dimension $d \geq 1$. Since the exact computation of sHD is slow especially for higher dimensions, we also propose two fast approximate algorithms. All our programs are freely available in the R package scatterdepth.