K-Volume Clustering Algorithms for scRNA-Seq Data Analysis

TL;DR

This paper introduces a new clustering algorithm for single-cell RNA sequencing data that improves performance by using geometric volume as a clustering criterion.

Contribution

The novel K-volume clustering algorithm uses convex volume as a criterion to optimize hierarchical structure and cluster count simultaneously.

Findings

K-volume clustering outperforms traditional methods in various biological applications.

The algorithm optimizes hierarchical structure and cluster number through nonlinear optimization.

Validation on real datasets confirms the algorithm's effectiveness.

Abstract

Clustering high-dimensional and structural data remains a significant challenge in computational biology, particularly for complex single-cell and multi-omics datasets. In this work, we introduce a novel clustering algorithm that utilizes the total convex volume defined by points within a cluster as a biologically relevant and geometrically interpretable criterion. This approach simultaneously optimizes both the hierarchical structure and the number of clusters at each level through nonlinear optimization. We evaluate our algorithm against other clustering methods, and the results demonstrate that our approach outperforms traditional techniques across a variety of biological applications.

Clustering high-dimensional and structural data remains a key challenge in computational biology, especially for complex single-cell and multi-omics datasets. In this study, we present K-volume clustering, a novel algorithm that uses the total convex volume defined by points within a cluster as a biologically relevant and geometrically interpretable criterion. This method simultaneously optimizes both the hierarchical structure and the number of clusters at each level through nonlinear optimization. Validation on real datasets shows that K-volume clustering outperforms traditional methods across a range of biological applications. With its theoretical foundation and broad applicability, K-volume clustering holds great promise as a core tool for diverse data analysis tasks.

Figures

Full text

1. Introduction

Clustering is a fundamental task in data analysis, aiming to group similar data points based on their intrinsic properties and patterns, without relying on labeled data [1]. Among the various clustering methods, K-means and its variants are widely recognized for their simplicity, scalability, and efficiency [2,3]. The K-means algorithm partitions a dataset into K clusters by iteratively minimizing intra-cluster variance and maximizing inter-cluster separation. Variants such as K-center [4], K-median [5] and K-density [6] have been introduced to better handle diverse data distributions, improve robustness to outliers, and capture nonlinear relationships. However, these methods typically require the number of clusters, K, to be predefined, which presents a significant challenge, especially with complex datasets.

Clustering methods have become essential tools in biological data analysis, particularly with the rise of high-throughput, high-dimensional datasets like single-cell RNA sequencing (scRNA-seq) [7]. scRNA-seq provides a powerful means of analyzing gene expression at the single-cell level, enabling the identification of diverse cell types, states, and lineages. In this context, clustering is critical for grouping cells with similar expression profiles to uncover cell heterogeneity and infer biological function [8,9]. Methods such as K-means, hierarchical clustering, and graph-based clustering have been adapted to handle the high dimensionality and sparsity of scRNA-seq data [7,10]. Advanced variations, including those that integrate dimensionality reduction techniques like Principal Component Analysis (PCA) or t-SNE, have further enhanced the interpretability of scRNA-seq clusters, facilitating the discovery of novel cell populations and their functional characteristics [11].

Identifying the optimal number of clusters (K) and defining hierarchical layers (H) remain significant challenges in both theoretical analysis and practical applications [3]. In K-means clustering, the requirement to predefine K can lead to complications; improper choices may result in over- or under-clustering, distorting the underlying data structure. Similarly, hierarchical clustering lacks a universal criterion for determining the appropriate number of layers or meaningful cut points in the dendrogram, often leading to subjective or inconsistent interpretations [12]. Techniques like the elbow method [13], silhouette analysis [14], and gap statistics [15] have been developed to estimate K, but these methods are frequently computationally expensive and highly sensitive to dataset characteristics. This issue is particularly pronounced in single-cell RNA sequencing (scRNA-seq) data, where biological processes span multiple scales and resolutions [16]. Additionally, the number of clusters can vary depending on the specific objectives, such as identifying common cell populations or detecting rare cell types [7,9]. The hierarchical relationships between cell populations, shaped by lineage progression and differentiation pathways, further complicate the analysis [17,18]. Therefore, there is a critical need for robust, automated approaches that can simultaneously determine the number of clusters and uncover hierarchical structures, especially given the complex, high-dimensional nature of scRNA-seq data.

In this study, we introduce K-volume clustering, a novel and foundational algorithm for optimal hierarchical clustering. This algorithm simultaneously optimizes both the number of hierarchical layers H and the number of clusters K within each layer using nonlinear optimization principles (see Figure 1). At each level, it maximizes the difference between the area of the convex hull encompassing all sample points and the cumulative areas of the K sub-convex hulls (Figure 1a). By jointly optimizing H and K, the method determines the optimal hierarchical structure, resulting in a robust clustering solution (Figure 1b). This new approach is particularly useful for revealing novel insights from complex biological datasets, enabling a deeper understanding of biological architectures and potentially uncovering hidden functional or evolutionary relationships.

1.1. Related Work

Clustering partitions data points into subsets, where each subset forms a cluster based on a specified criterion [1]. Clustering algorithms group ‘similar’ data points to uncover the relationships among them. Common similarity measures include Euclidean distance [19], squared Euclidean distance [20,21], and Hamming distance [22].

Clustering algorithms are generally categorized into center-based and density-based methods. Specifically, the K-volume clustering algorithm defines the “center” as an inner point of the convex hull formed by all the points in the cluster, which aligns it with the characteristics of a center-based clustering method. At the same time, the size of the convex hull defined by the K-volume algorithm corresponds to the “density” of the cluster. Thus, our K-volume clustering method effectively integrates aspects of both approaches.

Table 1 summarizes widely used clustering algorithms, with entries marked by ∗ denoting our contributions in this paper. Notably, center-based clustering problems are NP-hard, and finding optimal K-means clusters remains NP-hard, even in two-dimensional space [23]. Table 1 also lists the most recent approximation algorithms for these clustering problems.

In the work of Crescenzi and Giuliani [26], the authors identified mutual similarities from cluster analysis among differently labeled statistical units, generating new classifications at varying levels of detail that can be interpreted in biological terms. Similar to this work, our work offers another effort in the clustering analysis of complex scRNA-seq data at multiple levels of detail.

1.2. Our Contribution and Paper Organization

In this paper, we develop algorithms specifically for clustering structural data, such as biological and single-cell omics. In particular, we propose a novel clustering algorithm that uses the total convex volume enclosed by points within the same cluster as a clustering measure.

We present an intuitive and efficient algorithm along with its theoretical analysis in Section 2. In Section 3, we evaluate its performance against well-established clustering algorithms on biological data. Finally, we conclude our paper in Section 5.

2. A Greedy K-Volume Clustering Algorithm

We consider a clustering problem with the following input: a set of N data points, $[eqn]$ , where each point is represented as a D-dimensional real vector in $[eqn]$ . The distance between any two points, $[eqn]$ and $[eqn]$ , is denoted by $[eqn]$ . We introduce a new measure based on the convex volume. Specifically, the clustering cost is defined as the minimal convex volume enclosing all points within a cluster. Under this definition, a single data point has a volume of 0, and any set of collinear points also has a volume of 0.

We define the problem of partitioning data into K clusters within a hierarchical structure while minimizing the total convex volume as the K-volume clustering problem. When the dimensionality is $[eqn]$ , we refer to it as the K-area clustering problem.

2.1. The Algorithm’s Idea

We propose a simple yet elegant algorithm based on a greedy approach. In this work, we illustrate the idea for the case where $[eqn]$ , with the approach being naturally extendable to higher dimensions ( $[eqn]$ ). The algorithm starts with the convex hull of the data points. The key idea is to partition a cluster into two in such a way that the convex volume of the remaining clusters is minimized. This process is repeated iteratively until the total number of clusters reaches the specified value K or the given maximal hierarchical level H.

Given a set of points, finding the initial convex hull takes time $[eqn]$  [27]. The main algorithmic challenges in finding K-volume clusters are $[eqn]$ how to partition the data points into K clusters, and $[eqn]$ how to calculate the total volume of these clusters.

Consider a set of N points, $[eqn]$ , where each point $[eqn]$ is described by its coordinates $[eqn]$ in 2D space. Without loss of generality, we assume all x-axis and y-axis values are non-negative. Our goal is to partition the points into K clusters, i.e., K convex areas on the 2D plane, such that the total area is minimized. Let S represent a set of points, and $[eqn]$ denote the area of the convex polygon enclosing S.

The algorithm proceeds as follows: Initially, using Graham’s scan [28], we calculate the convex polygon for all N points. If $[eqn]$ , this convex polygon is the optimal solution, and all points belong to the same cluster. If $[eqn]$ , we proceed to build a tree that represents the set of clusters created during the execution of the algorithm. The tree structure allows the algorithm to efficiently organize and partition the data into meaningful groups. Since each layer in the tree represents a finer partitioning, the algorithm benefits from the hierarchical structure when it comes to refining clusters, especially as the number of clusters increases. The tree has the following properties:

1. The root node represents the initial convex polygon for all the data points.
2. Each node in the tree corresponds to a cluster. The root of each subtree represents the cluster for all the data points clustered in its subtree.
3. The leaf nodes represent the set of clusters that we currently have at any point in time.

For $[eqn]$ , we iteratively partition one of the current clusters, say C, into two subclusters, $[eqn]$ and $[eqn]$ . These subclusters, $[eqn]$ and $[eqn]$ , become the children of node C. This partitioning process continues until we have exactly K clusters or reach the maximal hierarchical level H. When partitioning a cluster C into two subclusters $[eqn]$ and $[eqn]$ , we use a greedy approach that maximizes the size of the removed area, ensuring an optimal division.

$$[eqn]$$

2.2. The Algorithm’s Description

The algorithm is described in detail as follows. Consider the set C. If C is partitioned into two convex polygons, $[eqn]$ and $[eqn]$ , we observe that $[eqn]$ and $[eqn]$ do not overlap, and there must exist a straight line $[eqn]$ , defined by two points a and b, that separates $[eqn]$ and $[eqn]$ . We maintain the points of C in two lists: $[eqn]$ , the set of points defining the convex polygon, and $[eqn]$ , the set of points in the interior of the convex polygon. The algorithm consists of three steps:

1. Identify the straight lines $[eqn]$ that can separate a cluster C.Given the convexity requirements for the clusters after any partitioning, any partition of the set C into two convex polygons $[eqn]$ and $[eqn]$ can be achieved by introducing a straight line that crosses two points from the set $[eqn]$ . Given two points $[eqn]$ and $[eqn]$ , the line $[eqn]$ is defined by the equation

$$[eqn]$$

which simplifies to

$$[eqn]$$

In total, there are $[eqn]$ such straight lines that can partition the cluster C. 2. Given a convex polygon C and a straight line $[eqn]$ , calculate the two convex polygons $[eqn]$ and $[eqn]$ separately by $[eqn]$ using Graham’s scan algorithm [28].The line $[eqn]$ partitions C into two sets $[eqn]$ and $[eqn]$ , where $[eqn]$ and $[eqn]$ . These sets are defined as follows:

$$[eqn]$$ $$[eqn]$$

Due to the convexity of C, $[eqn]$ , and $[eqn]$ , we also have the following relationship:

$$[eqn]$$

After partitioning the set C using the line $[eqn]$ into the point sets $[eqn]$ and $[eqn]$ , we construct convex polygons for each set using Graham’s scan algorithm [28]. This process results in sets $[eqn]$ , $[eqn]$ , $[eqn]$ , and $[eqn]$ . Additionally, we index the points in these sets in a clockwise order for the next step. 3. Calculate the area $[eqn]$ of a convex polygon $[eqn]$ . Then, determine the maximum area that can be removed by a single partition of the cluster.Consider the point $[eqn]$ and label the points clockwise on the convex polygon C. Using the Shoelace formula [29], the area of the convex polygon defined by the points $[eqn]$ , $[eqn]$ , …, $[eqn]$ is given by

$$[eqn]$$

For each line $[eqn]$ , we calculate the convex polygons resulting from the partition of the convex polygon C. Then, we compute the value $[eqn]$ . The best straight line $[eqn]$ is chosen to minimize the maximum area removed by partitioning the set C.

In Algorithm 1, we present the K-area clustering algorithm, denoted as K-AC, for clustering biological data. This algorithm partitions N data points into at most K clusters, C, or organizes them into at most H hierarchical levels of clusters. **Algorithm 1 **K-AC: K-area clustering algorithm (N, C, K, H)

• 1:Identify the convex hull C for the set of N data points.
• 2:Maintain a set S of convex polygons, initialized as $[eqn]$ .
• 3:Set $[eqn]$ and $[eqn]$ .
• 4:while  $[eqn]$ or $[eqn]$  do
• 5:    define a variable $[eqn]$ to represent the cost reduction due to partitioning, initialized to 0;
• 6:    for each $[eqn]$ of the k convex polygons $[eqn]$ in S do
• 7:        identify the line $[eqn]$ that partitions $[eqn]$ into two convex polygons, $[eqn]$ and $[eqn]$
• 8:        calculate the reduction in area, $[eqn]$ , resulting from partitioning the set $[eqn]$ ;
• 9:    end for
• 10:     $[eqn]$ ;
• 11:     $[eqn]$ ;
• 12:     $[eqn]$ ;
• 13:    update $[eqn]$ if the tree of $[eqn]$ increases its depth by 1 due to the addition of $[eqn]$ and $[eqn]$ ;
• 14:    update S as $[eqn]$ ;
• 15:end while
• 16:return S

Below, we present the running time complexity of K-AC in Theorem 1.

Theorem 1. The K-AC algorithm runs in $[eqn]$ time.

Proof. Using Graham’s scan algorithm, line 1 of K-AC runs in $[eqn]$ time. Lines 2 and 3 take constant time. The WHILE loop runs for K rounds. In each round, we identify $[eqn]$ straight lines, where $[eqn]$ . For each straight line, it takes linear time $[eqn]$ to locate $[eqn]$ and $[eqn]$ using Equations (1)–(3). The area calculations require $[eqn]$ time. Identifying the best partition to reduce the total area size takes $[eqn]$ time on lines 10 and 11. Line 12 takes constant time, while Line 13 takes $[eqn]$ time. Updating the clusters in line 14 takes $[eqn]$ time, using Graham’s scan algorithm. Thus, the total running time of the WHILE loop is

$$[eqn]$$ $$[eqn]$$

Inequality (5) holds by Jensen’s inequality, and Equation (6) holds under the assumption that $[eqn]$ . □

From Theorem 1, we observe that the K-AC algorithm scales efficiently with the total sample size N, highlighting both its theoretical and practical efficiency.

3. Experiments

In this section, we conduct experiments to evaluate our K-AC algorithm. Using real biological data, we compare K-area clustering with other well-known clustering algorithms, including K-center, K-median, and K-means clustering. The experiments are conducted on a MacBook Air featuring a 2.3 GHz Dual-Core Intel Core i5 processor and 16 GB of 2133 MHz LPDDR3 memory. This MacBook Air Apple M1 was manufactured by Apple in partnership with Foxconn and Pegatron. Apple designs the processors, but Taiwan Semiconductor Manufacturing Co. (TSMC) (Hsinchu, Taiwan) builds them.

3.1. Datasets

To evaluate the clustering performance of the K-AC algorithm, we use real scRNA-seq datasets from human and mouse cells (also listed in Table 2). These datasets are downloaded from an online repository (https://hemberg-lab.github.io/scRNA.seq.datasets/, accessed on 5 January 2025) and were previously used to assess scRNA-seq clustering tools [10]. The optimal number of clusters is determined based on the true cell labels provided by the original authors’ annotations.

A pre-processing step involves applying the PCA analysis to reduce the dimensionality of the data and generate values representing them in a 2D space. The use of PCA to highlight information relevant to the classification of genes or cell lines has been studied by Crescenzi and Giuliani [26] in the context of microarray analysis. Their work uncovers mutual similarities from cluster analysis among differently labeled statistical units, generating new classifications at varying levels of detail that can be interpreted in biological terms. The PCA results are shown in Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6.

3.2. Metrics

We set $[eqn]$ for the Biase, Deng, Goolam, Ting, and Yan datasets, respectively. We then measure the cost of the following clustering algorithms: the K-center clustering algorithm [29], the K-median clustering algorithm [30], the K-means Lloyd’s clustering algorithm [21] and our proposed K-area clustering algorithm K-AC. To evaluate the performance of these clustering algorithms, we use Normalized Mutual Information (NMI) [31] to assess the quality of the generated clusters.

NMI is a widely used metric for evaluating clustering performance by measuring the mutual dependence between the true labels and the predicted clusters. The mutual information is normalized to ensure values between 0 and 1, where 1 indicates perfect clustering and 0 implies no mutual information between the clustering and ground truth. The NMI between two cluster assignments C (ground truth) and $[eqn]$ (predicted clusters) is computed as

$$[eqn]$$

where $[eqn]$ is the mutual information between C and $[eqn]$ , and $[eqn]$ and $[eqn]$ are the entropy of C and $[eqn]$ , respectively.

3.3. K-Center Clustering Algorithm

The K-center clustering algorithm aims to find a partition $[eqn]$ of the data points into K clusters, with corresponding centers $[eqn]$ such that the maximum distance between any data point and the center of its assigned cluster is minimized. Specifically, the goal is to minimize

$$[eqn]$$

When K is not a fixed input or may vary as a function of the number of data points N, this K-center clustering problem becomes NP-hard. In our experiments, we apply the farthest-traversal algorithm to solve the K-center problem and generate clusters. The resulting clustering is shown in Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11.

3.4. K-Median Clustering Algorithm

The K-median clustering algorithm seeks to find a partition $[eqn]$ of the data points into K clusters, with corresponding centers $[eqn]$ to minimize the total distance between each data point and the center of its assigned cluster. Specifically, the objective is to minimize

$$[eqn]$$

When K is not part of an input or may be a function of n, then this K-median clustering problem is NP-hard. In our experiments, we run the Lloyd-style iteration algorithm [32] for the K-median clustering problem to generate clusters. The resulting clustering is shown in Figure 12, Figure 13, Figure 14, Figure 15 and Figure 16.

3.5. K-Means Clustering Algorithm

The K-means clustering algorithm aims to find a partition $[eqn]$ of the data points into K clusters, with corresponding centers $[eqn]$ such that the sum of squared distances between each data point and the center of its assigned cluster is minimized. Specifically, the objective is to minimize

$$[eqn]$$

When K is not part of an input or may be a function of n, then this K-means clustering problem is NP-hard. In our experiments, we run Llyod’s algorithm [21] for the K-means clustering problem to generate clusters. The resulting clustering is shown in Figure 17, Figure 18, Figure 19, Figure 20 and Figure 21.

3.6. K-Area Clustering Algorithm

The K-area clustering algorithm aims to find a partition $[eqn]$ of the data points into K clusters such that the sum of convex hull areas is minimized. Specifically, the objective is to minimize

$$[eqn]$$

where $[eqn]$ denotes the area of the region defined by the polygon $[eqn]$ .

In our experiments, we run our algorithm K-AC for the K-area clustering problem to generate clusters. The resulting clustering is shown in Figure 22, Figure 23, Figure 24, Figure 25 and Figure 26.

3.7. Performance Comparison

We summarize the NMI scores for these algorithms in Table 3. Our algorithm, K-AC, outperforms all other algorithms in three cases but is slightly less effective in two: the K-median algorithm on the Deng dataset and on the Ting dataset. We also evaluate the results from a biological perspective by utilizing the experimental labels provided by the original authors’ analysis. Across all datasets, the K-AC algorithm demonstrates a remarkable ability to accurately group outlier cells, which exhibits significant distances from the majority of cells within the same cell type. For example, in the Biase dataset [33], K-AC correctly clusters one outlier cell (a zygotic cell highlighted in the red dotted circle in Figure 22), while the other three methods incorrectly group it with its neighboring cluster (see Figure 7, Figure 12 and Figure 17). By identifying these outlier cells, the K-AC algorithm provides valuable insights into cellular heterogeneity, potentially leading to novel biological discoveries. These results underscore that K-AC offers a biologically relevant and geometrically interpretable clustering criterion, enabling more meaningful clusterings of high-dimensional biological data compared to the other three algorithms.

We report the algorithms’ running time and the convex areas produced by these four clustering algorithms in Table 4 and Table 5, respectively. For some instances, the K-AC algorithm runs faster than the other clustering algorithms; however, when the value of K is large, the K-AC algorithm runs slower. It is evident that the K area algorithm results in the smallest convex area, and in some cases, such as the Goolam dataset, it reduces the space by up to 60%. The comparison results are visualized in Figure 27 and Figure 28.

3.8. Hierarchical Clustering

To demonstrate the effectiveness of the K-area algorithm in interpreting the hierarchical structure of biological data—particularly in identifying the number of natural clusters—we conduct experiments to explore the construction of the cluster tree. We use the ratio of the convex hull areas between two neighboring rounds as a criterion for determining when to halt the partitioning of existing clusters.

The dataset Yan used in our experiments has an optimal number of clusters, $[eqn]$ , as provided by the authors’ original annotation [34]. We calculate the total cluster areas for $[eqn]$ and summarize the results in Table 6, with visualizations in Figure 29. The area ratio between two neighboring rounds increases from 0 to 1. As the ratio approaches 1, the benefit of further partitioning into additional clusters diminishes. Notably, after $[eqn]$ , the area ratio stabilizes near 1, indicating that further partitioning yields diminishing returns. This feature of the K-area algorithm aids in identifying an appropriate value for K, the number of clusters.

We construct the cluster tree generated by the K-area algorithm to visually represent its progression for the dataset Yan from Figure 30, Figure 31, Figure 32, Figure 33, Figure 34, Figure 35 and Figure 36. The tree structure of the partitions is shown in Figure 37. In this tree, the root node represents the initial convex hull, with each parent node (a cluster) being divided into two child nodes (subclusters). The leaves of the tree correspond to the final set of clusters. The height of the tree is constrained by the given hierarchical level H. By examining the authors’ labels for the cell types [34], we observe that the hierarchical tree organization aligns well with biological relationships. Specifically, C1 (4-cell embryo) is closely related to C2 (8-cell embryo), and they are initially merged, reflecting their sequential developmental stages. The combined group of C1 and C2 then merges with C3, which consists of a heterogeneous mixture of Oocyte, Zygote, and 2-cell embryo cells, indicating their shared early embryonic origin. C4 and C5 represent two distinct groups of Morulae cells, maintaining their biological distinction. Additionally, C6 (Late blastocyst) and C7 (hESC passage) cluster together, consistent with their advanced differentiation stages. This biologically meaningful clustering validates the hierarchical structure inferred by the K-AC algorithm.

4. Summary and Discussion

A key innovation of the K-volume algorithm is its ability to simultaneously optimize both the hierarchical structure and the number of clusters at each level, leveraging nonlinear optimization. Unlike traditional clustering methods, which require a predefined number of clusters, our method dynamically determines the optimal number of clusters and hierarchical layers based on the convex volume distribution. This adaptive nature makes it particularly well suited for complex biological datasets, such as scRNA-seq and spatial transcriptomics, where cell populations exhibit multi-scale hierarchical relationships.

Moreover, the flexibility of our method extends beyond biological applications, making it a powerful tool for a wide range of data-intensive domains, including finance, social network analysis and image segmentation, where hierarchical patterns naturally emerge. The method holds significant potential as a core component in modern data analysis pipelines. By integrating this method into widely used computational frameworks, we envision it becoming a standard clustering technique for high-dimensional and hierarchical data, offering a versatile and scalable solution for data-driven discovery.

Looking ahead, we will investigate incorporating performance metrics into the clustering of datasets. This consideration arises from the observation that the K-area algorithm performs slightly worse than the K-median algorithm on two datasets (Table 1). We believe this is due to the metric used for degenerate polygons. We define the area of a degenerate polygon as 0, which may result in multiple lines with zero area, where the line endpoints are close to each other. In such cases, we may consider incorporating the K-median algorithm’s metric. The optimal combination of performance metrics for various data points is currently under investigation.

Our future research will focus on extending this framework by incorporating advanced high-dimensional data processing techniques, such as dimensionality reduction methods and deep learning-based feature extraction, to further enhance clustering performance [35,36]. We also plan to integrate multi-omics data, including single-cell spatial transcriptomics [37], epigenomics [38], and proteomics [39], to create a more holistic representation of cellular states and regulatory mechanisms. By combining these advancements, we aim to develop a powerful, scalable clustering framework capable of uncovering complex biological patterns across multiple layers of molecular information, ultimately contributing to a deeper understanding of cellular functions and disease mechanisms. Note that in [40], it was reported that cultured cells exhibit high correlations between the values of certain shape descriptors, with the geometrical features postulated to be responsible for each correlation. Additionally, these correlations show a complex dependence on the size, shape, and number of invaginations. We plan to further investigate the geometrical characteristics of various cells and provide more biological interpretations for the clusters generated based on our geometrical approach.

5. Conclusions

In this paper, we address the challenge of clustering biological data by developing a novel algorithm based on the total convex volume encompassed by points within the same cluster. This approach offers a biologically relevant and geometrically interpretable clustering criterion, enabling more meaningful groupings of high-dimensional biological data. We introduce the K-volume algorithm, which effectively balances accuracy and computational feasibility. Experimental evaluations demonstrate the algorithm’s effectiveness compared to well-established clustering methods, showing promising results in capturing biological structures and relationships within the data.