arxiv: v2 [econ.gn] 12 Oct 2018

Transcription

1 A new and stable estimation method of ountry eonomi fitness and rodut omlexity Vito D. P. Servedio 1, Paolo Buttà 2, Dario Mazzilli 3, Andrea Tahella 4, Luiano Pietronero 3 arxiv: v2 [eon.gn] 12 Ot Comlexity Siene Hub Vienna, Josefstätter-Strasse 39, A-80 Vienna, Austria 2 Deartment of Mathematis, Saienza University of Rome, P.le Aldo Moro 5, Roma, Italy 3 Physis Deartment, Saienza University of Rome, P.le Aldo Moro 5, Roma, Italy 4 Institute for Comlex Systems, CNR, Via dei Taurini 19, Rome, Italy Corresondene: v.servedio@gmail.om Otober 15, 2018 Abstrat We resent a new metri estimating fitness of ountries and omlexity of roduts by exloiting a non-linear non-homogeneous ma alied to the ublily available information on the goods exorted by a ountry. The non homogeneous terms guarantee both onvergene and stability. After a suitable resaling of the relevant quantities, the non homogeneous terms are eventually set to zero so that this new metri is arameter free. This new ma almost rerodues the results of the original homogeneous metris already defined in literature and allows for an aroximate analyti solution in ase of atual binarized matries based on the Revealed Comarative Advantage (RCA indiator. This solution is onneted with a new quantity desribing the neighborhood of nodes in biartite grahs, reresenting in this work the relations between ountries and exorted roduts. Moreover, we define the new indiator of ountry net-effiieny quantifying how a ountry effiiently invests in aabilities able to generate innovative omlex high quality roduts. Eventually, we demonstrate analytially the loal onvergene of the algorithm involved. Keywords: Eonomi Comlexity; Non-linear ma; Biartite networks 1 Introdution In the last deade a new aroah to maroeonomis has been develoed to better understand the growth of ountries [1, 2]. The key idea is to onsider the international trade of ountries as a roxy of their internal rodution system. By desribing the international trade as a biartite network, where ountries and roduts are sites of the two layers, new metris for the eonomy of ountries and the quality of roduts an be onstruted by leveraging the network strutures only [1]. These metris quantify the fitness of ountries, the quality of their industrial system, and the omlexity of ommodities by indiretly inferring the tehnologial requirements needed to rodue them. The mathematial roerties of the algorithm involved in the evaluation of the metris, as well as the eonomi meaning of the metris and ossible aliations, have been disussed in revious aers [3, 4, 5]. Moreover, these two new metris have been suessfully used to develo state-of-the-art foreasting aroahes for eonomi growth [6, 7, 8]. The very same aroah has been alied to different soial and eologial systems resenting a biartite network struture and a ometition between the omonents of the system [9, ]. Thus, it is natural to interret fitness and omlexity as roerties of the network underlying those systems. The revised version of the fitness-omlexity estimation metri that we show here, results in a lear 1

2 and natural interretation in terms of network roerties and hels to better understand the different omonents that ontribute to the fitness of ountries. In the following, we first desribe the original metris of fitness and omlexity and their roerties, underlining some ritial issues that we solve with this new version. Then, we define the new roedure ste by ste and highlight its advantages in the ase of ountries-roduts biartite networks. Finally, we devise an aroximated solution and disuss its interretation. In Aendix B we list the main quantities aearing in the text. 2 Metri Definition 2.1 The Original Metri Objet of this work is the network of ountries and their exorted goods. This network is of biartite tye (ountries and roduts are mutually linked, but no link exists between ountries as well as between roduts and weighted (links arry a weight s, i.e., the exorted volume of rodut of ountry, measured in US$. Data ranging from year 1995 to year 2015 an be freely retrieved from the Web [11], though we use them after a roedure to enhane their onsisteny [8]. Eventually, we ome u with data about 161 ountries and more than 4000 roduts, whih were ategorized aording to the Harmonized System 2007 oding system, at 6 digits level of oarse-graining. The weighted biartite network of ountries and roduts an be rojeted onto an unweighted network desribed solely by the M matrix with elements set to unity when a given ountry meaningfully exorts a good and zero otherwise (See Methods. The original metri estimating the fitness of ountries and omlexity of roduts was defined by the following non-linear iterative ma: F (n Q (n = = M Q (n 1 with 1 C ( 1 M /F (n 1 with 1 P, with initial values F (0 = Q (0 = 1,,. In the revious exression, F and Q stand for the fitness of a ountry and quality (omlexity of a rodut ; C and P are the total number of ountries and exorted roduts, resetively; and from the dataset we have that C P. By multilying all F and Q by the same numerial fator k, the ma remains unaltered, so that the fixed oint of the ma (as n is defined u to a normalization onstant. In the original method this onstant is hosen at eah iteration n suh that fitness and omlexity are onstrained to lie on the double simlex defined by: F (n = C and (1 Q (n = P. (2 The metri defined in Equations (1 and (2 suessfully ranks the ountries of our world aording to their otential tehnologial develoment and, when alied to different yearly time intervals, an be used to suggest reise strategies to imrove ountry eonomies. It has also been roved to give the orret ranking of imortane of seies in a omlex eologial system [9]. Desite its suess, some oints an still be imroved: i. Convergene issues: As stated in a reent aer dealing with the stability of alulating this metri [12]: If the belly of the matrix [M ] is outward, all the fitnesses and omlexities onverge to numbers greater than zero. If the belly is inward, some of the fitnesses will onverge to zero. Sine an inward belly is the rule rather than the exetion, some ountries will have zero fitness and as a result all the roduts exorted by them get zero omlexity (quality. This is mathematially aetable, though it heavily underestimates the quality of suh roduts: Even natural resoures need the right know-how to be extrated so that their quality would be better reresented by a ositive quantity. To ure this issue one has to introdue the notion of rank onvergene rather than absolute onvergene, i.e., the fixed oint is onsidered ahieved when the ranking of ountries stays unaltered ste by ste. 2

3 ii. Zero exorts: The ountries that do not exort any good do have zero fitness indeendently from their finite aabilities. iii. Seialized world: In an hyothetial world where eah ountry would exort only one rodut, different from all other roduts exorted by other ountries, this metri would assign a unity fitness and quality to all ountries and roduts. Though mathematially aetable, this solution does not take into aount the intrinsi omlexity of roduts. iv. Equation symmetry: This is rather an aestheti oint, in that Equation (1 are not ast in a symmetri form. 2.2 The New Metri First, we reshae Equation (1 in a symmetri form by introduing the variable P = Q 1, i.e.: { F (n P (n = M /P (n 1 with 1 C = M /F (n 1 with 1 P. (3 Now the quality of roduts are given by the quantities P 1 and the metri is trivially equivalent to the original one rovided one uses the normalization onditions F (n = C and (n (P 1 = P. Next, we introdue two set of quantities φ > 0 and π > 0 and onsider the inhomogeneous non-linear ma defined as: { F (n = φ + M /P (n 1 with 1 C P (n = π + M /F (n 1 (4 with 1 P. Sine the ma is no more defined u to a multiliative onstant, the normalization ondition is not required anymore, while the initial ondition an be set as in the original metri F (0 = P (0 = 1,,. The fixed oint of the transformation is now trivially haraterized by the onditions: F φ, P π, F P > M. (5 The arameters φ and π an be interreted as follows. The arameter φ reresents the intrinsi fitness of a ountry. In fat, for a ountry k that does not exort any good we have M k = 0 so that its fitness is simly equal to φ k. Irresetive of its exorts any ountry has a set of aabilities that haraterize it. The arameter π is more intriguing. If no ountry exorts it (robably beause no ountry rodues it, the rodut q has not been invented yet and its quality lies at its maximum value πq 1 sine M q = 0. Therefore, the inverse of π q may be interreted as a sort of innovation threshold: The smaller the arameter is, the higher is the quality of the rodut in his outset and more sohistiated aabilities are neessary to rodue it. On the other hand, roduts like natural resoures may be assoiated with a larger value of the arameter sine require less omlex aabilities for their extration. In order to kee the metri evaluating algorithm simle and arameter free as in the original ase, we set a ommon value φ = π = δ, then we study the deendene of the metris on δ, and finally we set δ = 0 (in fat renouning to ure the issue number iii. listed above. 3 Results 3.1 Deendene on the Non-Homogeneous Parameter We onsider φ = π = δ, and address the deendene of the fixed oint uon δ. To outline the deendene of F and P from the arameter δ, we use the relations defined in Equation (4 and introdue the resaled quantities P = P /δ and F = F δ. After some trivial algebra we get from Equation (4: { F (n P (n = δ 2 + (n 1 M / P with 1 C = 1 + (n 1 M / F with 1 P, (6 from whih we dedue that, as soon as the arameter δ 2 is muh smaller than the tyial value of M matrix elements, i.e., muh smaller than unity, the fixed oint in terms of F and P almost does not deend on δ (see Figure 1. It is worth noting that the values of fitness F and quality Q = P 1 of the 3

4 49.0 F * δ P / δ δ Figure 1: Deendene on the non-homogeneous arameter: Deendene of fitness and quality at the fixed oint on the arameter δ. One ountry (Afghanistan and one rodut (live horses were hosen arbitrarily from the samle of year original ma defined by Equations (1 and (2 annot be obtained in this new metri when the arameter δ tends to zero. In terms of F and P the fitness and quality obtained in the original metri an be exressed as F = F δ 1 1 and Q = P δ 1. Sine the new metri rovides finite non vanishing values of F and P, by taking the limit δ 0 would deliver infinite values of F and Q. We might think that the normalization roedure neessary in the old metri in order to fix the arbitrary onstant would get rid of the ommon fator δ 1 and deliver the same values of the new metri. Unfortunately, this is not the ase sine the new metri does not rely on a normalization roedure. Therefore, sine a self-onsistent roedure of normalization, i.e., a rojetion on the double simlex defined by Equation (2, is missing in the new metri, the results annot oinide. Sine the quantities F and P are well defined in the limit δ 0, we shall fous on them only, in the following. We remind that the omlexities of roduts and thus to the. In artiular, the 1 Q with the Q exressed as in the seond of Equation (1, but with the tilde quantities alulated in the new metri. Therefore, we shall assign to Q = ( P 1 1 the meaning of omlexity of roduts in our new metri. The differenes between the old and new metris are deited in Figure 2. delivered by the original metri are onneted to the set of P 1 seond of Equation (6 an be interreted at the fixed oint as: P = Analyti Aroximate Solution In this setion we shall rovide an aroximate analyti solution that an be used to estimate the values attained by the ma of Equation (6 at the fixed oint. Desite their symmetri shae, Equation (4 are not symmetri at all sine in ase of atual ountries and roduts, the matrix M is retangular with the number of its rows C being muh less than the number of its olumns P. To estimate the effet of this asymmetry, we first onsider Equation (4 in a mean field fashion, where eah element of M is set to the average value M =, M/CP, and write, at the fixed oint: { f = δ 2 + P M 1 = 1 + C M f 1 (7, with now all F and P set to be equal to their mean field value f and, resetively. By setting δ = 0, we find = 1/(1 C 1 + C and f = P C. Indeed, an aroximate exression for the P P fixed oint of Equation (6 in the regime δ 1 and C P an be derived also beyond the mean field aroximation. To this end, we set again δ = 0 and onsider the orresonding fixed oint equation assoiated to Equation (6, i.e.: { F = M / P with 1 C P = 1 + M / F (8 with 1 P. From the emirial struture of the matrix M, we observe that the quantity D = M,, reresenting the diversifiation of ountry, i.e., the number of different roduts exorted by, is of the order of P 1 4

5 original algorithm fitness (F Year 2007 Year original algorithm omlexity (Q ~ new algorithm fitness (F new algorithm omlexity (Q ~ Figure 2: Comarison between the original and the revised metri: Differenes in ountry fitness (left anel and rodut omlexity (right anel alulated in the original metri of Ref. [1] (vertial axes and new metri (horizontal axes as referred to year The green line in the left anel is the best least square aroximation of ower-law tye (orrelation oeffiient with exonent a The dark line in the right anel is the best ower-law aroximation (orrelation oeffiient resulting with an exonent of a The year 2007 was hosen randomly. Similar results aly to all the years onsidered. In artiular, the orrelation oeffiient and the exonent of the green line in the left anel lie between and 0.990, and 1.48 and 1.61 resetively throughout the years. For the blak line in the right anel we find a orrelation oeffiient between and 0.979, and an exonent between 1.34 and P, at least for the majority of ountries (as an average over all the years onsidered we find that 70% of the ountries have 0.1 D /P 1. Therefore, setting P = max P and F = min F, Equation (8 imlies: { F D / P onst P/ P with 1 C P 1 + C/ F. From the first estimate, F onst P/ P, and therefore, by the seond estimate, P 1+onst C P. As P P 1, we onlude that P = 1 + W with W in the order of magnitude of C/P, and, as a onsequene, F is of the order of magnitude of P. We next omute exliitly the values of F and P at the first order in this aroximation. The alulation of seond order terms an be found in Aendix A. By using the first order aroximation (1 + a 1 1 a twie, from Equation (8 we have: W M D D M W. Now let H be the square matrix of elements H = M T D 2 M. Letting D 1 be the olumn vetor with omonents 1/D and 1 the identity matrix, the last dislayed formula reads: (1 HW M T D 1. We now observe that: H 1/D 2 onst C/P2. Therefore, the matrix (1 H is lose to the identity (the orretion is of order C/P 2 and hene invertible (with also the inverse lose to the identity. In this aroximation, W = M T D 1, so that the resaled (reiroals of the qualities of roduts are given by: P = 1 + M T D 1. (9 In the same aroximation, we obtain the resaled fitnesses F ; sine: F = M 1 + W M (1 W, 5

6 ount % 0.5% 1% 1.5% 2% 2.5% relative differene F (fixed oint Figure 3: Numerial vs Analyti relative error: The histogram of the relative differene ( F (aroximated (fixed oint / F is lotted with the number of ountries on the vertial axis. The aroximated values are alulated using Equation (. we have: F = D KD 1, ( having introdued the o-rodution matrix K = MM T with elements K = M M, reresenting the number of the same roduts exorted by the two ountries and. It is interesting to note how, u to the first order aroximation, the values of the fitness of ountries are deending on the o-rodution matrix and diversifiation only. The goodness of the aroximations above an be areiated in Figure 3 that shows how the relative differene between the numerial values at the fixed oint and the aroximate solution of Equation ( is below 0.5% for more than 85% of the ountries. It is worth noting that in a reent work the Eonomi Comlexity Index (ECI defined in Ref. [2] has been onneted to the setral roerties of a weighted similarity matrix M resembling our o-rodution matrix K [13]. This similarity is only aarent sine in ECI the matrix M is defined as: M = M M /D U with U = M ubiquity of rodut, i.e., it ontains a further weighting term (the ubiquity in the sum defining it. Besides, the two metris of ECI and Fitness-Comlexity differ very muh from eah other: ECI relies on a linear homogeneous ma, while Fitness-Comlexity relies on a non-linear and in this work also non-homogeneous ma. 3.3 Country Ineffiieny and Net-Effiieny From Equation ( we dedue that the leading art of fitness F is given by the diversifiation D. The diversifiation of a ountry is indeed an imortant quantity, for the alulation of whih we do not need any omliated algorithm. On the other hand, what the non-linear ma roosed does, is to quantify how a ountry manages to suessfully differentiate its roduts, and indiretly offers an estimate of the aabilities of a nation. In fat, a ountry exorting mainly raw materials would be less effiient with reset to a ountry exorting high tehnologial goods, when they have the same diversifiation value. For this reason, we introdue the new quantity I = D F, ineffiieny of ountry : the smaller the value I the more effiient is the diversifiation it hooses. From the aroximate solution dislayed in Equation (, we get that I K /D, so that the ineffiieny of a ountry is a weighted average of its o-rodution matrix elements. The deendene of the ountry ineffiieny on the diversifiation is dislayed in Figure 4, while a visual reresentation of it is dislayed in Figure 5. It is interesting to notie how a lear ower-law deendene exists between the ineffiieny and the diversifiation of a ountry. The struture of the M matrix is suh that those ountries with high diversifiation also exort low quality goods in average. Therefore to a large diversifiation would statistially orresond a large ineffiieny, though the found ower-law is not trivial and deends on the struture of the M. A 6

7 0 Year 2007 Country ineffiieny (I net effiieny (N JPN KOR CHE D Country diversifiation (D Figure 4: Role of diversifiation: The ountry ineffiieny I = D F is lotted vs. the diversifiation D with the blak line reresenting the ower-law relation I D 0.75 (linear regression with orrelation oeffiient In the inset the net effiieny N, defined as the differene between the blak line and the ineffiieny of the main grah, is shown. Plotted data ertain to year We find a similar behaviour for all the years onsidered with the exonent of D between 0.73 and 0.76, and the orrelation oeffiient between and similar ower-law behaviour is found between the fitness alulated with the traditional metri and the diversifiation, but with a different exonent (from the left anel of Figure 2 we dedue that there is a ower-law relation between the fitnesses alulated with the original metri and this new metri, and the exonent is around 1.53; sine the fitness F alulated with the new metri goes as D at the first order, then the old fitnesses also go as D In order to better areiate the rodution strategies of ountries, we subtrated the ommon ower-law trend of the deendeny of the ineffiieny on the diversifiation for eah year, hanged its sign and lotted the result in the right anel of Figure 6, whih thus shows the time evolution of a quantity that we all ountry net-effiieny N (net in the sense oosed to gross over the years It interesting to note how ountries behave differently over the time lase onsidered. Some ountries dislay a dereasing net-effiieny, others an inreasing or a onstant one. What many of these urves have in ommon is the dereasing set u around year 2000, more ronouned in the ase of higher develoed ountries, whih lie at high net-effiieny in the grah. 3.4 Loal Convergene From the simulations it is lear that the fixed oint obtained by iterating Equation (4 is loally stable. We an also rove it by resorting to the Jaobian of the transformation, in the ase of ountries and roduts. First we reall that the sum over the indexes and of Equation (4 run from 1 to C and P, resetively, with usually C P. In the ase of ountries and roduts C/P 1. We also fix φ = π = δ 1, so that the fitnesses and the (reiroals of the qualities at the fixed oint are aroximately given by: F = F /δ and P = δ P with F and P the omonents of the vetors F and P given in Equations ( and (9 resetively. Next, we alulate the Jaobian of the transformation at the fixed oint, whih an be simly exressed as the blok anti-diagonal matrix: ( J = 0 M T F 2 MP 2 0, (11 having introdued the diagonal matries F = diag(f 1, F 2,..., F and P = diag(p 1, P 2,..., P, resetively. We laim that the setral radius ρ(j of the square matrix J is stritly smaller than one. Denoting by σ(j the setrum of J, this means that ρ(j := max{ λ : λ σ(j} < 1. From this it follows [14] that the fixed oint is asymtotially stable and the onvergene exonentially fast. To rove the laim we onsider the square of the Jaobian that an be written as a blok diagonal matrix, J 2 = ( M T F 2 MP MP 2 M T F 2, (12 7

8 Figure 5: Ineffiieny artoon: Large ovals reresent three ountries, while small irles reresent roduts. In this simle examle, the ineffiieny I 1 of ountry 1 is I 1 = K 12 /D 2 + K 13 /D 3. From the figure we get K 12 = 2 and K 13 = 4, i.e. the number of roduts exorted by both ountries (the ardinality of the intersetion sets, and the diversifiations D 1 = 17, D 2 = 5, D 3 = 20. Thus, I 1 = 2/5 + 4/20 = 0.6 and the aroximated fitness F and note that the traes of the two matries on the diagonal is the same by alying a yli ermutation. Notiing that F P = F P and using the aroximate solutions in Equations ( and (9, we find with simle algebra that: Tr(J 2 = 2, M 2, F 2 P 2 2, M 2, D 2 = 2 1 C < 1. (13 D P Moreover, we an write the two non trivial matries omosing J 2 as: and: M T F 2 MP 2 = P(P 1 M T F 1 (F 1 MP 1 P 1 = PA T AP 1 (14 MP 2 M T F 2 = F(F 1 MP 1 (P 1 M T F 1 F 1 = FAA T F 1, (15 with A = F 1 MP 1. The matries AA T and A T A are symmetri and ositive-semidefinite so that their eigenvalues are real and non negative, and the matries FAA T F 1 and PA T AP 1 have the same eigenvalues. Therefore, the eigenvalues of J 2 are real and non negative and we an write aording to Equation (13: Tr(J 2 = λ 2 i < 1, (16 i with λ i eigenvalues of J. Finally, from the reeding equation we have max λ 2 i < max λ i < 1 so that at the fixed oint ρ(j < Robustness to Noise Fitness and omlexity (quality values deend on the struture of the matrix M. Noise an affet its elements by fliing their value. Thus, we test the robustness of the new metri to noise as desribed in [15]. The idea is to introdue random noise by fliing eah single bit of the matrix with robability η, whih then is a arameter tuning the noise level. The rank of ountry fitnesses in resene of noise 8

9 500 CHN ITA DEU 20 F * δ USA IND BRA RUS PAN AFG net-effiieny (N BRA CHN DEU GRC IND ITA JAP RUS KOR CHE USA Year IRQ year Figure 6: Time evolution of fitness and net-effiieny: (Left anel Country fitness yearly evolution as estimated by the new metri. (Right anel yearly time evolution of ountry net effiieny. The net effiieny is a detrended version of the ineffiieny defined in the text and dislayed in the inset of Figure 4 for the year Curves were artifiially smoothed by a ubi sline for a better visual reresentation. 1.0 Searman orrelation (ρ s original new noise level (η Figure 7: Noise robustness: Searman orrelation between the ranking of ountries based on fitness at zero noise and at different noise levels η (see Setion 3.5 in the main text. The erformane of the two metris is ratially indistinguishable. Note that at η = 1 all the elements of matrix M are flied so that the erturbed system is erfetly anti-orrelated with the original one. R η is then omared with the rank obtained without noise R. 0 The Searman orrelation ρ s is then evaluated between these two sets and shown in Figure 7 as a funtion of η for both the original and the new metris: The new metris show a erfet stability to random noise as the original ones with an unavoidable transition around η 0.5, where noise is so strong to alter signifiantly the struture of the matrix M. 4 Disussion The roosed new inhomogeneous non-linear metri to estimate eonomi fitness and omlexity defined in Equations (4 and (6 arries many advantages with reset to the original one. The fitnesses and omlexities resulting from these two aroahes are not idential, but highly orrelated to eah other as witnessed by the lots in Figure 2. This high orrelation between the two metris ensures that all the studies arried on with the original metri so far, an be obtained by alying this new metri as well. Besides the stability of the metri and its robustness, one more advantage is that the fitness is well defined also for those ountries that have low exortation volumes and that in the original metri had 9

10 20 JAP DEU net effiieny 0 RUS KOR BRA CHE GBR USA ITA CHN ESP GRC IND ~ fitness (F Figure 8: Net-effiieny vs fitness: eah line orresonds to the time evolution of the onnetion between the fitness of a ountry and its net-effiieny in the eriod between 1995 and This figure onnets the quantities on the vertial axes of the lots dislayed in Figure 6. Lines start from a large irle (year 1995 and end with a small one (year their fitness tending to zero. For those ountries it is now ossible to undertake a omarative study based on hyothetial investments (hanging the elements of the M matrix so to make reditions on their eonomi imat. By first symmetrising the original equations, by adding an inhomogeneous arameter and by resaling the quantities, one obtains Equation (6, where the arameter an be safely set to zero. This ensures that this new metri is arameter free as the original one. As a leasant side effet, the fixed oint of the ma an be well aroximated analytially, with an error with reset to the iterative fixed oint of less than 3% (see Figure 3. The result is reresented by Equations (9 and ( at the first order (Equations (19 and (20 at the seond order, whih allow for a simle intuitive exlanation of the omlexity of roduts and fitness of ountries. Let us disuss Equation ( first. The result suggests that the fitness of a ountry is trivially related, at the first order, to its diversifiation: The more roduts a ountry exorts, the larger is its fitness, i.e., the more develoed its aabilities. This simle exliit deendene of the fitness on the diversifiation is also an advantage with reset to the original metri, where the deendene was not exliitly lear. The seond term of Equation (, whih we all ineffiieny, is also very interesting. If a ountry is the only one to exort a given rodut, the ontribution of this rodut to its fitness is a full one, or in other words, the ontribution to the ineffiieny is zero. This situation mimis a ondition of monooly on that rodut and it is logial that the exorting ountry has the full benefit of it. When a rodut is exorted by multile nations then it is ritial to assess whether those ountries exort few or many other roduts (see Figure 5. If a rodut is exorted by a ountry with low diversifiation (low aabilities, then that rodut is not suosed to be of high omlexity. The result is that the ratio K /D an be lose to one ( = 1, = 2 in the figure and the ineffiieny assoiated to the ommon roduts is high, resulting in a small ontribution to the fitness of. The ineffiieny an be interreted in terms of the biartite network of ountries and roduts: The K ounts the number of links that onnet ountries and to the same roduts, while the differentiation D is the node degree of ountry. In other words, for a ountry the ineffiieny ounts the links to ommon roduts of all other ountries and weights them aording to the degree of those. To our knowledge, this kind of measure has never been onsidered in omlex networks so far. Sine, statistially, ountries with an high diversifiation also exort many less omlex roduts, the ineffiieny is an inreasing funtion of the diversifiation (Figure 4, main grah. If we subtrat the general trend, whih stems from the struture of the matrix M, we an areiate the net effet of seleting the goods to exort. We all this new de-trended quantity net-effiieny. In this way we somehow remove the negative effet of less valuable roduts and highlight the ontribution of more sohistiated goods. In the inset of Figure 4

11 we show the net-effiieny as a funtion of diversifiation and underline the three nations (Jaan, Korea, and Switzerland that stand out among the others. The time evolution of this new quantity is shown in the right anel of Figure 6. We an ombine the time evolution of both fitness and net-effiieny for a given ountry to determine to what degree they are orrelated. Figure 8 shows the two quantities for seleted ountries. It is lear how these two quantities are not related to eah other and reresent two omlementary information. In fat, the fitness is mainly onneted to the rodut diversifiation of a ountry (Equation ( while the net-effiieny is onneted to how omlex are the exorted roduts. In the figure we see two extreme ases reresented by Switzerland (CHE and South Korea (KOR, whose lines are ratially orthogonal. While Switzerland have dereased the number of different exorted goods in the years but have still exorted omlex roduts, South Korea have ket the number of exorted roduts as almost onstant but have inreased their omlexity. The oosite situation of South Korea we notie for Germany (DEU, where the omlexity of the exorted goods has dereased in time. Interestingly, China (CHN has systematially inreased both the number of exorted goods and their omlexity, whih we interret as a symtom of a solid eonomy in exansion. The omlexity of roduts is estimated by Equation (9 as the reiroal of the seond term of the sum. Sine the diversifiation of a ountry D is a diret measure of its aabilities, we exet to find a simle relation between it and the omlexities of roduts Q. Indeed, if we indiate with i those ountries exorting the rodut, for whih obviously we have M i = 1, and with m = M, we an write: ( 1 Q D 1 D 2 D m from whih we orroborate the main idea that the omlexities of roduts are driven by the ountries with low diversifiation (aabilities that exort it. Just for amusement, we observe how the omlexity of roduts an be onsidered as the equivalent resistor of a arallel of resistors eah one with resistane D. Somehow, a high D reresents an effetive resistane to the reation of a rodut and its exort, so that if a ountry exists with a low diversifiation exorting it, the effort (resistane of roduing that rodut is also low. 5 Materials and Methods 5.1 Constrution of the M Matrix We exloit the UN-COMTRADE data set [11], where re-exort and re-imort fluxes are exliitly delared, allowing us to exlude them from the analysis. As reorted by UNSTAT in Ref. [16], the 81.8% of the whole data set (96.8% in ase of develoed ountries does not aount for goods in transit. Moreover, ommodities that do not ross borders are not inluded in the data. Given the exort volumes s of a ountry in a rodut one an evaluate the Revealed Comarative Advantage (RCA indiator [17] defined as the ratio: RCA = s s / s s (17 in this way one an filter out size effets. As desribed in the Sulementary information of [8], from the time series of the RCA we an evaluate the rodutive ometitiveness of eah ountry in eah rodut by assigning to it a rodutivity state from 1 to 4. State 1 means that the ountry does not rodue (or is very unometitive in roduing a rodut, state 4 means that it is one of the main roduer in the world. We an then rojet this states onto the binarized matrix M by simly setting its elements to unity whenever a state larger than 2 is enountered, and set them to null otherwise. Aknowledgments We aknowledge A. Zaaria and E. Pugliese for interesting disussions. We are grateful to an anonymous reviewer for suggesting to show Fig. 8. We aknowledge the EU FP7 Grant rojet GROWTH- COM and CNR PNR Projet CRISIS Lab and the Austrian Researh Promotion Ageny FFG under grant #857136, for finanial suort. FFG is also overing the osts to ublish in oen aess. 11

12 Seond order exansion of fitness and quali- Aendix A ties In this setion we omute exliitly the values of F and P for C P u to the seond order of magnitude of C/P. Letting ε = C/P, we exand W = εw (1 + ε 2 W (2 + O(ε 3. By assuming D of the order of P and by using the seond order aroximation (1 + a 1 1 a + a 2 twie, Eq. (8 imlies that ( F = D 1 ε M W (1 ε 2 M [ (2 W D (W (1 D 2] + O(ε 3, with 1 C, (18 and εw (1 + ε 2 W (2 = M + ε M D D, + ε 2 M D,, M D M W (1 + ε 2 M D D, M [ (2 W (W (1 D 2] M W (1 W (1 D + O(ε 3, with 1 P. By the assumtion on the magnitude of D, the first sum in the right-hand side is of the order of ε, the seond one is of the order of ε 2, while the last two sums are of the order ε 3. Therefore, εw (1 = M D, ε 2 W (2 =, M D M εw (1. D Realling H denotes the square matrix of elements H = M T D 2 M (hene H ε/p and D 1 the olumn vetor with omonents 1/D, we have just showed that W = M T D 1 + HM T D 1 + O(ε 3. Therefore, in the seond order aroximation, the resaled (reiroals of the qualities of roduts are given by P = 1 + M T D 1 + HM T D 1. (19 In the same aroximation, from Eq. (18 we finally alulate the resaled fitnesses F. Denoting by (M T D 1 2 the olumn vetor with omonents (M T D 1 2 we get F = D KD 1 + M(M T D 1 2 MHM T D 1, (20 where the o-rodution matrix K = MM T has been introdued just below Eq. (. Aendix B text Imortant quantities defined throughout the C, P: Total number of ountries and roduts M: Binary matrix with element M = 1 if ountry is a ometitive ountry in exorting rodut ; M = 0 otherwise; exort ometitiveness is estimated by means of exort volumes F, Q : Fitness of ountry and quality (omlexity of rodut at the fixed oint P : Inverse of the quality of rodut ; it is a sort of rodut simliity (P = (Q 1 δ: Inhomogeneous arameter; this arameter is ruial in ahieving a stable algorithm to evaluate the metris; it will be eventually let go to 0 to get a arameter free metri F, P : Resaled versions of the orresonding un-tilded quantities: F = F δ, P = P /δ; these quantities do not deend on δ as soon as δ 0 and are better suited to reresent fitness and omlexity rather than the un-tilded ones Q : Similar to the omlexity Q above, but for the new metris alulated with the inhomogeneous algorithm: Q = ( P 1 1 K: Corodution matrix with element K equal to the number of the same roduts exorted by ountries and ; K = MM T D : Diversifiation of ountry, i.e., the number of roduts the ountry is ometitive in exorting I : Ineffiieny of ountry defined as I = D F ; it reresents the fitness enalty resulting from exorting goods that are also exorted by other ountries N : Net-effiieny of ountry ; it is a de-trended version of the ineffiieny; in the dataset onsidered N D 0.75 I ; it reresents how effetively a ountry diversifies its exorted goods by fousing on roduts not exorted by others, whih are usually among the most omlex 12

13 Referenes [1] Andrea Tahella, Matthieu Cristelli, Guido Caldarelli, Andrea Gabrielli, and Luiano Pietronero. A new metris for ountries fitness and roduts omlexity. Sientifi Reorts, 2:723, [2] César A Hidalgo and Riardo Hausmann. The building bloks of eonomi omlexity. Proeedings of the National Aademy of Sienes, 6(26: , [3] Andrea Tahella, Matthieu Cristelli, Guido Caldarelli, Andrea Gabrielli, and Luiano Pietronero. Eonomi omlexity: onetual grounding of a new metris for global ometitiveness. Journal of Eonomi Dynamis and Control, 37(8: , [4] Matthieu Cristelli, Andrea Gabrielli, Andrea Tahella, Guido Caldarelli, and Luiano Pietronero. Measuring the intangibles: A metris for the eonomi omlexity of ountries and roduts. PloS one, 8(8:e70726, [5] Andrea Zaaria, Matthieu Cristelli, Roland Kuers, Andrea Tahella, and Luiano Pietronero. A ase study for a new metris for eonomi omlexity: The netherlands. Journal of Eonomi Interation and Coordination, 11(1: , [6] Fabio Sarao, Riardo Di Clemente, Andrea Gabrielli, and Tiziano Squartini. Deteting early signs of the risis in the world trade. Sientifi Reorts, 6:30286, Jul [7] Matthieu Cristelli, Andrea Tahella, Masud Cader, Kirstin Roster, and Luiano Pietronero. On the Preditability of Growth; The World Bank, Poliy Researh Working Paer , doi:.1596/ [8] Andrea Tahella, Dario Mazzilli, and Luiano Pietronero. A dynamial systems aroah to gross domesti rodut foreasting. Nature Physis, 14(8:861, [9] Virginia Domínguez-Garía and Miguel A. Muñoz. Ranking seies in mutualisti networks. Sientifi Reorts, 5:8182, Feb [] Giulio Cimini, Andrea Gabrielli, and Franeso Sylos Labini. The sientifi ometitiveness of nations. PLOS ONE, 9(12:1 11, [11] htt://omtrade.un.org/. UN Comtrade International Trade Statistis Database. Last aessed on Ot. 8th [12] Emanuele Pugliese, Andrea Zaaria, and Luiano Pietronero. On the onvergene of the fitnessomlexity algorithm. The Euroean Physial Journal Seial Tois, 225(: , Ot [13] P. Mealy, J. Doyne Farmer, and A. Teytelboym. A New Interretation of the Eonomi Comlexity Index. Available online: htt://dx.doi.org/.2139/ssrn (Aessed on 30 Setember [14] Boris Hasselblatt and Anatole Katok. A first ourse in dynamis. With a anorama of reent develoments. Cambridge University Press, New York, Proosition [15] Federio Battiston, Matthieu Cristelli, Andrea Tahella, and Luiano Pietronero. How metris for eonomi omlexity are affeted by noise. Comlexity Eonomis, 3(1:1 22, [16] United Nations, Statistial Division. Overview of national omilation and dissemination raties. In International Merhandise Trade Statistis: Sulement to the Comliers Manual; United Nations ubliations: New York, NY, USA, 2008; Chater 1, [17] Bela Balassa. Trade liberalisation and revealed omarative advantage. The Manhester Shool, 33(2:99 123,