The Influence Function of Transport-based Quantiles

技術要約：輸送ベースの分位数の影響関数

問題提起
輸送ベースの分位数（Transport-based quantiles）は、最適輸送（OT）を通じて、単変量分位数を多変量分布へと正準に拡張したものである。近年の文献では、適切な条件下においてこれらの推定量の高分解点（high breakdown points）の存在が確立されているが、微小な摂動、具体的には影響関数（influence function; IF）に関する頑健性については、体系的な調査が行われてこなかった。本論文は、ターゲット分布 $P$ がハバ型（Huber-type）の汚染 $P_t = (1-t)P + t\delta_{x_0}$ を受ける際の、輸送写像 $Q_P$（固定された参照測度 $\mu$ をターゲット分布 $P$ へ押し出すOT写像）の安定性に関する理解の空白を埋めるものである。

手法
著者らは、最適輸送写像の偏微分方程式（PDE）による特徴付けに基づいた厳密な解析を用いている。ここでは、点質量による摂動が非平滑であるために、古典的なアプローチである陰関数定理は適用できないため、以下の多段階的な解析戦略を採用している。

1. $L^2$ 安定性評価: 解析は、摂動された輸送写像 $Q_{P_t}$ と未摂動の写像 $Q_P$ との間の $L^2(\mu)$ 偏差に関する $L^2$ 境界を確立することから始まる。著者らは、$d \ge 3$ の場合で平方 $L^2(\mu)$ 変位が $O(t^{1+2/d})$ となり、$d=2$ の場合には $O(t^2 |\log t|)$ となることを証明している。
2. 線形化と $L^q$ 境界: 1次極限（影響関数）を解析するために、著者らは数列 $\xi_t = (Q_{P_t} - Q_P)/t$ を研究する。この数列は $L^2$ において発散するものの、特定の範囲 $q \in (1, \frac{d(d+2)}{d^2+2d-4})$ において $L^q(\Omega_\mu)$ で一様に有界であることを示している。
3. PDE による特徴付け: $\xi_t$ に付随するポテンシャル関数の極限が、ディラック源項と均質ノイマン境界条件を持つ一様楕円型PDEの唯一の弱解であることが特定される。
4. 局所化と特異性解析: 摂動点を除外し、局所最大原理（Moser 反復法による）および Schauder 推定を利用することで、特異点の近傍における $C^{2,\alpha}$ 境界を導出する。これにより、特異点付近での挙動の精密な特徴付けが可能となる。

主要な貢献と結果

• 存在性と表現: 本論文は、すべての $x_0 \neq Q_P(z)$ に対して影響関数 $I(x_0; Q_P(z))$ の存在と一意性を確立している。$I(x_0; Q_P(z)) = \nabla G_{x_0}(z)$ であることが示されており、ここで $G_{x_0}$ は以下の楕円型PDEを解くポテンシャルである：
  $$\text{div}(f_\mu [\nabla Q_P]^{-1} \nabla G_{x_0}) = \mu - \delta_{z_0}$$
  ただし、$z_0 = F_P(x_0)$ は摂動点における輸送分布関数（$F_P = Q_P^{-1}$）の評価値である。

• 極型の特異性（Pole-Type Singularity）: 中心的な発見は、$d \ge 2$ の場合、影響関数が極型の特異性を示すことである。具体的には、摂動点 $x_0$ が分位数 $Q_P(z)$ に近づくとき（あるいは $F_P(x_0) \to z$ となるとき）、影響関数のノルムは発散する：
  $$\|I(x_0; Q_P(z))\| \asymp \frac{1}{\|z - F_P(x_0)\|^{d-1}}$$
  これは、単変量分位数（$d=1$）の有界な影響関数とは対照的である。

• 逆像の幾何学的特徴付け: 証明戦略の副産物として、論文は $K_{x_0, t} = Q_{P_t}^{-1}(x_0)$（摂動位置に写される参照領域内の点）の集合を特徴付けている。$d \ge 3$ の場合、この集合は $t \to 0$ のとき、半径 $t^{1/d}$ の球体として漸近的に振る舞う。

• 経験的分位数への含意: 影響関数が非有界であることは、確率変数 $I(X; Q_P(z))$ （ただし $X \sim P$）が無限の2次モーメントを持つ（$I \notin L^2(P)$）ことを意味する。その結果、経験的輸送分位数 $Q_{P_n}(z)$ は、平方可積分な影響関数を持つ標準的な $\sqrt{n}$-漸近線形表現を持たない。

• 重い裾を持つゆらぎに関する予想: 影響関数の裾の挙動（$\lambda^{-d/(d-1)}$ として減衰する）に基づき、著者らは経験的輸送分位数が非パラメトリックな速度 $r_n$（$d \ge 3$ の場合 $r_n \sim n^{1/d}$）で収束し、安定な非ガウス限界分布に従うという予想を立てている。次元 $d=3$ における数値実験は、経験分布が影響関数から導かれる線形統計量と一致し、ガウス分布から逸脱した重い裾を持つことを示しており、この予想を支持している。

意義
本論文は、輸送ベースの分位数の頑健性に関する理解を根本的に変えるものである。これらの推定量は高い分解点（粗大な誤差に対する耐性）を持つ一方で、影響関数における極型の特異性の発見は、「インライア（inliers）」（関心のある分位数に近い点）に対する独特の感受性を明らかにしている。この特異な振る舞いは、なぜ標準的な漸近正規性が成立しないのかを説明しており、経験的輸送分位数の極限分布が重い裾を持ち、潜在的に安定分布であることを示唆している。本研究は、最適輸送理論とロバスト統計学を多変量設定で結びつけるために必要な理論的枠組み（PDEによる特徴付けを通じて）を提供している。