key: cord-0753583-1yodfhgm authors: Huang, Xinyu; Chen, Dongming; Wang, Dongqi; Ren, Tao title: Identifying Influencers in Social Networks date: 2020-04-15 journal: Entropy (Basel) DOI: 10.3390/e22040450 sha: dcf9d7f00c45be4c107878ec26845ad14fc20bb5 doc_id: 753583 cord_uid: 1yodfhgm Social network analysis is a multidisciplinary research covering informatics, mathematics, sociology, management, psychology, etc. In the last decade, the development of online social media has provided individuals with a fascinating platform of sharing knowledge and interests. The emergence of various social networks has greatly enriched our daily life, and simultaneously, it brings a challenging task to identify influencers among multiple social networks. The key problem lies in the various interactions among individuals and huge data scale. Aiming at solving the problem, this paper employs a general multilayer network model to represent the multiple social networks, and then proposes the node influence indicator merely based on the local neighboring information. Extensive experiments on 21 real-world datasets are conducted to verify the performance of the proposed method, which shows superiority to the competitors. It is of remarkable significance in revealing the evolutions in social networks and we hope this work will shed light for more and more forthcoming researchers to further explore the uncharted part of this promising field. The research of network science is experiencing a blossom in the last decade, which provides profound implications in very different fields, from finance to social and biological networks [1] . Considering the enormous data scale, most studies merely focus on a small group of influential nodes rather than the whole network. Take social networks for instance, influential nodes are those that have the most spreading ability, or playing a predominant role in the network evolution. Notably, a popular star in online social media may remarkably accelerate the spreading of rumors, and a few super spreaders [2] could largely expand the epidemic prevalence of a disease (e.g., [3] . The research of influencer identification is beneficial to understanding and controlling the spreading dynamics in social networks with diverse applications such as epidemiology, collective dynamics and viral marketing [4, 5] . Nowadays, individuals interact with each other in more complicated patterns than ever. It is a challenging task to identify influencers in social networks for the various kinds of interactions. As we have known, the graph model is widely utilized to represent social networks, however, it is incapable of dealing with the multiple social links. For example, people use Facebook or WeChat to keep communication with family members or friends, use Twitter to post news, use LinkedIn to search for jobs, and use TikTok to create and share short videos [6] . It is easy to represent each social scenario via a graph model separately, in spite of they are belonging to the same group of individuals. The neglect of the multiple relationships between social actors may lead to an incorrect result of the most versatile users [7] . With the proposal of multilayer networks [8, 9] , we are able to encode the various interactions, which is of great importance and necessity of identifying influencers in multiple social networks. In this paper, we design a novel node centrality measure for monolayer network, and then apply it to multilayer networks to identify influencers in multiple social networks. This method is solely based on the local knowledge of a network's topology in order to be fast and scalable due to the huge size of networks, and thus suitable for both real-time applications and offline mining. The rest of this paper is organized as follows. Section 2 introduces the related works on influencers identification in monolayer network and multilayer networks. Section 3 presents the mathematical model and the method for detecting influencers. Section 4 exhibits the experiments and analysis, including comparison experiments on twenty-one real-world datasets, which verifies the feasibility and veracity of the proposed method. Section 5 summarizes the whole paper and provides concluding remarks. The initial research on influencers identification may date back to the study of node centrality, which means to measure how "central" a focal node is [10] . A plethora of methods for influencers identification are proposed in the past 40 years, which can be mainly classified into centrality measures, link topological ranking measures, entropy measures, and node embedding measures [11, 12] . Some of these measures take only the local information into account, while others even employ machine learning methods. Nowadays, it has been one of the most popular research topics and yielded a variety of applications [7] such as identifying essential proteins and potential drug targets for the survival of the cell [13] , controlling the outbreak of epidemics [14] , preventing catastrophic outages in power grids [15] , driving the network toward a desired state [16] , improving transport capacity [17] , promoting cooperation in evolutionary games [18] , etc. This paper investigates the problem of identifying influencers in social networks, by introducing a family of centrality-like measures and gives a brief comparison in Table 1 . Degree Centrality (DC) [19] is the simplest centrality measure, which merely counts how many social connections (i.e., the number of neighbors) a focal node has, defined as where N is the total number of nodes, a ij is the weight of edge (i, j) if i is connected to j, and 0 otherwise. The degree centrality is simple and merely considers the local structure around a focal node [20] . However, this method is probably mistaken for the negligence of global information, i.e., a node might be in a central position to reach others quickly although it is not holding a large number of neighbors [21] . Thus, Betweenness Centrality (BC) [22] is proposed to assess the degree to which a node lies on the shortest path between two other nodes, defined as where g st is the total number of shortest paths, g st (i) is the shortest path between s and t that pass through node i. The betweenness centrality considers global information and can be applied to networks with disconnected components. However, there is a great proportion of nodes that do not lie on the shortest path between any two other nodes, thereby the computational result receives the same score of 0. Besides, high computational complexity is also a limitation of applying for large-scale networks. Analogously, Closeness Centrality (CC) [23] is proposed to represent the inverse sum of shortest distances to all other nodes from a focal node, defined as where N is the total number of nodes, d ij is the shortest path length from node i to node j. The closeness centrality is capable of measuring the core position of a focal node via the utilization of global shortest path length, while it suffers from the lack of applicability to networks with disconnected components, e.g., if two nodes that belong to different components do not have a finite distance between them, it will be unavailable. Besides, it is also criticized by high computational complexity. Eigenvector Centrality [24] (EC) is a positive multiple of the sum of adjacent centralities. Relative scores are assigned to all nodes in a network based on an assumption that connections to high-scoring nodes contribute more to the score of the node than connections to low-scoring nodes, defined as where k i depicts the eigenvalue of adjacency matrix A, x = k 1 Ax depicts the eigenvector stable state of interactions with eigvenvalue k −1 1 . This measure considers the number of neighbors and the centrality of neighbors simultaneously, however, it is incapable of dealing with non-cyclical graphs. In 1998, Brin and Page developed the PageRank algorithm [25] , which is the fundamental search engine mechanism of Google. PageRank (PR) is a positive multiple of the sum of adjacent centralities, defined as where N depicts the total number of nodes, ∑ N i=1 PR 0 (i) = 0, k out j is the number of edges from node j point to i. Likewise, this method is efficient but also criticized by non-convergence in cyclical structures. As we have known, the clustering coefficient [26, 27] is a measure of the degree to which nodes in a graph tend to cluster together, defined as It is widely considered that a node with a higher clustering coefficient may benefit forming communities and enhancing local information spreading. However, Chen et al. expressed contrary views that the local clustering has negative impacts on information spreading. They proposed a ClusterRank algorithm for ranking nodes in large-scale directed networks and verified its superiority to PageRank and LeaderRank [28] . Therefore, the effect of clustering coefficient on information spreading is uncertain, which may benefit local information spreading but prohibit global (especially directional network) information spreading. In 2016, Ma et al. proposed a gravity centrality [29] (GR) by considering the interactions comes from the neighbors within three steps, defined as where k s (i) and k s (j) are the k-shell index of i and j, respectively. ψ i is the neighborhood set whose distance to node i is less than or equal to 3, d ij is the shortest path length between i and j. These methods consider semi-local knowledge of a focal node, i.e., the neighboring nodes within three steps, which are successful in many real-world datasets, such as Jazz [30] , NS [31] and USAir network [32] , etc. However, they are also with high computational complexity by globally calculating k-shell. In 2019, Li et al. improved the gravity centrality and proposed a Local-Gravity centrality (LGR) [33] by replacing k-shell computing and merely considering the neighbors within R steps, defined as where k i and k j are the degrees of i and j, respectively, d ij is the shortest path length between i and j. This method had been extremely successful in a variety of real-world datasets, however, the parameter R requires the calculating of network diameter, which is also a time-consuming process. LGR [33] Semi-local O(n 2 ) simple and capable in most cases additional parameters R determination The above-mentioned centrality measures have been utilized to rank nodes' spreading abilities in monolayer networks. The ranking of nodes in multilayer networks is a more challenging task and is still an open issue. The information propagation process over multiple social networks is more complicated, and conventional models are incapable without any modifications. Zhuang and Yaǧan [36] proposed a clustered multilayer network model, where all constituent layers are random networks with high clustering to simulate the information propagation process in multiple social networks. Likewise, Basaras et al. [37] proposed an improved susceptible-infected-recovered (SIR) model with information propagation probability parameters (i.e., λ ii for intralayer connections and λ ij for interlayer connections). Most of the recent endeavors concentrated on the multiplex networks, (e.g., clustering coefficient in multiplex networks [38] ), where all layers share the identical set of nodes but may have multiple types of interactions. Rahmede et al. proposed a MultiRank algorithm [39] for the weighted ranking of nodes and layers in large multiplex networks. The basic idea is to assign more centrality to nodes that are linked to central nodes in highly influential layers. The layers are more influential if highly central nodes are active in them. Wang et al. proposed a tensor decomposition method (i.e., EDCPTD centrality) [7] , which utilize the fourth-order tensor to represent multilayer networks and identify essential nodes based on CANDECOMP/PARAFAC (CP) tensor decomposition. They also exhibited the superiority to traditional solutions by comparing the performance of the proposed method with the aggregated monolayer networks. In a word, it is of great significance in identifying influencers in multiplex networks. Our purpose in this work is to devise a measure that can accurately detect influential nodes in a general multilayer network. The problem of finding influential nodes is described as extracting a small set of nodes that can bring the greatest influence on the network dynamics. With a given network model is the edge set. The identification of influential nodes is to pick a minimum of nodes as the initial seeds, which can achieve the maximum influenced scope, described as where A is the initially infected nodes, σ(A) denotes the final influenced node set. This problem is simplified as top-k influencers identification by additional setting |A| = k, which has recently attracted great research interests [40] [41] [42] . A variety of real-world social networks are, in fact, interconnected by different types of interactions between nodes, forming what is known as multilayer networks. In this paper, we employ a multilayer network model [9] , which can represent nodes sharing links in different layers. The multilayer network model is defined as where G = {G α ; α ∈ {1, . . . , L}} is a family of (directed or undirected, weighted or unweighted) graphs G α = (V α , E α ), which represents layers of M and C depicts the interactions between nodes of any two different layer, given by The corresponding supra-adjacency matrix can be represented as where A 1 , A 2 , ..., A L are the adjacency matrix of layer 1, 2, ..., L, respectively. N is the total number of the nodes, which can be calculated by N = ∑ 1≤l≤L |V l |. The non-diagonal block I αβ represents the inter-layer edges of layer α and layer β. Thus, the interlayer edges can be represented as Take the 9/11 terrorists network [43] for instance, the edges are classified into three categories (i.e., layers) according to the observed interactions which are plotted in Figure 1 . We employ the susceptible-infected-recovered (SIR) spreading model [44] as the influence analysis model. It has three possible states: • Susceptible (S) state, where a node is vulnerable to infection. • Infectious (I) state, where a node tries to infect its susceptible neighbors. • Recovered (R) state, where a node has recovered (or isolated) and can no longer infect others. In a network, if two nodes are connected then they are considered to have "contact". If one node is "infected", and the other is susceptible, then with a certain probability the latter may become infected through contact [45] . A node is considered to be recovered if it is isolated or immune to the disease. In detail, to check the spreading influence of one given node, we set this node as an infected node and the other nodes are susceptible nodes. At each time step, each infected node can infect its susceptible neighbors with infection probability β, and then it recovered from the diseases with probability γ, the differential equations are shown in Figure 2 . For simplicity, we set γ = 1. The process of the SIR model is plotted in Figures 3 and 4 with the famous Krackhardt's Kite network [46] . , all the nodes are in Susceptible state; while we select one node to be infected, and the neighbors will be infected soon, as shown in panel (b); Finally, the network will reach a stable state, i.e., the number of recovered nodes will reach a maximum, as shown in panel (c). 1XPEHURI6XVFHSWLEOHV,QIHFWLRXVDQG5HFRYHUHGV 6XVFHSWLEOHV ,QIHFWLRXV 5HFRYHUHGV In this paper, we define the node influence (INF, for short) as the energy derived from the neighbors, given as where R is the truncation radius, Γ(i) is the set of neighbors of node i, d ij is the shortest path length between node i and node j, k j is the degree of node j, w ij is the weight of edge e ij . For unweighted networks, w ij = 1. Analogously, we apply the proposed INF measure to multilayer networks (represented as I NF M R ) by the following modifications where R is the truncation radius, Γ α (i) is the set of neighbors of node i at layer α, k α j is the degree of node j at layer α, d ij is the shortest path length between node i and node j. For simplicity, we choose R = 1, thus d ij = 1 if node i and node j is connected through an intralayer edge or interlayer edge, and 0 otherwise. To explain the effect, we take the above-mentioned Krackhardt's Kite network (as plotted in Figure 5 ) and the 9/11 terrorists network (as plotted in Figure 1 ) for examples. The nodes centralities in Krackhardt's Kite network are shown in Table 2 . As shown in Table 2 , Node 4 is considered to be the most important node under the Degree, Katz and the proposed INF measure, while Node 8 has greater Betweenness, Node 6 and node 7 has greater Closeness or (Eigenvector). Thus, the node list (i.e., [4, 6, 7, 8] ) is considered to be the influencers. Furthermore, to evaluate the nodes' influence, we set each node as the initially infected and recorded the final recovered nodes, respectively. This process is repeated for 10,000 times and the results are shown in Table 3 . Table 3 . Averaging recovered nodes and iterations times of each node as initially infected spreaders under 10,000 times SIR stimulations with parameters setting β = 0.35, γ = 1. As shown in Table 3 , Node 4 (i.e., Diane), which is considered to be more influential under Degree, Katz, and INF centrality, shows more recovered nodes (i.e., 5.3182) after 10,000 times SIR stimulations. This experiment is available at https://neusncp.com/api/sir. Analogously, we conduct experiments on the three-layer 9/11 terrorists network. Particularly, we set the infected probability between intralayer edges as β and the probability between interlayer edges as β M = w ij β. The experimental results are plotted in Figure 6 . By conducting SIR simulations on the three-layer 9/11 terrorists network, we can obtain the influential nodes of each layer by calculating the number of finally recovered nodes. Afterward, we sort the nodes by the averaging recovered nodes, and compare the order with the results computed from the proposed INF indicator. It is shown in Figure 6 that the compared values (i.e., recovered nodes and INF) are in the same tendency, which verifies the feasibility of the proposed INF measure. Notably, several influential nodes, such as "Essid Sami Ben Khemais", "Mohamed Atta", and "Marwan Al-Shehhi" are also in the central position of the network, as shown in Figure 1a . The experimental results on the two sample networks show the feasibility of the proposed measure on monolayer and multilayer networks, respectively. Experiments on more real-world networks will be given in Section 4. Suppose m and n are the numbers of edges and nodes, respectively, L is the number of layers, the average degree of nodes is d, R is the truncation radius (commonly setting as R = 1). The complexity of INF for monolayer network is O(n + d R ). As for multilayer networks, the computational complexity is O(n + Ld R ), where L is also a small positive integer. Thus, the time complexity is acceptable as O(n + Ld). Overall, the proposed measure considers more neighboring information than the degree centrality and has a lower computational complexity than betweenness centrality and closeness centrality (i.e., O(nm + n 2 log n)). The experimental environment was with Intel(R) Core (TM) i5-7200U CPU @ 2.50 GHz (4 CPUs), 2.7 GHz, the memory was 8 GB DDR3. The operating system was Windows 10 64 bit, the programming language was Python 3.7.1, and the relevant libs were NetworkX 2.2 and Multinetx. The goal of the experiments was to compare the performance of the proposed INF measure with competitive indicators. In this paper, 21 real-world datasets were employed to verify the performance of the proposed method, which were classified into two groups. The first group covered 12 monolayer networks, which comprised four social networks (i.e., Club, Dolphins, 911 and Lesmis), three biological networks (i.e., Escherichia, C.elegans and DMLC), collaboration networks (i.e., Jazz and NS), a communication network (i.e., Eron), a power network (i.e., Power) and a transport network (i.e., USAir), as shown in Table 4 . Note: |V| and |E| denotes the number of nodes and edges, respectively. is the average degree; is the average shortest path length; |C| is the average clustering index; is the assortativity coefficient; |H| is the degree heterogeneity and β c represents the epidemic threshold of the SIR model. Club contains the friendships between the 34 members of a karate club at a US university. Dolphins dataset is a animals social network. 911 represents a monolayer terrorist network of September 11 attacks. Lesmis is the coappearance network of characters in the novel Les Miserables. Escherichia represetns transcriptional regulation networks in cells orchestrate gene expression, where nodes are operons, and each edge is directed from an operon that encodes a transcription factor to an operon that it directly regulates (an operon is one or more genes transcribed on the same mRNA). Eron is a email network collected from Eron company. Jazz lists the collaboration patterns of jazz musicians. USAir is an undirected weighted network as obtained by considering the 500 US airports with the largest amount of traffic from publicly available data. Nodes represent US airports and edges represent air travel connections among them. NS represents coauthorships between 379 scientists whose research centers on the properties of networks of one kind or another. C.elegans represents the edges of the metabolic network of C.elegans. DMLC represents the inferred Links by small/medium-scale rotein-protein interactions (collected from protein-protein interaction data bases). Power is a power grid of the western United States. The second group covered nine multilayer networks, which comprised six social networks (i.e., Padgett, Krackhardt, Vickers, Kapferer, Lazega and CS-Aarhus), two transport networks (i.e., LondonTransport and EUAirTransportation) and a biological network (i.e., humanHIV), as shown in Table 5 . Data availability: http://www.neusncp.com/user/file?id=12&code=data. To verify the performance of the proposed node influence in networks, this paper carries out a comparison experiment on the above-mentioned datasets: The nodes were removed by a certain indicator in descending order, and the number of subgraphs was recorded. This process repeated until there were not any nodes left. The varying tendency of the subgraphs' number exhibited the influence of a focal centrality. The experimental results are shown in Figure 7 . As shown in Figure 7 , with the nodes removing, the number of subgraphs was increasing and reached a maximum when the network was totally broken up, i.e., there were no edges at this moment. Afterward, the number of subgraphs (i.e., the number of nodes) was decreasing and finally reached zero when all the nodes were removed. The maximum numbers of subgraphs were obtained by the proposed INF measure on all the datasets except C.Elegans. However, the result of C.Elegans obtained by INF was very close to the best situation of BC, which suggests the feasibility of the proposed INF measure. We applied the SIR model to compare the rankings of influences calculated by each indicator among the above-mentioned networks. Initially, one node was set as "infected" state to infect its neighbors with probability β. Afterward, the infected nodes were recovered and never be infected again with probability γ. This spreading process repeated until there were no more infected nodes in the network. The influence of any node i can be estimated by where N R is the number of recovered nodes after the spreading process, and N is the total number of nodes in the network. For simplicity, we set γ = 1 and the epidemic threshold was After having obtained the standard nodes' influence sequence via SIR model simulations, we employed the Kendall's Tau coefficient [65] to compare the performance of each indicator. The Kendall's Tau coefficient is an index measuring the correlation strength between two sequences. Suppose given the standard sequence X = (x 1 , x 2 , . . . , x N ), and we obtained the computational sequence Y = (y 1 , y 2 , . . . , y N ) by a certain indicator. Any pair of two-tuples (x i , y i ) and (x j , y j ) (x = j) are concordant if both x i > x j and y i > y j or x i < x j and y i < y j . Meanwhile, they are considered as discordant, if x i > x j and y i < y j or x i < x j and y i > y j . If x i = x j or y i = y j , pairs are neither concordant nor discordant. Therefore, Kendall's Tau coefficient is defined as where N c and N d indicate the number of concordant and discordant pairs, respectively. The range of τ is [−1, 1]. Table 6 shows the computational Tau results with the comparison of standard sequence from SIR model simulations. q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q qq q qqq qq q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q C.Elegans As shown in Table 6 , the proposed measure outperformed the competitors in most cases, even in the Escherichia network, the computed Tau result of INF (0.0692) was close to that of CC (0.0971). Thus, it was also competitive in this network. If setting the limitation of identifying k influencers, we conducted experiments on the real-world datasets with top-k nodes by computational centrality nodes and compared the recovered nodes (i.e., the final number of nodes with recovered states). To compare the varying parameter k with the obtained τ, we conducted experiments on the above-mentioned datasets and set the ratio of β/β c , as shown in Figure 8 . Note: Given a network, the parameters of SIR model are given with the transmission probability β = 0.35 and recovering probability µ = 1 for simplicity. To obtain the standard ranking of nodes' influences, we conducted 1000 independent simulations, in each process every node is selected once as the infect seed once. The best perfromed indicator for each network is emphasized by bold. As shown in Figure 8 , the proposed node influence method is quite competitive in most of the datasets, although second to the performance of betweenness indicator in DMLC and Jazz datasets. Analogously, we conducted experiments on the nine multilayer networks by removing nodes with maximum centralities; the results are plotted in Figure 9 . q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q CS−Aarhus LondonTransport EUAirTransportation As shown in Figure 10 , the runtime accumulated from either group indicated that the proposed INF measure was efficient, which was close to that of DC and superior to BC, CC and LGR. Influencers identification is a fundamental issue with wide applications in different fields of reality, such as epidemic control, information diffusion, viral marketing, etc. Currently, degree centrality [19] is the simplest method, which considers nodes with larger degrees are more influential. However, for the lack of global information, a node lying in a "bridge" position might be neglected for holding a small degree. The betweenness [22] and closeness [23] centrality consider global information, but they are holding a high complexity, which are not suitable for applications in large-scale networks. Local gravity is a balanced method, however, the determination of parameter R requires computing network diameter, which is also time-consuming. Thus, a novel node influence measure is proposed in this paper, which merely considers the local neighboring information of a focal node with the complexity of O(n + Ld). Experimental results on 21 real-world datasets indicate the feasibility of the proposed measure. Firstly, the experiments of counting subgraphs with removing influential nodes show that the capability of the proposed INF measure. By removing the nodes according to the INF indicator, the networks are more easily broken up, as shown in Figures 7 and 9 . Secondly, we apply the SIR model to evaluate the node influence, which suggested the proposed INF measure is competitive to other indicators in most cases. Although inferior to BC on Jazz and DMLC networks, it is also competitive. By analyzing the structures of these two networks, we find that the nodes of Jazz network are densely connected (i.e., the average degree of 27.6970) and most of the nodes are holding the same number of neighbors (approximately 28 neighbors), which brings difficulties to identify which node is more influential. On the contrary, there is only one node (i.e., Node 2) holding a large number of neighbors (i.e., 439 neighbors) and the others only holding few neighbors (approximately four neighbors) in DMLC network, which is also difficult to identify influencers. Overall, the proposed method outperforms the other indicators in most cases. Finally, we compare the running time of each indicator on the 21 real-world datasets. Experimental results show the efficiency of the proposed measure. Aiming at solving the problem of identifying influencers in social networks, this paper proposes a novel node influence indicator. This method merely considers the local neighboring information in order to be fast and suitable for applications in large-scale networks. Extensive experiments on 21 real-world datasets are conducted, and the experimental results show that the proposed method outperforms competitors. Afterwards, the time complexity is compared, and we verify the efficiency of the proposed indicator. Overall, the proposed node influence indicator is capable of identifying influencers in social networks. The contribution of this work is likely to benefit many real-world social applications, such as promoting network evolutions, preventing the spreading of rumors, etc. As part of future works, the influencers in dynamic networks can be further studied by applying the proposed INF measure into a multilayer network model with numerous ordinal layers. The node's influence can be calculated by accumulating the local neighbors across all the layers. Besides, the effect of layers needs to be taken into consideration. In a word, we hope the findings in this work will help to improve the researches in this promising field. Surprise maximization reveals the community structure of complex networks Searching for superspreaders of information in real-world social media Superspreaders and superblockers based community evolution tracking in dynamic social networks Locating influential nodes in complex networks BridgeRank: A novel fast centrality measure based on local structure of the network Social network coalescence based on multilayer network model Identifying key nodes in multilayer networks based on tensor decomposition Multilayer networks The structure and dynamics of multilayer networks Network centrality: an introduction. In A Mathematical Modeling Approach from Nonlinear Dynamics to Complex Systems Influence analysis in social networks: A survey Scalable feature learning for networks Lethality and centrality in protein networks Identification of influential spreaders in complex networks Structural vulnerability of the North American power grid Controllability of complex networks Efficient routing on complex networks Diversity-maintained differential evolution embedded with gradient-based local search. Soft Comput Factoring and weighting approaches to status scores and clique identification Birds of a feather: Homophily in social networks A set of measures of centrality based on betweenness Power and centrality: A family of measures The anatomy of a large-scale hypertextual web search engine Collective dynamics of 'small-world'networks The architecture of complex weighted networks Identifying influential nodes in large-scale directed networks: The role of clustering Identifying influential spreaders in complex networks based on gravity formula Community structure in jazz Modularity and community structure in networks Reaction-diffusion processes and metapopulation models in heterogeneous networks Identifying influential spreaders by gravity model An index to quantify an individual's scientific research output Leaders in social networks, the delicious case Information propagation in clustered multilayer networks Identifying influential spreaders in complex multilayer networks: A centrality perspective Clustering coefficients in multiplex networks. arXiv 2013 Centralities of nodes and influences of layers in large multiplex networks Identifying Top-k Most Influential Nodes by using the Topological Diffusion Models in the Complex Networks Resampling-based predictive simulation framework of stochastic diffusion model for identifying top-K influential nodes Identification of top-k influential nodes based on enhanced discrete particle swarm optimization for influence maximization The Al Qaeda Factor: Plots against the West The mathematics of infectious diseases Introduction to Complex Networks: Models, Structures and Dynamics Assessing the political landscape: Structure, cognition, and power in organizations An information flow model for conflict and fission in small groups The bottlenose dolphin community of Doubtful Sound features a large proportion of long-lasting associations Social Network Analysis for Startups: Finding Connections on the Social Web The Stanford GraphBase: A Platform for Combinatorial Algorithms Network motifs in the transcriptional regulation network of Escherichia coli The network data repository with interactive graph analytics and visualization Community detection in complex networks using extremal optimization WormNet v3: A network-assisted hypothesis-generating server for Caenorhabditis elegans The Rise of the Medici Strategy and Transaction in an African Factory: African Workers and Indian Management in a Zambian Town The Collegial Phenomenon: The Social Mechanisms of Cooperation Among Peers in a Corporate Law Partnership New specifications for exponential random graph models BioGRID: a general repository for interaction datasets Combinatorial analysis of multiple networks. arXiv 2013 Navigability of interconnected networks under random failures Emergence of network features from multiplexity A new measure of rank correlation We would like to thank the anonymous reviewers for their careful reading and useful comments that helped us to improve the final version of this paper. The authors declare that there are no conflicts of interest regarding the publication of this paper.