A Riemannian View on Active Subspaces

技術要約：能動的部分空間（Active Subspaces）に関するリーマン幾何学的視点

問題提起
能動的部分空間は、高次元のユークリッド・ドメインにおける次元削減のための、固有値に基づいた原理的な手法を提供する。これらは、スカラー量 $f: \mathbb{R}^n \to \mathbb{R}$ が平均的に最も大きく変化する低次元の部分空間を特定し、その直交補空間では値がほぼ一定であることを示す。しかし、現代の多くのアプリケーションでは、パラメータがユークリッド空間ではなく、リーマン多様体（例：統計的形状解析におけるプレシェイプ空間である超球面）に自然に存在する。標準的な能動的部分空間の定義は、ユークリッド空間の標準的な同型性に依存しており、これにより勾配ベクトルの大域的な平均化が可能となっている。リーマン多様体上の曲がった多様体では、異なる点における接ベクトルは異なる接空間に属するため、比較のためのメカニズムなしには直接的な大域的平均を定義することができない。本論文は、平均二乗方向微分としての固有値の解釈性を維持しつつ、リーマン多様体上の滑らかな関数に対して能動的部分空間の枠組みを拡張するという課題に取り組んでいる。

手法
本論文は、勾配の外積を一般化することにより、リーマン多様体 $(M, g)$ 上の能動的部分空間の固有的な定式化を開発している。核心となる手法は、以下の3つの主要なステップで構成される。

1. 平行移動による固有の平均化： 勾配のグローバルな平均を定義するために、著者らは各リーマン勾配 $\nabla_M f(p)$ を、点 $p$ と中心点 $p_0$ を結ぶ一意な測地線に沿って、局所的な接空間 $T_p M$ から中心的な接空間 $T_{p_0} M$ へ輸送する。これは平行移動演算子 $P_p$ を用いて達成される。結果として得られる「勾配のリーマン平均テンソル積」は次のように定義される：
   $$G_0 = \int_X P_p [\nabla_M f(p)] \otimes P_p [\nabla_M f(p)] \, d\mu(p)$$
   ここで $\mu$ は多様体上の確率測度（通常は測地球内に支持される）である。この構成により、単一のベクトル空間 $T_{p_0} M$ 上の対称正定値テンソルが得られる。

2. 能動的多様体測地線 (Active Manifold-Geodesics, AMG)： $G_0$ の固有値分解により、固有値 $\lambda_i$ と $T_{p_0} M$ 内の固有ベクトル $w_i$ が得られる。著者らは「能動的多様体測地線」を、指数写像の下でのこれらの固有ベクトルによって張られる直線像として定義する：$W_{1,i} = \{ \exp_{p_0}(t w_i) \}$。これらの測地線は、多様体上の能動的な方向として機能する。

3. 固有的アプローチと外延的アプローチの比較： 本論文は、多様体がより高次元のユークリッド空間 $\mathbb{R}^m$ に埋め込まれている場合の外延的（extrinsic）なアプローチと、この固有的（intrinsic）なアプローチを対比させる。外延的な手法は、周囲の勾配を平均化し、それを中心の接空間に投影する。著者らは、投影の部分等長性（partial isometry）により、外延的なスペクトルが固有のスペクトルの「混同された（conflated）」バージョンになることを示し、両者の間の正確な関係を導出している。

主な貢献
本論文には、主に4つの貢献がある。

• 固有の一般化： 平行移動されたフレーム場に沿った平均二乗方向微分としての解釈を厳密に与える座標不変なテンソルとして、平行移動勾配外積（式12）を導入した（補題1）。
• 能動的多様体測地線 (AMG)： 順序付けられた測地線方向の集合を定義した（定義3）。全非線形な次元削減（例：能動的多様体やニューラルネットワーク）とは異なり、AMGは探索範囲を単一の接空間における直線の指数像に制限することで、低い複雑性と識別性を維持する。
• 定量的における固有・外延的一致： 固有のアプローチと外延的なアプローチは根本的に異なるものの、サンプリングされた測地球の半径 $R$ に対して、2次オーダー ($O(R^2)$) で一致することを確立した。具体的には、固有値は $O(R^2)$ で一致し、支配的な固有空間は $O(R^2/\eta)$ （$\eta$ はスペクトルギャップ）で一致する（命題3）。決定的なのは、支配的な周囲の固有ベクトルを投影するという一般的な外延的慣習が、その方向が投影の核（null space）にある場合に失敗し得ることを示し、固有的手法が堅牢であることを示した点である。
• 形状解析への適用： 理論は、統計的形状解析におけるハイパー球面であるプレシェイプ空間によって明確に動機付けられ、適用されている。その結果は、形状関数における応答駆動型の次元削減のための直接的な枠組みを提供する。

結果と数値的例証
著者らは、2次元球面 ($S^2$) 上の数値例を用いて定式化を検証している。

• リッジの回復： 球面に制限された周囲の線形関数（リッジ関数）に対して、固有のAMGは能動的な測地線方向を正常に特定する。論文では、リッジ構造が正規座標において、曲率に制限された二次的な収束率 ($O(R^2)$) で誤差が減少しながら回復されることが示されている。これは命題3と一致する。
• 退化したケース： 中心点が周囲のリッジの勾配方向と一致するシナリオを検討する。この場合、外延的な手法（支配的な周囲の固有ベクトルを投影する方法）は失敗し、自明なゼロベクトルを返す。対照的に、固有の手法は単一の方向が支配的ではないことを正しく識別し、2つの等しい固有値を返し、「投影された優位性」という失敗モードを回避する。
• 非リッジ関数： リッジではない二次関数に対して、固有の手法は関数の変動のみに基づいて方向の忠実なランキングを提供するが、外延的な手法は、この変動を埋め込みの幾何学と絡み合わせてしまう。

意義と主張
本論文は、固有の定式化が、特定の埋め込みに依存しない、ドメインの幾何学をより忠実に表現することを提供すると主張している。外延的な手法は計算上は便利であるが、活動を中心点を超えて拡張する場合に、重要度の順序を歪める可能性がある「部分等長性」を導入することを論じている。

著者らは、提案された手法が**測地線的に局所的（geodesically local）**であることを強調している。固有値が不活性な部分多様体の指標として解釈されるのは、中心点 $p_0$ を中心とする正規近傍（測地球）内においてのみ厳密に有効である。これをグローバルなドメインに拡張するには、局所的な解析のアトラス（地図）またはより複雑な葉層構造が必要となるが、これは今後の課題として留めている。本論文は、任意の多様体に対するグローバルな次元削減問題を解決すると主張するものではなく、滑らかな関数上のリーマン多様体に対する、説明可能で厳密かつ局所的に最適な、能動的部分空間の拡張を提示するものである。これは、形状解析におけるハイパー球面データへの即時的な適用可能性を持つ。