The Influence Function of Transport-based Quantiles

Technische Zusammenfassung: Die Einflussfunktion von Transport-basierten Quantilen

Problemstellung
Transport-basierte Quantile bieten eine kanonische Erweiterung der univariaten Quantile auf multivariate Verteilungen mittels optimaler Transporttheorie (OT). Während die jüngere Literatur etabliert hat, dass diese Schätzer unter geeigneten Bedingungen hohe Breakdown-Points besitzen, wurden ihre Robustheitseigenschaften hinsichtlich infinitesimaler Störungen – spezifisch die Einflussfunktion (Influence Function, IF) – bisher nicht systematisch untersucht. Diese Arbeit adressiert die Wissenslücke im Verständnis der Stabilität der Transport-Quantil-Abbildung $Q_P$ (die OT-Abbildung, die ein festes Referenzmaß $\mu$ auf eine Zielverteilung $P$ überträgt), wenn $P$ einer Huber-Typ-Kontamination unterliegt, definiert als $P_t = (1-t)P + t\delta_{x_0}$.

Methodik
Die Autoren verwenden eine rigorose Analyse, die auf der durch partielle Differentialgleichungen (PDE) charakterisierten Struktur von optimalen Transport-Abbildungen basiert. Im Gegensatz zu klassischen Ansätzen, die auf dem impliziten Funktionalen Theorem beruhen – welches hier aufgrund der nicht-glatten Punktmassen-Störung versagt –, nutzt das Paper eine mehrstufige analytische Strategie:

1. $L^2$-Stabilitätsabschätzungen: Die Analyse beginnt mit der Etablierung von $L^2$-Schranken für die Abweichung der gestörten Transport-Abbildung $Q_{P_t}$ von der ungestörten Abbildung $Q_P$. Die Autoren beweisen, dass die quadratische $L^2(\mu)$-Dislokation als $O(t^{1+2/d})$ für $d \ge 3$ und als $O(t^2 |\log t|)$ für $d=2$ skaliert.
2. Linearisierung und $L^q$-Schranken: Um den ersten Ordnungslimit (die Einflussfunktion) zu analysieren, untersuchen die Autoren die Folge $\xi_t = (Q_{P_t} - Q_P)/t$. Sie zeigen, dass diese Folge zwar in $L^2$ divergiert, aber in einem spezifischen Bereich $q \in (1, \frac{d(d+2)}{d^2+2d-4})$ gleichmäßig beschränkt bleibt in $L^q(\Omega_\mu)$.
3. PDE-Charakterisierung: Das Limit der Potenzialfunktion, die mit $\xi_t$ assoziiert ist, wird als die eindeutige schwache Lösung einer gleichmäßig elliptischen PDE mit einem Dirac-Quellenterm und homogenen Neumann-Randbedingungen identifiziert.
4. Lokalisierung und Singularitätsanalyse: Durch das Ausschließen des Störungspunktes und die Nutzung lokaler Maximumprinzipien (via Moser-Iteration) sowie Schauder-Abschätzungen leiten die Autoren lokale $C^{2,\alpha}$-Schranken fernab der Singularität ab. Dies ermöglicht eine präzise Charakterisierung des Verhaltens in der Nähe der Singularität.

Zentrale Beiträge und Ergebnisse

• Existenz und Repräsentation: Das Paper etabliert die Existenz und Eindeutigkeit der Einflussfunktion $I(x_0; Q_P(z))$ für alle $x_0 \neq Q_P(z)$. Es wird gezeigt, dass $I(x_0; Q_P(z)) = \nabla G_{x_0}(z)$ gilt, wobei $G_{x_0}$ das Potenzial ist, welches die elliptische PDE löst:
  $$\text{div}(f_\mu [\nabla Q_P]^{-1} \nabla G_{x_0}) = \mu - \delta_{z_0}$$
  wobei $z_0 = F_P(x_0)$ die Transport-Verteilungsfunktion am Störungspunkt ist und $F_P = Q_P^{-1}$.

• Pol-artige Singularität: Eine zentrale Erkenntnis ist, dass für Dimensionen $d \ge 2$ die Einflussfunktion eine Pol-artige Singularität aufweist. Speziell, wenn der Störungspunkt $x_0$ sich dem Quantil $Q_P(z)$ nähert (äquivalenterweise, wenn $F_P(x_0) \to z$), divergiert die Norm der Einflussfunktion:
  $$\|I(x_0; Q_P(z))\| \asymp \frac{1}{\|z - F_P(x_0)\|^{d-1}}$$
  Dies steht in starkem Kontrast zur beschränkten Einflussfunktion univariater Quantile ($d=1$).

• Geometrische Charakterisierung der Präbilder: Als Nebenprodukt der Beweisstrategie charakterisiert das Paper die Menge $K_{x_0, t} = Q_{P_t}^{-1}(x_0)$, welche die Punkte im Referenzbereich darstellt, die auf die Störungsstelle abgebildet werden. Für $d \ge 3$ verhält sich diese Menge asymptotisch wie eine Kugel mit Radius $t^{1/d}$ für $t \to 0$.

• Implikationen für empirische Quantile: Die unbeschränkte Natur der Einflussfunktion impliziert, dass die Zufallsvariable $I(X; Q_P(z))$ für $X \sim P$ einen unendlichen zweiten Moment besitzt ($I \notin L^2(P)$). Folglich besitzt das empirische Transport-Quantil $Q_{P_n}(z)$ keine Standard-$\sqrt{n}$-asymptotisch lineare Repräsentation mit einer quadratisch integrierbaren Einflussfunktion.

• Vermutung über Heavy-Tailed-Fluktuationen: Basierend auf dem Tail-Verhalten der Einflussfunktion (welches als $\lambda^{-d/(d-1)}$ abfällt), vermuten die Autoren, dass das empirische Transport-Quantil mit einer nicht-parametrischen Rate $r_n$ konvergiert (wobei $r_n \sim n^{1/d}$ für $d \ge 3$) und einer stabilen, nicht-gaußschen Grenzverteilung folgt. Numerische Experimente in Dimension $d=3$ stützen diese Vermutung und zeigen, dass die empirische Verteilung mit der aus der Einflussfunktion abgeleiteten linearen Statistik übereinstimmt und schwere Tails (heavy tails) aufweist, was selbst bei großen Stichprobenumfängen von der Gaußschen Verteilung abweicht.

Bedeutung
Das Paper verändert grundlegend das Verständnis der Robustheit von Transport-basierten Quantilen. Während diese Schätzer hohe Breakdown-Points besitzen (Resistenz gegen grobe Fehler), offenbart die Entdeckung einer Pol-artigen Singularität in ihrer Einflussfunktion eine distinkte Sensitivität gegenüber „Inliern“ (Punkten nahe dem Quantil von Interesse). Dieses singuläre Verhalten erklärt, warum die Standard-Asymptotische Normalität fehlschlägt, und legt nahe, dass die Grenzverteilung empirischer Transport-Quantile schwer tail-verteilt und potenziell stabil ist. Die Arbeit liefert den notwendigen theoretischen Rahmen (via PDE-Charakterisierung), um diese nicht-standardmäßigen asymptotischen Verhaltensweisen zu analysieren und schlägt so eine Brücke zwischen der Optimalen Transporttheorie und der robusten Statistik im multivariaten Setting.