Degree-preserving randomization explained

Background

Cataloged as early as 1996,^[1] the simplest implementation of degree preserving randomization relies on a Monte Carlo algorithm that rearranges, or "rewires" the network at random such that, with a sufficient number of rewires, the network's degree distribution is identical to the initial degree distribution of the network, though the topological structure of the network has become completely distinct from the original network.

Degree preserving randomization, while it has many different forms, typically takes on the form of a relatively simple approach: for any network consisting of

nodes with

edges, select two dyadically tied nodes. For each of these dyadic pairs, switch the edges such that the new dyadic pairs are mismatched. After a sufficient number of these mismatches, the network increasingly loses its original observed topography.

As is common with algorithms based on Markov chains, the number of iterations, or individual rewires, that must occur on a given graph such that the graph is sufficiently random and distinct from the original graph is unknown, though Espinoza^[2] asserts that a safe minimum threshold is

Q*E

, where

"is at least 100" (Espinoza). Others have provided input for this issue, including one author who states that a safe minimum may instead be at least

	E
	2

*ln(1/\epsilon)

, where epsilon lies in a range between

10^-6

and

10^-7

, though ultimately the correct number is not presently known.

Uses

There are several cases in which published research have explicitly employed degree preserving randomization in order to analyze network properties. Dekker used rewiring in order to more accurately model observed social networks by adding a secondary variable,

\pi

, which introduces a high-degree attachment bias. Liu et al. have additionally employed degree preserving randomization to assert that the Control Centrality, a metric they identify, alters little when compared to the Control Centrality of an Erdős–Rényi model containing the same number of

nodes in their simulations - Liu et al. have also used degree preserving randomization models in subsequent work exploring network controllability.

Additionally, some work has been done in investigating how Degree Preserving Randomization may be used in addressing considerations of anonymity in networked data research, which has been shown to be a cause for concern in Social Network Analysis, as in the case of a study by Lewis et al. Ultimately the work conducted by Ying and Wu, starting from a foundation of Degree Preserving Randomization, and then forwarding several modifications, has showed moderate advances in protecting anonymity without compromising the integrity of the underlying utility of the observed network.

Additionally, the method is similar in nature to the broadly used Exponential random graph models popularized in social science, and indeed the various forms of modeling networks against observed networks in order to identify and theorize about the differences expressed in real networks. Importantly, Degree Preserving Randomization provides a simple algorithmic design for those familiar with programming to apply a model to an available observed network.

Example

What follows is a small example showing how one may apply Degree Preserving Randomization to an observed network in an effort to understand the network against otherwise random variation while maintaining the degree distributional aspect of the network. The Association of Internet Researchers has a Listserv that constitutes the majority of discussion threads surrounding their work. On it, members post updates about their own research, upcoming conferences, calls for papers and also engage one another in substantive discussions in their field. These emails can in turn constitute a directed and temporal network graph, where nodes are individual e-mail accounts belonging to the Listserv and edges are cases in which one e-mail address responds to another e-mail address on the Listserv.

In this observed network, the properties of the Listserv are relatively simple to calculate - for the network of 3,235 individual e-mail accounts and 9,824 exchanges in total, the observed reciprocity of the network is about 0.074, and the [Average path length|average path length] is about 4.46. Could these values be arrived at simply through the nature of the network's inherent structure?

Applying the

	E
	2

*ln(1/\epsilon)

rule, this network would require around 67,861 individual edge rewires to construct a likely sufficiently random degree-preserved graph. If we construct many random, degree preserving graphs from the real graph, we can then create a probability space for characteristics, such as reciprocity and average path length, and assess the degree to which the network could have expressed these characteristics at random. 534 networks were generated using Degree Preserving Randomization. As both reciprocity and average path length in this graph are normally distributed, and as the standard deviation for both reciprocity and average path length are far too narrow to include the observed case, we can reasonably posit that this network is expressing characteristics that are non-random (and thus open for further theory and modeling).

External links

Dataset for example provided

Notes and References

Rao . A Ramachandra . Jana . Rabindranath . Bandyopadhyay . Suraj . 1996 . A Markov chain Monte Carlo method for generating random (0, 1)-matrices with given marginals . Indian Journal of Statistics Series A . November 5, 2014.
Espinoza . Max . On Network Randomization Methods: A Negative Control Study . 2014-11-06 . 2016-03-04 . https://web.archive.org/web/20160304114446/http://biogrid.engr.uconn.edu/REU/Reports_12/Final_Reports/Max.pdf . dead .