---
title: Skew-Adjusted Projection Depth (SPD)
url: https://www.emergentmind.com/topics/skew-adjusted-projection-depth-spd
type: topic
---

# Skew-Adjusted Projection Depth (SPD)

Skew-Adjusted Projection Depth (SPD) is a robust, affine-invariant data depth function for multivariate and functional data that generalizes classical projection depth by adapting to skewness in the underlying distribution. SPD delivers asymmetrically scaled outlyingness scores via a skew-adjusted univariate density estimation in projected subspaces, correcting a key limitation of symmetric scale estimators in the presence of skewed classes. The construct emerges from replacing the median absolute deviation (MAD)-based scaling in the Stahel–Donoho outlyingness (SDO) with asymmetric "fences" parameterized by the medcouple, leading to substantially improved classification performance under heavy-tailed or skewed populations [1504.01128].

## 1. Projection Depth and Skew-Adjustment

Classical Stahel–Donoho projection depth quantifies the “outlyingness” of a point $x$ with respect to a multivariate distribution $P_Y$ by considering all one-dimensional projections $v^\top x$ for $v\in\mathbb R^p$ with $\|v\|=1$, then normalizing the distance from the projected median by the MAD:
$$
\mathrm{SDO}(x; P_Y) = \sup_{\|v\| = 1} \frac{\left|v^\top x - \mathrm{med}(v^\top Y)\right|}{\mathrm{MAD}(v^\top Y)}
$$
The associated projection depth is
$$
\mathrm{PD}(x;P_Y) = \frac{1}{1+\mathrm{SDO}(x;P_Y)}
$$
When $P_Y$ is skewed, the symmetric MAD- or IQR-based scale is not optimal. Skew-adjusted projection depth (SPD) remodels outlyingness via
$$
\mathrm{AO}(x; P_Y) = \sup_{\|v\|=1} \mathrm{AO}_1(v^\top x; P_{v^\top Y})
$$
where the univariate adjusted outlyingness is
$$
\mathrm{AO}_1(z;P_Z)=
\begin{cases}
\frac{z-\mathrm{med}(Z)}{w_2(Z)-\mathrm{med}(Z)} & z>\mathrm{med}(Z)\\
\frac{\mathrm{med}(Z)-z}{\mathrm{med}(Z)-w_1(Z)} & z\le\mathrm{med}(Z)
\end{cases}
$$
The fences $(w_1, w_2)$ are skew-adapted:
$$
w_1 = Q_1 - 1.5 \exp(-4\,\mathrm{MC}(Z))\,\mathrm{IQR}(Z)
$$
$$
w_2 = Q_3 + 1.5 \exp(+3\,\mathrm{MC}(Z))\,\mathrm{IQR}(Z)
$$
where $Q_1, Q_3$ are quartiles, $\mathrm{MC}(Z)$ is the medcouple statistic—the median of all pairwise slopes in $Z$—serving as a robust skewness measure, and $\mathrm{IQR}=Q_3 - Q_1$. If $\mathrm{MC}(Z)<0$ one replaces $(Z,z)$ by $(-Z,-z)$ to work in the right-skewed regime. Then
$$
\mathrm{SPD}(x;P_Y) = \frac{1}{1+\mathrm{AO}(x;P_Y)}
$$
SPD thus replaces symmetric scaling with one-sided, skew-adaptive scaling, allowing the depth to respond asymmetrically to heavy-tailed data [1504.01128].

## 2. Efficient Estimation and Computational Complexity

For finite-sample estimation, consider $X_1,\dots,X_n\in\mathbb R^p$ and query $x$. The affine-invariant procedure:

1. Generate $m$ directions by sampling $p$ points and choosing normals to their affine hulls, repeated to obtain $\{v_j\}_{j=1}^m$.
2. For each $v_j$:
   - Project: $Z_i \gets v_j^\top X_i$, $z \gets v_j^\top x$.
   - Compute $\mathrm{med}(Z)$, $Q_1(Z)$, $Q_3(Z)$, $\mathrm{IQR}(Z)$, $\mathrm{MC}(Z)$.
   - Set $w_1$, $w_2$ as above.
   - Compute $\mathrm{AO}_j = \mathrm{AO}_1(z;P_Z)$.
3. Set $\mathrm{AO}(x) = \max_{j} \mathrm{AO}_j$.
4. Return $\mathrm{SPD}(x) = 1/(1+\mathrm{AO}(x))$.

The dominant cost is $O(m n p)$: with $m\approx 250p$, the effort per query $x$ is $O(p^2 n)$, accounting for projection, order statistics computation ($O(n\log n)$ for median/quartiles, $O(n\log n)$ for medcouple using available optimized algorithms) [1504.01128].

## 3. Theoretical Properties of SPD

SPD possesses the following properties:

- **Affine invariance**: SPD is preserved under all invertible linear transformations, as each projection and all order statistics are affine-invariant univariate functionals. The supremum over directions inherits this property.
- **Robustness**: SPD is bounded within $(0,1]$, with its denominator (AO) growing roughly linearly in $\lVert x\rVert$. As $x\to\infty$, $\mathrm{SPD}(x)\to 0$. High breakdown point ($\approx 50\%$) is inherited from the robust medcouple and quartiles. The bounded influence of medcouple and order statistics ensures resistance to outliers.
- **Sensitivity to skewness**: Exponentially adjusted fences ($\exp(\pm 3,\pm 4\mathrm{MC})$) cause the upper and lower fences to expand differently, adapting the depth to asymmetry and long tails. This mechanism prevents heavy-side values in skewed distributions from being spuriously classified as extreme.
- **Depth-like behavior**: $\mathrm{SPD}(x;P_Y)$ attains maximum 1 at points of zero adjusted outlyingness (typically the multivariate “median”) and decays to zero as $x$ departs from the bulk of $P_Y$ [1504.01128].

## 4. Algorithmic and Implementation Choices

Several practical choices influence SPD's estimation:

| Parameter                   | Default/Recommended          | Trade-offs and Context                 |
|-----------------------------|------------------------------|----------------------------------------|
| Number of directions $m$    | $250p$                       | Fewer $m$ speeds up computation; more directions yield a closer approximation to the supremum. |
| Direction generation        | Normals to hyperplanes of $p$ samples; random Gaussians also viable | Normals ensure exploration of convex hull facets; random often needs larger $m$.                 |
| Medcouple computation       | $O(n\log n)$ algorithm, can subsample for large $n$ | Subsampling provides scalable proxy to exact $\mathrm{MC}$.             |
| Fence multipliers           | $1.5$ (boxplot literature), exponents $\pm3$, $\pm4$ | Default values robust in experiments; tuning impacts fence aggressiveness.                        |

These choices underpin the scalability and practical adaptability of SPD to high dimensions and large sample sizes [1504.01128].

## 5. Comparative Performance in Classification

SPD is particularly effective in classification tasks where class distributions exhibit skewness. In Simulation Setting 2 of [1504.01128], with one Normal and one highly right-skewed exponential class (both 6-variate), the DistSpace transform with $k$-nearest neighbor classification was compared across three depth-based outlyingness measures:

| Method                        | Misclassification Rate (%) |
|-------------------------------|--------------------------|
| DistSpace + kNN (SPD)         | $\approx 2.1$            |
| DistSpace + kNN (PD)          | $\approx 5.4$            |
| DistSpace + kNN (SDO)         | $\approx 4.8$            |

SPD halved the error of PD. This enhancement is attributed to the skew-adaptive AO$_1$ fences, which avoid excessive contraction of the heavy-tailed region, ensuring effective discrimination under skew and heavy tails. The performance gain demonstrates that SPD’s asymmetric treatment of outlyingness is fundamental for depth-based learning in non-symmetric populations [1504.01128].

## 6. Summary and Significance

Skew-adjusted projection depth merges projection-based multivariate depth geometry with univariate, medcouple-driven skew adjustment, providing a robust, affine-invariant, and computationally feasible depth measure. SPD is explicitly constructed for scenarios where classical symmetric data depth fails due to skewness or heavy-tailedness. It is straightforward to estimate ($O(p^2 n)$ per query), deploys robust statistics at every stage, and can significantly reduce misclassification error in statistical learning pipelines. SPD exemplifies how robust univariate summaries combined with directional analysis yield substantial practical and theoretical gains in multivariate and functional data analysis [1504.01128].

Source: https://www.emergentmind.com/topics/skew-adjusted-projection-depth-spd