---
title: Tukey Median in Multivariate Analysis
url: https://www.emergentmind.com/topics/location-halfspace-median
type: topic
---

# Tukey Median in Multivariate Analysis

The location halfspace median, universally referred to as the Tukey median, is a canonical multivariate center estimator grounded in the concept of statistical depth. Defined via the minimum probability mass over all closed halfspaces containing a given point, it generalizes the univariate median to arbitrary dimension. The Tukey median is the maximizer(s) of halfspace depth for a given measure or data set, yielding a robust, affine-equivariant center with rich theoretical and computational properties. This article develops the definition, key structural results, robustness metrics, algorithmic techniques, and statistical guarantees associated with the location halfspace median, integrating contemporary advances for high-dimensional inference.

## 1. Formal Definition and Basic Structural Properties

Let $X_1, …, X_n \subset \mathbb{R}^d$ be a data sample, with empirical measure $\mu(A) = \frac{1}{n} \# \{ i : X_i \in A \}$. The halfspace depth of $x \in \mathbb{R}^d$ is defined as
$$
D(x) = \min \{ \mu(H) : x \in H, H \text{ a closed halfspace} \}.
$$
Equivalently, in the empirical case,
$$
D(x) = \min_{u \in S^{d-1}} \frac{1}{n} \# \{ i : u^\top X_i \leq u^\top x \}.
$$
The location halfspace median set is the set of maximizers
$$
C = \{ x \in \mathbb{R}^d : D(x) = \max_{y \in \mathbb{R}^d} D(y) \},
$$
often referred to as $M(\mu)$ or $\mathcal{M}(\mu)$. In general, $C$ may not be a singleton; it is always a closed (possibly lower-dimensional) convex polytope [1201.1171][1207.4988].

Central (level) regions $C_\alpha = \{x : D(x) \geq \alpha\}$ are convex and nested in $\alpha$, with the deepest (non-empty) region serving as the Tukey median set.

## 2. Geometric, Convex Analytic, and Covering Characterizations

Halfspace depth regions admit multiple geometric representations. For the population measure, $D_\alpha(\mu)$ coincides with
- Intersection of all closed halfspaces of mass at least $\alpha$;
- “Floating body” at level $\alpha$: intersection of all halfspaces whose *complements* carry mass at most $1-\alpha$ [2106.00616];
- For empirical measures (samples not in general position), depth regions are finite intersections of “irrotatable” halfspaces: each facet lies in a hyperplane spanned by at least $d$ data points, with rotation about any $(d-2)$-face increasing the cut-off count [1604.07104].

A salient result is the “ray basis theorem,” stating that the union of minimal-mass halfspaces at the median covers the entire space if and only if the point is indeed a median. For smooth measures, this gives a geometric covering property for every median point [2106.00616].

Recent advances introduced flag halfspaces—nested, open-in-face sets—showing that, even for atomic or non-smooth measures, the mass minimizing halfspace in the definition of depth is always attained in a flag halfspace, unifying the theory for both continuous and discrete cases [2209.11445].

## 3. Robustness: Breakdown Point and Affine Equivariance

The finite-sample breakdown point is the maximal fraction of contamination that an estimator tolerates before it can be made arbitrarily large. For the Tukey median, the breakdown point under the additive contamination (Huber) model satisfies
$$
\varepsilon(T^*, X_n) = \frac{\lambda^*}{1+\lambda^*},
$$
where $\lambda^*$ is the deepest sample depth [1604.07104][1604.07039]. Under halfspace symmetry and weak smoothness, the asymptotic breakdown approaches $1/3$:
$$
\varepsilon(T^*, X_n) \rightarrow \frac{1}{3}, \quad n \to \infty.
$$
For data in general position, the breakdown is at least $\lfloor n/(d+1) \rfloor / n$, which is maximal among affine-equivariant location estimators [1207.4988]. Under total variation contamination, the breakdown point is $1/4$ in high dimensions but remains $1/3$ in $d=2$ [2001.07805].

Affine equivariance is intrinsic: any invertible affine map $A$ and $b$ carries medians to medians, with
$$
D(Ax + b; P_{AX + b}) = D(x; P_X),
$$
and thus
$$
M(P_{AX + b}) = A M(P_X) + b,
$$
guaranteeing coordinate-free centrality [1207.4988][1201.1171][2109.00493].

## 4. Computational Methodologies

Computation of the Tukey median is NP-hard in general for $d \ge 3$, but tractable and efficient algorithms exist for small and moderate dimensions:
- **Exact bivariate algorithm:** $O(n \log n)$ via rotating-calipers or “circular sequence” on pairwise lines [1207.4988].
- **Exact general-$d$ algorithms:** Polytope intersection via enumeration of candidate halfspace facets, $O(n^{d-1}\log n)$ [1207.4988][1604.05927].
- **Approximate algorithms:**
  - **Random Tukey:** Sample $K$ random directions, estimate depth as minimum among observed halfspace fractions, $O(K n d)$ [1207.4988].
  - **ABCDepth:** Ball-intersection representation (Merkle equivalence): approximate level sets as intersections of balls containing at least $n-k+1$ data points, yielding median computation complexity linear in $d$ and quadratic in $n$ [1812.03174].
  - **Metric halfspace depth:** Generic anchor set $(i,j)$, compute $p_{ij}$, approximate depth at candidate $y$ by the minimal $p_{ij}$ where $\|y-X_i\| \leq \|y-X_j\|$; total complexity $O(n^3)$ for $n$ candidates [2109.00493].

Modern improvements leverage sharp upper bounds on maximum sample depth, such as $\lambda^* \leq \lfloor (n-p+2)/2 \rfloor / n$ for data in general position ($p$ dimension), avoiding redundant depth region searches and reducing computational burden in large $p/n$ regimes [1604.05927].

## 5. Statistical Properties: Consistency, Efficiency, Convergence Rates

Halfspace median estimators inherit uniform convergence from the finite VC-dimension of halfspace families:
$$
\sup_{x \in \mathbb{R}^d} |D_n(x) - D(x)| = O_p(n^{-1/2}),
$$
yielding
$$
C_\alpha(P_n) \rightarrow C_\alpha(P)
$$
in Hausdorff metric and, if the population median is unique, sample medians converge in probability [1207.4988][1201.1171].

In elliptical and $\alpha$-symmetric models (characteristic function depends on $\|t\|_\alpha$), the Tukey median achieves minimax-optimal convergence rates under contamination:
$$
\|\widehat{\mu}_n - \mu^*\|_\beta \lesssim \varepsilon + \sqrt{\frac{p}{n}} + \sqrt{\frac{\log(1/\delta)}{n}}
$$
in the dual norm $\beta = \frac{\alpha}{\alpha-1}$, for contamination level $\varepsilon < 1/3$ and $n$ large [2512.07634]. No moment or bounded support assumptions are required. For sub-Gaussian or heavy-tailed models, estimation error bounds adapt via the appropriate decay function $h(t)$.

Breakdown point coincides with theoretical upper limits: $1/3$ under symmetry for additive contamination, $1/4$ for total variation. The halfspace-metric projection estimator can push breakdown to $1/2$ if the depth decay function is known, but forfeits affine invariance [2001.07805].

## 6. Topological, Dimensional, and General Position Results

Sample median regions are closed convex polytopes; under general position they cannot be $(d-1)$-dimensional except in degenerate settings ($n=d$) [2209.11445]. In $d=2$, the median set is either a polygon (full-dimensional) or a singleton coinciding with a sample point.

If the data are not in general position, depth regions are still finite intersections of “irrotatable” halfspaces; non-uniqueness of the median is possible, but the average of deepest points provides a canonical center estimate [1604.07104].

In infinite-dimensional Banach or Hilbert spaces, halfspace depth degenerates: for most $x$, $\mathrm{HD}(x;P) = 0$, and the median at the center remains at depth $1/2$ [1201.1171].

## 7. Practical Recommendations and Applications

- For robust multivariate location, the Tukey median provides a nonparametric, affine-equivariant center with breakdown point optimal among affine-equivariant estimators [1207.4988][1201.1171].
- In high dimensions, employ approximate algorithms (e.g., ABCDepth) that scale linearly with $d$ and quadratically with $n$ [1812.03174].
- In $d=2$, exploit the fact that TukeyRegion package computes exact medians in all empirical cases due to dimensionality results [2209.11445].
- For robust inference under contamination, utilize depth-based central regions; in $\alpha$-symmetric models, finite-sample bounds follow the minimax rate [2512.07634].
- Outlier detection, depth-based classification, symmetry diagnostics, and nonparametric confidence regions are all supported foundationally by the halfspace median and depth central regions [1201.1171].

**Summary Table: Tukey Median Core Properties**

| Property                 | Mathematical Formulation                                    | Reference          |
|--------------------------|-------------------------------------------------------------|--------------------|
| Definition               | $D(x) = \min_{u} \frac{1}{n} \#\{i : u^\top X_i \leq u^\top x\}$ | [1207.4988][1201.1171] |
| Median region            | $C = \{x : D(x) = \max D(y)\}$                              | [1207.4988]        |
| Breakdown point          | $1/(d+1)$, up to $1/3$ asymptotically under symmetry        | [1604.07104][1207.4988] |
| Approximate algorithm    | Ball intersection, $O(d n^2)$ complexity                    | [1812.03174]       |
| Statistical consistency  | $O_p(n^{-1/2})$ uniform convergence under VC theory         | [1207.4988][2512.07634] |

In contemporary multivariate analysis, the location halfspace median stands as a foundational tool for robust, affine-invariant estimation, undergirded by rigorous depth theory, convex geometry, and scalable computation.

Source: https://www.emergentmind.com/topics/location-halfspace-median