Scale-Adjusted Metrics for Predicting the Evolution of Urban Indicators and Quantifying the Performance of Cities

More than a half of world population is now living in cities and this number is expected to be two-thirds by 2050. Fostered by the relevancy of a scientific characterization of cities and for the availability of an unprecedented amount of data, academics have recently immersed in this topic and one of the most striking and universal finding was the discovery of robust allometric scaling laws between several urban indicators and the population size. Despite that, most governmental reports and several academic works still ignore these nonlinearities by often analyzing the raw or the per capita value of urban indicators, a practice that actually makes the urban metrics biased towards small or large cities depending on whether we have super or sublinear allometries. By following the ideas of Bettencourt et al. [PLoS ONE 5 (2010) e13541], we account for this bias by evaluating the difference between the actual value of an urban indicator and the value expected by the allometry with the population size. We show that this scale-adjusted metric provides a more appropriate/informative summary of the evolution of urban indicators and reveals patterns that do not appear in the evolution of per capita values of indicators obtained from Brazilian cities. We also show that these scale-adjusted metrics are strongly correlated with their past values by a linear correspondence and that they also display crosscorrelations among themselves. Simple linear models account for 31%–97% of the observed variance in data and correctly reproduce the average of the scale-adjusted metric when grouping the cities in above and below the allometric laws. We further employ these models to forecast future values of urban indicators and, by visualizing the predicted changes, we verify the emergence of spatial clusters characterized by regions of the Brazilian territory where we expect an increase or a decrease in the values of urban indicators.

Introduction the world urban population has surpassed the rural one and, if this process persists, two-thirds of the world population are expected to be living in urban areas by 2050. On one hand, cities are usually associated with higher levels of literacy, health care and better opportunities; on the other hand, unplanned urbanization and bad political decisions also lead to pollution, environmental degradation, growth of crime, unequal opportunities and the increase in the number of people living in substandard conditions. In this sense, there is a vast necessity for finding patterns, quantifying and predicting the evolution of urban indicators, since these investigations may provide guidance for better political decisions and resources allocation.
Fostered by this need and also due the availability of an unprecedented amount of data at city level, several researchers have recently promoted an impressive progress into what has been called Science of Cities [2]. These new data allowed researchers to probe patterns of the cities to a degree not before possible and one of the most striking and universal finding of these studies was the discovery of robust allometric scaling laws between several urban indicators and the population size. Patents, gasoline stations, gross domestic product [3][4][5][6], crime [6][7][8][9][10][11], indicators of education [11], suicides [12], number of election candidates [13,14], transportation networks [15,16], employees from several sectors [17], measures of social interaction [18] are just a few examples where robust scaling laws have been found. Despite these intrinsic nonlinear relationships it is common to find works that try unveiling relationships between population size and urban indicators by employing linear regressions in the raw data; also, per capita indicators (that is, divided by population size) are ubiquitous in reports of government agencies and are often used as a guide for public policies when analyzing the temporal evolution of a given city or for comparing/ranking the performance of cities with different population sizes. However, these linear regressions may result in misguided/controversial relationships [9,19] and per capita indicators are oblivious to the allometric scaling laws that make city a complex agglomeration that cannot be modeled as a linear combination of its individual components [6,20,21]. In this sense, there is a paucity of alternative metrics for urban indicators that may overcome the foregoing problems and provide a fairer comparison between cities of different sizes as well as a better understanding of the city evolution. Bettencourt et al. [6] have recently proposed a simple alternative to overcome these nonlinearities by evaluating the difference between the actual value of the urban indicators and the value expected by the allometries with the population size (that is, the residuals in the allometric relationships). This scale-adjusted metric explicitly considers the allometric relationships and already have proved to be useful in the economic context [22,23] and for unveiling relationships between crime and urban metrics that are not properly carried out by regression analysis [9].
In this article, we follow the evolution of this relative metric for eight urban indicators from Brazilian cities in three years (1991, 2000 and 2010) in which the national census took place. By grouping cities in above and below the allometric laws, we argue that the average of this scale-adjusted metric provides a more appropriate/informative summary of the evolution of the urban indicators when compared with the per capita values; it also reveals patterns that do not appear when analyzing only the evolution of the per capita values. For instance, while the per capita values of homicides have systematically increased over the last three decades, both the average of the scale-adjusted metric for cities above and below the allometric law are approaching zero, that is, cities where the number of homicides is above the expected by the allometry have managed to reduce this crime, whereas this crime has been increased in those cities where number of homicides is below the allometry. We argue that the nonlinearities may affect the per capita indicators by creating a bias towards large cities for superlinear allometries and a bias towards small cities for sublinear allometries. We further show that these scale-adjusted metrics are strongly correlated with their past values by a linear correspondence, making them particularly good for predicting future values of urban indicators through linear regressions. We have tested this hypothesis via a linear model where the scale-adjusted metric for one indicator in a given census was predicted by a linear combination of all eight metrics evaluated from the preceding census. These simple models account for 31%-97% of the observed variance in data and correctly reproduce the average value of the scaled-adjusted metric when grouping the cities in above and below the allometric laws. Motivated by these good agreements, we present a prediction for the values of urban indicators in the year of 2020 by assuming the linear coefficients constants over time. By visualizing the predicted changes, we verify the emergence of spatial clusters characterized by regions of the Brazilian territory where the model predicts an increase or a decrease in the values of urban indicators. We further report a list containing all the scale-adjusted metrics as well as the predictions for each city in the hope that government agencies find these informations useful (S1 Dataset).

Results and Discussion
We start by considering the average of the per capita values for eight urban indicators described in the Methods Section. This is a common practice of government agencies for tracking the evolution of a particular city or for comparing a group of cities with different populations. We observe in Fig. 1 that almost all per capita indicators show a clear temporal trend: elderly population, female population, homicides and family income have increased over the years; whereas child labor, illiteracy and male population have decreased (the unemployment rates have evolved in a more complex manner, exhibiting no clear tendency). We could also list the cities in which these indicators have considerably changed or rank the cities that have made more progress in reducing, for instance, the homicide or illiteracy rates.
One of the main problems with this analysis is that it completely ignores the hypothesis that most urban indicators Y i (t) displays allometric relationships (or scale invariance) with the population N (t), that is, where A i is a constant, β i is the allometric (or scaling) exponent and t stands for time. This simple relationship summarizes the (average) effects of increasing the population size on the urban indicators; it states that cities are self-similar in terms of their population, in the sense that average properties of a given city can be inferred only by knowledge of its population. Urban indicators are thus expected to display a deterministic component emerging from very few and general properties of the urban networks related to social and infrastructure aspects of cities [21]. When analyzing per capita values, it is implicitly assumed that the value of an urban indicator is proportional to the population size (β i = 1) or, in other words, that cities are extensive systems. This idea is opposed to the complex systems approach of cities: complex systems are non-extensive (β i ≠ 1), meaning that its isolated parts do not behave in the same manner as when they are interacting. Cities have similar properties and only make sense as an entire "organism" and, in fact, there is robust evidence favoring the non-extensive and universal nature of cities (that is, the urban scaling hypothesis) across different cultures and historical periods [21,24,25]. Thus, several properties of a city of a given size cannot be linearly scaled for another city with larger or smaller population size. The dynamical processes mediated by the urban networks make the scale operation a nonlinear transformation for several urban indicators, often resulting in per capita savings of material infrastructure and in gain of socio-economic  The tiny error bars are 95% bootstrap confidence intervals for the average values. We note that elderly population, female population, homicides and family income display an increasing tendency; whereas child labor, illiteracy and male population have decreased over time. Despite having increased, the unemployment rates show a not very clear trend.
productivity [3]. From a more technical point of view, whenever the allometric exponent β i is different from one, there is a remaining component related to the population size when evaluating the per capita values of these urban indicators, that is, y i = Y i N ∼ N (βi−1) , which creates a bias towards large cities for β i > 1 and towards small cities for β i < 1; per capita measures are only efficient in correctly removing the effect of the population size in an urban indicator for β i = 1. For our data, a complete description of the allometric relationships between the eight urban indicators and the population is presented in the Methods Section, where we have confirmed the presence of allometries in our data as summarized in Table 1 (see also S1 and S2 Figs.).
The problem with per capita indicators thus prompts the question on how we can account for these allometries and correctly remove the effect of the population size on the urban indicators. Bettencourt et al. [6] have proposed a simple and efficient procedure for overcoming this problem by defining the so-called scale-adjusted urban indicator. The approach consists in evaluating the logarithmic difference between the actual value of an urban indicator Y i (t) and the value expected by allometric relationship with the population N (t) (that is, the residuals in the allometric relationships) in given year t, namely, The previous quantity explicitly considers the allometry between an urban indicator and the population size, creating a relative measure that is not biased by the population size (size-independent) for any value of β i . The scale-adjusted metric D Yi (t) captures the exceptionality (either good or bad) of a city, which somehow is the result of the nonlinear agglomeration process related to the socio-economic choices and historical path of a city [6]. Furthermore, D Yi establish a more "natural" scale for ranking cities by identifying whether an urban indicator of a given city is above (D Yi > 0) or below (D Yi < 0) the expected value from cities of similar sizes. This approach already have proved to be useful in the economic context [22,23] and was also employed for unveiling relationships between crime and urban metrics that are not properly carried out by linear regression analysis [9]. Figure 2 illustrates the definition of D Yi (t) and shows an example of allometry between homicides and population size (see also S1 and S2 Figs. for all urban indicators). In order to show the informations that this scale-adjusted metric provides, we study the evolution of the average D Yi (t) by creating two groups of cities: those whose the urban indicator Y i (t) was above (D Yi (t) > 0) and those whose Y i (t) was below (D Yi (t) < 0) the allometric law in the year t = 1991. Figure 3 shows these averages for the eight urban indicators over the three years in our database. For the majority of the urban indicators, the average D Yi (t) displays a statistically significant decreasing tendency for cities initially above the allometric law; whereas an increasing tendency for the average D Yi (t) is observed for those cities that were initially below the allometric law. For instance, cities where the number of children laboring, homicides and unemployment were above the value expected by the allometric law have been successful (on average) in reducing them; in contrast, cities where these indicators were below the allometric laws proved unable (on average) to improve or maintain this situation. The case of illiteracy is an exception to the previous pattern, since the average D Yi (t) has increased for cities initially above the allometric law and it is almost a constant for those cities initially below. Thus, cities where illiteracy was initially above the allometric law have failed (on average) in increasing the number of literates; on the other hand, those cities initially below the allometric law have not only kept this feature (on average) but also managed to further improve the literacy levels. In the case of population metrics, particularly regarding female and male populations, the approaching to the allometric laws (together with the decrease in proportion of males in the population) may represent a good aspect with respect to reduction of violence, since an excessive contingent of men can drive an increase in antisocial behavior [26]; moreover, both male population above and female population below the allometric law correlate with number of homicides above the allometric law [9]. Similarly, family income above the allometric law correlates with homicides above, while family income below the allometric law correlates with homicides below [9]. Remarkably, these informations remain hidden or distorted when we look only at the per capita values of these urban indicators.
The trends in the average values of the scale-adjusted metrics prompt the question of how the values of this relative metric are affected by their previous values, that is, are the values of D Yi (t + ∆t) and D Yi (t) correlated in some particular fashion? To answer this question, we analyze the scatter plots of the scale-adjusted metric in a given year versus its past values for each urban indicator for all cities.  (1991). We observe that, despite the different scattering degrees (Pearson correlation ranging from 0.47 to 0.98), linear functions are good approximations for the average tendency of these relationships. We thus adjust the linear model to each urban indicator via least square method and the best fitting parameter α i (Pearson correlations as well) is shown in Table 2 for the two combinations of years. We have omitted the values of A i because they are very small (∼ 10 −6 ). Despite the increasing tendency observed for the values of α i over the years (expect for homicides and family income), we observe that α i < 1 for almost all urban indicators except for illiteracy (α i = 1.01 ± 0.01 for 2000 versus 1991 and α i = 1.05 ± 0.01 for 2010 versus 2000). These results agree with the evolution of the average D Yi (t), that is, indicators characterized by α i < 1 present a tendency of approaching the allometric laws, while for α i ≥ 1 there should be a departing tendency from the allometric laws.
In order to better understand the role of the parameter α i , we consider the limit   The bar plots show the evolution of the average D Yi (t) when grouping the cities whose urban indicator was (A) above and (B) below the allometric law in the year t = 1991. The colorful bars are empirical data (red for t = 1991, blue for t = 2000 and green for t = 2010) and gray bars represent the predictions obtained from the linear model of Eq. 6 (see discussion in main text). The error bars are 95% bootstrap confidence intervals for the average values. Note that the average D Yi (t) displays a statistically significant decreasing tendency for cities whose urban indicator was initially above the allometric law and an increasing tendency is observed for cities whose urban indicator was initially below the allometric law. The only exception to this pattern is the case of illiteracy, where the average D Yi (t) has increased for cities initially above the allometric law and it is almost a constant for those cities initially below. We further observe that the predicted values (gray bars) keep this main tendency.
where ∆t is small to rewrite Eq. 3 as the following differential equation     Table 2. We note that α i < 1 for almost all indicators, a fact that is in agreement with the approach to allometric laws observed for the averages value of D Yi (t) (Fig. 3). Illiteracy is the only exception where α i ≳ 1. This result also agrees with the departing from the allometric law of the average D Yi (t) in the case of illiteracy. whose solution is where k is an integration constant. Thus, Eq. 5 predicts that D Yi (t) will exponentially approach the value A i (1 − α i ) for α i < 1 and that they will exponentially increase over time when α i > 1. For α i = 1 we have a linear behavior for the evolution of D Yi (t). It is worth remembering that A i is a very small number (∼ 10 −6 ) for all urban indicators and hence the values of D Yi (t) are actually approaching zero for α i < 1 (that is, the values of the urban indicators are getting closer to value expected by the allometric law). We further observe that 1 (1 − α i ) plays the role of a characteristic time and if we assume α i < 1, the smaller the value of α i is, the faster we expect to be changes in the urban indicator. For population-related indicators, we observe an increasing tendency in the values of α i and also that these values are among the largest ones, which agrees with  for socio-economic indicators, indicating that these latter are more susceptible to natural or public policies driven changes. The illiteracy is unique because its value of α i , that was very close to one (for D Yi (2000) versus D Yi (1991)), has increased, a result that agrees with the detachment from the allometry reported in Fig. 3 and suggests an acceleration in this process. In addition of being autocorrelated with their past values, D Yi (t + ∆t) for the urban indicator i also displays statistically significant crosscorrelations with D Yj (t) for other indicators j (S4 Fig.). These memory effects and also the fact that the residuals surrounding the relationships D Yi (t + ∆t) versus D Yi (t) are very close to Gaussian distributions (S5 Fig.) with standard deviations across windows practically constant (S6 Fig.) make these scale-adjusted metrics particularly good for being used in linear regressions aiming forecasts. We have thus adjusted the linear model (via ordinary least-squares method) by considering the relationships (t + ∆t) = 2000 versus t = 1991 and (t + ∆t) = 2010 versus t = 2000. In Eq. 6, C k is the linear coefficient quantifying the predictive power of D Y k (t) on D Yi (t + ∆t) (C 0 is the intercept coefficient) and η i (t) is the noise term accounting for the effect of unmeasurable factors. The results exhibiting the linear coefficients of each linear regression for the two combinations of years are shown in S1 Text. We note that these simple models account for 31%-97% of the observed variance in D Yi (t + ∆t) and that they correctly reproduce the average values of the scale-adjusted metric above and below the allometric laws for the years 2000 and 2010 only using data from the years 1991 and 2000, respectively (see S7 Fig.). We have further compared the distributions of the empirical values of D Yi with the predictions of these linear models and observed that the agreement is remarkable good for the indicators elderly, female and male population as well as for illiteracy and income (S8 and S9 Figs.). Motivated by these good agreements, we proposed to forecast the values of D Yi (t + ∆t) in the year of 2020 (next Brazilian national census). In order to do so, we have considered that the linear coefficients C k are constant over time and employed the average value of C k over the two combinations of years used in Eq. 6 for predicting the values of D Yi (t + ∆t). It is worth noting that by assuming C k constant, we are ignoring the evolution of socio-economic and policy factors. In an ideal scenario, one could track the evolution of the values of C k for achieving more reliable predictions. However, our data (that is, the two values for C k ) do not enable us to probe possible evolutionary behaviors in the values of C k . Even so, as pointed by Bettencourt [6], the dynamics of the urban metrics seems to be dominated by long timescales (∼30 years), and thus the approach of constant coefficients should be seen as a first approximation. The grays bars in Fig. 3 show the averages D Yi (2020) after grouping the cities with D Yi (1991) > 0 and D Yi (1991) < 0. We observe that predictions for the average values basically keep the trends presented in the previous years; for unemployment, in which the trend was not very clear, the predictions put this indicator together with most indicators, where the average D Yi (t) has been decreasing for cities initially above the allometric law and increasing for those initially below.
In order to gain further information on the predictions, we build a geographic visualization of the expected changes in D Yi (t) between the years t = 2010 and t = 2020. The circles over the maps in . Perhaps, the most striking feature of these visualizations is the fact that the predicted changes appear spatially clustered for almost all indicators, which somehow reflects the geographic inequalities existing in Brazil; however, some intriguing patterns are indicator-dependent.
For child labour, D Yi (t) is expected to increase around three of the most densely populated regions that contain the metropolitan areas of São Paulo, Rio de Janeiro and the metropolitan areas of almost all northeast capitals; we further observe that a decrease in child labour cases is expected in mostly of the inner and southern cities. For the indicators elderly, female and male populations the clustering of the changes in D Yi (t) is quite evident: elderly and female populations are foreseen to decrease in mostly of the northeast cities and display an increasing tendency in large part of the other regions (male population behaves anti-symmetrically to the female population). In the case of homicides, our model predicts a decrease in D Yi (t) for the vast majority of southern cities and we further observe a stripe near the east coast where D Yi (t) is expected to decrease for mostly of cities (both densely populated areas); on the other hand, inner cities (specially inner cities from the state of São Paulo and northeastern region) are expected to increase D Yi (t), suggesting that this violent crime may be "moving" towards less populated areas of the interior of Brazil.
For illiteracy, again, the clustering of the changes is easily perceptible: we expect an increase in D Yi (t) for mostly of the northeast and northernmost cities, while a decrease is predicted for the majority of the cities from other regions (excluding several inner cities of southernmost region). For family income, we also observe a clustering in the changes where most northern cities are expected to increase the value of D Yi (t) (specially the inner cities of these regions), while for most cites in the central part of Brazil are expected to decrease the value of D Yi (t); these expected changes may be (at least in part) related to the "bolsa família" (family allowance) program -a large scale social welfare program of the Brazilian government (more than 14 million families were beneficiaries in 2013) for providing financial aid to poor families via direct cash transfer -because large part of the families receiving this aid are from the north and northeast regions. It is worth noting that for participating in the program, families must insure that their children attend to school and thus one would expect a reduction in D Yi (t) for illiteracy in same regions that concentrate the beneficiaries of the program, which was not predicted by the model. This result suggest that simply enforce school attendance may not be efficient for reducing illiteracy, which can also be explained by common observed poor conditions of the public school system of those regions. Finally, for unemployment, there is no clear clustering of the changes in D Yi (t), but instead we expect a decrease in D Yi (t) widespread throughout the Brazilian territory (except in the southernmost region, where we note the prevalence of light shades of red).

Summary
Our article discusses that despite being widely used in government reports and also in academic works, per capita indicators completely ignore the universal allometric laws that appear ruling the growth dynamics of cities. We discuss that per capita indicators can be biased towards small or large cities depending on whether we have sub or superlinear allometries with the population size. We thus employed a scale-adjusted metric by evaluating the difference between the actual values of urban indicators and the ones that are expected by the allometric relationships. When investigating the evolution of the scale-adjusted metrics, we have reported patterns that do not appear in the per capita indicators. The scale-adjusted metrics also display a linear correspondence with their past values, a feature that facilitates the use of linear regressions for modeling the urban indicators. By employing simple linear models for describing the scale-adjusted metrics based on their past values, we verified that these models account for 31%-97% of the observed variance and correctly reproduce the average values of the scale-adjusted metric. Assuming the linear coefficients constant over time, we present predictions for the values of the scale-adjusted metrics in year of 2020, when the next Brazilian census will happen. We observe that the predicted changes for the urban indicators appear (for most cases) spatially clustered, that is, forming regions where most cities are expected to increase or decrease the value of the scale-adjusted metric. Apart from this visualization of the predicted changes, we also provide a table (S1 Dataset) with the values of the scale-adjusted metrics for the three past census data as well as the predictions for the next one in the supplementary materials. We believe that our analysis may find potential applications on development of new policies and resources allocation in the context of urban planing. Finally, we note that the methods worked out here can also be directly applied in other contexts where allometries are present such as for economic indexes and biological quantities.

Data presentation
The data we analyzed consist of the population size N (t) and eight urban indicators Y i (t) for each Brazilian city in the years t = 1991, 2000 and 2010 in which the national census took place. We filter these data by selecting the 1605 cities for which all the eight urban indicators were available, this corresponds to 28.8% of the total number of Brazilian cities but account for 76.5% of the total population of Brazil. These data are maintained and made freely available by the Department of Informatics of the Brazilian Public Health System -DATASUS [27]. The eight urban metrics are defined as follows: Child labour: the proportion of the population aged 10 to 15 years who is working or looking for work during the reference week, in a given geographic area, in the current year; Elderly population: the number of inhabitant of a given city aged 60 years or older; Female population: the number of inhabitant of a given city that is female; Homicides: injuries inflicted by another person with intent to injure or kill, by any means. Illiteracy: it gives the number of inhabitants in a given geographic area, in the current year, aged 15 years or older, who cannot read and write at least a single ticket in the language they know; Family income: this indicator gives the average household incomes of residents in a given geographic area, in the current year. It was considered as family income the sum of the monthly income of the household, in Reals (Brazilian currency) divided by the number of its residents; Male population: the number of inhabitant of a given city that is male; Unemployment: it gives the number of inhabitant aged 16 years or older who is without working or looking for work during the reference week, in a given geographic area, in the current year.
Despite there being other definitions [28], the results presented here have been obtained by considering that cities are the smallest administrative units with a local government (municipalities or municípios). The other commonly employed definition is the metropolitan area, which is composed of more than one municipality and its is usually associated with the coalescence of several municipalities. As discussed by Bettencourt et al. [24], the choice of the "unit of analysis" is crucial when studying properties of cities. Regarding the scaling analysis: on one hand, the disaggregation of the correct urban definition can introduce a bias in the value of the scaling exponent (either by reducing or increasing its expectation value); on the other hand, the aggregation of the correct urban definition usually make the allometry more linear [24]. In fact, changes in the scaling exponents have been reported when choosing different definitions of city [29][30][31]. However, there is no fail-safe procedure for defining the correct boundaries of a city, also some urban indicators are actually more spatially restricted than others (for instance, homicides versus family income). We have also analyzed our data after considering simultaneously the municipalities that do not belong to any metropolitan area and by aggregating the municipalities of the 39 metropolitan areas existing in Brazil. Despite the observation of relatively small changes in the scaling exponents, our conclusions remain unaltered under this scenario.

Fitting allometric laws between urban indicators and population
As we previously mentioned, several works have reported the existence of robust allometric relationships between several urban indicators and the population size. Regarding Brazilian cities, scaling laws already have been identified for several urban indicators [5,[7][8][9][10][11][12][13][14], mainly because of the existence of reliable data made freely available by Brazilian agencies (DATASUS and IBGE). Here, we want to confirm that the urban scaling hypothesis holds for all our urban indicators and if these allometries have been changed over time. Specifically, we test the hypothesis that an urban indicator can be described by a power-law function of the population size, that is, Y i (t) = 10 Ai N (t) βi , where β i is the allometric (or scaling) exponent and A i is constant. In order to do so, we have plotted the logarithm of each urban indicator against the logarithm of the population and adjusted a linear model via orthogonal distance regression (as implement in the package scipy.odr of the Python library SciPy [32]) to all these relationships. Although the empirical relationships present different scattering degrees, they all display good quality linear relationships (Pearson correlation ranging from 0.44 to 1.00 -see Table 1) which are well described by linear models in log-log scale (see Fig 2 and S1 and S2 Figs.). From Table 1, we further note that the values of β i for illiteracy, family income and unemployment display a weak decreasing tendency over the years, while the other indicators show only small fluctuations (no clear evolutive tendency). A weak decreasing tendency for unemployment also appears in the work of Ignazzi [11] on the same data from the years of 2000 and 2010. The values of β i thus classify our indicators in two groups: female population, homicides and unemployment have superlinear relationships with the population (β i > 0); while child labor, elderly population, illiteracy, family income and male population have sublinear ones (β i < 0). It is worth noting that despite the allometric exponents β i be close to one for elderly, female and male populations, the allometries between these indicators and total population are almost perfectly correlated, producing values of β i very close but statistically different from one. We further observe that the values of the allometric exponents β i reported here may slightly differ from previous-reported one due the different fitting procedures as well as different urban definitions; however, these discrepancies are often very small (for instance, by considering generalized least squares via the Cochrane-Orcutt procedure and another definition of city, Ignazzi [11] have found β = 1.23 and β = 1.19 for unemployment, respectively in years of 2000 and 2010).