How to apply distance-based clustering or dimensionality reduction for too many samples

I have a dataset with 200K samples (cases) and 30 variables. Every distance-based method for clustering or dimension reduction technique that I use, such as DBSCAN, Hierarchical Clustering, LLE, Isomap and … fail to run on my machine (normally I get R Session Terminated error) due to the large distance file being generated. (Distance calculation requires o(n^2) time and space)

Is there any solution to handle this problem? Is there any good package for the mentioned clustering or dimensionality reduction in R or Matlab that is suitable ?

Cross Validated Asked by Matin Kh on December 27, 2020

1 Answers

One Answer

Maybe you could try Mini-Batch K-Means. I have Matlab code for it:

function [c,counts,idx] = mbkmeans(x,k,c,counts)
    [N,D] = size(x);
    if ~exist('c','var') || isempty(c)
        c = x(1:min([k N]),:) + bsxfun(@times,randn(min([k N]),D)*0.001,std(x));
        if N < k
            c(N+1:k,:) = bsxfun(@plus,mean(x),bsxfun(@times,randn(k-N,D),std(x)));
    if ~exist('counts','var') || isempty(counts)
        counts = zeros(k,1);
    idx = knnsearch(c,x,'k',1);
    add = full(sparse(idx,1,1));
    counts(idx) = counts(idx) + add(idx);
    lr = 1 ./ counts(idx);   
    for i = 1:N
        c(idx(i),:) = (1 - lr(i)) * c(idx(i),:) + lr(i) * x(i,:);


clusters = mbkmeans(yourdata,numberofclusters);

You may feed it your entire dataset at once and you're done. Or you may feed it smaller subsets. In this case, use it like this:

[c1, counts1] = mbkmeans(subset1,numberofclusters);
[c2, counts2] = mbkmeans(subset2,numberofclusters, c1, counts1); %start clustering using previously created clusters
[c3, counts3] = mbkmeans(subset3,numberofclusters, c2, counts2);
[cn, countsn, indices] = mbkmeans(subsetn,numberofclusters, c(n-1), counts(n-1));

For R, there is the stream package (explanation here). You may also take a look at this, this and this.

Answered by rcpinto on December 27, 2020

Add your own answers!

Related Questions

linear causal model

1  Asked on November 29, 2021 by markowitz


What is the point of test set in ML?

4  Asked on November 29, 2021 by lelouche-lamperouge


Ask a Question

Get help from others!

© 2022 All rights reserved.