www.gusucode.com > emg_functions > emg_functions/distance_measures/remove_data_mahal.m

     % returns sorted data with the farthest 'removeQty' datapoints removed based
 % on Mahalanobis distance. A is m x n where m is number of features.
function B = remove_data_mahal(A, removeQty)
    distances = mahal(A',A');
    [distances, ind] = sort(distances);
    B = A(:,ind);
    B = B(:,1:length(distances) - removeQty);
end