Clustering high-dimensional data

Example 2D space with subspace clusters

The adjacent image shows a mere two-dimensional space where a number of clusters can be identified. In the one-dimensional subspaces, the clusters ${\displaystyle c_{a}}$ (in subspace ${\displaystyle \{x\}}$) and ${\displaystyle c_{b}}$, ${\displaystyle c_{c}}$, ${\displaystyle c_{d}}$ (in subspace ${\displaystyle \{y\}}$) can be found. ${\displaystyle c_{c}}$ cannot be considered a cluster in a two-dimensional (sub-)space, since it is too sparsely distributed in the ${\displaystyle x}$ axis. In two dimensions, the two clusters ${\displaystyle c_{ab}}$ and ${\displaystyle c_{ad}}$ can be identified.

The problem of subspace clustering is given by the fact that there are ${\displaystyle 2^{d}}$ different subspaces of a space with ${\displaystyle d}$ dimensions. If the subspaces are not axis-parallel, an infinite number of subspaces is possible. Hence, subspace clustering algorithms utilize some kind of heuristic to remain computationally feasible, at the risk of producing inferior results. For example, the downward-closure property (cf. association rules) can be used to build higher-dimensional subspaces only by combining lower-dimensional ones, as any subspace T containing a cluster, will result in a full space S also to contain that cluster (i.e. S ⊆ T), an approach taken by most of the traditional algorithms such as CLIQUE,[3] SUBCLU.[4] It is also possible to define a subspace using different degrees of relevance for each dimension, an approach taken by iMWK-Means,[5] EBK-Modes[6] and CBK-Modes.[7]

Projected clustering seeks to assign each point to a unique cluster, but clusters may exist in different subspaces. The general approach is to use a special distance function together with a regular clustering algorithm.

For example, the PreDeCon algorithm checks which attributes seem to support a clustering for each point, and adjusts the distance function such that dimensions with low variance are amplified in the distance function.[8] In the figure above, the cluster ${\displaystyle c_{c}}$ might be found using DBSCAN with a distance function that places less emphasis on the ${\displaystyle x}$-axis and thus exaggerates the low difference in the ${\displaystyle y}$-axis sufficiently enough to group the points into a cluster.

PROCLUS uses a similar approach with a k-medoid clustering.[9] Initial medoids are guessed, and for each medoid the subspace spanned by attributes with low variance is determined. Points are assigned to the medoid closest, considering only the subspace of that medoid in determining the distance. The algorithm then proceeds as the regular PAM algorithm.

If the distance function weights attributes differently, but never with 0 (and hence never drops irrelevant attributes), the algorithm is called a "soft"-projected clustering algorithm.

Projection-based clustering is based on a nonlinear projection of high-dimensional data into a two-dimensional space.[10] Typical projection-methods like t-distributed stochastic neighbor embedding (t-SNE),[11] or neighbor retrieval visualizer (NerV) [12] are used to project data explicitly into two dimensions disregarding the subspaces of higher dimension than two and preserving only relevant neighborhoods in high-dimensional data. In the next step, the Delaunay graph[13] between the projected points is calculated, and each vertex between two projected points is weighted with the high-dimensional distance between the corresponding high-dimensional data points. Thereafter the shortest path between every pair of points is computed using the Dijkstra algorithm.[14] The shortest paths are then used in the clustering process, which involves two choices depending on the structure type in the high-dimensional data.[10] This Boolean choice can be decided by looking at the topographic map of high-dimensional structures.[15] In a benchmarking of 34 comparable clustering methods, projection-based clustering was the only algorithm that always was able to find the high-dimensional distance or density-based structure of the dataset.[10] Projection-based clustering is accessible in the open-source R package "ProjectionBasedClustering" on CRAN.[16]

Bootstrap aggregation (bagging) can be used to create multiple clusters and aggregate the findings. This is done by taking random subsamples of the data, performing a cluster analysis on each of them and then aggregating the results of the clusterings to generate a dissimilarity measure which can then be used to explore and cluster the original data.[17][18] Since high-dimensional data are likely to have many non-informative features, weights can be used during the bagging process to increase the impact of the more informative aspects. This produces "ABC dissimilarities" which can then be used to explore and cluster the original data and also to assess which features appear to be more impactful in defining the clusters. [19] [20] [21]

Not all algorithms try to either find a unique cluster assignment for each point or all clusters in all subspaces; many settle for a result in between, where a number of possibly overlapping, but not necessarily exhaustive set of clusters are found. An example is FIRES, which is from its basic approach a subspace clustering algorithm, but uses a heuristic too aggressive to credibly produce all subspace clusters.[22] Another hybrid approach is to include a human-into-the-algorithmic-loop: Human domain expertise can help to reduce an exponential search space through heuristic selection of samples. This can be beneficial in the health domain where, e.g., medical doctors are confronted with high-dimensional descriptions of patient conditions and measurements on the success of certain therapies. An important question in such data is to compare and correlate patient conditions and therapy results along with combinations of dimensions. The number of dimensions is often very large, consequently one needs to map them to a smaller number of relevant dimensions to be more amenable for expert analysis. This is because irrelevant, redundant, and conflicting dimensions can negatively affect effectiveness and efficiency of the whole analytic process.[23]

Another type of subspaces is considered in Correlation clustering (Data Mining).