Next Article in Journal
Calculation of Entropy for a Sinusoid with Beta-Distributed Phase
Next Article in Special Issue
Dispersal (Entropy) and Recognition (Information) as Foundations of Emergence and Dissolvence
Previous Article in Journal
A Story and a Recommendation about the Principle of Maximum Entropy Production
Previous Article in Special Issue
A Lower-Bound for the Maximin Redundancy in Pattern Coding
 
 
Comment published on 22 December 2009, see Entropy 2009, 11(4), 1121-1122.
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Equiprobability, Entropy, Gamma Distributions and Other Geometrical Questions in Multi-Agent Systems

1
Department of Computer Science, Universidad de Zaragoza, Pedro Cerbuna 12, Zaragoza 50009, Spain
2
Department of Physics, Universidad de Extremadura, Avda. de Elvas s/n, Badajoz 06071, Spain
3
BIFI, Universidad de Zaragoza, Corona de Aragón 42, Zaragoza 50009, Spain
*
Author to whom correspondence should be addressed.
Submission received: 3 November 2009 / Accepted: 30 November 2009 / Published: 2 December 2009
(This article belongs to the Special Issue Information and Entropy)

Abstract

:
A set of many identical interacting agents obeying a global additive constraint is considered. Under the hypothesis of equiprobability in the high-dimensional volume delimited in phase space by the constraint, the statistical behavior of a generic agent over the ensemble is worked out. The asymptotic distribution of that statistical behavior is derived from geometrical arguments. This distribution is related with the Gamma distributions found in several multi-agent economy models. The parallelism with all these systems is established. Also, as a collateral result, a formula for the volume of high-dimensional symmetrical bodies is proposed.
Classification:
PACS 87.23.Ge; 05.90.+m; 89.90.+n

1. Introduction

Nowadays different approaches to study economic systems can be taken within a multi-agent framework. Econophysics applies the techniques and tools of statistical physics to the understanding of this type of systems [1]. From this point of view, an economic system is regarded as a set of many agents exchanging money or trading with products of the most diverse origins. This exchange, which can be done in a deterministic or in a random way, finally generates an asymptotic wealth distribution. One of the goals of econophysics is to propose different models for the agents’ interactions and, by making the respective simulation experiments, to try to reproduce in a statistical way the real asymptotic wealth distribution in society [2].
In this paper, a general multi-agent system evolving under an additive constraint is considered. The geometrical properties of this system in phase space are exploited in order to obtain its statistical behavior. A striking coincidence is observed. Some dynamical mechanisms that have been proposed in the literature to model the interaction among agents in economic systems [3,4,5,6,7,8,9] provoke the same statistical results in the asymptotic wealth distribution than those derived from the geometrical properties in our multi-agent system under the assumption of equiprobability. This fact bring us to speculate on the possible close relationship that there exists between the local interactions among the agents in the former economic systems and the global geometrical conformation in phase space of our general system. Hence, with an adequate change of coordinates, we suggest that those economic systems asymptotically evolve in an equiprobable way over the volume of accessible states in the transformed phase space.
We start in Section 2 by recalling the derivation done for some particular geometries for the cases in which the constraint has a linear or quadratic dependence on the variables defining the agents. Then, in Section 3, the statistical behavior for a more general constraint is obtained. In Section 4, we speculate on the possible relationship with some economic systems [3,4,5,6,7,8,9] in which the Gamma distributions are also obtained. A formula for the volume of high-dimensional symmetrical bodies is proposed in Section 5. The last Section 6 contains our conclusions.

2. Recalling Some Results

2.1. Linear constraint

Let us assume N agents interacting in an open economy, each one with coordinate x i , i = 1 , , N , with x i 0 representing the wealth or money of the agent i, and a total available amount of money E. The additive linear constraint reads:
x 1 + x 2 + + x N - 1 + x N E
The result here explained for an open economy is also obtained when constraint (1) is an equality, i.e., when the economy is closed [10]. Under random evolution rules for the exchanging of money among agents [3], let us suppose that this system evolves in the interior of the N-dimensional pyramid given by Equation (1). We can suppose that the state or the bank system of western societies plays in this model the role of a heat reservoir that supplies money instead of energy. The formula for the volume V N ( E ) of an equilateral N-dimensional pyramid formed by N + 1 vertices linked by N perpendicular sides of length E is
V N ( E ) = E N N !
If each point on the N-dimensional pyramid is equiprobable, then the probability f ( x i ) d x i of finding the agent i with money x i , with normalization condition 0 E f ( x i ) d x i = 1 , is proportional to the volume formed by all the points into the ( N - 1 ) -dimensional pyramid having the ith-coordinate equal to x i . Then f ( x i ) verifies
f ( x i ) = V N - 1 ( E - x i ) V N ( E )
If we call ϵ the mean wealth per agent, E = N ϵ , then in the limit of large N ( N ), we have
lim N 1 V N - 1 ( E - x ) V N ( E ) = 1 ϵ e - x / ϵ
where the index i has been removed because the distribution is the same for each agent, and thus the wealth distribution can be obtained by averaging over all the agents,
f ( x ) d x = ϵ - 1 e - x / ϵ d x
This Boltzmann-Gibbs distribution has been found to fit the real distribution of incomes in western societies [3].

2.2. Quadratic constraint

Now let us suppose a one-dimensional ideal gas of N non-identical classical particles with masses m i , with i = 1 , , N , and total maximum energy E. If particle i has a momentum m i v i , we define a kinetic energy:
K p i 2 1 2 m i v i 2
where p i is the square root of the kinetic energy. Then the quadratic constraint reads:
p 1 2 + p 2 2 + + p N - 1 2 + p N 2 E
The distribution for p i here derived for an open system is also obtained when constraint (7) is an equality, i.e., when the energy is fixed [11]. When a finite number of particles are present in the system the asymptotic distribution maximizes the Tsallis entropy [12]. In the case with constraint (7) the system has accessible states with different energy, which is supposed to be supplied by a heat reservoir. These states are all those enclosed into the volume of the N-sphere given by Equation (7), with radius E 1 / 2 . The formula for the volume V N ( R ) of an N-sphere of radius R is
V N ( R ) = π N 2 Γ ( N 2 + 1 ) R N
where Γ ( · ) is the Gamma function. If we suppose that each point into the N-sphere is equiprobable, then the probability f ( p i ) d p i of finding the particle i, with coordinate p i (energy p i 2 ) and normalization condition - R R f ( p i ) d p i = 1 , is proportional to the volume formed by all the points into the ( N - 1 ) -sphere having the ith-coordinate equal to p i . Then f ( p i ) verifies
f ( p i ) = V N - 1 ( ( E - p i 2 ) 1 / 2 ) V N ( E 1 / 2 )
If we call ϵ the mean energy per particle, E = N ϵ , then in the limit of large N ( N ), we have
lim N 1 V N - 1 ( ( E - p 2 ) 1 / 2 ) V N ( E 1 / 2 ) = 1 2 π ϵ - 1 / 2 e - p 2 / 2 ϵ
where the index i has been removed because the distribution is the same for each particle. Thus the asymptotic distribution
f ( p ) = 1 2 π ϵ - 1 / 2 e - p 2 / 2 ϵ
can be obtained by averaging over all the particles. If the change of variables p = m 2 v is performed, with v the generic velocity of a particle, then the Maxwellian distribution is just derived from geometrical arguments.

3. Multi-Agent Systems and Equiprobability: General Derivation of the Asymptotic Distribution

In this section, we address the same problem above presented but in a general way. Let b be a positive real constant (cases b = 1 , 2 have been indicated in the former section). If we have a set of positive variables ( x 1 , x 2 , , x N ) verifying the constraint
x 1 b + x 2 b + + x N - 1 b + x N b E
with an adequate mechanism assuring the equiprobability of all the possible states ( x 1 , x 2 , , x N ) into the volume given by expression (12), will we have for the generic variable x the distribution
f ( x ) d x ϵ - 1 / b e - x b / b ϵ d x
when we average over the ensemble in the limit N , E , with E = N ϵ , and constant ϵ?. Now it is shown that the answer is affirmative. Similarly, we claim that if the weak inequality (12) is transformed in equality the result will be the same, as it has been proved for the cases b = 1 , 2 in References [10,11].
From the cases b = 1 , 2 , (see Equations (3) and (9)), we can extrapolate the general formula that will give us the statistical behavior f ( x ) of the generic variable x, when the system runs equiprobably into the volume defined by a constraint of type (12). The probability f ( x ) d x of finding an agent with generic coordinate x is proportional to the volume V N - 1 ( ( E - x b ) 1 / b ) formed by all the points into the ( N - 1 ) -dimensional symmetrical body limited by the constraint ( E - x b ) . Thus, the N-dimensional volume can be written as
V N ( E 1 / b ) = 0 E 1 / b V N - 1 ( ( E - x b ) 1 / b ) d x
Taking into account the normalization condition 0 E 1 / b f ( x ) d x = 1 , the expression for f ( x ) is obtained:
f ( x ) = V N - 1 ( ( E - x b ) 1 / b ) V N ( E 1 / b )
The N-dimensional volume, V N ( b , ρ ) , of a b-symmetrical body with side of length ρ is proportional to the term ρ N and to a coefficient g b ( N ) that depends on N:
V N ( b , ρ ) = g b ( N ) ρ N
The parameter b indicates the original Equation (12) that defines the boundaries of the volume V N ( b , ρ ) . Thus, for instance, from Equation (2), we have g b = 1 ( N ) = 1 / N ! .
Coming back to Equation (15), we can manipulate V N ( ( E - x b ) 1 / b ) to obtain (the index b is omitted in the formula of V N ):
V N ( ( E - x b ) 1 / b ) = g b ( N ) ( E - x b ) 1 / b N = g b ( N ) E N b 1 - x b E N b
If we suppose E = N ϵ , then ϵ represents the mean value of x b in the collectivity, that is, ϵ = < x b > . If N tends toward infinity, it results:
lim N 1 1 - x b E N b = e - x b / b ϵ
Thus,
V N ( ( E - x b ) 1 / b ) = V N ( E 1 / b ) e - x b / b ϵ
Substituting this last expression in formula (15), the exact form for f ( x ) is found in the thermodynamic limit ( N , E ):
f ( x ) d x = c b ϵ - 1 / b e - x b / b ϵ d x
with c b given by
c b = g b ( N - 1 ) g b ( N ) N 1 / b
Hence, the conjecture (13) is proved.
Doing a thermodinamical simile, we can calculate the dependence of ϵ on the temperature by differentiating the entropy with respect to the energy. The entropy can be written as S = - k N 0 f ( x ) ln f ( x ) d x , where f ( x ) is given by Equation (20) and k is the Boltzmann constant. If we recall that ϵ = E / N , we obtain
S ( E ) = k N b ln E N + k N b ( 1 - b ln c b )
where it has been used that ϵ = < x b > = 0 x b f ( x ) d x
The calculation of the temperature T gives
T - 1 = S E N = k N b E = k b ϵ
Thus ϵ = k T / b , a result that recovers the theorem of equipartition of energy for the quadratic case b = 2 . The distribution for all b is finally obtained:
f ( x ) d x = c b b k T 1 / b e - x b / k T d x

4. Gamma Distributions, Economic Gas Models and Geometry: A Speculation

If we perform the change of variables y = ϵ - 1 / b x in the normalization condition of f ( x ) , 0 f ( x ) d x = 1 , where f ( x ) is expressed in (20), we find that
c b = 0 e - y b / b d y - 1
If we introduce the new variable z = y b / b , the distribution f ( x ) as function of z reads:
f ( z ) d z = c b b 1 - 1 b z 1 b - 1 e - z d z
Let us observe that the Gamma function appears in the normalization condition,
0 f ( z ) d z = c b b 1 - 1 b 0 z 1 b - 1 e - z d z = c b b 1 - 1 b Γ 1 b = 1
This implies that
c b = b 1 - 1 b Γ 1 b
By using Mathematica the positive constant c b is plotted versus b in Figure 1. We see that lim b 0 c b = , and that lim b c b = 1 . The minimum of c b is reached for b = 3 . 1605 , taking the value c b = 0 . 7762 . Still further, we can calculate from Equation (28) the asymptotic dependence of c b on b:
lim b 0 c b = 1 2 π b e 1 / b 1 - b 12 +
lim b c b = b - 1 / b 1 + γ b +
where γ is the Euler constant, γ = 0 . 5772 . The asymptotic function (29) is obtained after substituting in (28) the value of Γ ( 1 / b ) by ( 1 / b - 1 ) ! , and performing the Stirling approximation on this last expression, knowing that 1 / b . The function (30) is found after looking for the first Taylor expansion terms of the Gamma function around the origin x = 0 . They can be derived from the Euler’s reflection formula, Γ ( x ) Γ ( 1 - x ) = π / sin ( π x ) . We obtain Γ ( x 0 ) = x - 1 + Γ ( 1 ) + . From here, recalling that Γ ( 1 ) = - γ , we get Γ ( 1 / b ) = b - γ + , when b . Although this last term of the Taylor expansion, - γ , is negligible we maintain it in expression (30). The only minimum of c b is reached for the solution b = 3 . 1605 of the equation ψ ( 1 / b ) + log b + b - 1 = 0 , where ψ ( · ) is the digamma function (see Figure 1).
Let us now recall two interesting statistical economic models that display a statistical behavior given by distributions of the form (26), that is, the standard Gamma distributions with shape parameter 1 / b ,
f ( z ) d z = 1 Γ ( 1 b ) z 1 b - 1 e - z d z
Figure 1. Normalization constant c b versus b, calculated from Equation (28). The asymptotic behavior is: lim b 0 c b = , and lim b c b = 1 . This last asymptote is represented by the dotted line. The minimum of c b is reached for b = 3 . 1605 , taking the value c b = 0 . 7762 .
Figure 1. Normalization constant c b versus b, calculated from Equation (28). The asymptotic behavior is: lim b 0 c b = , and lim b c b = 1 . This last asymptote is represented by the dotted line. The minimum of c b is reached for b = 3 . 1605 , taking the value c b = 0 . 7762 .
Entropy 11 00959 g001
ECONOMIC MODEL A: The first one is the saving propensity model introduced by Chakraborti and Chakrabarti [4]. In this model a set of N economic agents, having each agent i (with i = 1 , 2 , , N ) an amount of money, u i , exchanges it under random binary ( i , j ) interactions, ( u i , u j ) ( u i , u j ) , by the following the exchange rule:
u i = λ u i + ϵ ( 1 - λ ) ( u i + u j )
u j = λ u j + ϵ ¯ ( 1 - λ ) ( u i + u j )
with ϵ ¯ = ( 1 - ϵ ) , and ϵ a random number in the interval ( 0 , 1 ) . The parameter λ, with 0 < λ < 1 , is fixed, and represents the fraction of money saved before carrying out the transaction. Let us observe that money is conserved, i.e., u i + u j = u i + u j , hence in this model the economy is closed. Defining the parameter n ( λ ) as
n ( λ ) = 1 + 2 λ 1 - λ
and scaling the wealth of the agents as z ¯ = n u / < u > , with < u > representing the average money over the ensemble of agents, it is found that the asymptotic wealth distribution in this system obeys the standard Gamma distribution [5]
f ( z ¯ ) d z ¯ = 1 Γ ( n ) z ¯ n - 1 e - z ¯ d z ¯
The case n = 1 , which means a null saving propensity, λ = 0 , recovers the model of Dragulescu and Yakovenko [3] in which the Gibbs distribution is observed. If we compare Equations (35) and (31), a close relationship between this economic model and the geometrical problem solved in the former section can be established. It is enough to make
n = 1 / b
z ¯ = z
to have two equivalent systems. This means that, from Equation (36), we can calculate b from the saving parameter λ with the formula
b = 1 - λ 1 + 2 λ
As λ takes its values in the interval ( 0 , 1 ) , then the parameter b also runs in the same interval ( 0 , 1 ) . On the other hand, recalling that z = x b / b ϵ , we can get the equivalent variable x from Equation (37),
x = ϵ < u > u 1 / b
where ϵ is a free parameter that determines the mean value of x b in the equivalent geometrical system. Formula (39) means to perform the change of variables u i x i , with i = 1 , 2 , , N , for all the particles/agents of the ensemble. Then, we conjecture that the economic system represented by the generic pair ( λ , u ) , when it is transformed in the geometrical system given by the generic pair ( b , x ) , as indicated by the rules (38) and (39), runs in an equiprobable form on the surface defined by the relationship (12), where the inequality has been transformed in equality. This last detail is due to the fact the economic system is closed, and then it conserves the total money, whose equivalent quantity in the geometrical problem is E. If the economic system were open, with an upper limit in the wealth, then the transformed system would evolve in an equiprobable way over the volume defined by the inequality (12), although its statistical behavior would continue to be the same as it has been proved for the cases b = 1 , 2 in [10,11].
ECONOMIC MODEL B: The second one is a model introduced in [6]. In this model a set of N economic agents, having each agent i (with i = 1 , 2 , , N ) an amount of money, u i , exchanges it under random binary ( i , j ) interactions, ( u i , u j ) ( u i , u j ) , by the following the exchange rule:
u i = u i - Δ u
u j = u j + Δ u
where
Δ u = η ( x i - x j ) ϵ ω x i - [ 1 - η ( x i - x j ) ] ϵ ω x j
with ϵ a continuous uniform random number in the interval ( 0 , 1 ) . When this variable is transformed in a Bernouilli variable, i.e., a discrete uniform random variable taking on the values 0 or 1, we have the model studied by Angle [7], that gives very different asymptotic results. The exchange parameter, ω, represents the maximum fraction of wealth lost by one of the two interacting agents ( 0 < ω < 1 ). Whether the agent who is going to loose part of the money is the i-th or the j-th agent, depends nonlinearly on ( x i - x j ) , and this is decided by the random dichotomous function η ( t ) : η ( t > 0 ) = 1 (with additional probability 1 / 2 ) and η ( t < 0 ) = 0 (with additional probability 1 / 2 ). Hence, when x i > x j , the value η = 1 produces a wealth transfer from agent i to agent j with probability 1 / 2 , and when x i < x j , the value η = 0 produces a wealth transfer from agent j to agent i with probability 1 / 2 . Defining in this case the parameter n ( ω ) as
n ( ω ) = 3 - 2 ω 2 ω
and scaling the wealth of the agents as z ¯ = n u / < u > , with < u > representing the average money over the ensemble of agents, it is found that the asymptotic wealth distribution in this system obeys the standard Gamma distribution [6]
f ( z ¯ ) d z ¯ = 1 Γ ( n ) z ¯ n - 1 e - z ¯ d z ¯
The case n = 1 , which means an exchange parameter ω = 3 / 4 , recovers the model of Dragulescu and Yakovenko [3] in which the Gibbs distribution is observed. If we compare Equations (44) and (31), a close relationship between this economic model and the geometrical problem solved in the last section can be established. It is enough to make
n = 1 / b
z ¯ = z
to have two equivalent systems. This means that, from Equation (45), we can calculate b from the exchange parameter ω with the formula
b = 2 ω 3 - 2 ω
As ω takes its values in the interval ( 0 , 1 ) , then the parameter b runs in the interval ( 0 , 2 ) . It is curious to observe that in this model the interval ω ( 3 / 4 , 1 ) maps on b ( 1 , 2 ) , a fact that does not occur in MODEL A. On the other hand, recalling that z = x b / b ϵ , we can get the equivalent variable x from Equation (46),
x = ϵ < u > u 1 / b
where ϵ is a free parameter that determines the mean value of x b in the equivalent geometrical system. Formula (48) means to perform the change of variables u i x i , with i = 1 , 2 , , N , for all the particles/agents of the ensemble. Then, we conjecture that the economic system represented by the generic pair ( λ , u ) , when it is transformed in the geometrical system given by the generic pair ( b , x ) , as indicated by the rules (47) and (48), runs in an equiprobable form on the surface defined by the relationship (12), where the inequality has been transformed in equality. As explained above, this last detail is due to the fact the economic system is closed, and then it conserves the total money, whose equivalent quantity in the geometrical problem is E. If the economic system were open, with an upper limit in the wealth, then the transformed system would evolve in an equiprobable way over the volume defined by the inequality (12), although its statistical behavior would continue to be the same as it has been proved for the cases b = 1 , 2 in References [10,11].

5. Other Geometrical Questions

We shall proceed now to derive an asymptotic formula ( N ) for the volume of the N-dimensional symmetrical body enclosed by the surface
x 1 b + x 2 b + + x N - 1 b + x N b = E
The linear dimension ρ of this volume, i.e., the length of one of its sides verifies ρ E 1 / b . As argued in Equation (16), the N-dimensional volume, V N ( b , ρ ) , is proportional to the term ρ N and to a coefficient g b ( N ) that depends on N. Thus,
V N ( b , ρ ) = g b ( N ) ρ N
where the characteristic b indicates the particular boundary given by Equation (49).
For instance, from Equation (2), we can write in a formal way:
g b = 1 ( N ) = 1 N 1 Γ ( N 1 + 1 )
From Equation (8), if we take the diameter, ρ = 2 R , as the linear dimension of the N-sphere, we obtain:
g b = 2 ( N ) = π 4 N 2 Γ N 2 + 1
These expressions (51) and (52) suggest a possible general formula for the factor g b ( N ) , let us say
g b ( N ) = a N b Γ N b + 1
where a is a b-dependent constant to be determined. For example, a = 1 for b = 1 and a = π / 4 for b = 2 .
In order to find the dependence of a on the parameter b, the regime N is supposed. Applying Stirling approximation for the factorial ( N b ) ! in the denominator of expression (53), and inserting it in expression (21), it is straightforward to find out the relationship:
c b = ( a b ) - 1 / b
From here and formula (28), we get:
a = Γ 1 b + 1 b
that recovers the exact results for b = 1 , 2 . The behavior of a is monotonous decreasing when b is varied from b = 0 , where a diverges as a 1 / b + , up to the limit b , where a decays asymptotically toward the value a = e - γ = 0 . 5614 .
Hence, the formula for g b ( N ) is obtained:
g b ( N ) = Γ 1 b + 1 N Γ N b + 1
It would be also possible to multiply this last expression (56) by a general polynomial K ( N ) in the variable N, and all the derivation done from Equation (53) would continue to be correct. We omit this possibility in our calculations. For a fixed N, we have that g b ( N ) increases monotonously from g b ( N ) = 0 , for b = 0 , up to g b ( N ) = 1 , in the limit b (see Figure 2). For a fixed b, we have that g b ( N ) decreases monotonously from g b ( N ) = 1 , for N = 1 , up to g b ( N ) = 0 , in the limit N (see Figure 3).
Figure 2. The factor g b ( N ) versus b for N = 10 , 40 , 100 , calculated from Equation (56). Observe that g b ( N ) = 0 for b = 0 , and lim b g b ( N ) = 1 .
Figure 2. The factor g b ( N ) versus b for N = 10 , 40 , 100 , calculated from Equation (56). Observe that g b ( N ) = 0 for b = 0 , and lim b g b ( N ) = 1 .
Entropy 11 00959 g002
Figure 3. The factor g b ( N ) versus N for b = 10 , 40 , 100 , calculated from Equation (56). Observe that g b ( N ) = 1 for N = 1 , and lim N g b ( N ) = 0 .
Figure 3. The factor g b ( N ) versus N for b = 10 , 40 , 100 , calculated from Equation (56). Observe that g b ( N ) = 1 for N = 1 , and lim N g b ( N ) = 0 .
Entropy 11 00959 g003
The final result for the volume of an N-dimensional symmetrical body of characteristic b given by the boundary (49) reads:
V N ( b , ρ ) = Γ 1 b + 1 N Γ N b + 1 ρ N
with ρ E 1 / b .

6. Conclusions

In this work, we have considered a general multi-agent open system verifying an additive constraint. Its statistical behavior has been derived from geometrical arguments. The Maxwellian and the Boltzmann-Gibbs distributions are particular cases of this type of systems. Also, other multi-agent economy models, such as the Dragalescu and Yakovenko’s model [3], the Chakraborti and Chakrabarti’s model [4] and the modified Angle’s model [6], show similar statistical behaviors compared to our general system. This fact fosters a geometrical interpretation of all those models. This geometrical speculation allows us to suggest the equivalence with the Chakraborti and Chakrabarti’s model when the geometrical characteristic b of our model runs in the interval ( 0 , 1 ) . The equivalence with the modified Angle’s model is suggested when b varies in the interval ( 0 , 2 ) . As a particular case of both types of models, the Dragulescu and Yakovenko’s model is obtained for b = 1 .
We have not found in the literature other multi-agent models to establish an equivalence with our system in the range b ( 2 , ) . This point remains an open question and a challenge that will probably trigger other works in this direction.

Acknowledgements

The authors acknowledge some financial support from Spanish DGICYT Projects FIS2006-12781-C02-01.

References

  1. Mantegna, R.; Stanley, H.E. An Introduction to Econophysics: Correlations and Complexity in Finance; Cambridge University Press: Cambridge, UK, 1999. [Google Scholar]
  2. Yakovenko, V.M. Econophysics, Statistical Mechanics Approach to. In Encyclopedia of Complexity and System Science; Meyers, R.A., Ed.; Springer: Berlin, Germany, 2009; pp. 2800–2826. [Google Scholar]
  3. Dragulescu, A.; Yakovenko, V.M. Statistical mechanics of money. Eur. Phys. J. 2000, B17, 723–729. [Google Scholar] [CrossRef]
  4. Chakraborti, A.; Chakrabarti, B.K. Statistical mechanics of money: How saving propensity affects its distribution. Eur. Phys. J. B 2000, 17, 167–170. [Google Scholar] [CrossRef]
  5. Patriarca, M.; Chakraborti, A.; Kaski, K. Statistical model with a standard Gamma distribution. Phys. Rev. E 2004, 70, 016104–016105. [Google Scholar] [CrossRef]
  6. Patriarca, M.; Heinsalu, E.; Chakraborti, A. The ABCD’s of statistical many-agent economy models. 2006; arXiv:physics/0611245. [Google Scholar]
  7. Angle, J. The inequality process as a wealth maximizing process. Physica A 2006, 367, 388–414, (references therein). [Google Scholar] [CrossRef]
  8. Gonzalez-Estevez, J.; Cosenza, M.G.; Lopez-Ruiz, R.; Sanchez, J.R. Pareto and Boltzmann-Gibbs behaviors in a deterministic multi-agent system. Physica A 2008, 387, 4637–4642. [Google Scholar] [CrossRef]
  9. Pellicer-Lostao, C.; Lopez-Ruiz, R. Economic models with chaotic money exchange. In Proceedings of the ICCS 2009, Baton Rouge, LA, USA, 2009; Part I. pp. 43–52.
  10. Lopez-Ruiz, R.; Sañudo, J.; Calbet, X. Geometrical derivation of the Boltzmann factor. Am. J. Phys. 2008, 76, 780–781. [Google Scholar] [CrossRef]
  11. Lopez-Ruiz, R.; Calbet, X. Derivation of the Maxwellian distribution from the microcanonical ensemble. Am. J. Phys. 2007, 75, 752–753. [Google Scholar] [CrossRef]
  12. Harremoës, P. Projections maximizing Tsallis entropy. AIP Conf. Proc. 2007, 965, 90–95. [Google Scholar]

Share and Cite

MDPI and ACS Style

López-Ruiz, R.; Sañudo, J.; Calbet, X. Equiprobability, Entropy, Gamma Distributions and Other Geometrical Questions in Multi-Agent Systems. Entropy 2009, 11, 959-971. https://0-doi-org.brum.beds.ac.uk/10.3390/e11040959

AMA Style

López-Ruiz R, Sañudo J, Calbet X. Equiprobability, Entropy, Gamma Distributions and Other Geometrical Questions in Multi-Agent Systems. Entropy. 2009; 11(4):959-971. https://0-doi-org.brum.beds.ac.uk/10.3390/e11040959

Chicago/Turabian Style

López-Ruiz, Ricardo, Jaime Sañudo, and Xavier Calbet. 2009. "Equiprobability, Entropy, Gamma Distributions and Other Geometrical Questions in Multi-Agent Systems" Entropy 11, no. 4: 959-971. https://0-doi-org.brum.beds.ac.uk/10.3390/e11040959

Article Metrics

Back to TopTop