# Robust network community detection using balanced propagation

###### Abstract

Label propagation has proven to be an extremely fast method for detecting communities in large complex networks. Furthermore, due to its simplicity, it is also currently one of the most commonly adopted algorithms in the literature. Despite various subsequent advances, an important issue of the algorithm has not yet been properly addressed. Random (node) update orders within the algorithm severely hamper its robustness, and consequently also the stability of the identified community structure. We note that an update order can be seen as increasing propagation preferences from certain nodes, and propose a balanced propagation that counteracts for the introduced randomness by utilizing node balancers. We have evaluated the proposed approach on synthetic networks with planted partition, and on several real-world networks with community structure. The results confirm that balanced propagation is significantly more robust than label propagation, when the performance of community detection is even improved. Thus, balanced propagation retains high scalability and algorithmic simplicity of label propagation, but improves on its stability and performance.

## 1 Introduction

Complex real-world networks can comprise local structural modules (i.e., communities GN02 ()) that are groups of nodes densely connected within and only loosely connected with the rest of the network. Communities may play important roles in different real-world systems – they can be related to functional modules in biochemical networks PDFV05 () or individuals with common interests in social networks GN02 (). Moreover, community structure also has a strong impact on dynamic processes taking place on such networks ADP06 () and can thus provide an important insight into not only structural organization but also functional behavior of various real-world systems.

As a consequence, analysis of network community structure has been the focus of recent endeavor in different fields of science. There has also been a substantial number of community detection algorithms proposed in the literature over the last years CNM04 (); WH04 (); PDFV05 (); SJN06 (); RAK07 (); AK08 (); BGLL08 (); RB08 (); Liu10 (); RN10 (); SB11 () (for a comprehensive survey see For10 ()). Nevertheless, due to scalability issues, only a small minority of these algorithms can be applied to large real-world networks with several millions, billions of nodes, edges respectively.

A notable step towards this end was made by RAK07 ()Raghavanet al., who employed a simple label propagation to reveal significant communities in large real-world networks. Communities are identified by propagating (community) labels among nodes, thus, each node is assigned the label shared by most of its neighbors. Due to very fast structural inference of label propagation, densely connected sets of nodes form a consensus on some particular label after only a few iterations RAK07 (); SB11 (). The algorithm thus exhibits near linear complexity, which makes it applicable on networks with millions of nodes in a matter of minutes SB11 (). The basic algorithm was further analyzed and refined by various authors HCZLDF08 (); TK08 (); BC09 (); LHLC09 (); LM09b (); LM09c (); PCW09 (); PSSL09 (); Gre10 (); LM10 (); SB11 (); SB10a (); YWGW10 (), when, due to its simplicity, label propagation is also currently one of the most commonly adopted algorithms in the literature.

Despite the above efforts, an important issue of label propagation has not yet been properly addressed. To overcome convergence problems in some types of networks, RAK07 ()Raghavanet al. have proposed propagating labels among nodes (i.e., updating nodes’ labels) in a random order. Although this updating strategy solves the aforementioned problem, introduction of randomness severely hampers the robustness of the algorithm, and consequently also the stability of the identified community structure. It has been noted that the algorithm reveals a large number of distinct community structures even in smaller networks RAK07 (); TK08 (); LM09b (); SB11 (), when these structures are also relatively different among themselves TK08 (); SB11 (). Still, the robustness of the algorithm can also be related to the significance of community structure in a network SB11 ().

We argue that updating the nodes in some particular order can be seen as placing higher propagation preference LHLC09 () to the nodes that are updated at the beginning, and lower propagation preference to the nodes that are updated towards the end (and updating the nodes in a random order). The order of node updates thus governs the dynamics of the algorithm in a similar manner as (corresponding) node propagation preferences. This observation allows us to stabilize the label propagation algorithm by utilizing node preferences to counteract (i.e., balance) the randomness introduced by random node updates. The resulting algorithm is denoted balanced propagation and differs from label propagation merely in the introduction of node balancers.

We have evaluated the proposed algorithm on synthetic benchmark networks with planted partition, and on various real-world networks with community structure. The results confirm that balanced propagation is significantly more robust than simple label propagation, when the performance of community detection is even improved (in most cases). We also apply the algorithm to an entire European road network, which is not considered to reveal clear community structure. Nevertheless, the algorithm accurately identifies communities that correspond to different (geographical) regions of Europe, without any serious issues with stability.

The rest of the article is organized as follows. In Section 2 we formally present label propagation, and review issues and advances relevant for this research. Section 3 introduces balanced propagation and discusses the main rationale behind it. Empirical evaluation with discussion is given in Section 4 and conclusion in Section 5.

## 2 Label propagation

Let the network be represented by a simple undirected graph , where is the set of nodes and is the set of edges^{1}^{1}1In directed networks, each edge is treated as undirected, and in multi-networks, multiple edges among nodes are encoded into edge weights.. Furthermore, let be the weight of the edge incident to nodes . Moreover, let denote the community (label) of node and let denote the set of its neighbors.

Basic label propagation algorithm (LPA) RAK07 () reveals network communities by exploiting the following simple procedure. At first, each node is labeled with an unique label, . Then, at each iteration, each node adopts the label shared by most of its neighbors (considering also edge weights). Hence,

(1) |

where is the set of neighbors of that share label (ties are broken uniformly at random). Due to the existence of many intra-community edges, relative to the number of inter-community edges, densely connected sets of nodes form a consensus on some particular label after a few iterations. Thus, when the algorithm converges (i.e., equilibrium is reached), disconnected sets of nodes sharing the same label are classified into the same community. Due to extremely fast structural inference of label propagation, the algorithm exhibits near linear time complexity RAK07 (); SB11 () (in the number of edges of the network) and can easily scale to networks with millions, or even billions, of nodes and edges SB11 (); SB10a ().

LHLC09 ()Leunget al. have first noticed that label propagation can be substantially improved by increasing propagation preference (i.e., propagation strength) from certain nodes. The updating rule of the algorithm (i.e., equation (1)) is thus rewritten into

(2) |

where is the preference of node . Adequate node preferences can alter the dynamics of label propagation, in order to guide the algorithm towards a more significant community structure SB11 (). For the analysis and comparison of different node preference strategies, and corresponding algorithms, see LHLC09 (); SB11 (); SB10a ().

Next, we also discuss two main issues of label propagation and its advances. First, consider a bipartite network with two sets of nodes, denoted red and green nodes. Further assume that, at some point of the algorithm, all red nodes share label , and all green nodes share label . Due to bipartite structure, at the next iteration, all red nodes will adopt label , and all green nodes will adopt label . Moreover, at next iteration, all nodes will recover their initial labels, failing the algorithm to converge. It should be noted that such oscillations of labels are not limited to bipartite networks, but occur in various real-world networks that are commonly analyzed in the literature.

To ensure convergence, RAK07 ()Raghavanet al. have proposed asynchronous updating of nodes. Hence, nodes are no longer updated all together, but sequentially, in some random order. Thus, when node’s label is updated, possibly already updated labels of its neighbors are considered (in contrast to synchronous updating, where only labels from the previous iteration are considered). Although asynchronous updating eliminates aforementioned oscillations of labels, introduction of randomness severely disturbs the robustness of the algorithm, and consequently also the stability of the identified community structure. The stability of label propagation presents a severe issue for the algorithm, however, it has not yet been properly addressed in the past (to the best of our knowledge).

Second, consider a network with overlapping communities PDFV05 () and let be a node that has equally strong connections with two or more such communities. As ties are broken uniformly at random (see equation (1)), label would then, in general, constantly change. Furthermore, when many of such nodes exist, the algorithm would obviously never converge. Again, the issue is not limited to networks with overlapping communities.

Two possible solutions have been proposed in the literature. LHLC09 ()Leunget al. suggested including label into the maximal label consideration (besides merely neighbors’ labels), when RAK07 ()Raghavanet al. proposed a slightly modified approach. When there are multiple maximal labels (among neighbors’ labels), and one of them equals the concerned label , the node retains its label. In contrast to the former, the latter approach considers concerned label only when there indeed exist multiple maximal labels. Although both presented approaches work well for simple label propagation (i.e., equation (1)), this is not necessarily the case for different advances of the algorithm (e.g., equation (2)). Still, for the analysis in this article we adopt the approach proposed by RAK07 ()Raghavanet al..

In the proceeding section we revisit both issues discussed above, and propose solutions to overcome them.

## 3 Balanced propagation

Label propagation with asynchronous updating accesses the nodes in a random order. In particular, nodes are (re)shuffled before each iteration, in order to address convergence issues in some networks. However, as already discussed in Section 2, this incorporation of randomness severely hampers the robustness of the algorithm.

The issue can be addressed in an ad hoc fashion by simply accessing the nodes in some predefined (deterministic) order. This would clearly stabilize the algorithm, and possibly also perform well on real-world networks. We have conducted several experiments with different update orders, based on various node statistics (i.e., degree and eigenvector centrality Fre77 (); Fre79 (), clustering coefficient WS98 ()). Exact results are omitted, however, they indicate that, although none of these deterministic orders performs well in all networks, best order commonly corresponds to node preference strategy that also performs well. For instance, when ordering the nodes based on their degrees (decreasingly) gives good results, setting propagation preferences to the degrees of the nodes (and updating them in a random order) also performs well (and vice-versa).

Based on the above discussion we pose a hypothesis that the order of node updates within asynchronous label propagation governs algorithm’s dynamics in a similar manner as the corresponding node propagation preferences. Intuitively, nodes that are updated at the end of some iteration cannot efficiently propagate their final labels onward, as (most of) their neighbors have already been updated. On the other hand, a node that is considered first can possibly propagate its label to all of its neighbors, and thus form a community. Hence, nodes updated at the beginning exhibit higher propagation strength than those that are considered towards the end.

We further study the proposed hypothesis on a toy example network in Figure 1. The network consists of two communities, namely and , that are defined in a strong sense RCCLP04 () (i.e., each node has more intra-community than inter-community edges). Further assume that, at some point of the algorithm, nodes in , namely , and , are labeled with unique (community) labels, when all nodes in have already been classified to their right community (see Figure 1).

We first analyze how different orders of node updates affect the final outcome of the algorithm. When node is considered first, it will adopt the label of either or . Due to symmetry, we can assume that it adopts the label of node . No matter which of the nodes or is updated next, at the end of this iteration, all nodes in community will be labeled with the same label (that initially belongs to node ). The outcome thus corresponds to the natural community structure of the network.

On the other hand, when node is updated last, the results can differ. Again, we can assume that node is considered before node . If node adopts the label of either or , the algorithm proceeds similar as above. However, node can also adopt the label of the second community (with some probability). In that case, it is straightforward to see that nodes and will also adopt the same label, thus, at the end, all nodes in the network will be classified to the same community .

To summarize, if we first consider the core of community (i.e., node ), the label propagation will inevitably lead to the natural community structure of the network. However, if we access the border of community first (i.e., nodes and ), the algorithm could potentially classify all nodes into the same community (mainly due to the fact that community is already established). The example shows that even in such simple network, label propagation is extremely sensitive to the order of node updates.

Similar behavior as above can be observed, when we set higher propagation preference to either core or border of community (and update the nodes in a random order). When core node has the highest preference in the network, nodes and would obviously adopt the label of node . This would unavoidably lead to identification of the natural community structure, no matter the order of updates. However, when higher preference is given to border nodes and (i.e., lowest preference is given to node ), outcome of the algorithm can again correspond to the trivial community structure, where all nodes are classified into the same community (depends on the preference of other nodes and the order of updates). We thus conclude that, at least for this toy example, order of node updates can be seen as placing higher propagation preference to the nodes that are updated first, and lower propagation preference to the nodes that are updated last.

The latter enables us to stabilize the basic label propagation algorithm. As random node updates cannot be avoided (Section 2), node propagation preferences can be utilized to counteract the randomness introduced by random updates. Node preferences are thus employed to balance the algorithm (i.e., node balancers) and are set according to the reverse order in which the nodes are assessed. This retains the dynamics of the basic algorithm, but greatly improves its robustness and the stability of the identified community structure.

Let nodes be ordered in some random way, and let denote the normalized position of node in this order. Hence,

(3) |

where . Assuming linearity, we introduce node balancers as

(4) |

where is the preference of node (see equation (2)). Note that node balancers have to be recomputed at the beginning of each iteration (i.e., after each random shuffling of nodes). The resulting algorithm is else identical to the basic label propagation (with node preferences) and is denoted balanced propagation algorithm (BPA). Empirical evaluation in Section 4 shows that balanced propagation is not only more stable than basic label propagation, but also improves its community detection. Note also that the revealed community structure could be even further stabilized by, e.g., combining multiple network partitions SG02 ().

We also analyze a variant of the algorithm, where logistic function is used to model the relation between update orders and propagation preferences (the algorithm is denoted BPA). Hence, node balancers are set due to

(5) |

where and are parameters of the algorithm. We fix and based on some preliminary experiments. Empirical analysis reveals that BPA usually performs slightly better than BPA (Section 4).

Last, we also briefly consider the second main issue of label propagation. As already discussed in Section 2, nodes having equally strong connections with several (overlapping) communities might prevent the algorithm from converging. The problem is even enhanced in the case of balanced propagation, as random node preferences, introduced through random update orders, can extend the issue to cases, where node has only similarly strong connections with different communities. Consequently, solutions proposed in the literature RAK07 (); LHLC09 () do not necessarily overcome the problem in the case of balanced propagation.

Still, the true reason behind these convergence problems is the existence of overlapping communities in real-world networks. However, the purpose of this research is to address issues with random update orders, and not to extend balanced propagation to overlapping communities (see, e.g., Gre10 ()). Thus, for the sake of the empirical analysis, we adopt the following simple approach (and limit the analysis to non-overlapping communities).

As the discussed problems of balanced propagation (i.e., BPA and BPA algorithms) are actually an artifact of node balancers, we simply discard their use, when the algorithm does not converge after at most some maximal number of iterations. Note that this is in fact identical to applying the basic label propagation (i.e., LPA algorithm) afterwards, which obviously ensures the algorithm’s convergence. We fix the maximum number of iterations to , what should suffice for networks with almost a billion edges SB11 ().

## 4 Experiments and discussion

First, balanced propagation was analyzed, and compared against label propagation, on synthetic benchmark networks with planted partition and on several real-world networks with community structure (sections 4.1, 4.2 respectively). We address the stability of the algorithms and also the accuracy of community detection. Next, the proposed algorithm was further applied to a complete European road network, when the results are analyzed and discussed in Section 4.3.

Due to generality, results in the following sections are assessed in terms of different measures of community structure significance. Earlier work commonly reported the modularity NG04 () of the identified community structure. Modularity measures the significance of communities due to some null model (which is considered to be without community structure). Commonly, a random graph with the same degree sequence is selected for the null model. Hence,

(6) |

where is the adjacency matrix of the network, is degree of node and is the Kronecker delta. Higher values represent more significant community structure (), however, recent work shows that modularity has a number of severe deficiencies FB07 (); KSKK07 (); GMC10 () and should not be considered as a reliable indicator of community structure.

For a more adequate assessment of the significance of revealed communities we also adopt the conductance Bol98 (). Let be some community in the network thus . Conductance of a set of nodes is then defined as

(7) |

where is the complement of and is the cumulative degree of (i.e., ). Conductance thus measures the goodness of community , or equivalently, the quality of corresponding network cut . Lower values represent more significant communities (). Nevertheless, conductance cannot be easily extended to an entire community structure of a network. Thus, results are commonly assessed at different scales separately, in the form of network community profile (NCP) LLDM09 () plots. Still, due to simplicity, we also define as the average conductance over all communities in a network.

For networks with known community structure, identified communities are also compared against the true ones. We adopt two measures from the field of information theory Mac03 (). First, normalized mutual information (NMI) DDDA05 (), has become a de facto standard in the community detection literature. Let be a partition (i.e., communities) extracted by some algorithm, and let be the known partition for some network (corresponding random variables are and respectively). NMI of and is then

(8) |

where is the mutual information of the partitions (i.e., ), and , and are standard and conditional entropies. NMI of identical partitions equals , and is for independent partitions ().

Second, variation of information (VOI) Mei07 (), has several desirable properties with respect to NMI. In particular, it is symmetric local measure that also has the properties of a distance in the space of partitions. VOI of and is defined as

(9) |

thus, lower values represent better correlation between partitions. The maximum value of VOI depends on the size of the network (), therefore, for meaningful comparisons, we divide the obtained values with KLN08 ().

### 4.1 Synthetic networks with planted partition

We have first analyzed the balanced propagation on a class of synthetic benchmark networks with planted partition LFR08 (). The significance of community structure is controlled by a mixing parameter , where smaller values give clearer community structure. Networks exhibit power-law degree and community size distributions, as commonly observed in real-world networks BA99 (); New04 (). Power-law exponents are set to and respectively (i.e., ). Moreover, we fix the number of nodes to and vary the sizes of communities between and nodes. Results are assessed in terms of NMI and are shown in Figure 2.

Considering only the average performance (Figure 2, top), no clear difference between balanced propagation (i.e., BPA and BPA algorithms) and label propagation (i.e., LPA algorithm) is observed. However, scatter plots showing individual runs (Figure 2, bottom) reveal that there is actually a significant disparity between the approaches. When community structure is only roughly defined (i.e., for ), balanced propagation either relatively accurately identifies communities in the network (i.e., ) or classifies all nodes into a single community (i.e., ). On the other hand, label propagation also commonly reports community structures, whose correspondence to the actual communities is only marginal (i.e., , respectively). The latter is particularly apparent in the case of larger communities (note also the difference in error bars).

The results thus confirm that balanced propagation is much more robust than simple label propagation, when the community detection strength of the basic algorithm is largely retained in the refined versions (on average). Still, to obtain results comparable with current state-of-the-art community detection algorithms (see LF09a ()), different advances of the basic approach have to be employed SB11 (); SB10a ().

To further address the validity of balanced propagation, we have also applied the algorithms to a random graph à la Erdös-Rényi ER59 () that (presumably) has no community structure. The number of nodes is again fixed to , when we vary the average degree between and . Both balanced propagation algorithms reveal no community structure in these networks – all nodes are classified into a single community (or multiple communities in the case of disconnected networks) in all realizations of random networks. On the other hand, label propagation also partitions the networks into non-trivial communities, when the average degree is small enough (i.e., for ).

### 4.2 Real-world networks with community structure

Balanced propagation was further analyzed on eight real-world networks with community structure (Table 1). All these network are commonly employed in the community detection literature, and include different social, biological and technological networks. Due to simplicity, all networks were treated as unweighed and undirected.

Network | Description | Nodes | Edges |
---|---|---|---|

karate | Zachary’s karate club. Zac77 () | ||

dolphins | Lusseau’s dolphins. LSBHSD03 () | ||

books | Political books. Kre08 () | ||

football | American football. GN02 () | ||

jazz | Jazz musicians. GD03 () | ||

elegans | Nematode C. elegans. JTAOB00 () | ||

netsci | Network scientists. New06a () | ||

power | U.S. power grid. WS98 () |

We first directly compare the stability of the revealed community structures for balanced and label propagation (i.e., BPA and BPA, and LPA algorithms respectively). We apply the algorithms to each network times and count the number of distinct community structures obtained. We also measure the pairwise VOI of the partitions, to further evaluate the robustness of the algorithms. Due to space complexity, analysis is reduced to smaller networks (with at most hundreds of nodes). Results can be seen in Table 2.

Network | Distinct | Pairwise | ||||
---|---|---|---|---|---|---|

LPA | BPA | BPA | LPA | BPA | BPA | |

karate | ||||||

dolphins | ||||||

books | ||||||

football | ||||||

jazz | ||||||

elegans |

The analysis confirms earlier observations that basic label propagation is relatively unstable, even on smaller networks RAK07 (); TK08 (); LM09b (); SB11 (). However, the latter does not hold for balanced propagation that reveals only a small number of distinct community structures in each network. In most cases, this number is for a scale smaller than in the case of label propagation. Moreover, the pairwise similarity between the structures is also significantly improved, when the same trend is observed if we measure similarity only among distinct structures (e.g., for elegans network, average pairwise VOI equals , and for LPA, BPA and BPA algorithms respectively).

We conclude that balanced propagation is significantly more robust than label propagation, and can be, despite its randomized nature, considered as fairly stable. Note also that balanced propagation with logistic model (i.e., BPA algorithm) performs slightly better than the basic algorithm with a linear model (i.e., BPA algorithm).

Three of the networks in Table 1, namely karate, dolphins and football, have known natural partitions into communities (that result from earlier studies). To analyze also the community detection strength of balanced propagation, we measure the VOI between the natural partitions and those identified by different algorithms. The results appear in Table 3, when we also report the results for a classical modularity optimization algorithm (MO) proposed by CNM04 ()Clausetet al. (for reference).

Network | Number | ||||
---|---|---|---|---|---|

LPA | BPA | BPA | MO | ||

karate | |||||

dolphins | |||||

football |

Note that, in the case of karate and dolphins networks, balanced propagation performs significantly better than label propagation (and modularity optimization), when in the case of football network, the obtained VOI is roughly the same. Thus, despite relatively similar performance on synthetic benchmark networks (Section 4.1), balanced propagation more accurately identifies the true communities within these real-world networks than label propagation (and also modularity optimization).

For a better comprehension, the fraction of correctly classified GN02 () nodes for BPA algorithm equals , and for karate, dolphins and football networks respectively (on average).

In Table 4 we also report average conductance and modularity of the revealed community structures for all networks in Table 1 (mainly to enable comparison with earlier work). Balanced propagation also performs better in terms of conductance. Still, results should be taken with caution as BPA and BPA algorithms commonly return larger communities than LPA algorithm, which implies lower average conductance (see below). On the other hand, according to modularity, performance depends on the size of the network. We argue that this is an artifact of an intrinsic scale incorporated into the measure of modularity (i.e., resolution limit FB07 (); GMC10 ()), thus, lower values of modularity obtained by balanced propagation on smaller networks should not be attributed to weaker community structure (see Table 3).

Again, a general pattern can be observed between both balanced propagation algorithms.

Net. | ||||||
---|---|---|---|---|---|---|

LPA | BPA | BPA | LPA | BPA | BPA | |

kara. | ||||||

dolph. | ||||||

books | ||||||

foot. | ||||||

jazz | ||||||

eleg. | ||||||

netsci | ||||||

power |

Next, we further analyze the larger two networks in Table 1, namely, netsci and power. We apply each algorithm times and analyze the conductance of obtained communities at different scales. The results are reported in the form of network community profile (NCP) LLDM09 () plots, and are shown in Figure 3. NCP plots measure the quality of the best community (due to conductance) as a function of its size (Figure 3, below). Social and information, and also technological, networks commonly reveal rather characteristic structure of NCP plots, with initial decreasing and subsequent increasing trend (for more see LLDM09 ()).

Observe that balanced propagation identifies communities on a much wider scale, including also larger communities. The structure of NCP plots thus better coincides with the analysis of LLDM09 ()Leskovecet al., where a natural (i.e., best) community size was estimated to a round nodes. In other words, basic label propagation finds best communities at much smaller scale than balanced propagation (i.e., at a round nodes), when the conductance is also significantly higher on average (Table 4). Note also that label propagation reveals a number of communities with very high conductance (i.e., (black) circles in the uppermost part of Figure 3, top), which can be directly related to the issues of the algorithm discussed in Section 2.

We conclude that, at least for the networks analyzed, balanced propagation is indeed more stable than basic label propagation, when the quality of the identified community structure is also improved in most cases.

Last, we also briefly analyze the scalability of the proposed balanced propagation. In Table 5 we report the average number of iterations^{2}^{2}2Each iteration has linear time complexity . made by the algorithms over runs. As discussed in Section 3, we do not directly address the issues with overlapping communities. Therefore, nodes, having strong connections with different communities, can prevent basic balanced propagation from converging. The results in Table 5 thus include only the runs where the algorithms converged in a fixed (maximal) number of iterations (this includes at least of runs in each case). For the same reason, netsci and power networks were not included in the analysis.

Network | Iterations | ||
---|---|---|---|

LPA | BPA | BPA | |

karate | |||

dolphins | |||

books | |||

football | |||

jazz | |||

elegans |

The complexity of label propagation is quite lower compared to balanced propagation. Still, all algorithms reveal communities in a relatively small number of iterations and can be easily scaled to larger networks (exhibit near linear time complexity ). It should also be noted that extremely fast convergence of label propagation can be somewhat related to random node updates (Section 2). Random update order can be seen as increasing propagation strength from certain nodes (Section 3), which limits the dynamics of the algorithm, and instantly leads it towards some stable, probably suboptimal (i.e., random), partition. The convergence of the algorithm is thus indeed fast, still, the identified community structure is extremely unstable and often suboptimal (as also observed by previous work RAK07 (); TK08 (); LM09b (); SB11 ()).

### 4.3 European road network

Road networks are not considered to convey a clear community structure, consisting of densely connected modules (due to sparsity of such networks). However, the network can still contain groups of nodes that are well isolated from others (i.e., connected through only few edges) and community detection algorithms can be employed to reveal such partition of the network. Communities should in this case largely relate to the properties of the road transport within the region, and also coincide with the geographical characteristics of the area.

We have constructed a network of all roads included in the International E-road Network (Figure 4). Nodes thus correspond to European cities and edges represent direct (class A, B) road connections among them. We limit the analysis to the main component of the network that consists of nodes and edges (a complete network has nodes and edges). Note that the network is neither scale-free BA99 () (i.e., maximum degree equals , when the degree distribution is, e.g., log-normal) nor small-world WS98 () (i.e., average distance among nodes is and the clustering coefficient WS98 () equals ).

Due to long average distances among different parts of the network, road networks are particularly hard to partition with standard community detection algorithms. Furthermore, as the network has almost tree-like structure, it is often hard to decide where to split long paths of nodes. Indeed, if we apply the basic label propagation (i.e., LPA algorithm) we obtain communities with and (on average over runs). Hence, communities consist of only nodes on average, thus, they can only hardly be considered as meaningful.

On the other hand, balanced propagation (i.e., BPA algorithm) partitions the network into communities with and (on average over runs). In Figure 4 we show the community structure that obtained minimum average conductance . Note how the largest communities quite accurately coincide with different (geographical) regions of Europe. In particular, from left to right (top to bottom), communities represent cities of Iberian Peninsula (e.g., Madrid), eastern Central Europe (e.g., Berlin), western Central Europe (e.g., Paris), Apennine Peninsula (e.g., Rome), eastern Russia, western Russia and Finland (e.g., Moscow), northern East Europe (e.g., Bratislava), southern East Europe (e.g., Bucharest), Balkan Peninsula (e.g., Skopje), Scandinavian Peninsula (e.g., Stockholm), etc. It is ought to be mentioned that, although community structures revealed by the algorithm through different runs indeed differ, in most cases, largest communities correspond to the same regions as discussed above. The latter thus further confirms the robustness of the balanced propagation.

## 5 Conclusions

The article addresses one of the main issues of label propagation algorithm for community detection – the stability of the identified community structure. We introduce balanced propagation that controls (i.e., stabilizes) the dynamics of basic label propagation through utilization of node balancers. The resulting approach is significantly more robust than its label propagation counterpart, when its community detection strength is even improved. Thus, balanced propagation retains high scalability and algorithmic simplicity of label propagation, but improves on its stability and performance. The proposition has been validated on synthetic networks with planted partition, and on several real-world networks with community structure. Moreover, the proposed algorithm was further applied to an entire European road network, where it accurately partitions the network with respect to (geographical) regions.

Due to its simplicity, balanced propagation can be easily incorporated into arbitrary (label) propagation algorithm, not limited to the field of community detection. Moreover, the work provides further comprehension of the propagation on networks, with different applications.

###### Acknowledgements.

The work has been supported by the Slovene Research Agency ARRS within the research program P2-0359.## References

- (1) M. Girvan, M.E.J. Newman, P. Natl. Acad. Sci. USA pp. 7821–7826 (2002)
- (2) G. Palla, I. Derényi, I. Farkas, T. Vicsek, Nature 435, 814 (2005)
- (3) A. Arenas, A. Díaz-Guilera, C.J. Pérez-Vicente, Phys. Rev. Lett. 96(11), 114102 (2006)
- (4) A. Clauset, M.E.J. Newman, C. Moore, Phys. Rev. E 70(6), 066111 (2004)
- (5) F. Wu, B.A. Huberman, Eur. Phys. J. B 38(2), 331 (2004)
- (6) S. Son, H. Jeong, J.D. Noh, Eur. Phys. J. B 50(3), 431 (2006)
- (7) U.N. Raghavan, R. Albert, S. Kumara, Phys. Rev. E 76(3), 036106 (2007)
- (8) G. Agarwal, D. Kempe, Eur. Phys. J. B 66(3), 409 (2008)
- (9) V.D. Blondel, J. Guillaume, R. Lambiotte, E. Lefebvre, J. Stat. Mech. P10008 (2008)
- (10) M. Rosvall, C.T. Bergstrom, P. Natl. Acad. Sci. USA 105, 1118 (2008)
- (11) J. Liu, Eur. Phys. J. B 77(4), 547 (2010)
- (12) P. Ronhovde, Z. Nussinov, Phys. Rev. E 81(4), 046114 (2010)
- (13) L. Subelj, M. Bajec, Phys. Rev. E (2011), accept.
- (14) S. Fortunato, Phys. Rep. 486(3-5), 75 (2010)
- (15) Y. Hu, H. Chen, P. Zhang, M. Li, Z. Di, Y. Fan, Phys. Rev. E 78(2), 026121 (2008)
- (16) G. Tibély, J. Kertész, Physica A 387(19-20), 4982 (2008)
- (17) M.J. Barber, J.W. Clark, Phys. Rev. E 80(2), 026129 (2009)
- (18) I.X.Y. Leung, P. Hui, P. Liò, J. Crowcroft, Phys. Rev. E 79(6), 066107 (2009)
- (19) X. Liu, T. Murata, Physica A 389(7), 1493 (2009)
- (20) X. Liu, T. Murata, Community detection in large-scale bipartite networks, in Proceedings of the International Conference on Web Intelligence and Intelligent Agent Technology (2009), Vol. 1, pp. 50–57
- (21) S. Pang, C. Chen, T. Wei, A realtime clique detection algorithm: Time-based incremental label propagation, in Proceedings of the International Conference on Intelligent Information Technology Application (2009), Vol. 3, pp. 459–462
- (22) C. Pang, F. Shao, R. Sun, S. Li, Detecting community structure in networks by propagating labels of nodes, in Proceedings of the International Symposium on Neural Networks (2009), pp. 839–846
- (23) S. Gregory, New J. Phys. 12, 103018 (2010)
- (24) X. Liu, T. Murata, Evaluating community structure in bipartite networks, in Proceedings of the IEEE International Conference on Social Computing (2010), pp. 576–581
- (25) L. Subelj, M. Bajec, Unfolding network communities by combining defensive and offensive label propagation, in Proceedings of the ECML PKDD Workshop on the Analysis of Complex Networks (2010), pp. 87–104
- (26) Q. Ye, B. Wu, Y. Gao, B. Wang, Detecting communities in massive networks based on local community attractive force optimization, in Proceedings of the International Conference on Advances in Social Networks Analysis and Mining (2010), pp. 291–295
- (27) L. Freeman, Sociometry 40(1), 35 (1977)
- (28) L.C. Freeman, Soc. Networks 1(3), 215 (1979)
- (29) D.J. Watts, S.H. Strogatz, Nature 393(6684), 440 (1998)
- (30) F. Radicchi, C. Castellano, F. Cecconi, V. Loreto, D. Parisi, P. Natl. Acad. Sci. USA 101, 2658 (2004)
- (31) A. Strehl, J. Ghosh, J. Mach. Learn. Res. 3, 583–617 (2002)
- (32) M.E.J. Newman, M. Girvan, Phys. Rev. E 69(2), 026113 (2004)
- (33) S. Fortunato, M. Barthelemy, P. Natl. Acad. Sci. USA pp. 36–41 (2007)
- (34) J. Kumpula, J. Saramäki, K. Kaski, J. Kertész, Eur. Phys. J. B 56(1), 5 (2007)
- (35) B.H. Good, Y.A. de Montjoye, A. Clauset, Phys. Rev. E 81(4), 046106 (2010)
- (36) B. Bollobás, Modern graph theory (Springer, 1998)
- (37) J. Leskovec, K.J. Lang, A. Dasgupta, M.W. Mahoney, Internet Math. 6(1), 29 (2009)
- (38) D.J.C. MacKay, Information theory, inference, and learning algorithms (Cambridge University Press, 2003)
- (39) L. Danon, A. Díaz-Guilera, J. Duch, A. Arenas, J. Stat. Mech. P09008 (2005)
- (40) M. Meila, J. Multivariate Anal. 98(5), 873 (2007)
- (41) B. Karrer, E. Levina, M.E.J. Newman, Phys. Rev. E 77(4), 046119 (2008)
- (42) A. Lancichinetti, S. Fortunato, F. Radicchi, Phys. Rev. E 78(4), 046110 (2008)
- (43) A.L. Barabási, R. Albert, Science 286(5439), 509 (1999)
- (44) M.E.J. Newman, Eur. Phys. J. B 38(2), 321 (2004)
- (45) A. Lancichinetti, S. Fortunato, Phys. Rev. E 80(5), 056117 (2009)
- (46) P. Erdős, A. Rényi, Publ. Math. Debrecen 6, 290 (1959)
- (47) W.W. Zachary, J. Anthropol. Res. 33(4), 452 (1977)
- (48) D. Lusseau, K. Schneider, O.J. Boisseau, P. Haase, E. Slooten, S.M. Dawson, Behav. Ecol. Sociobiol. 54(4), 396 (2003)
- (49) V. Krebs, A network of co-purcheased books about U.S. politics, http://www.orgnet.com/ (2008), http://www.orgnet.com/
- (50) P. Gleiser, L. Danon, Adv. Complex Syst. 6(4), 565 (2003)
- (51) H. Jeong, B. Tombor, R. Albert, Z.N. Oltvai, A. Barabási, Nature 407, 651 (2000)
- (52) M.E.J. Newman, Phys. Rev. E 74(3), 036104 (2006)