Tomography of scaling.

Abstract

Scaling describes how a given quantity <i>Y</i> that characterizes a system varies with its size <i>P</i>. For most complex systems, it is of the form <math xmlns="http://www.w3.org/1998/Math/MathML"><mi>Y</mi><mo>∼</mo><msup><mi>P</mi><mrow><mi>β</mi></mrow></msup></math> with a non-trivial value of the exponent <i>β</i>, usually determined by regression methods. The presence of noise can make it difficult to conclude about the existence of a nonlinear behaviour with <i>β</i> ≠ 1 and we propose here to circumvent fitting problems by investigating how two different systems of sizes <i>P</i>₁ and <i>P</i>₂ are related to each other. This leads us to define a local scaling exponent <i>β</i>_loc that we study versus the ratio <i>P</i>₂/<i>P</i>₁ and provides some sort of 'tomography scan' of scaling across different values of the size ratio, allowing us to assess the relevance of nonlinearity in the system and to identify an effective exponent that minimizes the error for predicting the value of <i>Y</i>. We illustrate this method on various real-world datasets for cities and show that our method reinforces in some cases the standard analysis, but is also able to provide new insights in inconclusive cases and to detect problems in the scaling form such as the absence of a single scaling exponent or the presence of threshold effects.