A Node Influence Based Label Propagation Algorithm for Community Detection in Networks

Abstract

Label propagation algorithm (LPA) is an extremely fast community detection method and is widely used in large scale networks. In spite of the advantages of LPA, the issue of its poor stability has not yet been well addressed. We propose a novel node influence based label propagation algorithm for community detection (NIBLPA), which improves the performance of LPA by improving the node orders of label updating and the mechanism of label choosing when more than one label is contained by the maximum number of nodes. NIBLPA can get more stable results than LPA since it avoids the complete randomness of LPA. The experimental results on both synthetic and real networks demonstrate that NIBLPA maintains the efficiency of the traditional LPA algorithm, and, at the same time, it has a superior performance to some representative methods.

1. Introduction

In recent years, complex networks have been widely used in many fields, such as social networks, World Wide Web networks, scientist cooperation networks, literature networks, protein interaction networks, and communication networks [1, 2]. Extensive studies have shown that complex networks have the property of communities (modules or clusters), within which the interconnections are close, but between which the associations are sparse. This property reflects the extremely common and important topology structure of complex networks and it is very important for understanding the structure and function of complex networks.

A great number of community detection algorithms have been proposed in recent decades, including modularity optimization algorithms [3–5], spectral clustering algorithms [6–8], hierarchical partition algorithms [9, 10], label propagation algorithms (LPA) [11, 12], and information theory based algorithms [13]. Among them, LPA is by far one of the fastest community detection algorithms. The complexity of LPA algorithm is nearly linear time, and the design of the algorithm is simple, all of which make LPA algorithm receive quite a lot of attention from numerous scholars [14–17]. However, it still has a number of shortcomings; for example, the community detection results are unstable.

In this paper, we propose a novel node influence based label propagation algorithm for community detection in networks (NIBLPA), improving the performance of the traditional LPA algorithm by fixing the node sequence of label updating and changing the label choosing mechanism when more than one label is contained by the maximum number of nodes. Firstly, NIBLPA calculates the node influence value of each node as the importance measure of nodes on the networks and fixes the nodes updating sequence in the descending order of node influence value; secondly, NIBLPA processes the label propagation repeatedly until the community structure of networks is detected. During each label updating process, when more than one label returned with the maximum number of nodes, instead of randomly selecting one label, we introduce the label influence into label computing formula to reselect the label from the set of labels with the same maximum number of nodes to improve the stability. Finally, NIBLPA divides all nodes with the same label into a community. Extensive experimental studies, by using various networks, demonstrate that our algorithm NIBLPA can get better community detection results compared with the state-of-the-art methods.

The rest of this paper is organized as follows. Section 2 introduces the related works including the traditional label propagation algorithm and the 𝑘-shell decomposition method. In Section 3, we introduce the main idea and the detailed process of our algorithm. The experimental results on various networks in Section 4 confirm the effectiveness of the algorithm. The conclusion is given in Section 5.

2. Related Work

A complex network can be modeled as a graph 𝐺=(𝑉,𝐸), where 𝑉={𝑣₁,𝑣₂,…,𝑣_𝑛} is the set of nodes, 𝐸={𝑒₁,𝑒₂,…,𝑒_𝑚} represents the edges between nodes, and 𝑛 and 𝑚 represent the number of nodes and edges in the network, respectively. Each edge in 𝐸 has a pair of nodes in 𝑉 corresponding. The label of 𝑣_𝑖 is denoted as 𝑐_𝑖. 𝑁(𝑖) represents the neighborhood set of 𝑣_𝑖 and 𝑑_𝑖 is the degree of node 𝑖.

2.1. Label Propagation Algorithm for Community Detection in Networks

In 2007, Raghavan et al. [11] applied the label propagation algorithm (LPA) to community detection, and the main idea of LPA is to use the network structure as the guide to detect community structures. LPA starts by giving each node a unique label, such as integers and letters, and in every iteration, each node changes its label to the one carried by the largest number of its neighbors. If more than one label is contained by the same maximum number of its neighbors, then randomly select one from them. In this repeated process, the dense groups of nodes change their different labels into the same label and nodes with the same label will be grouped into the same community.

The following equation is the formula of label updating: 𝑐_𝑖=argmax_{𝑙𝑗∈𝑁𝑙(𝑖)}1, (1) where 𝑁^𝑙(𝑖) represents the set of neighbors of 𝑣_𝑖 with label 𝑙.

For a weighted graph 𝐺, the weight of the edge between 𝑣_𝑖 and 𝑣_𝑗 is denoted as 𝑤_𝑖𝑗 and the label updating formula is changed as follows: 𝑐_𝑖=argmax_{𝑙𝑗∈𝑁𝑙(𝑖)}𝑤_𝑖𝑗. (2)

However, the algorithm cannot guarantee the convergence after several iterations. When the algorithm takes synchronize updating of the node labels (during the 𝑡th iteration, the node 𝑥 adopts its label only based on the labels of its neighbors at the (𝑡−1)th iteration), oscillations will occur in bipartite or nearly bipartite graph. As shown in Figure 1, the labels on the nodes oscillate between 𝑎 and 𝑏 in a bipartite graph. Therefore, Raghavan et al. [11] proposed asynchronous updating where node 𝑥 in the 𝑡th iteration updates its label based on a portion of labels at the 𝑡th iteration of its neighbors which have already been updated in the current iteration and another part of labels at the (𝑡−1)th iteration which are not yet updated in the current iteration to avoid the oscillation of labels.

Figure 1 

The oscillation of labels in a bipartite graph.

The design of label propagation algorithm is simple and easy to be understood. The process of the algorithm is presented in Algorithm 1.

Input: 𝐺=(𝑉,𝐸). Output: the result of community detection. (1) Initialization: assign a unique label to each node in the network, 𝑐𝑖(0)=𝑖. (2) Iteration of label propagation:  (a) Set 𝑡=1;  (b) Arrange the nodes of the network in random order, and store the results in the vector 𝑋.  (c) For each node 𝑣𝑖∈𝑋,let𝑐𝑖(𝑡)=𝑓(𝑐𝑖1(𝑡),…,𝑐𝑖𝑚(𝑡),𝑐𝑖(𝑚+1)(𝑡−1),…,𝑐𝑖𝑘(𝑡−1)), where 𝑣𝑖1,…,𝑣𝑖𝑚 are neighbors of 𝑣𝑖 those have already been updated in the current iteration and 𝑣𝑖(𝑚+1),…,𝑣𝑖𝑘 are neighbors those are not yet updated in the current iteration. The function 𝑓 here returns the label that the maximum number of its neighbors has. If multiple labels simultaneously have the maximum number, then randomly select one of them to assign to the node.  (d) If the label of every node does not change anymore, then stop the algorithm. Else, set 𝑡=𝑡+1 and go to Step (b). (3) Community division: divide all nodes share the same label into a community; the type of labels indicates the number of communities.

Algorithm 1 

Label propagation algorithm for community detection in networks (LPA).

In large networks with a huge number of nodes, each time the network may have different divisions because of the randomness of LPA algorithm. Among the solutions, it is difficult to determine which is the optimal. So the stability issue of LPA is necessary to be settled.

2.2. The 𝑘-Shell Decomposition Method

There are many measures we usually use to calculate the node importance, such as degree centrality [21], clustering coefficient centrality [22], and betweenness centrality [23]. Degree and clustering coefficient of nodes can only characterize the local information of networks. The complexity of computing betweenness is very high due to the need to calculate the shortest path. Kitsak et al. [24] pointed out that nodes with large 𝑘-shell value are very important for spreading dynamics on networks.

A 𝑘-shell is a maximal connected subgraph of 𝐺 in which every vertex’s degree is at least 𝑘. The 𝑘-shell value of node 𝑖, denoted by 𝐾𝑠(𝑖), indicates that node 𝑖 belongs to a 𝑘-shell but not to any (𝑘+1)-shell. The 𝑘-shell decomposition method is often used to identify the core and periphery of networks. It starts by removing all nodes with only one link, until no such nodes remain and assigns them to the 1-shell. In the same manner, it recursively removes all nodes with degree 2 (or less), creating the 2-shell. The process continues, increasing 𝑘 until all nodes in the network have been assigned to a shell. The shells with higher indices lie in the network core. The 𝑘-shell decomposition method can be efficiently implemented with the linear time complexity of 𝑂(𝑚), where 𝑚 is the number of edges in the network.

The 𝑘-shell decomposition method is shown in Figure 2. It is a simple network which can be divided into three different shells.

(a) The simple network with three shells

(b) 1-shell

(c) 2-shell

(d) 3-shell

Figure 2 

The illustration of 𝑘-shell decomposition method.

3. Our Method

Although asynchronous updating method can avoid oscillation of labels, there still are many limitations. As nodes are not updated simultaneously, the updating order of nodes plays a crucial impact on the stability and the quality of the results. The randomness of LPA in selecting one label when more than one label contained by the maximum number of nodes also makes the results unstable.

We analyze traditional LPA on a toy sample network in Figure 3 [25]. There are two communities in the network, {𝑣₁,𝑣₂,𝑣₃} and {𝑣₄,𝑣₅,𝑣₆}. The numbers inside the nodes represent their labels. Assuming that 𝑣₁, 𝑣₂, and 𝑣₃ have already shared the same label 2, while 𝑣₄, 𝑣₅, and 𝑣₆ still have unique labels. If we update 𝑣₄ first and randomly choose label 2 as its new label, then update 𝑣₆ before 𝑣₅. As a consequence, all nodes are classified into the same community. On the other hand, if node 𝑣₄ chooses label 6 and then updates node 𝑣₅ before 𝑣₆, the output will correspond with the right communities.

Figure 3 

The toy sample network.

Seen from the above analysis, LPA is very sensitive to the node updating order and the label choosing method. In this section we propose solutions to overcome the issues discussed above to improve the traditional LPA algorithm.

3.1. The Basic Idea

In the new algorithm, we choose the asynchronous updating method to avoid oscillation of labels in Figure 1. But the randomly determined label updating order of nodes affects the stability of the algorithm. We should order the nodes based on their importance for the network and the more important nodes should be updated earlier.

A node with a big 𝑘-shell value indicates that it is located in the core of the network. However, in a network, there are too many nodes with the same 𝑘-shell value and we cannot rank the node effectively. In general, in a network a node with more connections to the neighbors located in the core of the network is more important for the network. Inspired by these previous studies, we propose a novel centrality measure by considering both the 𝑘-shell value and degree of node itself and its neighbor’s 𝑘-shell values. The node influence of node 𝑖 is defined as follows: NI(𝑖)=𝐾𝑠(𝑖)+𝛼∗_{(𝑖)𝑗∈𝑁}𝐾𝑠𝑗𝑑𝑗, (3) where 𝛼 is a tunable parameter from 0 to 1, which is used to adjust the effect of its neighbors on the centrality of node 𝑖.

We choose node influence value as the measure of node importance, so we arrange nodes in the descending order of node influence value. The fixed node updating sequence makes the algorithm more stable.

The other random factor causing the instability of LPA is that when the number of labels with maximum nodes is more than one, the algorithm randomly selects one of the labels to assign to the node. Instead of randomly selecting one of the labels contained by the maximum nodes, we improve the label updating formula using the information of the label influence.

The label influence of label 𝑙 on node 𝑖 is computed as follows: LI(𝑙)=_{𝑗∈𝑁𝑙(𝑖)}NI𝑗𝑑𝑗. (4)

The new formula of label updating is changed as follows: 𝑐_𝑖=argmax_{𝑙∈𝑙max}LI(𝑙), (5) where 𝑙max denotes the set of labels that are simultaneously contained by the maximum nodes.

When multiple labels are simultaneously contained by the maximum nodes, we recalculate the value of the labels contained by the greatest number of nodes according to (5) and choose the label with the maximum value to assign to node 𝑖.

3.2. The Steps of NIBLPA Algorithm

The main steps of NIBLPA include initialization, iteration, and community division. Then NIBLPA can be described as Algorithm 2.

Input: 𝐺=(𝑉,𝐸). Output: the result of community detection. (1) Initialization: assign a unique label to each node in the network, 𝑐𝑖(0)=𝑖. (2) Calculate the node influence value for each node and arrange nodes in descending order of NI storing the results in the vector 𝑋. (3) Iteration of label propagation:  (a) Set 𝑡=1;  (b) For each node 𝑣𝑖∈𝑋,let𝑐𝑖(𝑡)=𝑓(𝑐𝑖1(𝑡),…,𝑐𝑖𝑚(𝑡),𝑐𝑖(𝑚+1)(𝑡−1),…,𝑐𝑖𝑘(𝑡−1)), where 𝑣𝑖1,…,𝑣𝑖𝑚 are neighbors of 𝑣𝑖 that have already been updated in the current iteration and 𝑣𝑖(𝑚+1),…,𝑣𝑖𝑘 are neighbors that are not yet updated in the current iteration. The function 𝑓 here returns the label that the maximum number of their neighbors has. If multiple labels simultaneously contained by the greatest number of nodes, then we recalculate each of the values of labels contained by greatest number nodes according to (5) and choose the label with maximum value to assign to the node 𝑣𝑖.  (c) If the label of every node does not change anymore, then stop the algorithm. Else, set 𝑡=𝑡+1 and go to Step (b). (4) Community division: assign all nodes share the same label into a community; the type of labels indicates the number of communities.

Algorithm 2 

Node influence based label propagation algorithm for community detection in networks (NIBLPA).

We implement NIBLPA on the toy sample network in Figure 3 with 𝛼=1. The decimals outside the nodes are the node influence value. Using our method on this network, the node updating sequence is fixed as 𝑣₁-𝑣₃-𝑣₄-𝑣₆-𝑣₂-𝑣₅ (when NI(𝑖)=NI(𝑗), rank 𝑖 and 𝑗 by their node IDs). The label propagation process is shown in Figure 4.

(a) Updating the label of 𝑣1

(b) Updating the label of 𝑣3

(c) Updating the label of 𝑣2

Figure 4 

Label propagation process of NIBLPA.

Firstly, we update the label of node 𝑣₁. We label 𝑣₁ with a set of tuples (𝑙,𝑛,LI(𝑙)), where 𝑐 is a label contained by its neighbor, and 𝑛 represents the number of its neighbors having the label 𝑙, and LI(𝑙) is an optional value recalculated by (5) when multiple labels are contained by the maximum neighbors. As shown in Figure 4(a), 𝑣₁ has three neighbors and they all have different labels with each other, and the set of tuples is {(2,1,1.833),(3,1,1.667),(4,1,1.667)}. So we choose label 2 as its new label.

Then, node 𝑣₃ is the next. After the label updating of 𝑣₁, there are two neighbors of 𝑣₃  that share label 2 and only one contains label 6, so we relabel 𝑣₃ with label 2 as shown in Figure 4(b). The next label propagations of 𝑣₄ and 𝑣₆ are consistent with 𝑣₁ and 𝑣₃.

Now only 𝑣₂ and 𝑣₅ are not updated and, as shown in Figure 4(c), all of their neighbors contain the same labels with themselves, respectively, so we do not need to relabel them. After only one iteration using this method, we get the final solution that contains two communities exactly the same with the ground truth. Since there is no randomness, the outcome is deterministic and perfect.

3.3. Time Complexity

The time complexity of the algorithm is estimated below. 𝑛 is the number of nodes, and 𝑚 is the number of edges.(1)The time complexity of initialization for all nodes: 𝑂(𝑛).(2)The time complexity of calculating the node influence value of all nodes: 𝑂(𝑚). The time complexity of ranking the nodes in descending order of NI: 𝑂 (𝑛log(𝑛)).(3)Each iteration of label propagation consists of two parts:(1)the time complexity of normal label updating: 𝑂(𝑚);(2)the time complexity of recalculating the labels based on (5) if necessary: 𝑂(𝑚).(4)The time complexity of assigning the nodes with the same label to a community: 𝑂(𝑛).

Phases (3) are repeated, so the time complexity of the whole algorithm is 2×𝑂(𝑛)+(2×𝑡+1)×𝑂(𝑚)+𝑂(𝑛log(𝑛)), where 𝑡 is the number of iterations and it is a small integer.

4. Experimental Studies

This section evaluates the effectiveness and the efficiency of our algorithm. We compare the performance of NIBLPA with LPA, KBLPA, and CNM. Where KBLPA is an improved LPA algorithm changing the node updating sequence by descending order of 𝑘-shell value. All the simulations are carried out in a desktop PC with Pentium Core2 Duo 2.8 GHz processor and 3.25 GB memory under Windows 7 OS. We implement our algorithm in Microsoft Visual Studio 2008 environment.

4.1. Datasets

In this section, we choose two types of synthetic and eight real networks to make experiments.

4.1.1. Clique-Ring Networks

Clique-Ring networks [26]: each clique is a complete subgraph 𝐾_𝑛 consisting of 𝑛(𝑛≥3) nodes and 𝑛(𝑛−1)/2 edges. A Clique-Ring network composed of 𝑚 such as subgraphs (𝑚≥2) has 𝑁=𝑚𝑛 nodes and 𝑀=𝑚𝑛(𝑛−1)/2+𝑚 edges. The structure is shown in Figure 5.

Figure 5 

The structure of Clique-Ring networks.

According to the generation rules of Clique-Ring networks, we construct four different size Clique-Ring networks. The parameters are shown in Table 1.

Table 1 

The parameters of four Clique-Ring networks.

4.1.2. LFR Benchmark Networks

LFR benchmark networks [27, 28] are currently the most commonly used synthetic networks in community detection. It can generate networks based on users’ need by changing the following parameters in Table 2.

Table 2 

The main parameters of LFR.

We generate six groups of LFR benchmark networks and all the networks share the common parameters of max𝑘 = 50. Each group contains nine networks with mu ranging from 0.1 to 0.9 and they also share parameters 𝑁, 𝑘, min𝑐, and max𝑐, respectively. The other parameters are set to the default values. The details are shown in Table 3.

Table 3 

The parameters of six groups of LFR networks.

4.1.3. Real Networks

We also make experiments on eight well known real networks, including Zachary’s karate club networks, Dolphins social networks, and American College Football networks. The detailed information of each network is shown in Table 4.

Table 4 

The information of real networks.

4.2. Evaluation Criteria

In this paper, we use modularity (𝑄) [2], 𝐹-measure [29], and normalized mutual information (NMI) [30] as the evaluation criteria which are currently widely used in measuring the performance of network clustering algorithms. Computing 𝐹-measure and NMI needs to know the true community structure of the network, while the modularity does not. For synthetic networks, since the ground truth of the community structure has been known, we use both 𝐹-measure and NMI on Clique-Ring networks and LFR benchmark networks to evaluate the results of community detection. While the underlying class labels of most real networks are unknown, we can only adopt the modularity as the evaluation criteria on partial real networks and use both NMI and modularity on others with known community structure.

4.2.1. Modularity

Consider the following: 1𝑄=2𝑚_{𝑖,𝑗∈𝑉}𝐴_𝑖𝑗−𝑑_𝑖𝑑_𝑗𝑐2𝑚×𝛿_𝑖,𝑐_𝑗, (6) where 𝑚 represents the number of edges in the network; 𝐴 is the adjacency matrix of the network, if node 𝑖 and node 𝑗 are directly connected, 𝐴_𝑖𝑗 = 1; otherwise, 𝐴_𝑖𝑗 = 0; 𝑐_𝑖 and 𝑐_𝑗, respectively, denote the label of node 𝑖 and node 𝑗, if 𝑐_𝑖=𝑐_𝑗, then 𝛿(𝑐_𝑖, 𝑐_𝑗)=1, else 𝛿(𝑐_𝑖, 𝑐_𝑗)=0.

4.2.2. 𝐹-Measure

Consider the following: 𝐹-Measure=2×Precision×Recall,Precision+Recall (7) where precision and recall are written as (8), respectively, ||Precision=𝑆∩𝑇|𝑆|,||Recall=𝑆∩𝑇|𝑇|. (8)

𝑇 is the set of node pairs (𝑖,𝑗), where nodes 𝑖 and 𝑗 belong to the same classes in the ground truth, and 𝑆 is the set of node pairs that belong to the same clusters generated by the evaluated algorithm. Then 𝑆∩𝑇 represents the intersection of node pairs of the ground truth and the clustering result.

4.2.3. Normalized Mutual Information (NMI)

Consider the following: ()=NMI𝑋∣𝑌−2×^|𝑋|_𝑖=1^|𝑌|_𝑗=1𝑋_𝑖∩𝑌_𝑗𝑋×log𝑛×_𝑖∩𝑌_𝑗𝑋_𝑖×𝑌_𝑗×^|𝑋|_𝑖=1𝑋_𝑖𝐶log^𝑖𝑛+^|𝑌|_𝑗=1𝑌_𝑖𝐶log^𝑗𝑛⁻¹, (9) where 𝑛 represents the number of nodes in the network, 𝑋 represents a community detection result generated by the evaluated algorithm, and 𝑌 represents the ground truth community structure.

4.3. Experimental Results and Analysis

In this section, the synthetic and real networks are used to test the effectiveness of NIBLPA comparing with traditional LPA, KBLPA, and CNM. Where LPA and KBLPA are processed 100 times and the average value is used as the results because of the randomness of these algorithms. We compare the stability of the algorithms by analyzing the fluctuation range of all the results.

4.3.1. The Experiments on Clique-Ring Networks

Table 5 shows the comparative results of the four algorithms on four different Clique-Ring networks, and for each instance, the best results are presented in boldface. The 𝐹-measure and NMI of LPA and KBLPA are in the form of average value ± the maximum difference between one result and the average value.

Table 5 

The comparison of results on Clique-Ring networks.

It can be seen from Table 5 that in the Clique-Ring networks which have special structure, NIBLPA can exactly detect the correct communities and CNM gets the right community structure on the first three networks. But on network C4, the result of CNM is much worse than others because modularity has the resolution limit problem. While the average 𝐹-measure of KBLPA algorithm is the lowest among LPA, KBLPA, and NIBLPA on the four networks and the average NMI of KBLPA is the lowest on most of the four networks except C4. These results illustrate that the fixed node sequence descending by the 𝑘-shell value at each step of label propagation cannot get good results. The instability of KBLPA is caused by the randomness of selecting label when multiple labels are simultaneously contained by the greatest number of nodes.

4.3.2. The Experiments on LFR Benchmark Networks

The twelve figures in Figure 6 are the NMI and 𝐹-measure of the four algorithms on six groups of LFR benchmark networks (N1~N6). The abscissa represents the parameter 𝑚𝑢 from 0.1 to 0.9. The ordinate in the left figures is the NMI of the results and the ordinate in the right figures is the 𝐹-measure.

(a) The experimental results on networks N1

(b) The experimental results on networks N2

(c) The experimental results on networks N3

(d) The experimental results on networks N4

(e) The experimental results on networks N5

(f) The experimental results on networks N6

Figure 6 

The results of the four algorithms on LFR benchmark networks.

The twelve figures in Figure 6 show that with the increase of 𝑚𝑢, the network structure is more and more complex and the four algorithms cannot be effective to detect the community structure. When mu is especially larger than 0.5, the NMI and 𝐹-measure decrease quickly. But generally, the performance of NIBLPA is better than the other three algorithms. Although NIBLPA does not guarantee to get the best performance, it can return stable, unique, and satisfied results. It can also be seen in Figure 6 that the fluctuation range of NMI and 𝐹-measure of LPA algorithm is large. KBLPA is also relatively stable, but its results are worse than LPA and NIBLPA. On these complex networks, CNM algorithm cannot detect the network structure effectively and it generally gets less number of communities than the truth.

4.3.3. The Experiments on Different Sizes of Networks

In order to compare the time efficiency of the algorithms, we generate 10 LFR benchmark networks, the size of which is from 1,000 to 10,000, and the other parameters are the same (𝑘 = 10, max𝑘 = 50, min𝑐 = 10, max𝑐 = 50, and 𝑚𝑢 = 0.1).

The time consumption of the four algorithms on the 10 LFR benchmark networks is shown in Figure 7(a). Figure 7(b) is the enlarge line chart of LPA, KBLPA, and NIBLPA.

(a) The comparison of time efficiency

(b) The enlarge line chart

Figure 7 

The efficiency comparison of the four algorithms on different sizes of networks.

From Figure 7, it is observed that the four algorithms use more and more time with the increase of the size of networks and CNM uses the longest time. When the number of nodes 𝑁 is larger than 5000, CNM cannot get the community structure because of the limit of computer memory. From Figure 7(b), one can note that when the number of nodes is greater than 7000, the time consumption of NIBLPA is less than LPA. To some extent, we can say NIBLPA is more suitable for community detection on large scale networks.

4.3.4. The Experiments on Real Networks

The eight real-world networks shown in Table 4 are commonly employed in the community detection literature and the first four networks have known ground truth community structures. So we compare the modularity 𝑄 and normalized mutual information NMI on the first four networks and only compare the modularity 𝑄 on the last four networks.

Table 6 shows the experimental results on the eight real networks, and for each instance, the best 𝑄 and NMI are presented in boldface.

Table 6 

The comparison of results on real networks.

It can be seen from Table 6 that in all the real networks besides R7(Blog) and R8(PGP), the modularity of NIBLPA is higher than the other three algorithms. Simultaneously, the NMI of NIBLPA on the first four networks is the best. The stability of KBLPA is better than LPA, but the modularity and NMI of KBLPA are worse than LPA on almost all of the networks. On the large size of PGP-network, CNM cannot detect the community structure. In general, NIBLPA can get better and stable results than the other three algorithms.

4.3.5. Instance Analysis

We compare the community structure detected by NIBLPA when NMI achieves the maximum with the true community structure of Dolphins.

Figure 8(a) shows the real community structures of Dolphins and Figure 8(b) is a community detection result of NIBLPA on Dolphins. Comparing these two figures, the division of DN63 and SN90 based on NIBLPA is inconsistent with the real structure. From the topology structure of Dolphins, we can see that DN63 has two adjacent nodes and they, respectively, belong to the two communities; DN63 has five neighbors, NIBLPA algorithm assigns it to the community which its most neighbors belong to. The modularity of Dolphins real community structure is lower than the result of NIBLPA, which draws a conclusion that the community division of NIBLPA is a reasonable result.

(a) The real community structures of Dolphins

(b) The result of NIBLPA on Dolphins

Figure 8 

The community structures of Dolphins.

4.3.6. Parameter Selection

There is only one parameter in NIBLPA algorithm, tunable parameter 𝛼. In order to analyze the impact of the parameter, we run NIBLPA with different values of 𝛼 on synthetic networks and compare NMI to analyze the effect of the parameter on the algorithm. In this way, we can investigate that under which 𝛼 the NIBLPA can achieve the best results.

We generate five LFR benchmark networks with 𝑘 ranging from 10 to 50 and all the networks share the common parameters of 𝑁=1000, max𝑘 = 50, min𝑐 = 10, max𝑐 = 50, and 𝑚𝑢 = 0.1. Figure 9 shows the results of NIBLPA on these networks.

Figure 9 

The value of NMI under different tunable parameter 𝛼.

As it can be seen in Figure 9, under different parameter 𝛼, the value of NMI changed a lot. However, for each network, there is an optimal 𝛼 under which the NIBLPA method can achieve the largest NMI. Moreover, on each network, the first extreme large value is generally the best result.

5. Conclusion

This paper presents a node influence based label propagation algorithm for community detection in networks. The algorithm firstly calculates the node influence value for each node and ranks the node in the descending order of node influence value. During each label updating process, when more than one label is contained by the maximum number of nodes, we introduce the label influence value into the formula of label updating to improve the stability. After the algorithm converges, nodes with the same label are divided into a community. This algorithm maintains the advantages of the original LPA algorithm. Moreover, it can get the stable community detection results by avoiding the randomness of label propagation. By experimental studies on synthetic and real networks, we demonstrate that the proposed algorithm has better performance than some of the current representative algorithms.

Conflict of Interests

The authors declare that they have no financial and personal relationships with other people or organizations that can inappropriately influence their work; there is no professional or other personal interest of any nature or kind in any product, service, or company that could be construed as influencing the position presented in, or the review of, the paper entitled.

Acknowledgments

This work was supported by the National High Technology Research and Development Program of China (Grant no. 2012AA0622022 and Grant no. 2012AA011004), the Doctoral Fund of the Ministry of Education of China (Grant no. 20100095110003 and Grant no. 20110095110010), and the Fundamental Research Funds for the Central Universities under Grant (Grant no. 2013XK10).