Consensus Community Detection

Runs a stochastic community detection algorithm multiple times and finds consensus communities via co-occurrence matrix thresholding. This approach produces more robust and stable community assignments than single runs.

Usage

community_consensus(
  x,
  method = c("louvain", "leiden", "infomap", "label_propagation", "spinglass"),
  n_runs = 100,
  threshold = 0.5,
  seed = NULL,
  ...
)

com_consensus(
  x,
  method = c("louvain", "leiden", "infomap", "label_propagation", "spinglass"),
  n_runs = 100,
  threshold = 0.5,
  seed = NULL,
  ...
)

Arguments

x: Network input: matrix, igraph, network, cograph_network, or tna object
method: Community detection algorithm to use. Default "louvain". Must be a stochastic method (louvain, leiden, infomap, label_propagation, spinglass).
n_runs: Number of times to run the algorithm. Default 100.
threshold: Co-occurrence threshold for consensus. Default 0.5. Nodes that appear together in >= threshold proportion of runs are placed in the same community.
seed: Optional seed for reproducibility. If provided, seeds for individual runs are derived from this seed.
...: Additional arguments passed to the community detection method.

Value

A cograph_communities object with consensus membership.

Details

The algorithm works as follows:

Run the specified algorithm n_runs times (without seeds to allow variation)
Build a co-occurrence matrix counting how often each pair of nodes appears in the same community
Normalize to proportions (0-1)
Threshold to create a consensus graph (edge if co-occurrence >= threshold)
Run walktrap on the consensus graph to get final communities

References

Lancichinetti, A., & Fortunato, S. (2012). Consensus clustering in complex networks. Scientific Reports, 2, 336.

Examples

if (requireNamespace("igraph", quietly = TRUE)) {
  g <- igraph::make_graph("Zachary")

  # Consensus from 50 Louvain runs
  cc <- community_consensus(g, method = "louvain", n_runs = 50)
  print(cc)

  # Stricter threshold for more robust communities
  cc2 <- community_consensus(g, threshold = 0.7, n_runs = 100)
}
#> Community structure (consensus_louvain)
#>   Nodes: 34  | Communities: 4  | Modularity: 0.6385 
#>   Sizes: 6, 5, 12, 11 
#> 
#>  node community
#>     1         3
#>     2         3
#>     3         3
#>     4         3
#>     5         2
#>     6         2
#>     7         2
#>     8         3
#>     9         4
#>    10         3
#>    11         2
#>    12         3
#>    13         3
#>    14         3
#>    15         4
#>    16         4
#>    17         2
#>    18         3
#>    19         4
#>    20         3
#>    21         4
#>    22         3
#>    23         4
#>    24         1
#>    25         1
#>    26         1
#>    27         4
#>    28         1
#>    29         1
#>    30         4
#>    31         4
#>    32         1
#>    33         4
#>    34         4