# Measuring the modularity

In this example, we will use label_propagation to optimize the modularity of a network. The data are from the mcmullen dataset.

using EcologicalNetwork

# Get the data in an object
N = mcmullen();


The next step is to generate starting communities for every species. We will assign species to random initial modules:

L = rand(1:richness(N), richness(N));


We can now start a repeated number of attempts to find the best partition, here using brim. If julia is started in parallel, this will use all assigned CPUs.

M = modularity(N, L, brim, replicates=100);


Q: 0.43 |c|: 26