The Influence Function of Transport-based Quantiles

Technische Samenvatting: De invloedsfunctie van transport-gebaseerde kwantielen

Probleemstelling
Transport-gebaseerde kwantielen bieden een canonieke uitbreiding van univariate kwantielen naar multivariate verdelingen via optimale transport (OT). Hoewel recente literatuur heeft vastgesteld dat deze schatters onder geschikte omstandigheden beschikt over hoge breakdown points, zijn hun robuustheidseigenschappen met betrekking tot infinitesimale perturbaties — specifgens de invloedsfunctie (IF) — niet systematisch onderzocht. Dit artikel adresseert het gat in het begrip van de stabiliteit van de transport-kwantielafbeelding $Q_P$ (de OT-afbeelding die een vaste referentiemaat $\mu$ naar een doelverdeling $P$ pusht) wanneer $P$ wordt onderworpen aan Huber-type contaminatie, gedefinieerd als $P_t = (1-t)P + t\delta_{x_0}$.

Methodologie
De auteurs hanteren een rigoureuze analyse die geworteld is in de partiële differentiaalvergelijking (PDE) karakterisering van optimale transport-afbeeldingen. In tegenstelling tot klassieke benaderingen die vertrouwen op de impliciete functiestelling, die hier faalt vanwege de niet-gladde aard van de puntmassa-perturbatie, maakt het artikel gebruik van een meerstaps analytische strategie:

1. $L^2$ Stabiliteitsschattingen: De analyse begint met het vaststellen van $L^2$-bounds op de afwijking van de geperturbeerde transport-afbeelding $Q_{P_t}$ van de onverstoorde afbeelding $Q_P$. De auteurs bewijzen dat de gekwadrateerde $L^2(\mu)$-verplaatsing schaalt als $O(t^{1+2/d})$ voor $d \ge 3$ en $O(t^2 |\log t|)$ voor $d=2$.
2. Linearisatie en $L^q$ Bounds: Om de eerste-orde limiet (de invloedsfunctie) te analyseren, bestuderen de auteurs de sequentie $\xi_t = (Q_{P_t} - Q_P)/t$. Zij tonen aan dat hoewel deze sequentie divergeert in $L^2$, deze wel uniform begrensd blijft in $L^q(\Omega_\mu)$ voor een specifieke reeks $q \in (1, \frac{d(d+2)}{d^2+2d-4})$.
3. PDE Karakterisering: De limiet van de potentiaalfunctie geassocieerd met $\xi_t$ wordt geïdentificeerd als de unieke zwakke oplossing van een uniform elliptische PDE met een Dirac-bronterm en homogene Neumann-randvoorwaarden.
4. Lokalisatie en Singulariteitsanalyse: Door de perturbatiepunt uit te sluiten en gebruik te maken van lokale maximumprincipes (via Moser-iteratie) en Schauder-schattingen, leiden de auteurs lokale $C^{2,\alpha}$ bounds af weg van de singulariteit. Dit maakt een precieze karakterisering van het gedrag nabij de singulariteit mogelijk.

Kernbijdragen en Resultaten

• Bestaan en Representatie: Het artikel stelt het bestaan en de uniciteit van de invloedsfunctie $I(x_0; Q_P(z))$ vast voor alle $x_0 \neq Q_P(z)$. Er wordt aangetoond dat $I(x_0; Q_P(z)) = \nabla G_{x_0}(z)$, waarbij $G_{x_0}$ de potentiaal is die de elliptische PDE oplost:
  $$\text{div}(f_\mu [\nabla Q_P]^{-1} \nabla G_{x_0}) = \mu - \delta_{z_0}$$
  waarbij $z_0 = F_P(x_0)$ de transport-verdelingsfunctie is geëvalueerd in het perturbatiepunt, en $F_P = Q_P^{-1}$.

• Pool-type Singulariteit: Een centrale bevinding is dat voor dimensies $d \ge 2$, de invloedsfunctie een pool-type singulariteit vertoont. Specifiek, naarmate het perturbatiepunt $x_0$ de kwantiel $Q_P(z)$ nadert (oftewel, naarmate $F_P(x_0) \to z$), divergeert de norm van de invloedsfunctie:
  $$\|I(x_0; Q_P(z))\| \asymp \frac{1}{\|z - F_P(x_0)\|^{d-1}}$$
  Dit staat in scherp contrast met de begrensde invloedsfunctie van univariate kwantielen ($d=1$).

• Geometrische Karakterisering van Pre-images: Als bijproduct van de bewijsstrategie karakteriseert het artikel de verzameling $K_{x_0, t} = Q_{P_t}^{-1}(x_0)$, die de punten in het referentiedomein vertegenwoordigt die naar de perturbatielocatie worden afgebeeld. Voor $d \ge 3$ gedraagt deze verzameling zich asymptotisch als een bol met straal $t^{1/d}$ wanneer $t \to 0$.

• Implicaties voor Empirische Kwantielen: De onbegrensde aard van de invloedsfunctie impliceert dat de willekeurige variabele $I(X; Q_P(z))$ voor $X \sim P$ een oneindig tweede moment heeft ($I \notin L^2(P)$). Bijgevolg bezit het empirische transport-kwantiel $Q_{P_n}(z)$ geen standaard $\sqrt{n}$-asymptotisch lineaire representatie met een vierkant-integreerbare invloedsfunctie.

• Conjectuur over Heavy-Tailed Fluctuaties: Gebaseerd op de staartgedrag van de invloedsfunctie (die vervalt als $\lambda^{-d/(d-1)}$), conjectureren de auteurs dat het empirische transport-kwantiel convergeert met een niet-parametrische snelheid $r_n$ (waarbij $r_n \sim n^{1/d}$ voor $d \ge 3$) en een stabiele, niet-Gaussische limietverdeling volgt. Numerieke experimenten in dimensie $d=3$ ondersteunen deze conjectuur, waarbij ze laten zien dat de empirische verdeling overeenkomt met de lineaire statistiek afgeleid van de invloedsfunctie en zware staarten vertoont, wat afwijkt van het Gaussische gedrag, zelfs bij grote steekproefomvang.

Significantie
Dit artikel verandert fundamenteel het begrip van de robuustheid van transport-gebaseerde kwantielen. Hoewel deze schatters over hoge breakdown points beschikken (weerstand tegen grove fouten), onthult de ontdekking van een pool-type singulariteit in hun invloedsfunctie een duidelijke gevoeligheid voor "inliers" (punten dicht bij het kwantiel van belang). Dit singular gedrag verklaart waarom de standaard asymptotische normaliteit faalt en suggereert dat de limietverdeling van empirische transport-kwantielen heavy-tailed is en potentieel stabiel. Het werk biedt het noodzakelijke theoretische kader (via PDE-karakterisering) om deze niet-standaard asymptotische gedragingen te analyseren, waardoor de optimal transport-theorie wordt verbonden met robuuste statistiek in de multivariate setting.