[論文レビュー] Fréchet Means for Distributions of Persistence diagrams

Katharine Turner, Yuriy Mileyko+2 more|arXiv (Cornell University)|Jun 13, 2012

Topological and Geometric Data Analysis参考文献 19被引用数 11

ひとこと要約

本稿では、$L^2$-Wasserstein空間上の勾配降下法を用いて、パーシステンス図の経験的分布のフレシェ平均を計算するためのアルゴリズムを提案する。収束性を証明し、図が離散測度から独立同分布に従って抽出される場合に、経験的平均が弱大数の法則に従うことを確立することで、トポロジカルデータ解析における統計的推論を可能にする。

ABSTRACT

Given a distribution $ρ$ on persistence diagrams and observations $X_1,...X_n \stackrel{iid}{\sim} ρ$ we introduce an algorithm in this paper that estimates a Fréchet mean from the set of diagrams $X_1,...X_n$. If the underlying measure $ρ$ is a combination of Dirac masses $ρ= \frac{1}{m} \sum_{i=1}^m δ_{Z_i}$ then we prove the algorithm converges to a local minimum and a law of large numbers result for a Fréchet mean computed by the algorithm given observations drawn iid from $ρ$. We illustrate the convergence of an empirical mean computed by the algorithm to a population mean by simulations from Gaussian random fields.

研究の動機と目的

• トポロジカルデータ解析における統計的解析の中心的役割を果たすパーシステンス図のフレシェ平均を推定するための計算的に実行可能な手法を開発すること。
• 理論的には存在することが保証されているが、その計算的手続きを欠いているフレシェ平均の計算を補填すること。
• ディラック質量の有限和からなる離散確率測度からの独立同分布標本に対して、提案されたアルゴリズムの収束性を確立すること。
• パーシステンス図のフレシェ平均および分散を特徴づけることで、統計的推論の理論的基盤を提供すること。

提案手法

• パーシステンス図の空間に$L^2$-Wasserstein距離を備えた空間上で、勾配降下に基づくアプローチを用いる。
• パーシステンス図のアレクサンドロフ空間におけるフレシェ関数の半concavityを活用し、支援ベクトル（サブ勾配）を定義する。
• 点間の最適輸送双対写像を用いて、図の間の測地線経路を構築し、経路の$K$-リプシッツ連続性を保証する。
• 2つの図の間のすべての測地線経路の集合が前コンパクトであることを証明し、最小化列の収束を可能にする。
• 支援ベクトルを用いて候補平均を反復的に改善することで、経験的フレシェ関数の局所最小値を計算する。
• 標本の凸包内の図の集合が一様であることを確立し、コンパクト性および収束の議論を支援する。

実験結果

リサーチクエスチョン

• RQ1 理論的には存在することが保証されているパーシステンス図のフレシェ平均を、構成的アルゴリズムで計算可能か？
• RQ2 提案されたアルゴリズムは、パーシステンス図の経験的分布に適用した場合に、局所最小値に収束するか？
• RQ3 元の分布がディラック質量の有限混合である場合、経験的フレシェ平均の漸近的性質は何か？
• RQ4 $L^2$-Wasserstein距離がパーシステンス図上で統計的推論および収束保証をどのように支援するか？
• RQ5 パーシステンス図の空間は、下から有界な曲率を持つ測地線的アレクサンドロフ空間とみなせるか？これにより勾配に基づく最適化が可能になるか？

主な発見

• 任意の有限個のパーシステンス図の集合に対して、アルゴリズムは経験的フレシェ関数の局所最小値に収束する。
• 元の測度がディラック質量の有限和である場合、アルゴリズムによって計算された経験的フレシェ平均は弱大数の法則を満たす。
• $L^2$-Wasserstein距離を備えたパーシステンス図の空間は、測地線的空間であり、曲率が下から有界なアレクサンドロフ空間である。
• この空間におけるフレシェ関数は半concaveであり、最適化のためのサブ勾配として支援ベクトルを定義可能である。
• 2つの図の間のすべての測地線経路の集合は前コンパクトである。これにより最小化経路の存在が保証され、アルゴリズムの収束を支持する。
• 有限個の図の凸包の閉包はコンパクトである。これにより収束の証明においてアーヴェラ・アスコリの定理が適用可能になる。