Sinkhorn

Wasserstein Barycenter

メモPDF Wasserstein distanceが定義されて 運搬行列が推定できるとする それを使って、複数のデータを「重ね合わせる(分布の運搬コスト最小で重ねる、という意味の重ね合わせ)」としたとき、すべてのデータへのWasserstein distanceの総和が最小になるよう…

Sinkhorn's matrix balancingメモ

Barycenter package Subgradient()関数 Sinkhorn's algorithm python paper Matrix Scalling全般の長いレビュー python コードの説明書き .. math:: \gamma = arg\min_\gamma <\gamma,M>_F + reg\cdot\Omega(\gamma) s.t. \gamma 1 = a \gamma^T 1= b \gamma…

近似計算

参考(こちら) 正方行列を行と列の両方についてstochastic(足して1)になるように補正する処理 こちら これを使うと数独も解けるという(こちら) # sinkhorn sinkhorn<-function(A,epsilon=0.00001){ B<-A row_sum<-apply(A,1,sum) while(max(abs(t(row_sum)-1…