March 6, 2018

Download this article as a PDF file

In the past weeks, I have received several requests to
address the merits of the Anna D. Broido and Aaron Clauset
(BC) preprint [1] and their fruitless search for scale-free networks in nature. The preprint’s central claim is deceptively simple: It
starts from the textbook definition of a scale-free network as a network with a
power law in the degree distribution [2]. It then proceeds to fit a power law
to 927 networks, finding that only 4% are scale-free. The author's conclusion that ‘scale-free networks are rare,’ is turned into the title of the preprint, helping it to get maximal attention.
It worked—*Quanta* magazine accepted its
conclusions without reservations. After*The*
*Atlantic* carried the article, the
un-refereed preprint received a degree of media exposure that the original
discovery of scale-free networks never enjoyed.

While I saw the conceptual problems with the manuscript, I was convinced that the paper must be technically proficient. Yet, once I did dig into it, it was a real ride. If you have the patience to get to the end of this commentary, you will see where it fails at the conceptual level. But, we will learn that it also fails, repeatedly, at the technical level.

Let me start by summarizing an enormous body of work on scale-free networks as briefly as possible, essential to understand the moment when the BC paper arrives on the scene.

*Empirical
Discovery*. The scale-free property was observed in 1999,
independently in the WWW [3] and the
internet at the autonomous systems level [4]. The
empirical observation was simple: the degree distribution of these networks is well approximated with a power law,

P(k)∼k⁻ᵞ (1)

Back then this was an unexpected finding, given that the
prevailing model was the random network model of Erdős and Rényi [5] which predicted a Poisson degree
distribution, easily differentiable from a power law (Figure 1). We named these networks *scale-free*, inspired by the scale-free
nature of power laws observed in the vicinity of phase transitions.

*Mechanistic
Model*. In 1999 Réka Albert and I proposed a mechanism to
explain the origin of the observed power law [2],
finding that it is rooted in growth and preferential attachment. A simple model based on these two mechanisms, known today as the BA or the *scale-free
model*, generated a network whose degree distribution follows (1), with
degree exponent *γ**=3.
* The scale-free model is a mechanistic model,
meaning that it is not a model of the Internet, or the WWW or the cell—it only
aims to explain the mechanism behind the scale-free nature of a network. Real
networks are far more complicated, indicated by the fact that their degree exponent *γ* varies from 2 to 8. Growth and
preferential attachment *alone* cannot
explain that variation.

*Real
Networks*. Months after the 1999 paper several key discoveries
helped the community understand the origins of the different exponents. The
rate-equation based continuum theory, developed by Mendes,
Dorogovtsev, Redner, Krapivsky, and many others [6], helped incorporate into the scale-free model many effects that naturally occur in real
networks, like the disappearance of nodes, the addition of new links between
previously existing nodes, link deletion, aging, and so on. These papers have
shown analytically that the presence of the additional effects alter the degree distribution in two
ways. First, they change the degree exponent— successfully explaining the diverse exponents
observed in real networks. Equally important, *they induce corrections to the degree distribution, making P(**k*)* deviate in predictable ways from a pure power law*. Some of the
most common corrections include [7]:

*Small degree saturation:*In any system where two nodes that are already in the network can choose to connect to each other (internal links), the analytical form of the degree distribution becomes**P(k)=C(k+A)***⁻ᵞ,*leading to a small degree saturation. In social networks and the WWW, most links are such internal links. A similar small-degree saturation can also be induced by initial attractiveness, well documented in citation networks [8].*High degree cutoffs:*If preferential attachment is sublinear,. Similar high degree cutoffs also appear under node removal.**the degree distribution follows a stretched exponential, or a power law with an exponential cutoff***Fitness:*Nodes have different abilities to complete for the links, a diversity that can be modelled by giving each of them a unique and different fitness*η*. In the presence of fitness, the precise form of P(k) depends on the fitness distribution*ρ**(**η**)*. For example, a uniform fitness distribution induces a logarithmic correction in P(k) and other forms of*ρ**(**η**)*can lead to rather exotic forms for P(k).

In real networks many of the elementary
processes discussed above appear together, and so do many others. In other words, by 2001 it was pretty
clear *that there is no one-size-fits all
formula for the degree distribution for networks driven by the scale-free mechanism*. A *pure power law only emerges in simple
idealized models, driven by only growth and preferential attachment*, and
free of any additional effects. The theory was predicting that in real networks, if the scale-free is present, the power law tends to coexist with some corrective
function, expecting power laws with exponential cutoffs, stretched
exponentials, logarithmic corrections to the power law, and so on. These are
so important, that I devoted a full chapter in the *Network Science [9]*
textbook to them. We also learned that there are multiple ways of analyzing the presence of the scale-free property, as I explain in Box 1. The conclusion is simple, well understood
in the literature*: if we wish to obtain
an accurate fit to the degree distribution of a real network, we first must build
a generative model that analytically predicts the functional form of P(k)**.*

So it is time to go back to Ref [1] and
its key claim that “*Scale-free networks
are rare*.” How exactly did it arrive to this conclusion? By insisting to
fit a pure power law to every network, and ignoring what the theory predicts
for any of them. As it is difficult to find real systems that are free of
additional effects, *it makes no sense to
fit indiscriminately a power law to all of them. One must fit the distribution
that the theory predicts, which is predictably different for each system.*

Interestingly, the theory predicts that in many real networks driven by growth and preferential attachment, the degree distribution should follow a power law with an exponential cutoff. If you look at Table II of Ref [1], BC find that 51% of the networks they explored favor this distribution. In other words, the measurements of BC validate the theory, contradicting the authors' central claim.

Getting this far, you may ask yourself, if it is so difficult to fit the degree distribution, why do thousands of papers claim that a large number of real networks, from the Internet to protein interaction and social networks, are scale-free? The answer is simple—despite the many processes shaping their topology, for many real networks the fat tailed nature is so obvious that it’s hard to miss (see Figure 1). Indeed, Clauset’s most cited work found that many of the large networks that the community does consider scale-free, like the Internet, protein interaction network, citation network, co-authorship networks, do pass statistical significance without even considering the necessary corrections (see Table 6.1 in Ref [12]).

So the puzzling technical question is this: Why do the authors fail to find the scale-free networks, that everyone
else in the literature does, including Clauset himself in his earlier work?
This is where the first technical surprise comes: ** they fail to see the scale-free
property because they invent a new criterion of weak and strong scale-free
networks.** And the real surprise?

You would think that the nomenclature of 'weak' and 'strong' scale-free networks BC uses
has to do with statistical significance. Indeed, one could plausibly call some
of the large and well mapped real networks ‘strong scale-free’, implying that for
them the statistical significance is exceptional. And one could coin the term “weak
scale-free network” for a network for which a power law has lower statistical
significance. But that is not what BC
do. Instead, they take each network, and generate from the original *multiple* *synthetic networks*. This way their set
of 927 real networks turns into *23,999
synthetic networks*. Then they toss out 81% of them, deeming them ‘unlikely to be scale-free’ (pg 3, BC). They proceed with the
remaining 4,477 synthetic networks, which is the corpus they study.

For example, an unweighted directed network like the
WWW turns into three different degree sequences: the incoming, the
outgoing, and the total degree. Instead of asking which of these synthetic networks may be well fitted
with a power law, they ask, what *fraction
of these synthetic networks passes* the power law criteria.
Then they propose the following naming convention, taken from their paper:

*Super-Weak*: For at least 50% of graphs, none of the alternative distributions are favored over the power law.*Weakest*: For at least 50% of graphs, the power- law hypothesis cannot be rejected (p ≥ 0.1).*Weak*: The requirements of the Weakest set, and there are at least 50 nodes in the distribution’s tail (ntail > 50).*Strong*: The requirements of the Weak set, and that both 2 < αˆ < 3, and for at least 50% of graphs none of the alternative distributions are favored over the power-law.*Strongest*: The requirements of the Strong set for at least 90% of graphs, rather than 50%, and for at least 95% of graphs none of the alternative distributions are favored over the power-law.

To help understand what this is, let us take the citation graph, a well studied scale-free network. It is a directed network, whose incoming degrees (the number of citations a paper gets) is known since the 1960s to be well fitted with a power law [13]. However, the outgoing degrees follow a different distribution, as those degrees capture how many different papers a given publication cites, and that number is determined by journal policy, having artificial cutoffs. So, of the three graph they generate from a well-studied scale-free network, only one can pass their criteria, the one defined by the incoming degrees. But that is not enough to make it even super-weak scale-free in the language of BC.

In some cases, their methodology generates
as many as 900 synthetic networks from a single real system. If at least
50% of these synthetic networks pass the power law test they would call it super-weak. They require 90% to pass to place a network in the strongest category. Now, the truth is that typically ** only one of these networks matters: the one
that captures the purpose or the function of the original system.** But
they offer an equal vote to each fictional synthetic network, letting them decide whether the original system is scale-free or not (Box 3).

The classification BC choose to use has no precedent in network science, and has no relevance to the role of the network they study. They offer no explanation of the need for these made-up criteria, nor do they explan how they choose the 50% and 90% percentages. They are arbitrary.

But let's focus on what really
matters. And those are controls. That is, whatever criteria we choose, we
must test it on networks whose degree distribution is well understood. So if we feed into our methodology a scale-free network, like the one produced by the
scale-free model with a pure power law degree distribution, free of any corrections,
the method should predict that that network is in the strongest category.
Indeed, *if the gold standard fails the
strongest class, who would pass*? In
other words, the scale-free model should not be super-weak, weakest, or weak. It should be in the strongest class, correct? After all, we have a formal exact proof of the power law
nature of the scale-free model [14].

BC are truly aware of this crucial criteria, hence they do assure the reader on pg 5 that their method passes this obvious but essential test:

“Two mechanisms produce scale-free networks (a directed version of preferential attachment [42] and a directed vertex copy model [43]), and one does not (simple Erdős-Rényi random graphs). Applied to synthetic networks generated by these mechanisms, our methodology correctly placed the synthetic data sets into scale-free categories suitable for their generating parameters, i.e., preferential attachment and vertex-copy data sets were categorized as scale-free networks, while Erdős-Rényi random graphs were not (see Appendix E).”

So it’s all good. Yet, curious, I did look into Appendix E, which, I assumed that would simply offer the evidence. And it did. But not the evidence we would all expect based on the text above. Here is word by word what we find there:

“When we consider the plausibility of the power-law fit, we see fewer networks. 62% of the preferential attachment graphs fall into the Weakest and Weak categories,”

They also report for the Erdős-Rényi network, that

“51% and 50% of these networks fall into the Weakest and Weak categories, respectively."

In other words, according to the criteria
the authors invented, *38% of the time the scale-free model is not even ‘weak
scale-free.' In contrast, 51% of the time the ER model is classified as ‘weak
scale-free’!*

You may ask, how hard is to distinguish an Erdős-Rényi model from a scale-free one? In their test BC used networks of 5,000 nodes. In Figure 3 we can see how different are the degree distribution of a scale-free network and an Erdős-Rényi network at this size. It needs no sophisticated statistical test to capture the difference.

But wait, it gets even more interesting. If we read further (pg 7, BC), this is what we find about their gold standard, the scale-free model (they refer to it as ‘preferential attachment graph’):

“When we consider the plausibility of the power-law fit, we see fewer networks. 62% of the preferential attachment graphs fall into the Weakest and Weak categories, 60% in the Strong category,and 0 in the Strongest category. “

*That is, not a single realization
of the scale-free model is deemed strongly scale-free**.* **According to the criteria they invented, not even the scale-free model is scale-free any longer.**

At the end the preprint's main result, so prominently featured in the press, rests on a simple claim: the authors fail to find the numerous scale-free networks in nature. If we actually read the paper, we find that they refuse to see them: Table II in Ref [1] shows that a high fraction of networks have degree distributions in agreement with the predictions of the continuum theory describing scale-free networks. The true failure is their methodology: it fails to detect that the gold standard is scale-free. Now, if you invent an arbitrary criterion that not even the mathematically proven models can pass, what exactly do you expect to get?

What really puzzles me, after all this, is that 4% of real networks did pass their criteria, finding them to be strongly scale-free. Which are those networks that are more scale-free than the scale-free model? They must be walking on water.

It is crucial for network science to constantly test the solidity of its pillars, like the scale-free nature of real networks, and the mechanisms responsible for it. We must continue therefore to constantly rule out competing hypotheses, a process that is essential for science to advance. Cutting corners will not get us there, however.

So before accepting the attention-grabbing claims of
the BC paper, you must peek under its hood. And if you take the time to do
that, you will find that the study ** is oblivious to 18 years of knowledge accumulated in network science**. You will
also find a

[1] AD Broido, A Clauset. Scale-free networks are rare. arXiv:1801.03400

[2] A.-L. Barabási and R. Albert. Emergence of scaling in random networks. Science, 286:509-512, 1999.

[3] R. Kumar, P. Raghavan, S. Rajalopagan, and A.Tomkins. Extracting Large-Scale Knowledge Bases from the Web. Proceedings of the 25thVLDBConference, Edinburgh,Scotland,pp.639-650,1999; H. Jeong, R. Albert and A. L. Barabási. Internet: Diameter of the world-wide web. Nature, 401:130-131, 1999;

[4] M. Faloutsos, P. Faloutsos, and C. Faloutsos. On power-law relationships of the internet topology. Proceedings of SIGCOMM. Comput. Commun. Rev. 29: 251-262, 1999.

[5] P. Erdős and A. Rényi. On random graphs, I. Publicationes Mathematicae (Debrecen), 6:290-297, 1959.

[6] P. L. Krapivsky, S. Redner, and F. Leyvraz, Connectivity of Growing Random Networks”, Phys. Rev. Lett. 85, 4629-4632 (2000); S.N. Dorogovtsev, J.F.F. Mendes, and A.N. Samukhin. Structure of growing networks with preferential linking, Physical Review Letters, 85: 4633, 2000.

[7] For more details, see Chapter 6 of Network Science.

[8] Y.-H. Eom and S. Fortunato. Characterizing and Modeling Citation Dynamics. PLoS ONE, 6: e24926, 2011.

[9] A.-L. Barabási *Network
Science*. Cambridge University Press, 2017.

[10] Section 5.7 in Ref. [9].

[11] R. Cohen, K. Erez, D. ben-Avraham and S. Havlin. Resilience of the Internet to random breakdowns. Physical Review Letters, 85:4626, 2000; B. Bollobás and O. Riordan. Robustness and Vulnerability of ScaleFree Random Graphs. Internet Mathematics, 1, 2003; R. Pastor-Satorras and A.Vespignani. Epidemic spreading in scalefree networks. Physical Review Letters, 86:3200–3203, 2001.; R. Albert, H. Jeong, A.-L. Barabási. Error and attack tolerance of complex networks. Nature, 406: 378–482, 2000;

[12] A Clauset, CR Shalizi, MEJ Newman. Power-law distributions in empirical data. SIAM review 51 (4), 661-703, 5933 (2009).

[13] In some cases, depending how the data was collected, the log-normal also offers a reasonable fit to citation networks. For the power law, we have generative methods; the mechanistic origin of the log normal are less understood.

[14] B. Bollobás, O. Riordan, J. Spencer, and G. Tusnády. The degree sequence of a scale-free random graph process. Random Structures and Algorithms, 18:279-290, 2001.

[15] R. Cohen, K. Erez, D. ben-Avraham and S. Havlin. Resilience of the Internet to random breakdowns. Physical Review Letters, 85:4626, 2000; B. Bollobás and O. Riordan. Robustness and Vulnerability of ScaleFree Random Graphs. Internet Mathematics, 1, 2003; R. Pastor-Satorras and A.Vespignani. Epidemic spreading in scalefree networks. Physical Review Letters, 86:3200–3203, 2001.; R. Albert, H. Jeong, A.-L. Barabási. Error and attack tolerance of complex networks. Nature, 406: 378–482, 2000;

...

Recent posts

A study's failure to find scale-free networks where decades of research has documented their existence offers a cautionary tale on using search criteria that fails elementary tests.

continue readingA physicist probes an elegant law that governs cells, cities, and almost everything in between

continue reading