--- title: Multi-layer Stochastic Co-block Model url: https://www.emergentmind.com/topics/multi-layer-stochastic-co-block-model type: topic --- # Multi-layer Stochastic Co-block Model Searching arXiv for recent and directly relevant papers on multi-layer stochastic co-block models and neighboring multilayer block-model formulations. Multi-layer stochastic co-block model (ML-ScBM) denotes a class of latent block models for multi-layer directed or bipartite networks in which each layer has its own block connectivity matrix, but the latent structure is asymmetric: nodes are assigned separate sender/row and receiver/column communities, and the two community counts may differ. In the directed aligned-node setting, a layer-\(\ell\) edge satisfies \(\mathbb P(A^{(\ell)}_{ij}=1)=B^{(\ell)}(g^s(i),g^r(j))\); in the bipartite degree-corrected setting, \(\Omega_l=\Theta_r Z_r B_l Z_c^\top \Theta_c\). The model therefore differs from a standard multilayer SBM, which uses one common partition, by explicitly representing outgoing and incoming patterns or row and column communities separately [2602.21569][2307.10572][2405.04711]. ## 1. Core concept and asymmetric latent structure The defining feature of ML-ScBM is the separation between sending and receiving roles. In the formulation for multi-layer directed networks, one observes \(L\) directed layers on the same \(n\) nodes, with adjacency matrices \[ A^{(\ell)} \in \{0,1\}^{n\times n}, \qquad \ell=1,\dots,L,\quad A^{(\ell)}_{ii}=0, \] and \[ \mathbb{P}\!\left(A^{(\ell)}_{ij}=1\right)=B^{(\ell)}\!\bigl(g^s(i),g^r(j)\bigr), \qquad i\neq j. \] Here \(g^s\) gives sender community labels, \(g^r\) gives receiver community labels, and \(B^{(\ell)}\in[0,1]^{K_s\times K_r}\) is the layer-specific block probability matrix. The key structural point is that a node may send edges like one community but receive edges like another, so \(K_s\neq K_r\) is natural rather than exceptional [2602.21569]. In the directed co-clustering formulation of spectral community detection, the same asymmetry is written in row/column notation. There are \(K_y\) row clusters and \(K_z\) column clusters, membership matrices \(Y\in\{0,1\}^{n\times K_y}\) and \(Z\in\{0,1\}^{n\times K_z}\), and layer-specific block matrices \(B_l\in[0,1]^{K_y\times K_z}\). The model is \[ A_{l,ij} \sim \mathrm{Bernoulli}\big(\rho\, B_{l,g_i^y g_j^z}\big), \qquad i\neq j, \] with population matrix \[ \mathcal P_l = \rho Y B_l Z^\top. \] In this representation, row clusters encode who sends similarly, while column clusters encode who receives similarly [2307.10572]. For multi-layer bipartite networks, the asymmetry is attached to two distinct node sets rather than to two roles played by the same node set. The degree-corrected model uses \(n_r\) row nodes, \(n_c\) column nodes, row membership matrix \(Z_r\), column membership matrix \(Z_c\), degree parameters \(\theta_r,\theta_c\), and layer-specific block matrices \(B_l\). The model is \[ A_l(i_r,j_c)\sim \mathrm{Bernoulli}\!\big(\Omega_l(i_r,j_c)\big), \qquad \Omega_l:=\Theta_r Z_r B_l Z_c^\top \Theta_c. \] This formulation reduces to the degree-corrected stochastic co-block model when \(L=1\), to the stochastic co-block model when the degree parameters are constant, and to multilayer SBM when the network is undirected and \(Z_r=Z_c\) [2405.04711]. ## 2. Canonical multilayer formulations The directed ML-ScBM and the multilayer degree-corrected stochastic co-block model share the assumption that layers are conditionally independent given latent communities, but they organize asymmetry differently. In the aligned-node directed setting, the same node set appears on both sides of each adjacency matrix, and asymmetry comes from distinct sending and receiving clusterings. In the bipartite setting, each layer is an \(n_r\times n_c\) matrix, and asymmetry is anchored in separate row and column node types [2307.10572][2405.04711]. Population geometry is central to both formulations. For directed co-clustering, the row-side and column-side aggregate population matrices are \[ \sum_{l=1}^L \mathcal P_l\mathcal P_l^\top = \rho^2 Y \left(\sum_{l=1}^L B_l \Delta_z^2 B_l^\top\right)Y^\top, \] and \[ \sum_{l=1}^L \mathcal P_l^\top\mathcal P_l = \rho^2 Z \left(\sum_{l=1}^L B_l^\top \Delta_y^2 B_l\right)Z^\top. \] The leading eigenvector matrices of these population operators are \[ U = Y\Delta_y^{-1}Q^R, \qquad V = Z\Delta_z^{-1}Q^C. \] If the aggregated block matrix is full rank, distinct clusters correspond to distinct rows of \(U\) or \(V\); if it is rank deficient, cluster recovery still remains possible under explicit separation conditions. This relaxes the common assumption that the embedding rank must equal the number of clusters [2307.10572]. For the multilayer degree-corrected stochastic co-block model, the key population matrices are \[ \tilde S_r := \sum_{l\in[L]} \Omega_l\Omega_l^\top, \qquad \tilde S_c := \sum_{l\in[L]} \Omega_l^\top\Omega_l. \] A key lemma states that, under suitable rank conditions, the row-normalized eigenvectors of these matrices have exactly one distinct row per community. In particular, if \(\mathrm{rank}(\sum_l B_lB_l^\top)=K\le K_r\), then \[ U_r^*(i_r,:)=U_r^*(\tilde i_r,:)\iff Z_r(i_r,:)=Z_r(\tilde i_r,:), \] and when \(K=K_r\), different row communities are separated by \(\sqrt{2}\) in the normalized eigenspace [2405.04711]. These formulations show that ML-ScBM is not defined merely by multiple adjacency matrices. Its distinctive object is a multilayer collection of block operators with separate latent structures on the two sides of interaction. ## 3. Spectral co-clustering and debiased aggregation A major algorithmic theme in ML-ScBM is that naive layer aggregation can be misleading. In the directed setting, direct summation \(\sum_l A_l\) can suffer from community cancellation: different layers may offset each other and destroy cluster signal. To avoid this, the spectral method developed in [2307.10572] aggregates Gram matrices rather than adjacency matrices: \[ \sum_{l=1}^L A_lA_l^\top \quad \text{for row clusters}, \qquad \sum_{l=1}^L A_l^\top A_l \quad \text{for column clusters}. \] Because diagonal terms create bias, the method subtracts degree diagonals and forms the debiased sum-of-Gram matrices \[ S^R = \sum_{l=1}^L \big(A_lA_l^\top - D_l^{out}\big), \qquad S^C = \sum_{l=1}^L \big(A_l^\top A_l - D_l^{in}\big). \] The algorithm, called DSoG, computes leading eigenvectors of \(S^R\) and \(S^C\), then applies \(k\)-means to their rows. It also allows the embedding ranks \(K\le K_y\) and \(K'\le K_z\), so rank need not equal the number of sender or receiver communities [2307.10572]. The bipartite analogue is NcDSoS, normalized spectral co-clustering based on debiased sum of squared matrices. Its observed operators are \[ S_r=\sum_{l\in[L]}(A_lA_l^\top - D_l^r), \qquad S_c=\sum_{l\in[L]}(A_l^\top A_l - D_l^c), \] followed by spectral decomposition, row normalization, and \(K\)-means on the normalized eigenvector rows. The paper explicitly contrasts NcDSoS with non-debiased variants, arguing that the diagonal degree terms inflate the spectral norm error and that subtracting \(D_l^r\) and \(D_l^c\) gives a tighter approximation to the population matrices [2405.04711]. These constructions indicate that second-order aggregation is structurally important in multilayer asymmetric networks. This suggests that, in ML-ScBM, the informative signal often lives more naturally in layer-aggregated covariance-type operators than in raw adjacency sums. ## 4. Estimating the numbers of sender and receiver communities A distinctive inferential problem in ML-ScBM is that the number of sender communities and the number of receiver communities must be estimated separately. In [2602.21569], this is formulated as a sequence of goodness-of-fit tests for candidate pairs \((K_{s0},K_{r0})\): \[ H_0:\ (K_s,K_r)=(K_{s0},K_{r0}) \qquad\text{vs.}\qquad H_1:\ K_s>K_{s0}\ \text{or}\ K_r>K_{r0}. \] The alternative is underfitting: the candidate model merges true sender or receiver communities. The test is built from a normalized residual matrix. With estimated labels \(\hat g^s,\hat g^r\) and estimated block probabilities \(\hat B^{(\ell)}\), the empirical residual matrix is \[ \hat{R}(i,j)= \begin{cases} \dfrac{\sum_{\ell=1}^L\left(A^{(\ell)}(i,j)-\hat{\Omega}^{(\ell)}(i,j)\right)} {\sqrt{(n-1)\sum_{\ell=1}^L \hat{\Omega}^{(\ell)}(i,j)\bigl(1-\hat{\Omega}^{(\ell)}(i,j)\bigr)}}, & i\neq j,\\[1.2em] 0, & i=j, \end{cases} \] with \[ \hat T_n=\sigma_1(\hat R)-2. \] The largest singular value is used because it measures the strongest leftover low-rank structure after fitting the candidate model. Under the correct model, \(\hat T_n\) is near \(0\); under underfitting, \(\hat T_n\) diverges to infinity [2602.21569]. Two procedures are then built from this dichotomy. MLDiGoF searches candidate pairs in lexicographic order, first by increasing \(k_s+k_r\) and then by smaller \(k_s\), and stops at the first pair for which \(\hat T_n\tau_n\), with default \(\tau_n=8\log n\). Both procedures are proved to consistently recover \((K_s,K_r)\) under the multi-layer stochastic co-block model [2602.21569]. This testing framework addresses a limitation that remains explicit in several spectral co-clustering papers, namely that \(K_r\) and \(K_c\) are treated as known inputs [2405.04711]. ## 5. Statistical guarantees and the benefit of multiple layers A recurring theoretical conclusion is that additional layers improve recovery when the latent row/column structure is shared across layers and the block matrices vary by layer. In the directed spectral co-clustering analysis, the misclassification bounds for row and column clustering are \[ \mathcal L(Y,\widehat Y) \le \frac{c_4}{n\zeta_n^2} \left(\frac{1}{n^2}+\frac{\log(L+n)}{Ln^2\rho^2}\right), \] and \[ \mathcal L(Z,\widehat Z) \le \frac{c_{10}}{n\xi_n^2} \left(\frac{1}{n^2}+\frac{\log(L+n)}{Ln^2\rho^2}\right), \] under the stated assumptions. The bounds improve as \(L\) increases, as \(\rho\) increases, and as the eigenspace separation parameters \(\zeta_n,\xi_n\) increase. The paper emphasizes that multi-layers would bring benefits to the clustering performance, and that weak signal in one layer can be compensated by stronger signal in other layers [2307.10572]. For the multilayer degree-corrected stochastic co-block model, the central theorem yields consistency, and under balanced communities and homogeneous degree scaling the cleaner corollary gives \[ \hat f_r = O\!\left(\frac{\log(n_r+n_c+L)}{\rho^2 n_r n_c L}\right)+O\!\left(\frac1{n_r^2}\right), \] \[ \hat f_c = O\!\left(\frac{\log(n_r+n_c+L)}{\rho^2 n_r n_c L}\right)+O\!\left(\frac1{n_c^2}\right). \] If \(n_r=O(n_c)=O(n)\), then \[ \hat f_r=\hat f_c= O\!\left(\frac{\log(n+L)}{\rho^2 n^2 L}\right)+O\!\left(\frac1{n^2}\right). \] These rates make the role of \(L\) explicit: error decreases as the number of layers increases [2405.04711]. For model selection, [2602.21569] proves a sharp dichotomy for \(\hat T_n\): under the null, \[ \forall \epsilon>0,\qquad \mathbb{P}(\hat T_n<\epsilon)\to 1, \] while under underfitting, \[ \hat T_n \xrightarrow{P} \infty. \] The sequential procedures based on this dichotomy are consistent for the true \((K_s,K_r)\). Together, these results indicate that multilayer information improves inference in at least two distinct senses: it reduces clustering error for fixed community counts, and it supports separate estimation of sender and receiver community numbers. ## 6. Relation to neighboring multilayer block models ML-ScBM is often discussed alongside a broader multilayer block-model literature, but several neighboring models are not co-block models in the strict sense. The multilayer SBM of [1411.1098] treats an observed network as the aggregate of hidden interaction layers and allows different group assignments in different layers for the same node set, yet the asymmetry is across layers rather than across node types. The strata multilayer SBM of [1507.01826] clusters layers into strata that share a common SBM, coupling node-to-community and layer-to-stratum assignments, but it still uses SBM structure rather than separate row and column communities. The goodness-of-fit test in [2508.04957] and the connectivity-matrix inference theory in [2406.11152] concern standard multilayer SBM with common node memberships across layers, not a co-block model. The same distinction appears in mixed-membership and weighted multiplex settings. The multi-layer mixed membership stochastic block model of [2404.03916] uses a common mixed-membership matrix \(\Pi\) across layers and does not introduce distinct row/column memberships. The multiplex Dirichlet stochastic block model of [2412.11971] models directed compositional edge weights and has a clear co-block flavor through expected block-to-block shares \(v_{kh}^{(s)}\), but it uses a single shared partition of nodes across layers rather than separate sender and receiver clusterings. The partially exchangeable SBM of [2410.10619] is designed for node-colored multilayer networks and introduces a hierarchical partition prior coherent with within- and across-layer block-connectivity structures, yet its novelty lies in partial exchangeability rather than in stochastic co-blocking. This suggests a useful terminological boundary. “Multi-layer stochastic co-block model” is most precise when it refers to multilayer asymmetric block models with distinct latent structure on the two sides of interaction, as in directed sender/receiver models and bipartite row/column models [2602.21569][2307.10572][2405.04711]. By contrast, many multilayer SBMs are related in spirit because they share latent structure across layers, but they do not represent two-sided asymmetry and therefore occupy a different part of the multilayer block-model landscape. Source: https://www.emergentmind.com/topics/multi-layer-stochastic-co-block-model