Weighted Rich-club Effect

tnet » Two-mode Networks » Weighted Rich-club Effect

The weighted rich-club assesses whether, and the extent to which, of “prominent” nodes (e.g., the ones with high (out-)degree or node (out-)strength for (directed) networks) exchange among themselves the strongest ties in a network (for a review, see Weighted Rich-club in One-mode Networks). Since certain ways of defining the prominent nodes can be associated with having strong ties, we divided the measure calculated on the observed network by the measure calculated on an ensemble of random networks (similar to the topological rich-club measure that we built on; Colizza, 2006). To ensure that the random networks were comparable to the observed network, we constrained the random networks so that each node maintained its prominence (e.g., degree or node strength). When prominence was defined as degree, this was straight-forward by either globally reshuffling the weights in the network (the topology remained, and hence, the degree of nodes) or reshuffling ties while maintaining nodes’ degree using Molloy and Reed’s (1995) method.

Directed weight reshuffleHowever, when prominence was defined as node strength, we could not use either of these methods for undirected networks. (For directed networks, it is possible to redirect a node’s out-going ties randomly in the network, thus maintaining the node strength). Therefore, for undirected networks, we created directed networks by linking connected nodes with two directed ties. The weight of an undirected tie is duplicated to the two directed ones. A transformed network based on this process is is illustrated in the top panel of the diagram on the right. Then we reshuffled the weights attached to the outgoing ties of a node, thus, maintaining the sum of weights on outgoing ties. For example, in the bottom panel of the diagram (the weight of a directed tie is placed close to the origin of the tie in this figure) we have randomly the weights of node E’s outgoing ties, 1 and 5. This method has been refered to as the directed weight reshuffling method (Serrano et al., 2007). However, this method breaks the weight symmetry of ties.

The appropriateness of this method for undirected networks depends on the research setting and how tie weights are defined (I would like to thank Tom Snijders for highlighting this during my viva). For example, its applicability to undirected transportation networks is justified by the typically directed nature of traffic flows (although the US airport network displays a high symmetry; Barrat et al., 2004). Conversely, in an undirected collaboration network this might not be appropriate. More generally, for one-mode projections of two-mode networks, it might be more appropriate to reshuffle the two-mode network before projecting it onto a one-mode network (see randomisation without structural zeros in Rao et al., 1996; Snijders, 1991).

To test this method, I have applied it to the scientific collaboration network used in the Weighted Rich-club paper (Newman, 2001). This network is available on the Datasets-page (Thanks, Mark for sending me a copy and allowing me to post it!).

To calculate \phi_{\mathrm{null}}^w(s) using this null model, I first randomised the two-mode network by reshuffling ties while maintaining authors’ and papers’ degree (similarly to the method proposed by Molloy and Reed, 1995, for one-mode networks). One step of this method is exemplified in the diagram below. I then projected the randomised two-mode network onto a weighted one-mode network using the same method in the paper (for more details, see projecting two-mode networks). Finally, I calculated the weighted-rich club effect, \phi^w, on the one-mode projections. The \phi_{\mathrm{null}}^w(s) is the average \phi^w over many random networks.
Two-mode reshuffling

The result from this procedure is the displayed in the following diagram (the bars refer to the 95% confidence interval, see Section 3.2.2 of my thesis for more details):

Weighted rich-club coefficient for Newman's (2001) scientific collaboration network.

The diagram suggest that, unlike what we found using the Directed weight reshuffle, when prominence is defined as having co-authored more than 20 articles, there is a negative and significant weighted rich-club effect. Thus, this finding is a substantiation of the speculation that in certain social networks a negative effect exists among the very prominent people due to some form of competition. This might account for the reluctance of highly productive authors to establish strong ties among themselves, as is suggested by the lack of interaction among the two most productive authors in the Figure below: Barabasi and Newman are not connected.¹

Subset of the prominent nodes (s >= 5) in the network science collaboration network (Newman, 2006).

The null model used in this post could also be used in a number of other frameworks. For example, it could be used in the topological rich-club (Colizza, 2006; Zlatic et al., 2008). In particular, it might be more appropriate for Zlatic et al. (2008) as they defined prominence as node strength.
¹ Although Albert-László Barabási and Mark Newman have written a book together with Duncan Watts, it did not form part of this dataset. Nevertheless, this collaboration would add only a weak tie between these authors.

Want to test it with your data?

The following code was used to produced the above diagrams:

# Load tnet

# Download data

# Run function 


Barrat, A., Barthelemy, M., Pastor-Satorras, R., Vespignani, A., 2004. The architecture of complex weighted networks. Proceedings of the National Academy of Sciences 101 (11), 3747-3752. arXiv:cond-mat/0311416

Colizza, V., Flammini, A., Serrano, M. A., Vespignani, A., 2006. Detecting rich-club ordering in complex networks. Nature Physics 2, 110-115. arXiv:physics/0602134

Molloy, M., Reed, B., 1995. A critical point for random graphs with a given degree sequence. Random Structures and Algorithms 6, 161-180.

Newman, M. E. J., 2001. Scientific collaboration networks. II. Shortest paths, weighted networks, and centrality. Physical Review E 64, 016132.

Newman, M. E. J., 2006. Finding community structure in networks using the eigenvectors of matrices. Physical Review E 76 (036104). arXiv:physics/0605087

Opsahl, T., 2013. Triadic closure in two-mode networks: Redefining the global and local clustering coefficients. Social Networks 35, doi:10.1016/j.socnet.2011.07.001.

Opsahl, T., Colizza, V., Panzarasa, P., Ramasco, J. J., 2008. Prominence and control: The weighted rich-club effect. Physical Review Letters 101 (168702). arXiv:0804.0417.

Rao, A. R., Jana, R., Bandyopadhyay, S., 1996. A markov chain monte carlo method for generating random (0, 1)-matrices with given marginals. Sankhya A 58, 225- 242.

Serrano, M. A., Boguna, M., Vespignani, A., 2007. Patterns of dominant flows in the world trade web. Journal of Economic Interaction and Coordination 2, 111-124. arXiv:0704.1225.

Snijders, T. A. B., 1991. Enumeration and simulation methods for 0-1 matrices with given marginals. Psychometrika 56 (3), 397-417.

Zlatic, V., Bianconi, G., Diaz-Guilera, A., Garlaschelli, D., Rao, F., Caldarelli, G., 2008. On the rich-club effect in dense and weighted networks. arXiv:0807.0793.

If you use any of the information on this page, please cite: Opsahl, T., Colizza, V., Panzarasa, P., Ramasco, J. J., 2008. Prominence and control: The weighted rich-club effect. Physical Review Letters 101 (168702)

5 Comments Add your own

  • 1. Nico  |  December 12, 2017 at 12:57 pm

    Hey Tore,

    great blog, it’s been very interesting and usefull to me.

    Unfortunately I ran into some trouble replicating this entry. I’m not sure, if entirely understand the approach to get the diagram with the normalised phi coefficient. So far I tried the following:

    net <- as.tnet(edgelist, type="weighted two-mode tnet")
    random_net <- rg_reshuffling_tm(net, option = "links", seed = 1)
    random_one <- projecting_tm(random_net , method="sum")
    out <- weighted_richclub_w(random_one, rich="k", reshuffle="links", NR=1000, seed=1)

    Is this correct or should I rather reshuffle the two-mode-network 1000 times instead of the one mode rich club?

    I hope you find the time to respond!

    • 2. Tore Opsahl  |  December 20, 2017 at 5:13 pm

      Hi Nico,

      Glad you are interested in using the two-mode version of the weighted rich-club coefficient. I am not sure where that snippet of code is from; however, you should use the weighted_richclub_tm-function that computes the coefficient on an empirical network and compares it to ones found on an ensemble of reshuffled versions of the empirical network.


      • 3. Nico  |  January 3, 2018 at 8:40 am

        Thanks for your answer! I’ll try that :-)

  • 4. Nico  |  January 16, 2018 at 10:15 am

    Sorry to bother you again, but I was wondering, if there is a way to make it work with a weighted bipartite network?

    Because now I get the error message: “There are duplicated entries in the edgelist” and since `weighted_richclub_tm` takes binary TM networks only, I cannot enter a weighted TM net (right?).

    Hopefully you can help me out once more!

    • 5. Tore Opsahl  |  January 17, 2018 at 3:48 am

      Hi Nico,

      I haven’t written the code for a weighted bipartite network; however, there is nothing that would stop you from defining a rich club coefficient that incorporate weights, and then use the various reshuffling procedures for weighted bipartite networks.

      Let me know how you get along!


Leave a Reply to Nico Cancel reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s

Subscribe to the comments via RSS Feed

%d bloggers like this: