PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
 
Sci Rep. 2013; 3: 2171.
Published online 2013 July 10. doi:  10.1038/srep02171
PMCID: PMC3707349

Global efficiency of local immunization on complex networks

Abstract

Epidemics occur in all shapes and forms: infections propagating in our sparse sexual networks, rumours and diseases spreading through our much denser social interactions, or viruses circulating on the Internet. With the advent of large databases and efficient analysis algorithms, these processes can be better predicted and controlled. In this study, we use different characteristics of network organization to identify the influential spreaders in 17 empirical networks of diverse nature using 2 epidemic models. We find that a judicious choice of local measures, based either on the network's connectivity at a microscopic scale or on its community structure at a mesoscopic scale, compares favorably to global measures, such as betweenness centrality, in terms of efficiency, practicality and robustness. We also develop an analytical framework that highlights a transition in the characteristic scale of different epidemic regimes. This allows to decide which local measure should govern immunization in a given scenario.

Epidemics never occur randomly. Instead, they follow the structured pathways formed by the interactions and connections of the host population1,2. The spreading processes relevant to our everyday life take place on networks of all sorts: social (e.g. epidemics3,4), technological (e.g. computer viruses5,6) or ecological (cascading extinctions in food webs7). With a network representation, these completely different processes can be modelled as the propagation of a given agent on a set of nodes (the population) and links (the interactions). Different systems imply networks with different organizations, just as different agents require different epidemic models.

There has long been significant interest in identifying the influential spreaders in networks. Which nodes should be the target of immunization efforts in order to optimally protect the network against epidemics? Unfortunately, most studies feature two significant shortcomings. Firstly, the proposed methods are often based on optimization or heuristic algorithms requiring nearly perfect information on a static system8,9; this is rarely the case. Secondly, methods are usually tested on small numbers of real systems using a particular epidemic scenario10,11; this limits the scope of possible outcomes.

We first present a numerical study, perhaps the largest of its kind to date, where we argue that, depending on the nature of the network and of the disease, different immunization tactics have to be taken into consideration. In so doing, we formalize the notion of node influence and illustrate how local knowledge around a particular node is usually sufficient to estimate its role in an epidemic. We also show how, in certain cases, the influence of a node is not necessarily dictated by its number of connections, but rather by its role in the network's community structure (see Fig. 1). Far from trivial, it follows that an efficient immunization strategy can be obtained solely from local measures, which are easily estimated in practice and robust to noisy or incomplete information. We further develop an analytical formalism ideally suited to test the effects of local immunization on realistic network structures. Combining the insights gathered from the numerical study and this formalism, we finally formulate a readily applicable approach which can easily be implemented in practice.

Figure 1
Protein interactions of S. cerevisiae (subset)22.

Results

Models and measures

There exist two standard models emulating diverse types of epidemics: the susceptible-infectious-recovered (SIR) and susceptible-infectious-susceptible (SIS) dynamics. In both, an infectious node has a given probability of eventually infecting each of its susceptible neighbors during its infectious period, which is terminated by either death/immunity leading to the recovered state (SIR) or by returning to a susceptible state (SIS). In the SIR dynamics, for a given transmission probability T, the quantity of interest is the mean fraction Rf of recovered nodes once a disease, not subject to a stochastic extinction, has finished spreading (i.e. we focus on the giant component12). Since each edge can only be followed once, this dynamics investigates how a population is vulnerable to the invasion of a new pathogen. In the SIS dynamics, we are interested in the prevalence I* (fraction of infectious nodes) of the disease at equilibrium (equal amounts of infections and recoveries) as a function of the ratio λ = α/β of infection rate α and recovery rate β. This particular dynamics permits the study of how a given network structure can sustain an already established epidemic.

Should a fraction ε of the population be fully immunized, our objective is to identify the nodes whose absence would minimize Rf and I*. The epidemic influence of a node — that is the effect of its removal on Rf and I* — depends mainly on its role in the organization of the network. Hence to efficiently immunize a population, we must first understand its underlying structure.

Network organization can be characterized on different scales, each of which affect the dynamics of propagation. At the microscopic level, the most significant feature is the degree of a node (its number of links, noted k) which in turn defines the degree distribution of the network. The significance of the high-degree nodes (the hubs) for network structure in general13, for network robustness to random failure14 and for epidemic control15 has long been recognized.

At the macroscopic level, the role of a node can be described by its centrality, which may be defined in various ways. Frequently used in the social sciences is the betweenness centrality (b), quantifying the contributions of a given node to the shortest paths between every pair of nodes in the network16. Arguably, this method should be among the best estimate of a node's epidemic influence as it directly measures its role in the different pathways between all other individuals17, yet at a considerable computational cost. A simpler method, the k-core (or k-shell) decomposition18,19, assigns nodes to different layers (or coreness c) effectively defining the core and periphery of a network (high and low c respectively). It has recently been shown that coreness is well suited to identify nodes that are the most at risk of being infected during the course of an epidemic20. In light of our results, we will be able to discuss the distinction between a node's vulnerability to infection and its influence on the outcome of an epidemic.

The mesoscopic scale has recently been the subject of considerable attention. At this level of organization, the focus is on the redundancy of connections forming dense clusters referred to as the community structure of the network21,22. Nodes can be distinguished by their membership number m, i.e., the number of communities to which they belong. We will consider that two links of a given node are part of one community if the neighbours they reach lead to significantly overlapping neighbourhoods21. This definition is directly relevant to epidemic dynamics as links within communities do not lead to new potential infections. We call structural hubs the nodes connecting the largest number of different communities. These nodes act as bridges facilitating the propagation of the disease from one dense cluster to another. Targeting structural hubs to hinder propagation in structured populations has been previously proposed and investigated10,11, but has yet to be tested extensively.

Note that the microscopic and mesoscopic levels (as defined above) are characterized by local measures in the sense that they do not require a complete knowledge of the network, in contrast to global measures like the betweenness centrality. Moreover, as we will see, local measures are less sensitive to incomplete or incorrect information. Adding, removing or rewiring a link only affects the degree or membership of nodes directly in the neighbourhood of the modification; whereas the same alterations can potentially affect the centrality of nodes anywhere in the network through cascading effects. Furthermore, even if community detection often requires the tuning of a global resolution parameter, we will see that this additional step does not affect the identification of structural hubs, meaning that local information is sufficient to accurately determine a node's memberships.

In our numerical simulations we will have a perfect knowledge of static networks. This will allow us to use global measures as a reference to test the efficiency of local measures best suited in practice. We therefore ask without discrimination: which of the degree k, the coreness c, the betweenness centrality b or the membership number m is the best identifier of the most influential nodes on the outcome of an epidemic? To answer this question, we have simulated SIR and SIS dynamics with Monte Carlo calculations on 17 real-world networks. In each case, a fraction ε of the nodes was removed in decreasing order of the nodes' score for each of the four different measures. By comparing their efficiency to reduce Rf or I* as a function of ε, we are able to establish which measure is best suited for a given scenario characterized by a network structure, a propagation dynamics and a disease transmissibility (i.e. probability of transmission).

Case study: a data exchange network

We first illustrate our methods using the network of users of the Pretty-Good-Privacy algorithm for secure information interchange (hereafter, the PGP network)23, which could be the host of the propagation of computer viruses, rumors or viral marketing campaigns. Results for the 16 other networks are presented and discussed in the next section as well as in the Supporting Information (SI) document.

Communities in the network are extracted with the link community algorithm of Ahn et al.21. This algorithm groups links — and therefore the nodes they join — into communities based on the overlap of their respective neighbouring nodes. It is this overlap that reduces the number of new potential infections in a community structure, as opposed to a random network. This method thus reflects our understanding of how communities affect disease propagation. While it may not directly detect the social groups or functional modules of a network, it identifies significant clusters of redundant links. This redundancy or overlap is quantified through a Jaccard coefficient, and two links are grouped into the same community when their coefficient exceeds a certain threshold. The threshold value acts as a resolution, enabling to look at different levels of organization. As suggested21, the value of the threshold is chosen to maximize the average density ρ of the communities (see Methods). As this choice may seem arbitrary, Fig. 2 investigates the similarity between the nodes with the highest membership numbers, for different thresholds. It suggests that the membership number is fairly robust around the threshold. Moreover, Fig. 2 also demonstrates that the effect of the removal of the structural hubs on a SIS epidemics is very robust to the choice of the threshold. Thus, we will henceforth use the membership numbers obtained with the threshold value corresponding to the highest community density.

Figure 2
Robustness of structural hubs in the PGP network.

The differences, if any, between the efficiency of the different methods are due to the immunized nodes not being the same. Figure 3 (top) investigates the correlations between the different properties (k, b, c and m) of each node. Perhaps the most important result here is that nodes with a high membership number may have relatively small degree, coreness and betweenness centrality. Hence, we expect the immunizing method based on community structure to have a different influence on the outcome of epidemics. Figure 3 (bottom) shows the consistensy (or lack thereof) of a given measure, depending on the quality of the available data. The robustness of local (micro and meso) measures is of obvious practical advantage. Both robustness and correlations are further investigated in the SI.

Figure 3
Difference in immunization targets for the PGP network.

To study various epidemic scenarios, we consider both SIS and SIR dynamics (which may behave quite differently) with different values of the transmission probability (λ and T for SIS and SIR, respectively). In fact, each network features an epidemic threshold, i.e. critical values λc24 and Tc25, below which I* and Rf vanish to zero in an equivalent infinite network ensemble. As we will show, the observed behavior can differ significantly depending whether or not λ and T are close to their critical value.

Figure 4 presents results of different immunization methods against SIS dynamics for different values of λ. On the top figure, where λ is near λc, the most successful method of intervention is to target nodes according to their degree. At low transmissibility, the disease follows only a very small fraction of all links. The shortest paths are seldom used and the poor performance of betweenness centrality follows. Moreover, the disease will not be affected by the community structure, because even in dense neighbourhoods, most links will not be travelled. We then say that the disease, unaffected by link clustering, follows a tree-like structure (without loops), where community memberships are insignificant. It is therefore better to simply remove as many links as possible.

Figure 4
Efficiency of the immunization methods against an SIS epidemics on the PGP network.

As λ increases beyond λc, we see that immunization based on membership numbers quickly outperforms the other methods. As more links are travelled, the disease is more likely to follow superfluous links in already infected communities. Hubs sharing their many links within few communities are therefore not as efficient in causing secondary infections as one might expect. Similarly, targeting through betweenness centrality also performs better with higher λ, albeit not as well as membership-targeting in this case. For λ [dbl greater-than sign] λc, immunization based on membership numbers (local) and on betweenness centrality (global) converge toward similar efficiency, significantly outperforming degree-based immunization.

Another interesting feature of our results is the poor performance of immunization based on node coreness. A previous study had clearly shown that epidemics mostly flourished within the core of the network (see Fig. 5) because of its density20. Ironically, this density also implies redundancy. While the core nodes are highly at risk of being infected, their removal has a limited effect because there exist alternative paths within their neighbourhood: the core offers a perfect environment to the disease and is consequently robust to node removal. It is therefore more effective to stop the disease from reaching, or leaving, the core by removing the nodes bridging other neighbourhoods (i.e. the structural hubs).

Figure 5
k-core decomposition of the PGP network.

Similar conclusions are drawn for the SIR dynamics. As T moves away from Tc, the most significant level of organisation shifts from the degree (microscopic) to communities (mesoscopic) as membership-based immunization progressively outperforms the other strategies.

Results on networks of diverse nature

In this section, we highlight different behaviours observed in social, technological and communication networks using 7 other datasets (full results for the 17 datasets are available in the SI): subset of the World Wide Web (WWW)13, MathSciNet co-authorship network (MathSci)27, Western States Power Grid of the United States (Power Grid)28, Internet Movie Database since 2000 (IMDb)29, cond-mat arXiv co-authorship network (arXiv)22, e-mail interchanges between members of the University Rovira i Virgili (Email)30 and Gnutella peer-to-peer network (Gnutella)31.

The results for the WWW, MathSci and IMDb networks further support our previous conclusions, with the exception that membership-based immunization performs surprisingly better than the degree-based variant even near the epidemic threshold of the network (see WWW and MathSci). The betweenness-centrality-based immunization was not tested on IMDb because of computational constraints (its computation required over 800 hours with our available ressources and a standard algorithm32), which illustrates a significant limit of this measure. Approximations could have been used33, but the intricate (and mostly unknown) relationship between the efficiency of the measure and the accuracy of the approximation would have only caused additional uncertainties.

The results presented for the Power Grid network illustrate a fundamental difference between the SIS and the SIR dynamics: while we are interested in the fraction of the network sustaining an established epidemic in SIS, it is the fraction of nodes invaded by a new disease that is relevant in SIR. In fact, the structure of the Power Grid, a chain of small, easily disconnected modules, enhances the qualitative discrepancy between the epidemic influence of nodes subjected to these two dynamics. For the SIS dynamics, the membership-based intervention is the most efficient because it weakens all modules, limiting the prevalence of the disease. In distinction, targeting through betweenness centrality merely separates the modules, so that they indiviually remain infected. For the SIR dynamics, separating the modules is the best approach as it directly stops the infection from spreading; while weakened – but connected – modules still provide pathways. This effect is a direct consequence of the particular structure of the Power Grid and is insignificant on other networks.

Finally, the last set of results, on arXiv, Email and Gnutella, present the effect of the community density ρ on the performance of membership-based immunization. For very small ρ, the paths within communities do not qualitatively differ from the links bridging neighborhoods in their effect on the disease propagation. This targeting method is therefore expected to converge toward degree-based immunization if m and k are strongly correlated. However, as most tested networks had fairly dense communities, ρ ≥ 0.3, the relevance of memberships should not be understated.

Investigation of the epidemic regimes transition

The results of the previous sections suggest that local information (i.e., degree, membership) is often sufficient for a nearly optimal global immunization. More precisely, we found these two methods to outperform or to be as efficient as the betweenness centrality (the global method used for comparison) in 62 of the 68 studied scenarios (i.e., 17 networks / 2 dynamics / 2 transmissibility regimes). This implies that membership (e.g., on PGP), degree (e.g., Gnutella) or both (e.g. MathSci) lead to an immunization at least as efficient as global methods while having the noteworthy advantage of requiring much less information and of being less sensitive to incomplete information. This section focuses on the conditions guiding the choice between the degree-based or the membership-based immunization strategy. In this respect, Figs. 4 and and66 provide a useful hindsight: the membership-based strategy is more efficient than the degree-based one when transmissibility is high and/or when communities are dense. To further our understanding and test this hypothesis, we introduce a random network model featuring a community structure, and exactly solve its final state (Rf) under SIR dynamics using generating functions.

Figure 6
Efficiency of the immunization methods against SIS and SIR epidemics on several networks.

Our model is a slightly modified version of the configuration model12,34 where nodes are connected either through single links or through motifs (see Fig. 7 for an example). Motifs are used to simulate the effect of a community structure, that is the redundancy of the neighbourhoods of nodes. Our motifs are composed of M nodes, all connected to each other, and a node belongs to i motifs and has j single links with probability p(i,j). This node therefore has a degree (k) equal to (M − 1)i + j and a membership (m) equal to i + j. Networks are generated using a stub pairing scheme: a node belonging to i motifs and having j single links has i “motif stubs” and j “link stubs”. Groups and single links are then formed by randomly choosing M motif stubs and 2 link stubs, respectively, and then by linking the corresponding nodes to one another. This last step is repeated until none of the motif and link stubs remains. The distribution An external file that holds a picture, illustration, etc.
Object name is srep02171-m5.jpg therefore defines a maximally random network ensemble, and the results obtained are averaged over this ensemble.

Figure 7
Synthetic networks with tunable community structure.

Extending previous work35, we compute the expected value of Rf for the network ensemble just defined where nodes and links are randomly removed to simulate immunization and disease transmission (SIR dynamics), respectively. Full details are given in the SI. Using typical values for {p(i,j)}, our model illustrates and confirms our hypothesis by clearly showing in Fig. 8a transition of efficiency between the degree-based and the membership-based immunization strategy. Initially less efficient when the transmissibility is low (i.e., higher threshold, lower value of Rf), membership progressively outperforms degree as the transmissibility increases. As mentionned above, for lower values of T, the best option is therefore to immunize the hubs (high k) to shift the degree distribution towards lower degrees. For higher values of T, targeting structural hubs (high m) that act as bridges between “independent” neighbourhoods leads to a more efficient immunization as it reduces the number of paths between different regions of the network. Note that we do not explicitly model the effect of community density. This could have been done by letting links exist independently with a given probability η. This is however identical to letting the disease propagate with probability ηT. Thus, the value of T in Fig. 8 is related to the density of the communities, and our conclusions can therefore be extended to the cases of low/high community densities.

Figure 8
Results of local immunization methods on synthetic networks.

Discussion

One of the main contributions of this work is to offer a formal definition of the epidemic influence of nodes, i.e. the effect of its removal on I* of Rf, which is open to diverse methods of approximation. Our results confirm that standard measures such as the degree or betweenness centrality are not always the best indicators of a node's influence. Moreover, we have highlighted that the coreness, which has recently been proposed as an indicator of nodes' influence20, offers poor performances. This has brought us to distinguish between individual risk and global influence. We have also illustrated how a universal approach is still wanting, since different networks and different diseases require different methods of intervention.

Consequently, the fact that the numbers of links and/or communities to which a node belongs are excellent measure of its epidemic influence — at times better, at times equivalent, but never much worse than global centrality measures — is a particularly important result. The fact that they both are local measures is especially relevant considering that we rarely have access to the exact network structure of a system, either because it is simply too large (WWW), too dynamic (email networks) or because the links themselves are ill-defined (social networks). Not only are local measures computable from a limited subset of a network (which is often the only available information), but a coarse-grained measure like membership is even more interesting as it is easier to estimate than a node's actual degree. For instance, consider how much simpler it is to enumerate your social groups (work, family, etc.) than the totality of your acquaintances.

Finally, the existence of a transition between two epidemic regimes with different characteristic scales may well be the single most important conclusion of this work. In the first regime, for low transmissibility and sparse communities, the microscopic structural features (i.e. node connectivity or degree) offer the most relevant information; while for higher transmissibility and denser communities, mesoscopic features (i.e node communities or membership) appear more relevant. We expect to see an equivalent transition between any pair of measures which oppose the micro and meso scales (e.g. different range-limited measures of centrality36).

Based on our empirical and analytical results, we thus propose a simple procedure on how to judge which local measure can be expected to yield the best results in a given situation. From the available subset of a given network:

  1. Obtain the degree distribution to estimate the transmissibility of the disease in relation to the epidemic threshold λc24 or Tc25.
  2. If easily transmissible (λ [dbl greater-than sign] λc or T [dbl greater-than sign] Tc), evaluate the network's community structure; otherwise, go to 4.
  3. If the community density is high (An external file that holds a picture, illustration, etc.
Object name is srep02171-m6.jpg), immunize nodes according to their memberships; otherwise, go to 4.
  4. For a transmissibility near the epidemic threshold, or for sparse communities (low ρ), immunize according to the degree of the nodes.

The analytical and numerical frameworks used in this work are expected to guide immunization efforts toward simpler, more precise and efficient strategies. Likewise, the introduction of a node influence classification scheme opens a new avenue for finding better local estimates of a node's role in the global state of its system.

Methods

Betweenness centrality

For all pairs (a,b) of nodes excluding i, list the na,b shortest paths between a and b. Let na,b(i) be the number of these paths containing i. The betweenness centrality bi of node i is then given by:

An external file that holds a picture, illustration, etc.
Object name is srep02171-m1.jpg

Coreness

The coreness of node i is the highest integer ci such that the node is part of the set of all nodes with at least ci links within the set.

Community detection

Two links, eij and eik, from a given node i, are said to belong to the same community if their Jaccard coefficient J(eij,eik) (similarity measure) is above a given threshold Jc:

An external file that holds a picture, illustration, etc.
Object name is srep02171-m2.jpg

where n+(u) is the set containing the neighbors of u including u.

Community density

The density ρi of a community i of ni > 2 nodes and di links is the proportion of the possible redundant links that do exist; i.e., the fraction of existing links excluding the minimal ni – 1 links that are needed for this community to be connected:

An external file that holds a picture, illustration, etc.
Object name is srep02171-m3.jpg

The community density ρ is then calculated according to

An external file that holds a picture, illustration, etc.
Object name is srep02171-m4.jpg

where D is the total number of links not belonging to single link communities, for which ρi = 021.

Immunization

To perform the immunization of a fraction ε of the network according to a certain measure Γ, we remove the εN nodes with the highest Γ. When a choice must be made (nodes with equal Γ), all decisions are taken randomly and individually for each simulated epidemics.

Monte Carlo simulations

To investigate the fraction of a network which can sustain an epidemics, SIS simulations start with all nodes in an infectious state and are then relaxed until an equilibrium is reached. To investigate the mean fraction of a network which a disease can invade, SIR simulations start with a single randomly chosen infectious node and run until there are no more infectious nodes. Results shown in the figures are obtained by averaging over the outcome of several numerical simulations until the minimal possible standard deviation (limited by network structure and finite size) is obtained. For the SIR dynamics, only the simulations leading to largescale epidemics (at least 1% of the nodes) were considered. The complete procedure is given in the SI.

Author Contributions

L.H.-D. and A.A. designed the study. L.H.-D., A.A. and J.-G.Y. performed the computations. All authors have contributed to the analysis and wrote the manuscript.

Supplementary Material

Supplementary Information:

Global efficiency of local immunization on complex networks: Supplementary Information

Acknowledgments

The authors wish to thank Louis Roy for the development of a k-core visualization tool; Yong-Yeol Ahn et al. for their link community algorithm; all the authors of the cited papers for providing their network data; and Calcul Québec for computing facilities. This research was funded by CIHR, NSERC and FRQ-NT.

References

  • Caldarelli G. & Vespignani A. Large Scale Structure and Dynamics of Complex Networks. World Scientific Publishing Company, Singapore (2007).
  • Keeling M. J. & Rohani P. Modeling Infectious Diseases in Humans and Animals Princeton University Press, Princeton (2008).
  • Anderson R. M. & May R. M. Infectious Diseases of Humans: Dynamics and Control Oxford University Press, New York (1991).
  • Keeling M. J. & Eames K. T. D. Networks and epidemic models. J R Soc Interface 2, 295–307 (2005). [PMC free article] [PubMed]
  • Pastor-Satorras R. & Vespignani A. Epidemic Spreading in Scale-Free Networks. Phys. Rev. Lett. 86, 3200–3203 (2001). [PubMed]
  • Gómez-Gardeñes J., Echenique P. & Moreno Y. Immunization of real complex communication networks. Eur. Phys. J. B 49, 259–264 (2006).
  • Dunne J. A. & Williams R. J. Cascading extinctions and community collapse in model food webs. Philos Trans R Soc Lond B Biol Sci 364, 1711–1723 (2009). [PMC free article] [PubMed]
  • Gallos L. K., Liljeros F., Argyrakis P., Bunde A. & Havlin S. Improving immunization strategies. Phys. Rev. E 75, 045104(R) (2007). [PubMed]
  • Chen Y., Paul G., Havlin S., Liljeros F. & Stanley H. E. Finding a Better Immunization Strategy. Phys.Rev. Lett. 101, 058701 (2008). [PubMed]
  • Salathé M. & Jones J. H. Dynamics and Control of Diseases in Networks with Community Structure. PLoS comp. biol. 6, e1000736 (2010). [PMC free article] [PubMed]
  • Masuda N. Immunization of networks with community structure. New J. Phys 11, 123018 (2009).
  • Newman M. E. J., Strogatz S. H. & Watts D. J. Random graphs with arbitrary degree distributions and their applications. Phys. Rev. E 64, 026118 (2001). [PubMed]
  • Barabási A.-L. & Albert R. Emergence of scaling in random networks. Science 286, 509–512 (1999). [PubMed]
  • Albert R., Jeong H. & Barabási A.-L. Error and attack tolerance of complex networks. Nature 406, 378–382 (2000). [PubMed]
  • Pastor-Satorras R. & Vespignani A. Immunization of complex networks. Phys. Rev. E 65, 036104 (2002). [PubMed]
  • Freeman L. Centrality in social networks: Conceptual clarification. Social Networks 1, 215–239 (1979).
  • Barthélemy M. Betweenness centrality in large complex networks. Eur. Phys. J. B 38, 163–168 (2004).
  • Batagelj V. & Zaveršnik M. Generalized Cores.. arXiv:cs/0202039v1.
  • Batagelj V. & Zaveršnik M. An O(m) Algorithm for Cores Decomposition of Networks. arXiv:cs/0310049v1.
  • Kitsak M. et al. Identification of influential spreaders in complex networks. Nature Physics 6, 888–893 (2010).
  • Ahn Y.-Y., Bagrow J. P. & Lehmann S. Link communities reveal multiscale complexity in networks. Nature 466, 761–764 (2010). [PubMed]
  • Palla G., Derényi I., Farkas I. & Vicsek T. Uncovering the overlapping community structure of complex networks in nature and society. Nature 435, 814–818 (2005). [PubMed]
  • Boguñá M., Pastor-Satorras R., Díaz-Guilera A. & Arenas A. Models of social networks based on social distance attachment. Phys. Rev. E 70, 056122 (2004). [PubMed]
  • Hébert-Dufresne L., Noël P.-A., Allard A., Marceau V. & Dubé L. J. Propagation dynamics on networks featuring complex topologies. Phys. Rev. E 82, 036115 (2010). [PubMed]
  • Newman M. E. J. Spread of epidemic disease on networks. Phys. Rev. E 66, 016128 (2002). [PubMed]
  • Alvarez-Hamelin I., Dall'Asta L., Barrat A. & Vespignani A. k-core decomposition: A tool for the visualization of large scale networks. Advances in Neural Information Processing Systems 18, 41–50 (2006).
  • Palla G., Farkas I. J., Pollner P., Derényi I. & Vicsek T. Fundamental statistical features and self-similar properties of tagged networks. New J. Phys. 10, 123026 (2008).
  • Watts D. J. & Strogatz S. H. Collective dynamics of small-world networks. Nature 393, 440–442 (1998). [PubMed]
  • Hébert-Dufresne L., Allard A., Marceau V., Noël P.-A. & Dubé L. J. Structural Preferential Attachment: Network Organization beyond the Link. Phys. Rev. Lett. 107, 158702 (2011). [PubMed]
  • Guimera R., Danon L., Diaz-Guilera A., Giralt F. & Arenas A. Self-similar community structure in a network of human interactions. Phys. Rev. E 68, 065103(R) (2003). [PubMed]
  • Ripeanu M., Foster I. & Iamnitchi A. Mapping the Gnutella Network: Properties of Large-Scale Peer-to-Peer Systems and Implications for System Design. IEEE Internet Computing Journal 6, 50–57 (2002).
  • Brandes U. A Faster Algorithm for Betweenness Centrality. J. Math. Sociol. 25(2), 163–177 (2001).
  • Madduri K., Ediger D., Jiang K., Bader D. A. & Chavarría-Miranda D. G. A Faster Parallel Algorithm and Efficient Multithreaded Implementations for Evaluating Betweenness Centrality on Massive Datasets. Third Workshop MTAAP (2009).
  • Newman M. E. J. Properties of highly clustered networks. Phys. Rev. E 68, 026121 (2003). [PubMed]
  • Allard A., Hébert-Dufresne L., Noël P.-A., Marceau V. & Dubé L. J. Bond percolation on a class of correlated and clustered random graphs. J. Phys. A: Math. Theor. 45, 405005 (2012).
  • Ercsey-Ravasz M., Lichtenwalter R. N., Chawla N. V. & Toroczkai Z. Range-limited centrality measures in complex networks. Phys. Rev. E 85, 066103 (2012). [PubMed]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group