Compute Modularity

Computes the modularity Q of a partition with resolution parameter. Q = (1/2m) * sum_ij (A_ij - gamma * k_i * k_j / (2m)) * delta(c_i, c_j)

Usage

mskModularity(A, membership, gamma = 1)

Arguments

A

Adjacency matrix.

membership

Integer vector of community assignments.

gamma

Resolution parameter (default: 1.0).

Value

Numeric modularity value.

Examples

C <- matrix(c(1, 1, 0, 0,  0, 1, 1, 0,  0, 0, 1, 1,  1, 0, 0, 1),
            nrow = 4, dimnames = list(paste0("b", 1:4), paste0("m", 1:4)))
hg <- MSKHypergraph(C)
A <- projectMuscleGraph(hg)
cm <- mskCommunityDetect(hg, gamma = 1)
mskModularity(A, cm$membership, gamma = 1)
#> m1 
#>  0