Lancichinetti–Fortunato–Radicchi benchmark

Short description: Algorithm

Lancichinetti–Fortunato–Radicchi benchmark is an algorithm that generates benchmark networks (artificial networks that resemble real-world networks). They have a priori known communities and are used to compare different community detection methods. The advantage of the benchmark over other methods is that it accounts for the heterogeneity in the distributions of node degrees and of community sizes.

The algorithm

The node degrees and the community sizes are distributed according to a power law, with different exponents. The benchmark assumes that both the degree and the community size have power law distributions with different exponents, [math]\displaystyle{ \gamma }[/math] and [math]\displaystyle{ \beta }[/math], respectively. [math]\displaystyle{ N }[/math] is the number of nodes and the average degree is [math]\displaystyle{ \langle k \rangle }[/math]. There is a mixing parameter [math]\displaystyle{ \mu }[/math], which is the average fraction of neighboring nodes of a node that do not belong to any community that the benchmark node belongs to. This parameter controls the fraction of edges that are between communities. Thus, it reflects the amount of noise in the network. At the extremes, when [math]\displaystyle{ \mu = 0 }[/math] all links are within community links, if [math]\displaystyle{ \mu = 1 }[/math] all links are between nodes belonging to different communities.

One can generate the benchmark network using the following steps.

Step 1: Generate a network with nodes following a power law distribution with exponent [math]\displaystyle{ \gamma }[/math] and choose extremes of the distribution [math]\displaystyle{ k_{\min} }[/math] and [math]\displaystyle{ k_{\max} }[/math] to get desired average degree is [math]\displaystyle{ \langle k\rangle }[/math].

Step 2: [math]\displaystyle{ (1 - \mu) }[/math] fraction of links of every node is with nodes of the same community, while fraction [math]\displaystyle{ \mu }[/math] is with the other nodes.

Step 3: Generate community sizes from a power law distribution with exponent [math]\displaystyle{ \beta }[/math]. The sum of all sizes must be equal to [math]\displaystyle{ N }[/math]. The minimal and maximal community sizes [math]\displaystyle{ s_{\min} }[/math] and [math]\displaystyle{ s_{\max} }[/math] must satisfy the definition of community so that every non-isolated node is in at least in one community:

[math]\displaystyle{ s_{\min} \gt k_{\min} }[/math]

[math]\displaystyle{ s_{\max} \gt k_{\max} }[/math]

Step 4: Initially, no nodes are assigned to communities. Then, each node is randomly assigned to a community. As long as the number of neighboring nodes within the community does not exceed the community size a new node is added to the community, otherwise stays out. In the following iterations the “homeless” node is randomly assigned to some community. If that community is complete, i.e. the size is exhausted, a randomly selected node of that community must be unlinked. Stop the iteration when all the communities are complete and all the nodes belong to at least one community.

Step 5: Implement rewiring of nodes keeping the same node degrees but only affecting the fraction of internal and external links such that the number of links outside the community for each node is approximately equal to the mixing parameter [math]\displaystyle{ \mu }[/math].

Testing

Consider a partition into communities that do not overlap. The communities of randomly chosen nodes in each iteration follow a [math]\displaystyle{ p(C) }[/math] distribution that represents the probability that a randomly picked node is from the community [math]\displaystyle{ C }[/math]. Consider a partition of the same network that was predicted by some community finding algorithm and has [math]\displaystyle{ p(C_2) }[/math] distribution. The benchmark partition has [math]\displaystyle{ p(C_1) }[/math] distribution. The joint distribution is [math]\displaystyle{ p(C_1, C_2) }[/math]. The similarity of these two partitions is captured by the normalized mutual information.

[math]\displaystyle{ I_n = \frac{\sum_{C_1,C_2} p(C_1,C_2) \log_2 \frac{p(C_1,C_2)}{p(C_1)p(C_2)} }{\frac 1 2 H(\{p(C_1)\}) + \frac 1 2 H(\{p(C_2)\})} }[/math]

If [math]\displaystyle{ I_n=1 }[/math] the benchmark and the detected partitions are identical, and if [math]\displaystyle{ I_n=0 }[/math] then they are independent of each other.