Tutorials › Multi-Omics Data Integration › Shared Sample Allocations Across Omics Blocks

Bayesian Integrative Clustering · Tutorial 262 of 597

Shared Sample Allocations Across Omics Blocks

Learn how a single latent allocation can link matched samples across multiple omics data blocks without requiring their cluster-specific profiles to be identical.

Advanced 10 min read

What You'll Learn

  • Define a shared allocation for matched samples across several omics blocks
  • Distinguish a shared partition from identical measurements or identical profiles
  • Calculate a conditional allocation probability from prior weights and joint likelihoods
  • Recognize label switching and compare partitions rather than raw label numbers
  • Diagnose sample misalignment, weak evidence, and block-dominance risks
  • Compare shared allocations with separate clustering and simple concatenation

Why Share Sample Allocations?

A multi-omics study may measure RNA abundance, protein levels, and metabolites on the same biological samples. One scientific question is whether those samples can be grouped into a common set of latent subpopulations, with each sample's group membership linking its measurements across data types. The central modeling choice is then an allocation: one latent cluster label for each matched sample, used across all the blocks.

Sharing labels does not mean that the blocks contain the same features, have the same numerical scale, or show the same cluster-specific profiles. It means that a sample assigned to cluster \(k\) in one block is assigned to cluster \(k\) in every block. This is a statement about how sample memberships are organized, not a claim that the biological measurements are identical or that one block causes another.

Definition: Suppose \(M\) omics blocks are observed on the same \(n\) matched samples. Write \(X^{(m)}\) for the data in block \(m\), and \(x_i^{(m)}\) for the profile of sample \(i\) in that block. A shared allocation is a vector \(z=(z_1,\ldots,z_n)\), where \(z_i\in\{1,\ldots,K\}\) is the latent cluster label assigned to sample \(i\), with that same \(z_i\) applying to \(x_i^{(m)}\) for every block \(m\). The induced partition groups samples by their labels; the numeric names of the labels are arbitrary.

The matched-sample condition matters. Row \(i\) in every block must refer to the same biological unit under a compatible sampling definition. If one block has samples in a different order, the rows must be aligned before a shared allocation is used. A mistaken match can make unrelated profiles appear to support the same cluster assignment.

What the Shared Label Does—and Does Not—Constrain

A joint clustering model represents the observed data and the allocation together. In general terms, it has a distribution for the allocation and a joint data model conditional on that allocation and model parameters:

$$ p(z, X^{(1)},\ldots,X^{(M)}\mid\theta) = p(z\mid\theta)\, p(X^{(1)},\ldots,X^{(M)}\mid z,\theta). $$

Here \(\theta\) denotes model parameters. The shared allocation makes each sample's label common across blocks. It does not, by itself, say that the blocks are conditionally independent given \(z\), nor that their profiles have the same form. Those are additional modeling choices, not consequences of sharing labels. Block-specific parameters can describe different features and different cluster patterns within a common grouping of samples.

For a simple finite allocation model, one may assume \(z_i\) are conditionally independent given cluster probabilities \(\pi_1,\ldots,\pi_K\), with \(P(z_i=k\mid\pi)=\pi_k\) and \(\sum_{k=1}^K\pi_k=1\). This is one possible allocation model, not part of the definition of a shared allocation. Other allocation models can be used. In either case, labels summarize a model-based grouping; the data and scientific question determine whether that grouping is useful.

Worked Example: One Label for Two Matched Blocks

Consider four hypothetical samples measured in a continuous-expression block and a protein block. The entries below are simplified one-feature summaries for illustration; they are not real study data.

SampleExpression summaryProtein summary
A2.118
B2.420
C7.843
D8.246

Suppose a model assigns the shared allocation \(z=(1,1,2,2)\). Then samples A and B belong to cluster 1 in both blocks, while C and D belong to cluster 2 in both blocks. The protein measurements are not being equated with the expression measurements: their units and scales differ. The common claim is only that the two blocks use the same sample grouping.

Before interpreting the grouping, verify that the expression row for A and the protein row for A really refer to the same sample. If, for example, the protein rows were accidentally ordered B, A, C, D but not reordered, the shared labels would link the wrong profiles. A plausible-looking allocation cannot repair incorrect sample matching.

How Joint Evidence Informs an Allocation

The allocation is latent, so its uncertainty is represented by probabilities over possible labels or complete allocations. For a particular sample \(i\), consider its label while holding the other labels and model parameters fixed. The conditional probability of assigning it to cluster \(k\) is proportional to the prior allocation weight times the likelihood of the observed data under that proposed assignment:

$$ P(z_i=k\mid X^{(1)},\ldots,X^{(M)},z_{-i},\theta) \propto P(z_i=k\mid z_{-i},\theta)\, p(X^{(1)},\ldots,X^{(M)}\mid z_i=k,z_{-i},\theta). $$

The second term is the likelihood contribution of the observed data across the blocks under the proposed assignment. It is the joint likelihood unless the model makes and justifies a factorization into block-specific terms. Thus, a block can support or weaken an assignment through its contribution to the joint model, but the shared-label definition alone does not license multiplying separate block likelihoods.

Worked Example: Normalizing Conditional Allocation Weights

For a hypothetical sample, suppose there are two possible labels. Conditional on the other labels and parameters, let the prior allocation weights be \(0.60\) for cluster 1 and \(0.40\) for cluster 2. Suppose the model's joint likelihood for the observed profiles across all blocks is \(0.18\) under cluster 1 and \(0.06\) under cluster 2. These likelihood values are supplied by the model; no assumption that block likelihoods factor is needed here.

The unnormalized conditional weights are

$$ w_1=0.60(0.18)=0.108,\qquad w_2=0.40(0.06)=0.024. $$

Their sum is \(0.108+0.024=0.132\). Therefore,

$$ P(z_i=1\mid X,z_{-i},\theta) = \frac{0.108}{0.132} \approx 0.8182,\qquad P(z_i=2\mid X,z_{-i},\theta) = \frac{0.024}{0.132} \approx 0.1818. $$

Under these stated conditions, cluster 1 has the larger conditional allocation probability. The value \(0.8182\) is conditional on the specified model, the other labels, and the parameters; it is not a probability that the model or cluster 1 is true. It also does not show which block contributed more unless the joint likelihood is further examined under an appropriate model.

A single most probable label can conceal uncertainty. If two or more labels have similar posterior probabilities, reporting only the highest-probability label overstates how clearly the sample is assigned. Summaries should preserve uncertainty when it is material, for example by reporting posterior allocation probabilities or uncertainty in cluster membership, rather than treating a single estimated allocation as certain.

Partitions, Label Names, and Comparing Blocks

Cluster labels are names, not measurements. If every occurrence of label 1 is changed to 2 and every occurrence of 2 to 1, the grouping of samples is unchanged. A useful way to describe a partition is through pairwise co-clustering indicators:

$$ C(i,j)=\mathbf{1}(z_i=z_j), $$

where \(C(i,j)=1\) means samples \(i\) and \(j\) share a cluster, and \(C(i,j)=0\) means they do not. Under one shared allocation, the indicator is the same for every block because the labels \(z_i\) and \(z_j\) are shared. This gives a direct way to check whether separately displayed labels represent the same grouping.

1
Write each allocation against sample identities.
Use the same sample order in every block and verify the matching before comparing labels.
2
Compare which pairs are grouped.
For each pair, check whether its samples have equal labels; do not infer agreement merely because the numeric labels look similar.
3
Account for label permutations.
If two allocations group exactly the same samples but use different label numbers, they describe the same partition after relabeling.

Worked Example: Counting Pairwise Partition Disagreements

For four hypothetical samples, compare allocation A, \(z^{(A)}=(1,2,1,2)\), with allocation B, \(z^{(B)}=(1,2,2,1)\). There are \(\binom{4}{2}=6\) unordered sample pairs. In allocation A, the same-cluster pairs are \((1,3)\) and \((2,4)\). In allocation B, they are \((1,4)\) and \((2,3)\).

PairSame cluster in A?Same cluster in B?Agreement?
(1,2)NoNoYes
(1,3)YesNoNo
(1,4)NoYesNo
(2,3)NoYesNo
(2,4)YesNoNo
(3,4)NoNoYes

The two allocations disagree on four of the six pairs, so their pairwise disagreement proportion is \(4/6\approx0.6667\). They agree that pairs (1,2) and (3,4) are separated, but assign the other four pair relationships differently. This comparison is unaffected by a permutation of the label names within either allocation.

Justification: Shared Labels Imply a Shared Partition

The key implication follows directly from the definition. For any blocks \(m\) and \(r\), and any two samples \(i\) and \(j\), a shared allocation uses the same labels in both blocks. Therefore,

$$ \mathbf{1}(z_i^{(m)}=z_j^{(m)}) = \mathbf{1}(z_i=z_j) = \mathbf{1}(z_i^{(r)}=z_j^{(r)}). $$

The first and last indicators must agree because \(z_i^{(m)}=z_i=z_i^{(r)}\) and \(z_j^{(m)}=z_j=z_j^{(r)}\). Hence every block induces the same pairwise grouping. Conversely, if every block's allocation induces the same partition, its numeric labels can be permuted to represent a common allocation, provided the blocks contain the same matched samples. This is why partition comparisons, rather than raw label-number comparisons, are appropriate.

Alternatives, Diagnostics, and Limits

A shared allocation is useful when the scientific target is a common sample grouping supported by several data types. It is not the only reasonable strategy. The main alternatives answer somewhat different questions:

ApproachWhat it representsMain caution
Shared allocationOne sample partition used across all blocksCan be restrictive if meaningful groupings differ by data type
Separate clusteringA distinct partition estimated for each blockDoes not impose a common grouping; label numbers across fits do not automatically correspond
Concatenated featuresOne clustering of a combined feature matrixFeature counts, scales, and noise can let one block dominate unless handled appropriately

A shared allocation should be assessed for both data integrity and model adequacy. Check sample identifiers and row order first. Then examine whether assignments are uncertain, whether the proposed grouping is sensitive to modeling choices, and whether one block appears to dominate the joint evidence. Posterior predictive checks and residual diagnostics can help assess model fit, but a visually separated partition alone does not establish that the model describes each block adequately.

Several limitations remain even when the data are aligned. A common partition may be too restrictive if different assays reveal genuinely different subgroups. A shared label does not guarantee that every block contains information about every cluster, and weakly informative blocks may contribute little. Conversely, large or highly variable blocks may have disproportionate influence under an unsuitable model. Finally, the number of clusters and the model for block-specific profiles affect the inferred allocation; labels should be interpreted with those choices in view.

Common Mistakes

  • Reading label numbers as biological identities. Labels 1 and 2 can be swapped without changing a partition. Compare grouped sample pairs, not the numbers alone.
  • Assuming identical profiles from shared labels. A shared allocation links sample memberships; block-specific features, scales, and profiles can still differ.
  • Multiplying block evidence without a model justification. A product of block likelihoods requires a factorization assumption. Use the joint likelihood unless that assumption is stated and appropriate.
  • Ignoring sample matching. Shared labels are meaningful only when each row is correctly linked to the same sample across blocks.
  • Treating the most probable allocation as certain. A hard label hides uncertainty when posterior allocation probabilities are spread across clusters.
Key takeaway: A shared sample allocation assigns one latent label to each matched sample and uses that label in every omics block. It enforces a common partition, but does not by itself impose conditional independence, identical profiles, or certainty about membership.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions.

  1. Three blocks have rows for the same 20 samples. What must be true before one shared allocation can be interpreted across them?
  2. Two allocations have different numeric labels but exactly the same same-cluster sample pairs. Do they represent different partitions? Explain.
  3. A sample has prior allocation weights 0.50 and 0.50 for two clusters. Its joint likelihood is 0.12 under cluster 1 and 0.08 under cluster 2. Calculate the normalized conditional allocation probabilities.
  4. Does a shared allocation imply that measurements in the blocks are conditionally independent given the labels? Explain.
  5. What is one reason to prefer separate clustering over a shared allocation, and one caution when simply concatenating features?

Solutions. (1) The rows must be aligned so that the same row index identifies the same matched biological sample in all three blocks. (2) No. The partitions are the same; a permutation of label names does not change which sample pairs are grouped. (3) The unnormalized weights are \(0.50(0.12)=0.06\) and \(0.50(0.08)=0.04\); the total is \(0.10\). The probabilities are \(0.06/0.10=0.60\) for cluster 1 and \(0.04/0.10=0.40\) for cluster 2. (4) No. Sharing labels constrains sample memberships, not the conditional dependence structure of the measurements. (5) Separate clustering may be preferable when blocks are expected to reveal different groupings. With concatenation, blocks with more features or different scales may dominate the combined analysis.

Clinical Biostats