The Influence Function of Transport-based Quantiles

Resumo Técnico: A Função de Influência de Quantis Baseados em Transporte

Enunciado do Problema
Os quantis baseados em transporte oferecem uma extensão canônica dos quantis univariados para distribuições multivariadas via transporte ótimo (OT). Embora a literatura recente tenha estabelecido a existência de pontos de quebra (breakdown points) elevados para esses estimadores sob condições adequadas, suas propriedades de robustez em relação a perturbações infinitesimais — especificamente a função de influência (IF) — não foram sistematicamente investigadas. Este artigo aborda a lacuna na compreensão da estabilidade do mapa de quantil de transporte $Q_P$ (o mapa de OT que leva uma medida de referência fixa $\mu$ à distribuição alvo $P$) quando $P$ é submetida a uma contaminação do tipo Huber, definida como $P_t = (1-t)P + t\delta_{x_0}$.

Metodologia
Os autores empregam uma análise rigorosa fundamentada na caracterização por equações diferenciais parciais (EDPs) dos mapas de transporte ótimo. Diferente de abordagens clássicas que dependem do teorema da função implícita, que falha aqui devido à natureza não suave da perturbação de massa pontual, o artigo utiliza uma estratégia analítica de múltiplas etapas:

1. Estimativas de Estabilidade $L^2$: A análise começa estabelecendo limites $L^2$ sobre o desvio do mapa de transporte perturbado $Q_{P_t}$ em relação ao mapa não perturbado $Q_P$. Os autores provam que o deslocamento quadrático $L^2(\mu)$ escala como $O(t^{1+2/d})$ para $d \ge 3$ e $O(t^2 |\log t|)$ para $d=2$.
2. Linearização e Limites $L^q$: Para analisar o limite de primeira ordem (a função de influência), os autores estudam a sequência $\xi_t = (Q_{P_t} - Q_P)/t$. Eles demonstram que, embora essa sequência divirja em $L^2$, ela permanece uniformemente limitada em $L^q(\Omega_\mu)$ para um intervalo específico de $q \in (1, \frac{d(d+2)}{d^2+2d-4})$.
3. Caracterização por EDP: O limite da função potencial associada a $\xi_t$ é identificado como a única solução fraca de uma EDP uniformemente elíptica com um termo de fonte de Dirac e condições de Neumann homogêneas.
4. Localização e Análise de Singularidade: Ao excluir o ponto de perturbação e utilizar princípios máximos locais (via iteração de Moser) e estimativas de Schauder, os autores derivam limites $C^{2,\alpha}$ locais longe da singularidade. Isso permite uma caracterização precisa do comportamento próximo à singularidade.

Principais Contribuições e Resultados

• Existência e Representação: O artigo estabelece a existência e unicidade da função de influência $I(x_0; Q_P(z))$ para todo $x_0 \neq Q_P(z)$. Demonstra-se que $I(x_0; Q_P(z)) = \nabla G_{x_0}(z)$, onde $G_{x_0}$ é o potencial que resolve a EDP elíptica:
  $$\text{div}(f_\mu [\nabla Q_P]^{-1} \nabla G_{x_0}) = \mu - \delta_{z_0}$$
  onde $z_0 = F_P(x_0)$ é a função de distribuição de transporte avaliada no ponto de perturbação, e $F_P = Q_P^{-1}$.

• Singularidade do Tipo Polo: Uma descoberta central é que, para dimensões $d \ge 2$, a função de influência exibe uma singularidade do tipo polo. Especificamente, conforme o ponto de perturbação $x_0$ se aproxima do quantil $Q_P(z)$ (equivale a $F_P(x_0) \to z$), a norma da função de influência diverge:
  $$\|I(x_0; Q_P(z))\| \asymp \frac{1}{\|z - F_P(x_0)\|^{d-1}}$$
  Isso contrasta fortemente com a função de influência limitada dos quantis univariados ($d=1$).

• Caracterização Geométrica de Pré-imagens: Como um subproduto da estratégia de prova, o artigo caracteriza o conjunto $K_{x_0, t} = Q_{P_t}^{-1}(x_0)$, que representa os pontos no domínio de referência mapeados para a localização da perturbação. Para $d \ge 3$, este conjunto comporta-se assintoticamente como uma bola de raio $t^{1/d}$ conforme $t \to 0$.

• Implicações para Quantis Empíricos: A natureza não limitada da função de influência implica que a variável aleatória $I(X; Q_P(z))$ para $X \sim P$ possui um segundo momento infinito ($I \notin L^2(P)$). Consequentemente, o quantil de transporte empírico $Q_{P_n}(z)$ não admite uma representação assintoticamente linear padrão de $\sqrt{n}$ com uma função de influência integrável ao quadrado.

• Conjectura sobre Flutuações de Cauda Pesada: Com base no comportamento de cauda da função de influência (que decai como $\lambda^{-d/(d-1)}$), os autores conjecturam que o quantil empírico converge a uma taxa não paramétrica $r_n$ (onde $r_n \sim n^{1/d}$ para $d \ge 3$) e segue uma distribuição limite estável e não gaussiana. Experimentos numéricos em dimensão $d=3$ apoiam essa conjectura, mostrando que a distribuição empírica se alinha com a estatística linear derivada da função de influência e exibe caudas pesadas, desviando-se do comportamento gaussiano mesmo em grandes tamanhos de amostra.

Significância
Este trabalho altera fundamentalmente a compreensão da robustez dos quantis baseados em transporte. Embora esses estimadores possuam pontos de quebra elevados (resistência a erros grosseiros), a descoberta de uma singularidade do tipo polo em sua função de influência revela uma sensibilidade distinta a "inliers" (pontos próximos ao quantil de interesse). Esse comportamento singular explica por que a normalidade assintótica padrão falha e sugere que a distribuição limite dos quantis de transporte empíricos é de cauda pesada e potencialmente estável. O trabalho fornece o arcabouço teórico necessário (via caracterização por EDP) para analisar esses comportamentos assintóticos não padronizados, unindo a teoria do transporte ótimo à estatística robusta no cenário multivariado.