Preprint | arXiv:2504.04079v2

Scalable Robust Bayesian Co-Clustering with Compositional ELBOs

Ashwin Vinod, Chandrajit Bajaj

A fully variational co-clustering framework that learns row clusters, column clusters, and cell-level interaction structure in one noise-robust training pipeline.

Noisy input

Co-clustered blocks

Problem

Co-clustering rows and features under real-world noise

Co-clustering identifies homogeneous groups of both instances and features in a matrix. Traditional methods often depend on shallow factorizations, linear transformations, or assumptions that struggle with sparse, high-dimensional, corrupted, or missing data.

Approach

Variational row, column, and joint latent structure

SRVCC uses row-side and column-side VAEs with Gaussian mixture priors, adds a joint cell-level latent variable, and optimizes a compositional objective with mutual information coupling and doubly reparameterized gradients.

Core contributions

What the paper adds

Fully variational co-clustering

Learns row and column clusters directly in latent space, replacing separate post-hoc clustering with a generative, end-to-end objective.

Compositional ELBOs

Combines row-side, column-side, and joint cell-level variational objectives so each part of the matrix contributes structured evidence.

Noise and missing-input robustness

Regularized reconstruction, KL structure, and noise learning are designed to keep cluster assignments stable under corrupted or sparse data.

Coherent row-column partitions

A mutual-information cross-loss encourages row and column assignments to preserve dependence in the original matrix.

Method

Processing pipeline

The framework turns a raw matrix into soft row clusters, soft column clusters, and joint block labels while retaining uncertainty throughout the training objective.

  1. Input matrix

    Rows are instances; columns are features.

  2. Row and column VAEs

    Separate encoders learn probabilistic latent spaces for each axis.

  3. GMM priors

    Mixture components provide soft cluster anchors in both latent spaces.

  4. Joint cell latent space

    A third latent variable models row-column interaction and local deviations.

  5. MI cross-loss

    Soft row and column partitions are coupled through mutual information.

  6. Co-cluster structure

    Rows and columns reorder into interpretable blocks.

Reported results

Benchmarks across image, text, and web data

Selected values from the paper tables. ACC is clustering accuracy; NMI is normalized mutual information.

DatasetSRVCC ACCSRVCC NMINote
Fashion-MNIST-test68.2 +/- 1.865.0 +/- 1.6Strong image benchmark gain over DeepCC
WebKB483.2 +/- 1.642.3 +/- 1.2Large accuracy gain on sparse text/web data
Yale58.1 +/- 1.761.0 +/- 1.5Improved facial-image latent separation
WebKB wisconsin81.6 +/- 2.251.5 +/- 1.6Best reported result in both ACC and NMI table columns
IMDb movies actors26.2 +/- 2.419.4 +/- 1.8Modest gain on a harder sparse movie-feature split

Ablation

Why the two-stage ELBO matters

The paper compares feature-only clustering, a simple cascade, and the two-stage compositional ELBO design. The two-stage version with DREG reports the strongest accuracy and NMI on the shown ablations.

Fashion-MNIST-test

Feature only

62.1 / 58.2

DREG helps, but without the joint stage it remains weaker.

Fashion-MNIST-test

Simple cascade

64.4 / 59.6

Better than feature-only, still below two-stage ELBO.

Fashion-MNIST-test

Two-stage ELBO

68.2 / 65.0

Best reported ACC and NMI for this ablation.

WebKB wisconsin

Two-stage ELBO

81.6 / 51.5

Best reported ACC and NMI under DREG.

Biomedical appendix

PPMI co-clustering as a Parkinson's exploration

The paper appendix applies the co-clustering approach to PPMI clinical and imaging variables. This section is exploratory and should be read as biomarker-discovery analysis, not diagnostic deployment.

Mild outlier profile

One near-singleton group shows low motor burden, high MoCA, early Hoehn and Yahr stage, and mild or near-normal imaging values.

Older heterogeneous group

An older cluster displays wider motor severity and comparatively preserved DaTSCAN signals, suggesting non-linear relationships among age, imaging, and symptoms.

Moderate PD subgroups

Several broader groups share moderate severity but differ by cognition, imaging spread, age at diagnosis, and ventricular or striatal measures.

Candidate biomarker families

Putamen and caudate DaTSCAN ratios, ventricular and striatal MRI volumes, UPDRS-III, Hoehn and Yahr stage, and MoCA emerge as separating variables.

Resources

Paper and citation

@article{vinod2025scalable,
  title={Scalable Robust Bayesian Co-Clustering with Compositional ELBOs},
  author={Vinod, Ashwin and Bajaj, Chandrajit},
  journal={arXiv preprint arXiv:2504.04079},
  year={2025}
}
Download PDF