Вход в функцию earthmovdist::emdL1()
Я пытался использовать emdL1()
функция от earthmovdist
пакет для сравнения двух дистрибутивов. Входами в указанную функцию являются два вектора x
а также y
, Проблема в том, что я не знаю, должен ли я использовать исходный вектор значений как x
а также y
или значения функции распределения плотности (например, density(x)$y
) рассчитать это расстояние.
Я не смог найти больше информации по этой теме в интернете. В качестве примера рассмотрим эти два распределения:
Которые представляют собой сравнение двух сегментов по двум различным переменным. Результаты для EMD в первом случае являются 626
а также 673
соответственно, если код запускается emdL1(x,y)
а также 4522
а также 1491
если код запускается emdL1(density(x)$y,density(y)$y)
,
Ясно, что результат противоположен в зависимости от того, какой подход я выберу. Для меня было бы более разумно, чтобы первые распределения были более разными (парное t.test имеет значение <0,05), но я хотел бы понять, какой подход является правильным.