Sensitivity Analysis on the Sphere and a Spherical ANOVA Decomposition

技術概要：球面上の感度分析と球面 ANOVA 分解

問題定義
分散分析（ANOVA）分解は、多変量関数の構造を理解するための標準的な枠組みであり、特に感度分析と次元削減において重要な役割を果たす。古典的な設定では、定義域がテンソル積空間（例：$[0,1]^d$）である場合、この分解は自然な積構造に依存しており、これにより異なる変数部分集合に対応する成分関数の直交性が保証される。しかし、この枠組みを高次元球面 $S^d$ に拡張することは、理論的に大きな課題を伴う。球面にはテンソル積構造が存在せず、制約条件 $\|x\|=1$ が入力変数間の本質的な依存関係を誘起する。その結果、独立変数上で積分を行うことに依存する古典的な ANOVA 成分の定義は、球面上では一意性や直交性を持つ項を生成することができない。さらに、従属変数に対する既存の一般化は、しばしばグリッド閉領域や特定の座標変換を必要とし、これにより元の物理的変数の解釈可能性が損なわれる。

手法
本論文は、一様分布の下で関数 $f: S^d \to \mathbb{R}$ 向けに特化した新しい球面 ANOVA 分解を提案する。この手法は、以下の 3 つの主要なメカニズムを通じて球面の幾何学的制約に対処する。

1. パリティ分解: 古典的な場合とは異なり、著者はパリティベクトル $\xi \in \{0, 1\}^{d+1}$ を導入し、項を偶数部と奇数部に分解する。これは、球面上における変数の代数的依存関係（$x_i^2$ 項）が、パリティを区別しなければ解決できない冗長性を生み出すため必要である。分解は以下の形式をとる：
   $$f(x) = \sum_{u \in \mathcal{D}} \sum_{\xi \in \Upsilon_u} f_{u,\xi}(x_u)$$
   ここで、$u$ はインデックスの部分集合であり、$\xi$ は各変数に対する項のパリティを示す。

2. 反復射影演算子: 著者は、変数 $u$ の固定値に対応する球面の「ファイバー」（低次元球面）上で関数を積分する射影演算子 $P_u$ を定義する。古典的なテンソル積射影とは異なり、この演算子は因数分解しない。ANOVA 項は、古典的なホエディング分解と同様だが球面幾何学に適応されたこれらの射影を用いて反復的に定義される。

3. 冗長性の排除: 関係式 $\sum x_i^2 = 1$ により、変数 $x_{d+1}$ に依存する偶数項は、残りの変数の関数として表現できる。一意かつ非冗長な分解を確保するため、著者は最高インデックス変数が偶数パリティで現れる特定の項を省略することを提案する。これにより項の総数が削減され、基底における線形依存性が防止される。

4. 直交基底の構築: 分解を数値的に実装するために、本論文は項 $f_{u,\xi}$ が定義される低次元球 $B^{|u|}$ に対して縮小された直交基底を構築する。これらの基底は、1 次元項に対して Gegenbauer 多項式を、高次元項に対して Jacobi 多項式と球面調和関数を組み合わせて用い、球面上で誘起された重み関数に関する直交性を保証する。

5. 数値近似と感度分析: 本分解は、これらの基底関数の切り捨てられた級数を用いて近似される。係数は、ランダムなサンプリングノード上で最小二乗問題（LSQR アルゴリズムを使用）を解くことで決定される。これらの項の分散に基づき、著者は球面ソボル指数を定義して変数の重要性を定量化する。厳密な直交性が偶数項に対して失われることを認めつつ、得られる指数が変数の影響に対する安定した尺度を提供することを示す。

主な貢献

• 球面 ANOVA の定義: 本論文は、座標系変換（例：球座標）を必要とせず、物理的意味を曖昧にすることなく元の変数の解釈可能性を保持する、球面専用の最初の ANOVA 分解を確立する。
• パリティに基づく構造: 古典的 ANOVA には存在しないが、球面の幾何学を扱うために不可欠な、パリティベクトル $\xi$ による偶数項と奇数項の区別の必要性を導入する。
• 理論的性質: 著者は古典的性質の転用可能性を分析し、異なるパリティの項間では直交性が成り立つが、変数の依存性により同一パリティの項間では厳密には成り立たないことを示す。また、項が満たすべき特定の積分条件を提供する。
• 縮小基底と冗長性の処理: 球面の制約から生じる冗長な項を特定し省略する方法を定義し、線形独立な基底関数の辞書を作成する。
• 感度指数: 多様体上に制約された従属入力に対する大域感度分析の枠組みを提供する、球面領域に適応されたソボル指数を定義する。

結果
既知の相互作用構造を持つテスト関数（多項式、混合パリティ関数、および Ishigami 関数を含む）を用いて、10 次元球面（$S^{10}$）上で数値実験が行われた。

• 提案された手法は、テスト関数における既知の相互作用に対応する非ゼロ ANOVA 項を正常に特定した。
• 完全に奇数項で構成される関数については、分解は期待にほぼ一致し、非ゼロ項を正確に検出した。
• 偶数項を含む関数については、手法は主要な相互作用変数を正しく特定したが、本質的な依存性により、分散がより高次の項に分配された。著者は、これらの「誤解を招きうる」項の大きさは主要項と比較して小さいと指摘している。
• 項の次数を増加させる反復アプローチ（1 次、次に 2 次）は、すべての項を同時に解くことと同等の結果をもたらしており、この文脈における LSQR ソルバの安定性を検証した。

意義と主張
本論文は、提案された球面 ANOVA 分解が、球面上の関数に対する「意味のあるかつ解釈可能な表現」を提供し、専用の分解が存在しなかった文献のギャップを埋めると主張している。その意義は以下の点にある。

• 幾何学を考慮した分析: 古典的手法が解釈可能性を失わずに捉えられない、球面の本質的な曲率と変数の結合を尊重する。
• 次元削減: 球面上の高次元関数はしばしば低次相互作用によって支配されていることを特定することで、近似や積分における次元の呪いを緩和する道筋を提供する。
• 感度分析: 依存入力を持つ設定における変数の重要性の定量化を可能にし、正規化された特徴ベクトルを取り扱う分野（地球物理学、量子力学など）における特徴選択やノイズ低減のツールとして機能する。

著者は謙虚なトーンを維持し、この分解がすべての応用ベースの質問（特定の関数空間の特性記述など）に答えるものではなく、また偶数項に対する厳密な直交性の欠如により、ソボル指数は正確な分散分割ではなく相対的重要性の尺度として解釈されるべきであることを認めている。この研究は、球面領域上の関数を分析するための基礎的な理論的枠組みおよび実用的なアルゴリズムとして提示されている。