34

数理解析研究所講究録
第2057巻 2017年 34-50

位相的データ解析の現在
大林一平*

東北大学原子分子材料科学高等研究機構

Ippei Obayashi
WPI‐AIMR, Tohoku Umiversity

1

はじめに
本講究録では,位相的データ解析(Topological Data Analysis, TDA) について概説する.

TDA で重要なパーシステントホモロジーに注目し,その概要,数学的定義,現状の応用

例,将来の課題等について述べる.

位相的データ解析とは,トポロジーの概念を用いてデータを解析する手法の総称である.データ解析では統計や機械学習のように確率論に基づいた手法,またフーリエ解析のように解析の概念をベースにした手法などが既存の手法としてよく使われてきた.ボロノイ図のような計算幾何学の手法は以前から利用されてきた *1.

計算機によるホモロジーの計算が1990年代にある程度実用的になり,2000年代にパーシステントホモロジーの概念が発展したことにより TDA というアイデアが注目を集めるようになった [2, 3, 4, 5, 6].
パーシステントホモロジーはターゲッ トとなる図形データに関し1‐パラメータ族を考えることでより多くの情報を得ることができるようになったもので,例えばこの1‐パラメー
タ族に空間スケールを割り当てることでデータのマルチスケールな幾何的情報を得ることができるようになった.すでにアモルファスの構造解析 [7, 8] やタンパク質 [9] の解析,
ウイルスの遺伝的進化 [10] の解析,センサーネッ トワーク [11 】の解析,などへの応用例などがある.
**1〒980‐8577宮城県仙台市青葉区片平2‐1‐1, Email: ippei.obayashi.d8@tohoku.ac.jp
ボロノイ図の数学的な概念自体は非常に古いものらしい.実用的な計算アルゴリズムである Formne の走査線 (sweepl\underline{\mathrm{i}\mathrm{n}\mathrm{e}) アルゴリズムは1986年に発表された [1].

35

2

パーシステントホモロジーの概要
パーシステントホモロジーの応用で重要となるのが,パーシステント図である.典型的なパーシステント図の応用の場合,入カデータは2次元もしくは3次元の有限個の点の集合(ポイントクラウド) で,その出力がパーシステント図である.
パーシステント図の定義や計算方法はさておき,まず例題として [7] で述べられている

アモルファスの構造解析の例について見ていこう.図1の(a) と(c) は分子動力学 (MD)
シミュレーションで得た液体シリカとアモルファスシリカの原子配置を描画したものである.原子配置はどちらもランダムに見え,見た目で区別することは難しい.アモルファス

構造と液体の原子レベルでの構造の違いを捉えるのは難しい,というのは良く知られた問題で,この2つの特徴付けは材料科学の重要な問題である.しかし原子配置から計算した

パーシステント図(図1 (\mathrm{b})(\mathrm{d}) ) ははっきりと異なり,簡単にこの2つを区別できるこれはこの2つの原子配置が (見た目ではわからないものの) 決定的に異なる幾何構造を持っている,ということを意味している.図(b) と図(d) を比べると,(d) のほうが特徴的な筋状のものが見える.これは後で詳しく議論するが原子配置の典型的構造が低次元的制約を持
つということを意味している.
\l corner_{1}\mathrm{q}\mathrm{u}|\mathrm{d}

8

7

6

5

4

3

2

0_{0}

図1

1

2

3\mathfrak{g}|\mathrm{R}\mathrm{h}【 \mathrm{A}\wedge 2|56

78\ovalbx{tsmREJCT}_0^6infyo\almthr{B}w\ovalbx{tsmREJCT}_0^63no\iS

シリカの原子配置とそのパーシステント図([71のデータを元に計算)

36

それでは,この図はどのようにして計算されるのであろうか.まず数学的基礎となるのが,アルファ複体とそのフィルトレーションである.入力となるポイントクラウドの位相の情報は明らかに 「 n 個の連結成分」 だけである (図 2(\mathrm{a}) ). そのため,データの幾何的情報を見るために各点に半径 r の円を貼り付ける (図 2(\mathrm{b}) ). すると位相構造が新たに生成され,2つの穴ができる (すなわち, \dim H_{1} =2 である). こうすることによって図 (a) が「なんとなく」 持っている構造を定式化することができる.この図2(b) のような構造をアルファ複体と呼ぶ.
(a)\bullet\bullet\bullet\bullet\bullet\bullet\bullet\bullet\bullet\bullet\bullet

図2

ポイントクラウドと円する
小さい半径では穴が現れない図3

穴が発生する:半径 b

穴が消滅する:半径 d

半径を変えたときの穴の生成と消滅

さて,ここで問題になるのが円の半径 r をいかに選ぶか,ということである.

rによって穴が見えたり見えなかったりするのでこれは重要な問題である.ただ,この図2でいろんな半径を取ってみると,上側の穴は下側の穴より安定して存在する,ということがわかる.このような事実をうまく捉えるために,フィルトレーションの概念を用いる.半径 r

37を 0 から大きくしていき,対象の図形の増大列を作る.すると,図3のように r が増大す

るにつれ穴が発生したり消滅したりする.パーシステントホモロジーの技法を用いると,各穴が発生する半径 b とそれが埋められて消滅する半径 d を対で計算することができる.
この b を発生時刻 (birth time),

d を消滅時刻 (death fime),と呼び,対 (b,d) をパーシス

テンス対,生存対 (Uirth‐deaffi pair) と呼ぶ.最終的に得られるものは対の有限個の集合\{(b_{i}, d_{\dot{ $\iota$}})\}_{i=1}^{s}\cup\{(b_{i_{\rangle} \infty)\}_{i=s+1^{*2} ^{ $\epsilon$+t} であり,これを (x, y) ‐平面にプロットしたものがパーシステント図(persistence diagam) である.ホモロジーは空間が \mathbb{R}^{N} であれば 0 から N‐lの次数までは意味があるので,パーシステント図も普通 0 次から N—l次まで考える.

0

次が連結成分,1次が穴,2次が空洞,の情報を持っている.
実は図1は1次のパーシステント図である. 0 次や2次ではあまり特徴的な違いが見られないが,1次では非常に特徴的な違いが見える.これはシリカガラスの構造がネット
ワークガラスと呼ばれるもので,共有結合によるネットワ一ク構造が重要であるために穴,すなわちリング状の構造を調べられる1次のパーシステント図に特徴が現れるのである *3.
1.0

1.5

(a)
1.0

\bullet0.8\bullet:\bullet0. 5\bullet
.

0.6

.

.

\displaytefrc{B}\i$omega thr{n}

0.0

0.4

\bullet\bullet-0.5\bullet\bullet
.

..

.

0.2

.
\bullet‐1.0-1.\underline{5}_{15}0\mathrm{q}_{0}
-1.0

-0.5

0. 0

図4

0. 5

1. 0

1.5

0. 2

0. 4

0.6

O.S

1.0

Birth

円形の点集合データ (a) とそのパーシステンス図 (b)パーシステント図の意味をより深く考えるため,円形の点集合データ (図 4(\mathrm{a}) ) とその1
*20 次のホモロジーを考えると,どんなに半径を大きくしても連結成分が一つ残る.これは消滅時刻が \infty
であるとみなす*3

ネットワーク型のアモルファスのリング構造を調べるというのは 「リング統計」 と呼ばれる手法が良く使われている.ある意味でパーシステントホモロジーによる解析はリング解析の発展版と見なせる.

38

次のパーシステント図(図 4(\mathrm{b})) を調べよう.1次のパーシステント図には対角線から離れた点が1つある.この点が図 4(\mathrm{a}) の一番目立つリング構造を表現している.パーシステ
ンス対の意味から考えるとその発生時刻 (x 軸) はリングをなす点の間の距離の1/2に対応
し,その消滅時刻 (y 軸) はリングの半径に対応している.また,対角線のそばには数多くのパーシステンス対があるが,そのような対に対応するリング構造は発生してすぐに消滅するようなリングであり,すなわち点の位置をずらすとリングでなくなってしまうようなものが対応している.つまり対角線のそばのパーシステンス対はどちらかというとノイズ的なものであると言える.フィルトレーションを用いる,つまり半径パラメータ r を徐々に変化させることでこういったノイズ的なリングとしっかりした構造を持つリングを定量的に特定することが可能となる.また,対象となるデータに何らかの固有のスケール (例えば原子配置を考えるならば,原子間の距離の平均が重要なスケールになる) があるならば,それを目安にパーシステント図からノイズ的情報を除去することも可能である.

パーシステントホモロジーの数学的基礎

3

ここまでの議論は,数学的な厳密さや計算アルゴリズムの実現可能性などを避けてきた.この節では以下の内容について議論する.
\bullet

アルファ複体を単体複体として定義し,計算機上で実現しやすくする\bullet

単体複体の増大列からパーシステンス対を定義する\bulletパーシステンス対の計算アルゴリズム概要

アルファ複体

3.1

点集合の幾何構造の解析ツールとして,ボロノイ図が良く使われる.図5の8個の点に対し破線がその点集合のボロノイ図である.これは点集合 P=\{u_{i}\}_{i=1}^{N} に対し,空間を「それぞれの点に最も近い領域」 で分割したものである.そして 「点集合の各点を頂点とする」 「2つのボロノイ領域が接するなら対応する点の間に辺を張る」 「3つのボロノイ領

域が同時に接するなら対応する点の間に面を張る」 「4つ以上も同様」 として構成された単体複体はドロネー複体と呼ばれる (図5の実線と灰色の部分,Del(P) と書く).
半径パラメータ r に対応するアルファ複体 Alp (P, r) はDel (P) の部分複体として以下のように実現される [3].
\bulletドロネー複体の各頂点を頂点とする

39'\mathrm{c}\mathrm{c}\mathrm{c}図5\bulletnボロノイ図とドロネー複体

個のボロノイ領域が接するとき,領域に対応する n 頂点に半径 r の円盤を置いてその n 個の円盤が共通部分を持つならそこに (n-1) ‐単体を張る

図6では平面上の3点で r を変えた場合で,一方は面が張られ,もう一方は面が張られない例である.このようにして得られた単体複体は脈体定理 [4] より|\displaystyle \mathrm{A}\mathrm{l}\mathrm{p}(P, r)|\simeq\bigcup_{u_{i}\in P}B_{r}(u_{i})が任意の r>0 で成立する.ここで \simeq はホモトピー同値を意味し, |\cdot| は単体の幾何的実現を意味する.つまり半径 r の円を各点に置いたものとアルファ複体はホモトピー同値となり,ホモロジーの問題を考えるにあたってはAlp (P, r) を考えれば十分である.

図6

半径の大小によって面が張られない例と張られる例

実際に図6を良く見ると半径パラメータ r と3点の外接円半径との大小関係で中央の穴の存在/ 不存在が決まり,それとアルファ複体の定義が合致していることがわかる.ア40

ルファ複体は有限単体複体であるので,計算機で表現するには好都合である.最終的に半径 r を 0 から増大させると,単体複体の増大列が得られる.これはアルファフィルトレー

ションと呼ばれる.また,ドロネー複体上の各単体 $\sigma$ に対し,rを増大させたときにその単体が張られる瞬間の r を r_{ $\sigma$} と書くことにする.アルファフィルトレーションはcgal(http: // cgal. \mathrm{o}\mathrm{r}\mathrm{g} ) などで計算できる.

パーシステントホモロジー

3.2

実は,パーシステントホモロジーはアルファフィルトレ一ションに限らず一般の単体複体の増大列に対して定義できる.そこでまず単体 { $\sigma$_{1}

)

.

.

.,

$\sigma$_{K}} からなる単体複体を考え,以下の条件が成立していると仮定しよう.
Condition 1. \{$\sigma$_{1}, \cdots$\sigma$_{k}\} が任意の 1\leq k\leq K において部分複体となる

このとき, \{$\sigma$_{1},\cdots$\sigma$

緑の幾何的実現を X_{k} と置くと, X_{1} \hookrightarrow. ..\mapsto X_{K} という包含写像の列ができ,この列に射応する体 Q^{*4} を係数とするホモロジー加群の列 H_{*}(X_{1})\rightarrow.

.

.

\rightarrow H_{*}(X_{K}) ができる.これをパーシステントホモロジーという.するとこの列はパーシステントホモロジーの構造定理 [6] により,既約区間(intervaldecomposable)\mathrm{I}[i,j)=0\rightarrow\ldots\rightarrow 0\rightar ow Q\rightar ow 1\cdots\rightar ow 1Q\rightar ow 0\rightar ow\ldots\rightar ow 0
( i 番目から j-1 番目まで非零)の直和で一意に分解できる.例えば

H_{\ell}(X_{\mathrm{i} )11\rightarrow

恥 (X_{2})11\rightarrow

H_{\ell}(X_{3})\rightarrow

11

0 \rightarrow Q \rightarrow\sim Q\rightarrow

.

.

.

.

.

.

\oplus \oplus \oplus

0 \rightarrow Q \rightarrow 0\rightarrow\oplus \oplus \oplusのような分解が可能となる.この分解において 0\rightarrow Q はホモロジーの生成元の 「発生」
を,

Q\rightarrow 1Q はその 「持続」 を, Q\rightarrow 0 はその 「消滅」 を,それぞれ意味する.これが2節で説明した 「穴」 の生成と消滅の数学的に厳密な定義である.そして,各 \mathrm{I}[i,j ) に対して (i,j)*4という対を対応させたものがパーシステンス対である.

議論の簡単さ,計算機での実装の容易さ,のため \mathbb{Z}_{2} を使うことが多い.

41

アルファフィルトレーションは実数 r でパラメータ付けられているので,有限個での増大列の定義を直接適用することはできない.しかし,実質的に問題になるのは r が\{r_{ $\sigma$} | $\sigma$\in \mathrm{D}\mathrm{e}1(P)\} に含まれている場合のみであることがすぐにわかる.これは有限個であるのでドロネー複体の単体を r_{ $\sigma$} によって昇順に並べ,同じ発生半径を持つ単体がある
場合には適切に順序を決めると (この順で並べたものを \{$\sigma$_{1},\ldots,$\sigma$_{K}\} とする), 条件1を満たすようにできる.そして構造定理で得られる各既約区間 \mathbb{I}[i,j ) に対して,(r_{$\sigma$_{i} , r_{$\sigma$_{j} )を考えると,これがパーシステンス対となるのである.
この分解定理自体は体係数の多項式環が単項イデアル整域 (PI)) であること,また,有限生成 PID 係数の加群は有限生成加群の構造定理より一意に分解可能であること,など
から得られる帰結である *5.

ここではパーシステントホモロジーの問題をどのようにして多項式環の問題に還元するか\sear owという一般論は置いておいて,簡単な例題を代わりに用いることでこの分解について直感的に述べよう.ここでは, f_{1} : V_{1}\rightarrow V_{2}, f_{2}:V_{2}\rightarrow V_{3}(Vi, V_{2} V3:有限次元ベクトル空間, f_{1}, f_{2} :線形写像) という写像があるとしよう.すると,,以下のような分解が可能である.

V_{1} \rightarrow f_{1} V_{2} \rightarrow f_{2} V_{3}\mathrm{k}\mathrm{e}\mathrm{r}(f_{1})|| \rightarrow 0|| \rightarrow 0||\oplus\mathrm{k}\mathrm{e}\mathrm{r}(f_{2}\mathrm{o}f_{1})/\mathrm{k}\mathrm{e}\mathrm{r}(f_{1})\rightar ow\sim\rightar ow\sim\oplus\mathrm{k}\mathrm{e}\mathrm{r}(f_{2})\cap \mathrm{i}\mathrm{m}(f_{1})\oplusV_{1}/\mathrm{k}\mathrm{e}\mathrm{r}(f_{2}\mathrm{o}f_{1})\oplus\oplus\rightarrow\mathrm{i}\mathrm{m}(f_{1})/(\mathrm{i}\mathrm{m}(f_{1})\cap \mathrm{k}\mathrm{e}\mathrm{r}(f_{2})\oplus\oplus0\oplus\rightar ow\sim\mathrm{i}\mathrm{m}(f_{2}\mathrm{o}f_{1})\oplus0 \rightarrow \mathrm{k}\mathrm{e}\mathrm{r}(f_{2})/(\mathrm{i}\mathrm{m}(f_{1})\cap \mathrm{k}\mathrm{e}\mathrm{r}(f_{2})) \rightarrow 0\oplus\oplus\oplus0 \rightarrow V_{2}/(\mathrm{i}\mathrm{m}(f_{1})+\mathrm{k}\mathrm{e}\mathrm{r}(f_{2})) \rightarrow\sim \mathrm{i}\mathrm{m}(f_{2})/\mathrm{i}\mathrm{m}(f_{2}\mathrm{o}f_{1})\oplus\oplus\oplus0 \rightarrow 0 \rightarrow V_{3}/\mathrm{i}\mathrm{m}(f_{2})
この分解を見ればわかるように,分解定理は包含写像から導出されるホモロジー加群の列に限らず適用可能である.実際そういったアイデアによるパーシステントホモロジーの拡張(例えばジグザグパーシステンス [12, 13]) が考案されている.
*5体係数の多項式環が PID であることが本質的であるため,通常のようにホモロジーの係数は \mathb {Z} にはし
ない.

42

3.3

アルゴリズムの概要

この分解を計算するためには,上に挙げた分解に対応するようなうまいベクトル空間の基底を見つければよい.これは結局掃き出し法に帰結され,実際のアルゴリズムも掃き出し法に基いている.特に通常の単体複体の増大列に関するパーシステントホモロジーの分解は非常にシンプルかつ実用的なアルゴリズムが知られている.行列 B を基底\{$\sigma$_{1}, \cdots, $\sigma$_{K}\} に関する境界作用素 \partial の行列表示である,としよう.すると条件1より B
は上半三角行列となる.そこで 「左の列のスカラ倍を右の列に加える」 という操作のみで掃き出し法を進める.最終的に掃き出し法を進められなくなった段階での行列の各列が求める基底となり,求める分解が得られる.このアルゴリズムを [14] に述べられているように記述すると以下の通りである.

\displaystyle\frac{\mathrm{A}1\mathrm{g}\mathrm{o}\mathrm{n}^{\mathrm{i}\mathrm{t}\mathrm{h}\mathrm{m}1.\ovalbox{\t\smal REJ CT}.\text{き}\mathrm{f}\mathrm{f}\mathrm{i}\text{し}\mathrm{a}\mathrm{e}(\hslash^{$\theta$}\mathb {Z}_{2}\emptyset\ovalbox{\t\smal REJ CT}_{\mathrm{D}^{\mathrm{a}\mathrm{A}){\mathrm{f}\mathrm{o}\mathrm{r}j=1,.n\mathrm{d}\mathrm{o}while there exists i<j with low(i) =1\mathrm{o}\mathrm{w}(j) do
add column i to column j

ここでlow(i) は i 列で 0 でない要素の最大のインデックスであるとする (すべて 0 のときはlow(i) は未定義とする). たったこれだけのプログラムでパーシステンス対が計算できるのである.

このアルゴリズムは非常に簡単であるため,データ (行列) の表現法,並列化,といった様々な工夫が可能である.論文 [14] にはそういったアルゴリズムに関する話題や実際の計算ソフトウェアの性能評価などがなされている.

パーシステントホモロジーの応用例

4

この節では,パーシステントホモロジーの2つの応用例について詳しく紹介する.一つはウイルスの進化の解析 [10] で,もう一つはアモルファスの構造解析 [7, 8] である.これ

らの応用例によってパーシステントホモロジーの有用性を示す.
4.1

ウイルスの進化の解析

ここでは,[10] で述べられているウイルス進化のパーシステントホモロジーによる解析について紹介する.

43

ここで使われる道具は,ヴイートリス

リップス複体 (しばしばリップス複体と略される ) [4] である.これまで議論してきたアルファ複体は,点に半径 r の円を貼り付けること

で位相構造を構成している.しかし,ウイルスの遺伝子の研究をするにあたってはこの円という構造は意味を持たない.しかし遺伝子間の距離という概念は定義可能である.2つの遺伝子の距離をゲノム配列の違いの大きさで計測する,というのは遺伝子の研究でよく使われるアイデアである.例えばハミング距離のようなもので距離を計測する.そこで,点の間の距離だけから定義できるリップス複体を以下のように定義する.

\mathrm{V}\mathrm{R}(P,r)= { $\sigma$=\{v_{1_{\rangle}}\ldots, v_{k}\}\subset P|d(v_{i},v_{j})\leq r for all i,j }ここで P=\{u_{i}\} は点の集合で,d は P 上の距離, r>0 は距離パラメータである.

rを変化させることでアルファフィルトレーションの場合と同様にフィルトレーションを構築できる.

遺伝子の進化に関しては2つの重要な要素がある.一つは 「垂直伝播」 で,親から子への遺伝情報の伝播である.ウィルスの場合は増殖過程でこれが行われる.もう一つは 「水平伝播」 で,ある個体が他の固体の遺伝子を取り込むことである.ウイルスの場合は遺伝子再集合 (reassottment) と呼ばれる現象があり,複数のウイルスが同じ宿主に感染した場合にそのウイルスの遺伝子が混ざり合う.ウイルスの場合には遺伝情報を守る壁 (細胞) がないため,このような混合現象が比較的頻繁に生じる,ということである.
ここで垂直伝播とランダムな遺伝子の変異しかない世界を考えてみよう.すると元々同じ遺伝子を持っていた個体が変異によって異なる遺伝情報を持つようになり,それが子孫へ伝えられさらに変異する,となる.このようにして遺伝子の多様性が増大し,遺伝子の距離による構造は木構造になる.遠い祖先で異なる遺伝子を持つようになった2つの個体は距離が大きく,近い祖先で異なる遺伝子を持つようになった2つの個体は距離が近くなる.これによって木構造が実現される.そして木構造の1次のホモロジー加群は自明となるため,この世界では遺伝子距離による1次のホモロジー加群の次元は 0 に近くなるはずである.しかしここで水平伝播という概念が加わると事情は変化する.水平伝播が生じる

ことで現れる個体の遺伝子は,その元となった2つ(以上) の個体の遺伝子のどれともそれなりに近しいはずである.つまり水平伝播によって木構造の異なる枝がつながるという現象が生じる.これによって枝の間に橋ができることで1次のホモロジーが現れる.これをパーシステントホモロジーで解析するのである.

ここで紹介している論文ではシミュレーションでの遺伝子の変化の様子と実際のインフ44

ルエンザウイルスの遺伝子データをパーシステントバーコード *6 を用いて比較すること

で,例えば遺伝子再集合が生じる頻度を推定したり,どのようなパターンが特徴的に頻繁に生じるかを議論したりしている.

4.2

アモルファスの構造解析

次に述べるのは [7, 8] でのアモルファスの構造解析である.2節で少し紹介したが,この節でより詳しく解説する.

アモルファスは非晶質とも呼ばれ,結晶のような繰り返し構造 (長距離秩序,Long Range
Order と呼ばれる) は存在しないが,最近接原子や隣の隣の原子までの距離といった短距

離では構造を持つ (短距離秩序,Short Range Order と呼ばれる) ような固体である.代表的な物質としてシリカガラス (\mathrm{S}\mathrm{i}\mathrm{O}_{2}) がある.一方液体も短距離秩序がある一方で長距離構造を持たず,アモルファスと液体の違いは何なのかははっきりとはわかっていない.現在のところ,アモルファス特有の構造として中距離秩序(Midium Range Order) と呼べるものが存在するのではないかと考えられている.ここで紹介する論文では,この中距離秩序をパーシステント図で特徴付けている.

この論文の内容を紹介するために,まずopfimal cycle[15, 16] について紹介する.ホモロ
ジーはチェイン複体 \{C_{q}\}_{q} 上の境界作用素 \partial_{q} : C_{q}\rightarrow C_{q-1} によって H_{q}=Z_{q}/B_{q}, Z_{q}=\mathrm{k}\mathrm{e}\mathrm{r}\partial_{q}, B_{q}=\mathrm{i}\mathrm{m}\partial_{q+1} と定義される. H_{q} の基底ベクトルが図形上の穴や空洞などを表している.しかし, B_{\mathrm{q} で割っているため,この基底ベクトル自体は具体的な幾何的オブジェ

クトではない.例えば図7に対しては \dim H_{1}=1 であるが, z_{i}+B_{1} ( i=1,2 )3) はすべてこの H_{1} の基底である.しかしこのなかで掬がこの穴を一番よく表していると考えられる.この z_{3} の特徴付けとして何が良さそうかを考えると,「ループの長さが最も短い」のが良いのではないかと考えられる.このような z_{3} のことをoptimalcycleと呼ぶ.基本

的にこのようなoptimal cycle を探す問題は線形空間上のある種の最適化問題として定式化され,その計算アルゴリズムが良く調べられている.[15] は通常のホモロジー加群に関する optimal cycle についての論文であり,[16] によってパーシステントホモロジーに関する問題として拡張されている.このoptimal cycle を用いるとパーシステント図の上の各パーシステンス対が元データのどの部分にあたるのかを知ることができる.optimal
cycle はここで解説するアモルファスの解析に有効利用されている.
図8は[7] より引用した,液体シリカ/アモルファスシリカ/結晶シリカの原子配置より*6パーシステントホモロジーの可視化手法の一種.パーシステント図と同じ情報を別の形で可視化する.

45

図7

Optimal Cycle の問題計算された1次のパーシステント図である.結晶の場合は,パーシステンス対はいくつかの小さい島状に分布し,液体の場合は,2次元的に分布が広がっている.アモルファスの場合はどちらとも異なり,曲線的な分布をしている.パーシステント図を用いることでア
モルファスを結晶,液体の双方から区別することができたのである.アモルファスシリカ

は共有結合による乱雑なネットワーク構造が主要な構造である (このようなガラスをネッ

トワークガラスと呼ぶ) ので,1次のパーシステント図が有用なのである.
\mathfrak{l} \$}$\phi$\mathr{}_\overlin{\mathr{o}^$\delta$}3$\tau_{3}^frc\ovelin{}r $\alph}{simeq\athr{J}^2\hat{C_i}^{s 2}\wedg\wedg.
10Birth. $\Lambda$^{\el_{\mathrm{j}^ $\xi$}

図8\mathrm{B}\mathrm{i}\mathrm{r} $\iota$ t_{1.i}'\hat{1}^{\dot{k} |\mathrm{B}i_{1}\cdot\mathrm{t}\mathrm{k}_{\backsla h!^{:}^{:};.\mathrm{A}\sim|\backsla h

液体シリカ/ アモルファスシリカ/結晶シリカの1次のパーシステント図([7], Fig2より引用,カラーの図表をグレイスケール化している)

さらにこの分布の意味するところを調べるためにoptimal cycle が活用される.図8の中央のアモルファスに対する図には,パーシステント対に対するoptimal cycle を示している (大きい丸が酸素原子,小さい丸がケイ素原子を表している). C_{T} は最近接の3原子がなす構造と対応しており,これは上で言う短距離秩序に対応している. C_{O}, B_{O}, C_{P}
はより離れた原子からなる構造を表していることがわかる.これが中距離秩序に対応する構造であろうと考えられる.

曲線的な分布については,パーシステント図の安定性の概念と関連している.入力のポ
イントクラウドを連続的に動かしたとき,パーシステント図上の各パーシステンス対は連続的に変化することが知られている.図8の C_{O} は酸素原子3個からなる構造が典型的に原子配置上に現れていることを意味するが,さらにそれが曲線的分布をしているというこ46とは,この安定性の概念から典型的配置が低次元的な制約を受けているということを意味することがわかる.実際この3原子を取り出してきてそれがなす三角形の二辺の長さとその間の角度をプロットすると低次元的構造を持つことがこの論文では示されている.曲線的分布の幅はその低次元構造からのずれ具合を表している.
この論文ではパッキングガラスと呼ばれる金属原子によるガラスで典型的な構造も扱っている.半径のそろった原子を3次元空間につめこむと,典型的には体心立方格子 (BCC),
面心立方格子 (FCC), 六方最密充填 (HPC) といった結晶構造をなす.しかし大きさの異なる原子が3次元空間に詰め込まれることでランダムなばらつきが生じ,アモルファスとなる.この場合のアモルファスの特徴付けは1次のパーシステント図ではなく主に2次のパーシステント図でできることが示されている.球による空間充填構造が重要なので2次のパーシステント図が有効なのである.

このように,様々な種類のアモルファスがパーシステント図によって特徴付け可能である.ネットワークガラスとパッキングガラスの違いがパーシステント図の次数の違いに現れるというのも興味深い.

その他の問題など

5

ここでは,上の応用例では説明できなかったTDA の広がりについて簡単に解説する.

5.1

統計手法や機械学習との連携

この節では統計や機械学習の手法をいかにパーシステントホモロジーの技法と結び付けるかについての研究の現状について簡単に紹介する.

統計手法や機械学習は現在のデータ分析の主流の手法である.データから法則,パターンを見つけだし,新たなデータにそのパターンをあてはめることでその性質を推定したりすることができる.統計や機械学習では (例外もあるが) 通常入力としてベクトルを取る.

元データをなんらかの方法でベクトル化してから学習させることが多い.このベクトル化の善し悪しが機械学習の性能を決めることも多い *7.
パーシステント図を統計や機械学習で活用するためには,パーシステント図をベクトルに変換するのが最も自然である.元々のデータに含まれる幾何的特徴がパーシステント図に効率的に抽出されて含まれていると考えられるため,パーシステント図を経由することでより良い処理ができるはずである.現在,ベクトル化には以下のような様々な手法が提*7

このベクトルは 「特徴量」 「記述子」 「素性」 などと呼ばれる.分野によってどの用語を使うかが異なる.

47

案されている.
\bulletPersistence landscape[17]\bulletPersistence Scale Space Kemel (PSSK)[18]\bulletPersistence Weighted Gaussian Kernel (PWGK)[19]\bulletPersistence image[20]

今のところ決定版と呼べる手法はなく,様々な手法が試されている.基本的にはどれも適当な L^{2} 関数空間に埋め込む手法である.[18] では3次元形状データの分類問題への応用がなされ,[19] では上の応用例と同様のアモルファスの原子配置データを用いた液体からガラスへの転移点の検出などが応用としてなされている.
5.2

多次元パーシステンス

ここまでの解説では,パーシステントホモロジーは一列に並んだ単調増大列を考えてきた.しかし,二次元的な増大列を考えることができるならばより便利であろう.つまり,以下の図のような2軸の包含関係を考えるということである.

X_{1,1} \hookrightarrow X_{1,2} \hookrightarrow\upar ow\mathrm{J}^\cdot}

X_{2,1 ,\upar ow} \ho krightar ow X_{2,2 ,\upar ow} \ho krightar ow..

.

[21] で述べられている応用例としては,一つの軸はこれまで利用してきた,貼り付ける円を大きくする方向,もう一つはノイズ除去の強さによる方向,とするものである.パーシ

ステントホモロジーは点がずれるようなノイズには強いものの,外れ値には弱いという性質がある.そこで外れ値の検出技法によってそのような点を除去することを考える.しかしここでノイズ除去の強さパラメータを決めるという問題が新たに生じる.そこでこのパラメータもフィルトレーションに含めてしまおう,というアイデアである.
このためには,以下のようなホモロジーの可換図式の中にある構造を一次元パーシステ

ントホモロジーの構造定理のような方法で同定できればよいのである.

H_{\ell}(X_{1,1}) \rightarrow H_{l}(X_{1,2}) \rightarrow\upar ow

H_{\el }(X_{2,1})\upar ow\upar ow\rightarrowHゑ

(X_{2,2})\upar ow\rightarrow\cdots48

しかし,[21 】で述べられている事実は,そのような構造を一般的に見つけだすのは二次元以上の場合では不可能である,ということである.しかしそうはいっても多次元のパーシ
ステントホモロジーを考えることに需要はある.そこで考えられるアイデアとしては\bullet

この可換図式の構造を部分的に説明できるような道具を探す\bullet図式の形状を限定して考える\bulletよくあるデータ,典型的データでうまくいく道具立てを考える[21] では一つ目のアプローチにより rank invariant という概念を提案している.

RIVET[22] ではこのアイデアを使って多次元パーシステントホモロジーの見事な可視化
ソフトウェアを実現している.
[23] では二つ目のアプローチで comnutative ladderという名前で図式が小さな梯子型( 2\times 2, 2\times 3, 2\times 4) ならば構造を計算できる,ということを示している.この証明には表現論のツールなども使われている.

三つめのアプローチの成功例は寡聞して聞いたことがないが,データ処理の文脈から言えばこういったアイデアも何かあるのではないだろうか.例えば行列の対角化を考えると,一般的にはジョルダン標準形などを考えなければすべての場合にうまくいくとはいえないが,実際問題としては対象となる行列は対角化可能だろうと仮定して数値計算することが多い.また,対称行列に限定して対角化したい,といったときにはより効率的なアプローチもできるわけである.正直あまり良いアイデアがあるわけではないが,こういう方
向性も有望ではないだろうか.

6

おわりに
この講究録で述べられてきたように,パーシステントホモロジーは様々な応用がある.

さらにウイルスの遺伝子の解析でとアルファ複体ではなくリップス複体を用いたり,また

アモルファスの解析でoptimal cycle という概念が有効であったように,データからちょっとパーシステント図を計算すればおしまい,というものではなく,様々な工夫や新たな数学的概念が必要となる.多次元のパーシステントホモロジーのように数学的な困難を持った問題も多い.多次元のパーシステントホモロジーをやろうと思うと表現論が顔を出したり,ここでは説明しなかったが圏論が有効だったりもする.このようにパーシステントホ
モロジーと他の分野の数学との繋がりも見えてきた状況である.統計や機械学習との連携もまだ始まったばかりである.新たな応用に対しては新たな理論的道具が必要とされるはずである.パーシステントホモロジー,位相的データ解析はまだまだ理論的問題や応用の49

問題が山積みであり,これらは今後も発展しつづけると筆者は期待している.

参考文献[1] S. Fortune. A sweepline algorithm for Voromoi diagrams. Proceedings of the second
annual symposium on Computational geometry. 313‐322. (1986)[2] G. Carlsson. Topology and Data. Bun. Amer. Math. Soc. 46 (2009), 255‐308.
[3] H. Edelsbrunner and J. Harer. Computational Topology: An Introduction. AMS (2010).
[4] 平岡裕章,タンパク質構造とトポロジー ‐パーシステントホモロジー群入門‐, シリー

ズ現象を解明する数学 (三村昌泰・竹内康博・森田善久編集),共立出版 (2013).
[5] H. Edelsbrunner, D. Letscher, and A. Zomorodian. Topological Persistence and Simpli‐
ficauon. Discrete Comput. Geom. 28(4) (2002), 511‐533.

[6] A. Zomorodian and G. Carlsson. Computing Persistent Homology. Discrete Comput.
Geom. 33(2) (2005), 249‐274.

[7] Y. Hiraoka, T. Nakamura, A. Hirata, E. G. Escolar, K. Matsue, and Y. Nishiura. Hi‐
erarchical structures of amorphous solids characterized Uy persistent homology. PNAS

2016113 (26), 7035−7040[8] T. Nakamura, Y. Hiraoka, A. Hirata, E. G. Escolar, and Y. Nishiura. Persistent homol‐
ogy and many‐body atomic structure for medium‐range order in the glass. Accepted in

Nanotechnology.
[9] M. Gameiro, Y. Hiraoka, S. Izumi, M. Kramar, K. Mischaikow, and V. Nanda. Topo‐logical measurement of protein compressibility via persistent diagrams. Japan J. Indust.

Appl. Math. 32 (2015), 1‐17.
[10] J. M. Chan, G. Carlsson, and R. Rabadan. Topology of viral evolution. PNAS 110(46)(2013), 18566−18571.
[11] V. de Silva and R. Ghrist. Coverage in sensor networks via persistent homology. Alge‐
braic and Geometric Topology 7 (2007), 339‐358.

[12] G. Carlsson, V. de Silva, and D. Morozov. Zigzag Persistent Homology and Real‐valued
Functions. Proceedings of the Annual Symposium on Computational Geometry (2009)247‐256.

[13] G. Carlsson and V. de Silva. Zigzag Persistence. Foundations of Computational Mathe‐
matics 10(4) (2010), 367‐405.

[14] N. Otter, M. A. Porter, U. Tillmann, P. Grindrod, and H. A. Harrington. A roadmap

50

for the computation of persistent homology. http: // arxiv. \mathrm{o}\mathrm{r}\mathrm{g}/\mathrm{a}\mathrm{b}\mathrm{s}/1506.
08903.

[15] T. K. Dey, A. N. Hirani, and B. Krishnamoorthy. Optimal homologous cycles, total

unimodularity and linear programming. SIAM Joumal on Computing 40(4) (2011),
1026‐1044.

[16] E. G. Escolar, and Y. Hiraoka, Optimal Cycles for Persistent Homology Via Linear Pro‐
\cdot

gramming, optimization in the Real World: Toward Solving Real‐World Optim zation
Problems, Springer Japan (2016), 79‐96.
[17] P. Bubenik. Statistical topological data analysis using persistence landscapes. The Jour‐nal of Machine Learning Research archive 16(1) (2015), 77‐102.

[18] J. Reininghaus, S. Huber, U. Bauer, and R. Kwitt. A Stable Multi‐Scale Kemel for Topo‐logical Machine Learning. IEEE Conference on Computer Vision and Pattern Recogni‐tion (2015) 4741 −4748.

[19] K. Genki, K. Fukumizu, and Y. Hiraoka. Persistence weighted Gaussian kernel for topo‐

logical data analysis. Proceedings of the 33\mathrm{r}\mathrm{d} International Conference on Machine
Learning (ICML) (2016).
[20] H. Adams, S. Chepushtanova, T. Emerson, E. Hanson, M. Kirby, F. Motta, R. Neville,
C. Peterson, P. Shipman, and L. Ziegelmeier. Persistence Images: A Staule Vector Rep‐

resentation of Persistent Homology. http: // arxiv. \mathrm{o}\mathrm{r}\mathrm{g}/\mathrm{a}\mathrm{b}\mathrm{s}/1507.

06217.

[21] G. Carlsson and A. Zomorodian. Discrete and Computational Geometry 42 (2009)71‐93.

[22] M. Lesnick and M. Wright. RIVET. http: // rivet.online/.
[23] E. G. Emerson and Y. Hiraoka. Persistence Modules on Commutative Ladders of Fimte

Type. Discrete and Computational Geometry 55 (2016) 100‐157.