Расстояние Mahalanobis
d2 = mahal(Y,X)Y к ссылочным выборкам в X.
Сгенерируйте коррелированый набор данных двумерной выборки.
rng('default') % For reproducibility X = mvnrnd([0;0],[1 .9;.9 1],1000);
Задайте четыре наблюдения, которые являются равноотстоящими от среднего значения X в Евклидовом расстоянии.
Y = [1 1;1 -1;-1 1;-1 -1];
Вычислите расстояние Mahalanobis каждого наблюдения в Y к ссылочным выборкам в X.
d2_mahal = mahal(Y,X)
d2_mahal = 4×1
1.1095
20.3632
19.5939
1.0137
Вычислите Евклидово расстояние в квадрате каждого наблюдения в Y от среднего значения X .
d2_Euclidean = sum((Y-mean(X)).^2,2)
d2_Euclidean = 4×1
2.0931
2.0399
1.9625
1.9094
Постройте X и Y при помощи scatter и используйте цвет маркера, чтобы визуализировать расстояние Mahalanobis Y к ссылочным выборкам в X.
scatter(X(:,1),X(:,2),10,'.') % Scatter plot with points of size 10 hold on scatter(Y(:,1),Y(:,2),100,d2_mahal,'o','filled') hb = colorbar; ylabel(hb,'Mahalanobis Distance') legend('X','Y','Location','best')
Все наблюдения в Y([1,1] , [-1,-1,], [1,-1], и [-1,1]) являются равноотстоящими от среднего значения X в Евклидовом расстоянии. Однако [1,1] и [-1,-1] намного ближе к X, чем [1,-1] и [-1,1] в расстоянии Mahalanobis. Поскольку расстояние Mahalanobis рассматривает ковариацию данных и шкалы различных переменных, это полезно для обнаружения выбросов.