Source-linked AI summary
Entropy of stochastic blockmodel ensembles
Tiago P. Peixoto
TL;DR
Existing blockmodel log-likelihood expressions may not directly yield the ensemble entropy. The paper derives entropy formulas across blockmodel variants and uses them to define an inference likelihood, which works for intrinsic and extrinsic degree correlations in the reported settings.
Problem
Existing blockmodel log-likelihood expressions can be non-closed or omit terms, so they do not directly yield the ensemble entropy.
Source: §I · p. 1Method
The paper derives microcanonical entropy expressions for traditional and degree-corrected ensembles with soft or hard degree constraints, then uses entropy to define a block-inference log-likelihood.
Source: §I · p. 1, §II · p. 2, §VII · p. 14Results
The inference estimators work with intrinsic and extrinsic degree correlations, and perform better than the sparse approximation on networks with broad degree distributions.
Source: §I · p. 2, §VII · p. 14Takeaways & Limitations
The entropy-derived likelihood can be used to detect a likely block partition in observed network data, including settings with intrinsic or extrinsic degree correlations.
Source: §I · p. 2, §VII · p. 14The soft degree-corrected approximation assumes intrinsic degree correlations between block pairs can be neglected, yet may still omit entropic effects even when that assumption holds.
Source: B. Degree-corrected ensembles with… · p. 5, A. Intrinsic degree correlations · p. 12
Abstract
from arXiv · showhide
Stochastic blockmodels are generative network models where the vertices are separated into discrete groups, and the probability of an edge existing between two vertices is determined solely by their group membership. In this paper, we derive expressions for the entropy of stochastic blockmodel ensembles. We consider several ensemble variants, including the traditional model as well as the newly introduced degree-corrected version [Karrer et al. Phys. Rev. E 83, 016107 (2011)], which imposes a degree sequence on the vertices, in addition to the block structure. The imposed degree sequence is implemented both as "soft" constraints, where only the expected degrees are imposed, and as "hard" constraints, where they are required to be the same on all samples of the ensemble. We also consider generalizations to multigraphs and directed graphs. We illustrate one of many applications of this measure by directly deriving a log-likelihood function from the entropy expression, and using it to infer latent block structure in observed data. Due to the general nature of the ensembles considered, the method works well for ensembles with intrinsic degree correlations (i.e. with entropic origin) as well as extrinsic degree correlations, which go beyond the block structure.
I. INTRODUCTION
The paper derives entropy expressions for stochastic blockmodel ensembles and uses them to obtain a log-likelihood for block inference. It extends the model to degree variation and multiple graph and degree-constraint variants.
I. INTRODUCTION
- The standard model gives same-block vertices equal expected degrees, while degree correction accommodates observed degree variation within blocks.The paper derives entropy expressions for traditional and degree-corrected ensembles, including simple and multigraph, directed and undirected, and soft and hard degree constraints.Source: §I · pp. 1–2
- Earlier blockmodel log-likelihoods may be non-closed or omit terms, preventing direct recovery of the ensemble entropy.With equal probability for each graph, the log-likelihood is the negative entropy and can be used to infer a likely block structure.Source: §I · p. 1
II. TRADITIONAL AND DEGREE-CORRECTED BLOCKMODELS
- The traditional microcanonical model fixes block sizes and edge counts; degree-corrected variants additionally impose vertex degrees as averages or as an exact sequence.When nonzero block-edge counts are sufficiently large, the microcanonical and canonical formulations become equivalent because fluctuations around the mean can be neglected.Source: §II · p. 2
A. Standard stochastic blockmodel
- For simple graphs, sparse-limit entropy approximations can omit important terms when block pairs are locally dense, even if the whole graph is sparse.The paper says to use the fuller expressions instead in such cases.Source: A. Standard stochastic blockmodel · p. 3
1. Directed graphs
Directed blockmodels allow asymmetric edge counts between blocks and distinguish edges leaving from those entering a block. Their entropy follows the undirected counting expression with a factor of one-half.
1. Directed graphs
- Directed blockmodels allow an asymmetric block-edge matrix and distinguish edges leaving a block from those entering it.The number of edge choices is counted as in the undirected case, without a separate diagonal treatment.Source: §1 · p. 3
- The directed entropy is identical to the corresponding undirected expression except for a factor of 1/2.The paper also gives the corresponding alternative expression and sparse approximation with the same factor.Source: §1 · p. 3
- For directed graphs, the diagonal edge count represents internal edges, rather than twice that number as in the undirected convention.This convention is specified for the directed entropy expression.Source: §1 · p. 3
B. Degree-corrected ensembles with “soft” constraints
Soft degree constraints assign each vertex an expected degree, allowing degree variability within blocks. The resulting entropy depends on the full degree-class edge structure and can retain arbitrary degree correlations.
B. Degree-corrected ensembles with “soft” constraints
- Expected degrees label vertices within each block, with each degree imposed as an ensemble average rather than a fixed value on every graph.This soft constraint is implemented by splitting each block into degree-labeled classes.Source: B. Degree-corrected ensembles with… · pp. 3–4
- The full matrix of edge counts between degree-labeled classes accommodates arbitrary expected-degree sequences and degree correlations, generalizing the ensemble in.A less-constrained version fixes only total edges between blocks, irrespective of degree labels.Source: B. Degree-corrected ensembles with… · p. 4
- Maximizing entropy subject to block-pair edge totals yields a nonlinear system with no general closed-form solution; its sparse-limit approximation neglects intrinsic correlations between block pairs.The sparse or classical limit assumes exp(λ_rs + μ_rκ + μ_sκ′) ≫ 1.Source: B. Degree-corrected ensembles with… · p. 4
- Higher moments of the expected-degree sequence can contribute substantially to entropy, regardless of the blockwise mean degree.Their importance depends on the degree sequence; broad distributions can make these contributions non-negligible even for sparse graphs.Source: B. Degree-corrected ensembles with… · p. 4
- In the stated limit, block partition and expected-degree sequence contribute independent entropy terms, and the block term corresponds to mutual information between memberships at edge ends.The expected edge-end distribution is expressed using the total half-edges leaving each block.Source: B. Degree-corrected ensembles with… · p. 4
- The sparse approximation requires valid degree-class edge counts and can fail with strong intrinsic inter-block disassortative degree correlations.Even when the approximation is valid, entropy still reflects the cost of connecting to a block with a broad degree distribution.Source: B. Degree-corrected ensembles with… · p. 5
1. Directed graphs
Directed degree-corrected ensembles group vertices by block and expected in- and out-degree. Their classical-limit entropy expressions are derived analogously to the undirected case, with the same validity caveats.
1. Directed graphs
- Directed degree-corrected ensembles separate vertices by block and expected in- and out-degree, using labels (r, κ−, κ+).The labels are incorporated into the degree-corrected entropy expression.Source: §1 · p. 5
- Taking the classical limit for directed degree-corrected ensembles yields an entropy expression by inserting the directed edge-count approximation into the degree-corrected entropy.The expression is omitted in the passage.Source: §1 · p. 5
- The directed expression inherits the undirected case’s caveats about the suitability of the approximation and the resulting entropy formula.Source: §1 · p. 5
IV. MULTIGRAPH ENSEMBLES
The multigraph treatment counts edge choices when parallel edges are allowed, then derives entropy expressions for undirected and directed cases. In the sparse limit, the undirected multigraph entropy matches that of simple graphs.
IV. MULTIGRAPH ENSEMBLES
- For undirected multigraphs, allowing multiple edges changes the count of edge choices between blocks and leads to a corresponding entropy expression.The edge-choice count uses combinations with repetition.Source: §IV · p. 5
- In the sparse limit, the multigraph entropy matches the simple-graph entropy because the probability of multiple edges vanishes.The expanded expression differs by an alternating sign in its last term before that term is dropped.Source: §IV · p. 5
2. Directed graphs
- For directed multigraphs, graph counts multiply over directed block pairs, since the edge-count matrix is generally asymmetric.The diagonal block term is treated without the undirected-case adjustment.Source: §2 · p. 6
- Directed multigraph entropy is the corresponding undirected expression with a factor of 1/2, including in the stated approximation.For directed graphs, e_rr counts internal edges rather than twice that number.Source: §2 · p. 6
A. Degree-corrected ensembles with “soft” constraints
Soft degree-corrected ensembles allow arbitrary degree correlations, while hard degree constraints require sparse-limit assumptions for entropy counting. The section extends these analyses to directed graphs and discusses the conditions under which the resulting expressions apply.
A. Degree-corrected ensembles with “soft” constraints
- Soft degree-corrected ensembles allow arbitrary degree correlations; imposing an uncorrelated constraint yields a nonlinear system without a closed-form solution.In the classical limit, the entropy expression is approximated by the sparse-limit expression, where the simple-graph and multigraph ensembles become equivalent.Source: A. Degree-corrected ensembles with… · p. 6
1. Directed graphs
- The directed multigraph analysis uses block labels that include in- and out-degree, and its classical-limit entropy expression has caveats like the simple-graph case.The expression is similar to the simple-graph ensemble except for an alternating sign in the last term.Source: §1 · p. 6
V. DEGREE-CORRECTED ENSEMBLES WITH “HARD” CONSTRAINTS
- The hard-constraint counting approach starts from configurations obeying a degree sequence and assumes sufficiently sparse limits, neglecting interactions between edge placements.The assumption is described as non-rigorous but compatible with more rigorous analyses.Source: §V · p. 7
- For a single block, randomly matching half-edges generates configurations; rejecting multigraphs can yield uniformly generated simple graphs, but multigraphs are not equiprobable.The hard-constraint case has no universally valid asymptotic entropy expression, even for a single block.Source: §V · pp. 6–7
A. Configurations
The section counts degree-sequence-compatible half-edge configurations and relates those counts to graph entropies. It also identifies conditions under which hard- and soft-constraint ensembles become equivalent.
A. Configurations
- The configuration count corrects for equivalent permutations of out-neighbours, but additional corrections are needed to count simple graphs or multigraphs accurately.For multigraphs, repeated edges reduce the number of equivalent pairings by a factor of m! for m parallel edges.Source: A. Configurations · p. 7
- In the sparse limit with sufficiently large degrees, simple- and multigraph ensembles with soft constraints and the hard-constraint configuration ensemble become equivalent.The comparison uses the condition that individual degrees equal their expected values.Source: A. Configurations · pp. 7–8
1. Directed configurations
- Directed configuration counting accounts for the asymmetric block-edge matrix and corrects for permutations of both in- and out-degrees.Source: §1 · p. 8
B. Simple graphs
The simple-graph entropy calculation estimates the correction for configurations that produce repeated edges or self-loops. Its hard-constraint expression agrees with a known single-block result, but equivalence with soft constraints requires additional conditions.
- The calculation estimates the probability that configurations yield simple graphs by treating repeated-edge and self-loop probabilities as independent.It applies Stirling’s approximation and assumes vertex degrees are much smaller than the number of half-edges leaving each block.Source: B. Simple graphs · p. 8
- For one block, the hard-constraint expression matches a known degree-sequence entropy result in the limit k_i ≪ N.The expression is also compatible with a later result, apart from a trivial constant.Source: B. Simple graphs · p. 8
- Hard- and soft-constraint ensembles are equivalent only in sufficiently sparse cases when higher-order terms can be neglected and degrees are large enough or distributions broad enough.The self-loop term must also be negligible.Source: B. Simple graphs · pp. 8–9
1. Directed graphs
The section derives hard-constraint entropy expressions for directed simple graphs and multigraphs. The expressions recover known single-block cases under sparse degree limits, with multigraph counting requiring a correction for repeated edges.
1. Directed graphs
- For one block, the directed simple-graph expression recovers a known result for imposed in- and out-degree sequences when k−/+ ≪ N.Source: §1 · p. 9
C. Multigraphs
- Multigraph configurations require a correction because graphs with m parallel edges have fewer generating configurations by a factor of m!.The resulting entropy resembles the simple-graph expression except for the signs of its last two terms.Source: C. Multigraphs · p. 9
- For one block, the multigraph expression matches a known result in the limit k_i ≪ N; another expression omits the final term because it counts parallel self-edges differently.Source: C. Multigraphs · p. 9
1. Directed multigraphs
Entropy-derived likelihoods extend blockmodel detection beyond the sparse, narrow-degree regime, while analogous calculations cover directed multigraphs. In a broad-degree example, omitting higher-order entropic effects yields spurious block partitions.
1. Directed multigraphs
- For directed multigraphs, the self-loop correction is treated like other parallel edges, so there is no separate self-loop term.With B = 1, the expression recovers the known count for multigraphs with imposed in/out-degree sequence when k−/+ ≪ N.Source: §1 · p. 9
VI. BLOCKMODEL DETECTION
- Blockmodel inference maximizes the likelihood over vertex partitions to identify latent equivalence classes, unlike community detection’s focus on dense within-group connections.The dense-within-group case corresponds to blockmodels whose diagonal elements of ers are largest.Source: §VI · pp. 9–10
- The objective’s L parameter controls how many terms are included: L = 0 recovers, while L > 0 includes higher-order corrections that may matter for large within-block degree moments.The second term can be neglected entirely in the L = 0 case.Source: §VI · p. 10
- A greedy inference procedure repeatedly moves vertices to maximize L, but can get trapped in local maxima and require restarts from other random partitions.Updating one vertex has complexity O(B(B(L + 1) + ⟨k⟩)), independent of system size.Source: §VI · p. 10
- The entropy-derived log-likelihood generalizes, whose expression is recovered when graphs are sparse and degree distributions are not broad.The same entropy methods extend to multigraph and directed ensembles.Source: §VI · p. 10
- In a deliberately simple broad-degree example, failing to account for entropic effects produces spurious block partitions relative to the higher-order objective.The example is chosen to make block identification easy, not to test the approach in difficult cases.Source: §VI · p. 10
A. Intrinsic degree correlations
On a broad-degree sample with intrinsic degree correlations, the likelihood identifies the correct B = 4 partition when higher-order terms are included, while smaller L values can favor degree-based artifacts. The analysis also shows that likelihood growth alone does not reliably select the number of blocks.
A. Intrinsic degree correlations
- A typical generated sample has N = 103 vertices, B = 4, strong block structure with w = 0.99, and a broad degree distribution with γ = 1.1 and [kmin, kmax] = [30, 200], causing intrinsic dissortative degree correlations.Source: A. Intrinsic degree correlations · p. 11
- Because log-likelihood increases monotonically with B, choosing the partition with the largest likelihood can fail to identify the appropriate block count.As B grows, the number of suitable models grows while the data remain fixed.Source: A. Intrinsic degree correlations · p. 11
- An estimate of the expected true log-likelihood helps distinguish gains from discovering block structure from gains caused by limited data.The passage denotes this expected value by L* as sample size tends to infinity.Source: A. Intrinsic degree correlations · p. 11
- B = 4 yields the correct partition for every tested L, while B = 8 with L = 0 separates vertices by degree, an entropic artifact.For L = 2, the inferred partition is correct; larger values give the same result with negligible differences.Source: A. Intrinsic degree correlations · p. 11, Figure 4 · p. 12
- For L ≥ 2, the inferred partition's correlation with degree essentially vanishes, whereas lower L values remain significantly correlated.The right panel compares average normalized mutual information across B and L values.Source: A. Intrinsic degree correlations · p. 12, Figure 3 · p. 12
- The log-likelihood approximation neglects some entropic effects and may diverge for sufficiently large degrees and L.The exact entropy expression can address extreme cases, but requires more parameters and may overfit.Source: A. Intrinsic degree correlations · p. 12
B. Extrinsic degree correlations
Without degree information, increasing the block count can make inferred sub-blocks track degree differences rather than the intended structure; an auxiliary degree-bin partition helps distinguish these effects. The conclusion presents analytical entropies for several blockmodel variants and applications to inference and network-system analysis.
B. Extrinsic degree correlations
- Ignoring degree correlations recovers the correct B = 4 partition, but larger B creates degree-based sub-blocks and rising log-likelihood; full resolution may require B ∼ 4K.For the example, K = 168, making the full parameter matrix larger than the number of network edges.Source: B. Extrinsic degree correlations · p. 13
- An auxiliary partition bins vertices by degree and is held fixed while inferring B blocks, allowing degree correlations to be inferred separately from degree-independent block structure.Choose few bins with low within-bin degree variance and avoid bins with very few vertices; then swap labels within each degree bin to make the inferred block labels meaningful.Source: B. Extrinsic degree correlations · p. 14
- With the auxiliary partition, log-likelihood growth points more convincingly to B = 4, and the degree-partition NMI differs significantly from inference without it.For D = 8, a residual degree association remains.Source: Figure 7 · p. 13, B. Extrinsic degree correlations · p. 14
VII. CONCLUSION
- The paper derives analytical entropies for traditional and degree-corrected ensembles, including directed and undirected graphs and soft and hard degree constraints.Setting the number of blocks to one recovers known entropies for random graphs with arbitrary degree sequences.Source: §VII · p. 14
- Entropy-derived likelihoods support block inference with intrinsic degree correlations and can be adapted to arbitrary extrinsic degree correlations, generalizing the approach in.The approach is described as suitable beyond the sparse-graph, non-broad-degree-sequence setting in which is expected to work well.Source: §VII · p. 14
- Blockmodel entropy can help assess the relevance of empirical-network topological features and the statistical significance of modular partitions.Source: §VII · p. 14
- For network systems with an energy function depending on block structure, blockmodel entropy can be used to obtain equilibrium properties.A cited gene-regulatory model uses robustness as an energy function and reports a core-periphery transition at higher selective pressure.Source: §VII · p. 14