Wasserstein 距離(ワッサースタイン距離)は、2つの確率分布 P,Q の差を「確率質量を どれだけ・どれだけ遠くへ 移動させれば P を Q に変換できるか」
という 最適輸送(Optimal Transport) の観点で定義する距離。
直感的には、分布を「土の山」とみなし、別の分布の形にするために土を運ぶときの 最小の運搬コスト を測る量であり、 Earth Mover’s Distance (EMD) とも呼ばれる。
(自然言語処理における「編集距離」に近いかも?)
もっとも一般化された定義は p-Wasserstein
p-Wasserstein 距離¶
距離空間 (X,∥⋅∥) 上の確率分布 P,Q に対し、
p≥1 とすると、p-Wasserstein 距離は次で定義される。
Wp(P,Q)=(γ∈Π(P,Q)inf∫X×X∥x−y∥pdγ(x,y))1/p Π(P,Q):周辺分布がそれぞれ P,Q となる結合分布(輸送計画)
γ(x,y):点 x から y へ運ばれる確率質量
∥x−y∥:輸送コスト(通常はユークリッド距離)
1次 Wasserstein 距離¶
実務・理論の両面で最もよく使われるのが 1次 Wasserstein 距離 (Wasserstein-1 distance)。
W1(P,Q)=γ∈Π(P,Q)inf∫∥x−y∥dγ(x,y) from scipy.stats import wasserstein_distance
p = [0, 1, 1]
q = [1, 0, 0]
wasserstein_distance(p, q)
距離の公理を満たす¶
Wasserstein 距離は、次の性質をすべて満たす。
非負性:Wp(P,Q)≥0
同一性:Wp(P,Q)=0⟺P=Q
対称性:Wp(P,Q)=Wp(Q,P)
三角不等式:Wp(P,R)≤Wp(P,Q)+Wp(Q,R)
Kantorovich–Rubinstein 双対表現¶
W1 には重要な双対表現が存在する。
W1(P,Q)=∥f∥Lip≤1sup(EP[f(X)]−EQ[f(X)]) KL ダイバージェンスとの比較¶
KL ダイバージェンスは
DKL(P∥Q)=∫p(x)logq(x)p(x)dx で定義されるが、次のような違いがある。
| 観点 | KL | Wasserstein |
|---|
| 台の不一致 | 発散しやすい | 安定 |
| 幾何構造 | 反映しない | 反映する |
| 距離性 | なし | あり |