library(emdist) my.emd <- function(A,B,r=100){ tandem.AB <- c(A,B) # そのtop 100をカットオフとする s.t.AB <- sort(tandem.AB,decreasing = TRUE) #cutoff <- s.t.AB[floor(length(tandem.AB)*r)] cutoff <- s.t.AB[r] # カットオフ以上の値の番地をと…

2次元の分布としてのメタボロームデータのことをやっている EMdist,Wass...についてはこちらを そのほかとしてはTransportation problemであって、その解法はハンガリアンメソッドだという ハンガリアンメソッドは(重み付き)割り当てマッチング(こちら)でも…

Wasserstein metricは分布の差異を表す指標(Wikiの記事) 1st Mallows and Wasserstein distancesとしてEarth Mover's distance(EMD)がある(Wikiの記事) EMDの計算はRでは、emdistパッケージ EMDはTransportation problem(Wiki記事)を解くことで求められると…