Abstract:Advances in technology have made it convenient to obtain a large amount of single cell RNA sequencing (scRNA-seq) data. Since that clustering is a very important step in identifying or defining cellular phenotypes, many clustering approaches have been developed recently for these applications. The general methods can be roughly divided into normal clustering methods and integrated (ensemble) clustering methods which combine more than two normal clustering methods aiming to get much more informative performance. In order to make a contrast with the integrated clustering algorithm, the normal clustering method is often called individual or base clustering method. Note that the results of many individual clustering methods are often developed to capture one aspect of the data, and the results depend on the initial parameter settings, such as cluster number, distance metric and so on. Compared with individual clustering, although integrative clustering method may get much more accurate performance, the results depend on the base clustering results and integrated systems are often not self-regulation. Therefore, how to design a robust unsupervised clustering method is still a challenge. In order to tackle above limitations, we propose a novel Ensemble Clustering algorithm based on Probability Graphical Model with Graph Regularization, which is called EC-PGMGR for short. On one hand, we use parameter controlling in Probability Graphical Model (PGM) to automatically determine the cluster number without prior knowledge. On the other hand, we add a regularization term to reduce the effect deriving from some weak base clustering results. Particularly, the integrative results collected from base clustering methods can be assembled in the form of combination with self-regulation weights through a pre-learning process, which can efficiently enhance the effect of active clustering methods while weaken the effect of inactive clustering methods. Experiments are carried out on 7 data sets generated by different platforms with the number of single cells from 822 to 5,132. Results show that EC-PGMGR performs better than 4 alternative individual clustering methods and 2 ensemble methods in terms of accuracy including Adjusted Rand Index (ARI) and Normalized Mutual Information (NMI), robustness, effectiveness and so on. EC-PGMGR provides an effective way to integrate different clustering results for more accurate and reliable results in further biological analysis as well. It may provide some new insights to the other applications of clustering.

Graph-Aligned Random Partition Model (GARP)

Partitioning gene expression data by data-driven Markov chain Monte Carlo

Non-exchangeable random partition models for microclustering

Bayesian Analysis for Exponential Random Graph Models Using the Adaptive Exchange Sampler.

A bayesian multivariate mixture model for high throughput spatial transcriptomics

Bayesian Robust Learning in Chain Graph Models for Integrative Pharmacogenomics

Bayesian mixtures of common factor analyzers: Model, variational inference, and applications

Bayesian graphical models for computational network biology

Particle Gibbs for Bayesian Additive Regression Trees

Parallelizing MCMC with Random Partition Trees

Adaptive Variable Clustering in Gaussian Graphical Models.

Model-based Clustering with Sparse Covariance Matrices

Simultaneous Bayesian Clustering and Model Selection with Mixture of Robust Factor Analyzers

Model-based clustering based on sparse finite Gaussian mixtures

Incorporating Grouping Information into Bayesian Gaussian Graphical Model Selection

Network clustering analysis using mixture exponential-family random graph models and its application in genetic interaction data.

Bayesian inference of mixed Gaussian phylogenetic models

EC-PGMGR: Ensemble Clustering Based on Probability Graphical Model With Graph Regularization for Single-Cell RNA-seq Data

A Bayesian Approach To Graph Partitioning

Inferring Hierarchical Mixture Structures: A Bayesian Nonparametric Approach

Particle-Gibbs Sampling For Bayesian Feature Allocation Models