Econophysics: What Can Physicists Contribute To Economics?
Econophysics: What Can Physicists Contribute To Economics?
contribute to economics?
Parameswaran Gopikrishnan*, Luis A. Nunes Amaral*,
Yanhui Liu*, Martin Meyer*, Vasiliki Plerou1"*,
Bernd Rosenow**, and H. Eugene Stanley*.
* Center for Polymer Studies1 and Department of Physics,
Boston University, Boston MA 02215, USA.
t Department of Physics, Boston College, Chestnut Hill, MA 02617, USA.
Institut fur Theoretische Physik, Universitat zu Koln, D-50937 Koln, Germany.
The Center for Polymer Studies is supported by the National Science Foundation.
CP511, Unsolved Problems of Noise and Fluctuations, edited by D. Abbott and L. B. Kish
2000 American Institute of Physics 1-56396-826-6/007$ 17.00
233
Downloaded 13 Jan 2009 to [Link]. Redistribution subject to AIP license or copyright; see [Link]
INTRODUCTION
The analysis of financial data using concepts and methods developed for physical systems has a long tradition [1-4] and has recently attracted the interest of
physicists [5-9]. Possible reasons for this interest include the scientific challenge
of understanding the dynamics of a strongly fluctuating complex system with a
large number of interacting elements. Moreover economic fluctuations could have
many repercussions and understanding fluctuations in physical systems is a topic
where many physicists have contributed. In addition, it is possible that the experience gained by studying fluctuations in physical systems might yield new results
in economics.
One can ask how physicists can contribute to the search for solutions to the puzzles posed by modern economics that economists themselves have not yet solved?
One approachin the spirit of experimental physicsis to begin empirically, with
real data that one can analyze in some detail, but without prior models. In economic systems such as financial markets, one has available a great deal of real data.
Moreover, if one has at one's disposal the tools of statistical physics and the computing power to carry out any number of approaches, this abundance of data is to
great advantage. Thus, for many physicists, studying the economy means studying
a wealth of data on a strongly fluctuating complex system. Indeed, physicists in increasing numbers are finding problems posed by economics sufficiently challenging
to engage their attention [10-26].
Recent studies attempt to uncover and explain the peculiar statistical properties of financial time series such as stock prices, stock market indices or currency
exchange rates. The dynamics of financial markets is difficult to understand not
only because of the complexity of its internal elements but also due to the many
intractable external factors acting on it, which may differ from market to market.
Remarkably, the statistical properties of certain observables appear to be similar
for quite different markets [31], consistent with the possibility that there may exist
"universal" mechanisms.
The most challenging difficulty in the study of financial markets is that the
nature of the interactions between the different elements comprising the system
is unknown, as is the way in which external factors affect it. Therefore, as a
starting point, one may resort to empirical studies to help uncover the regularities or
"empirical laws" that may govern financial markets [27]. The interactions between
the different elements comprising financial markets generate many observables such
as the transaction price, the share volume traded, the trading frequency, and the
values of market indices. Recent empirical studies are based on the analysis of
price fluctuations. This talk reviews recent results on (a) the distribution of stock
price fluctuations and its scaling properties, (b) time-correlations in financial time
series, and (c) correlations among the price fluctuations of different stocks. Space
limitations restrict us to focusing mainly on our group's work; a more balanced
234
Downloaded 13 Jan 2009 to [Link]. Redistribution subject to AIP license or copyright; see [Link]
account can be found in two recent books [5,6], other articles in these proceedings,
and two other recent international conferences [7,9]. Recent work in this field
also focuses on applications such as risk control, derivative pricing, and portfolio
selection [28], which shall not be discussed in this talk. The interested reader should
consult, for example, Refs. [5,6,29].
10
CO
C
o
o
IT)
Q.
$
1Q2
A(i f
(a>:
^UWI^Y^
3 -is
r LJ'
0 v^i^
-15
0 h<^^^^^
Time
Year
FIGURE 1. (a) The daily records of the S&P 500 index for the 35-year period 1962-96 on a
linear-log scale. Note the large jump which occured during the market crash of October 19, 1987.
Sequence of (b) 10 min returns and (c) 1 month returns of the S&P 500 index, normalized to
unit variance, (d) Sequence of i.i.d. Gaussian random variables with unit variance, which was
proposed by Bachelier as a model for stock returns [1]. For all 3 panels, there are 850 events
i.e., in panel (b) 850 minutes and in panel (c) 850 months. Note that, in contrast to (b) and
(c), there are no large events in (d).
235
Downloaded 13 Jan 2009 to [Link]. Redistribution subject to AIP license or copyright; see [Link]
is apparent from Fig. Ib that when one analyzes returns on short time scales, large
events are much more likely to occur, in contrast to a sequence of Gaussian distributed random numbers of the same variance (Fig. Id). As one analyzes returns
on larger time scales, this difference is apparently much less pronounced (Fig. Ic).
In order to understand this process, one starts by analyzing the probability distribution of returns on a given time scale At, which in our study, varies from 1 min
up to a few months.
The nature of the distribution of price fluctuations in financial time series is a
long standing open problem in finance which dates back to the turn of the century.
In 1900, Bachelier proposed the first model for the stochastic process of returns
an uncorrelated random walk with independent, identically Gaussian distributed
(i.i.d) random variables [1]. This model is natural if one considers the return over
a time scale At to be the result of many independent "shocks", which then lead
by the central limit theorem to a Gaussian distribution of returns [1]. However,
empirical studies [4,18,19] show that the distribution of returns has pronounced
tails in striking contrast to that of a Gaussian. Despite this empirical fact, the
Gaussian assumption for the distribution of returns is widely used in theoretical
finance because of the simplifications it provides in analytical calculation; indeed,
it is one of the assumptions used in the classic Black-Scholes option pricing formula
[30].
In his pioneering analysis of cotton prices, Mandelbrot observed that in addition
to being non-Gaussian, the process of returns shows another interesting property:
"time scaling" that is, the distributions of returns for various choices of At,
ranging from 1 day up to 1 month have similar functional forms [4]. Motivated
by (i) pronounced tails, and (ii) a stable functional form for different time scales,
Mandelbrot [4] proposed that the distribution of returns is consistent with a Levy
stable distribution [2,3].
Conclusive results on the distribution of returns are difficult to obtain, and require a large amount of data to study the rare events that give rise to the tails.
More recently, the availability of high frequency data on financial market indices,
and the advent of improved computing capabilities, has facilitated the probing of
the asymptotic behavior of the distribution. For example, Mantegna and Stanley [18] analyzed approximately 1 million records of the S&P 500 index. They
report that the central part of the distribution of S&P 500 returns appears to be
well fit by a Levy distribution, but the asymptotic behavior of the distribution of
returns shows faster decay than predicted by a Levy distribution. Hence, Ref. [18]
proposed a truncated Levy distributiona Levy distribution in the central part followed by an approximately exponential truncationas a model for the distribution
of returns. The exponential truncation ensures the existence of a finite second moment, and hence the truncated Levy distribution is not a stable distribution [32,33].
The truncated Levy process with i.i.d. random variables has slow convergence to
Gaussian behavior due to the Levy distribution in the center, which could explain
the observed time scaling for a considerable range of time scales [18].
Recent studies [34,35] on considerably larger time series using larger databases
236
Downloaded 13 Jan 2009 to [Link]. Redistribution subject to AIP license or copyright; see [Link]
show quite different asymptotic behavior for the distribution of returns. Our recent
work [34] analyzed three different data bases covering securities from the three
major US stock markets. In total, we analyzed approximately 40 million records of
stock prices sampled at 5 min intervals for the 1000 leading US stocks for the 2-year
period 1994-95 and 35 million daily records for 16,000 US stocks for the 35-year
period 1962-96. We study the probability distribution of returns (Fig. 2(a,b,c)) for
individual stocks over a time interval At, where At varies approximately over a
10U
c
~ 10~1
Livy Regime ,
a; 17
Postive taila'3
o Negative tail
(a)
1
10
100
10
100
10
100
Normalized returns
FIGURE 2. (a) Log-log plot of the cumulative distribution of the normalized 1 min returns for
the S&P 500 index. Power-law regression fits in the region 3 < g < 50 yield a = 2.95 0.07
(positive tail), and a = 2.75 0.13 (negative tail). For the region 0.5 < g < 3, regression
fits give a = 1.6 0.1 (positive tail), and a = 1.7 0.1 (negative tail), (b) Log-log plot of
the cumulative distribution of normalized returns of the S&P 500 index. The positive tails are
shown for At = 16, 32,128, 512 mins. Power-law regression fits yield estimates of the asymptotic
power-law exponent a = 2.69 0.04, a = 2.53 0.06, a = 2.83 0.18 and a = 3.39 0.03 for
At = 16,32,128 and 512 mins, respectively, (c) The positive and negative tails of the cumulative
distribution of the normalized returns of the 1000 largest companies in the TAQ database for the
2-year period 1994-1995. The solid line is a power-law regression fit in the region 2 < x < 80.
237
Downloaded 13 Jan 2009 to [Link]. Redistribution subject to AIP license or copyright; see [Link]
238
Downloaded 13 Jan 2009 to [Link]. Redistribution subject to AIP license or copyright; see [Link]
(b)
10 .
one day
(a)
10
A A A A ^ A A A A A A A A .
one day
a2=0.93
|io-
O
0)
Q.
2 io~4
10
P2=0.90
210'
0)
I ,
Q. 10
10
10
10
10
10
10
10
102
103
104
105
t [min]
/[min' ]
FIGURE 3. Plot of (a) the power spectrum S(f) and (b) the detrended fluctuation analysis
F(t) of the absolute values of returns #(), after detrending the daily pattern [40,41] with the
sampling time interval A = 1 min. The lines show the best power law fits (R values are better
than 0.99) above and below the crossover frequency of /x = (1/570) min"1 in (a) and of the
crossover time, tx = 600min in (b). The triangles show the power spectrum and DFA results for
the "control", i.e., shuffled data.
volatility correlations show asymptotic I// behavior [40-42]. Using the same data
bases as above, Liu and his collaborators also study the cumulative distribution of
volatility [40,42] and find that it is consistent with a power-law asymptotic behavior,
characterized by an exponent JJL w 3, just the same as that for the distribution of
returns. For individual companies also, one finds a similar power law asymptotic
behavior [41]. In addition, it is also found that the volatility distribution scales for
a range of time intervals just as the distribution of returns.
239
Downloaded 13 Jan 2009 to [Link]. Redistribution subject to AIP license or copyright; see [Link]
240
Downloaded 13 Jan 2009 to [Link]. Redistribution subject to AIP license or copyright; see [Link]
regime, and yet preserves its shapescalesfor a range of time scales and (ii) the
long range correlations in the amplitude of price fluctuations. How are the two
related ?
Previous explanations of scaling relied on Levy stable [4] and exponentiallytruncated Levy processes [6,18]. However, the empirical data that we analyze
are not consistent with either of these two processes. In order to confirm that
the scaling is not due to a stable distribution, one can randomize the time series
of 1 min returns, thereby creating a new time series which contains statisticallyindependent returns. By adding up n consecutive returns of the shuffled series, one
can construct the nmin returns. Both the distribution and its moments show a
rapid convergence to Gaussian behavior with increasing n, showing that the time
dependencies, specifically volatility correlations are intimately connected to the
observed scaling behavior [34].
Using the statistical properties summarized above, can we attempt to deduce a
statistical description of the process which gives rise to this output? For example,
the standard ARCH model [27,44] reproduces the power-law distribution of returns;
however it assumes finite memory on past events and hence is not consistent with
long-range correlations in volatility. On the other hand, the distribution of volatility and that of returns which have similar asymptotic behavior, however support
the central ARCH hypothesis that g(t) = tv(t), where e is an i.i.d. Gaussian
random variable independent of the volatility v ( t ] , and g(t) denotes the returns.
A consistent statistical description may involve extending the traditional ARCH
model to include long-range volatility correlations [45].
A more fundamental question would be to understand the above results starting
from a microscopic setting. Researchers have also studied microscopic models that
might give rise to the empirically observed statistical properties of returns [5,10].
For example, Lux and Marchesi [10] recently simulated a microscopic model of
financial markets with two types of traders, what they refer to as 'fundamentalist'
and 'noise' traders. Their results reproduce the power-law tail for the distribution
of returns and also the long range correlations in volatility.
In the last section, we found evidence for different modes of correlations between
different companies. For example, the largest eigenvalue of the cross-correlation
matrix showed correlations that pervade the entire market. Could it be that the
above observed scaling properties are related to how correlations propogate from
one unit to the other such as occur in critical phenomena? Researchers have studied
economic data from the physics perspective of a complex system with each unit
depending on the other. Specifically, the possibility that all the companies in a
given economy might interact, more or less, like a spin glass. In a spin glass, each
spin interacts with every other spinbut not with the same coupling and not even
with the same sign. For example, if the stock price of a given business firm A
decrease by, e.g., 10% , this will have an impact in the economy. Some of these will
be favorablefirm B, which competes with A, may experience an increase in market
share. Others will be negativeservice industries that provide personal services for
firm A employees may experience a drop-off in sales as employee salaries will surely
241
Downloaded 13 Jan 2009 to [Link]. Redistribution subject to AIP license or copyright; see [Link]
decline. There must be positive and negative correlations for almost any economic
change. Can we view the economy as a complicated spin glass?
To approch this problem, M. H. R. Stanley and M. A. Salinger first located and
secured a databasecalled COMPUSTATthat lists the annual sales of every firm
in the United States. With this information, M. H. R. Stanley and co-workers calculated histograms of how firm sizes change from one year to the next [52]. They
find that the distribution of growth rates of firm sales has the same functional
proximately 1/6. Recently, similar statistical properties were found for the GDP of
countries [53] and for university research fundings [54]. Hence, it is not impossible
to imagine that there are some very general principles of complex organizations
at work here, because similar empirical laws appear to hold for data on a range of
systems that at first sight might not seem to be so closely related. Buldyrev models
this firm structure as an approximate Cayley tree, in which each subunit of a firm
reacts to its directives from above with a certain probability distribution [55]. More
recently, Amaral et al. [56] have proposed a microscopic model that reproduces both
the exponent and the distribution function. Takayasu and Okuyama [57] extended
Buldyrev, D. Canning, P. Cizeau, X. Gabaix, S. Havlin, P. Ch. Ivanov, R. N. Mantegna, C.-K. Peng, M. A. Salinger, and M. H. R. Stanley. We also thank M.
Barthelemy, J.-P. Bouchaud, D. Sornette, D. Stauffer, S. Solomon, and J. Voit for
helpful dicussions and comments.
REFERENCES
1. L. Bachelier, Ann. Sci. Ecole Norm. Sup. 3, 21 (1900).
2. V. Pareto, Cours d'Economie Politique, Lausanne and Paris, 1897.
3. P. Levy, Theorie de I'Addition des Variables Aleatoires, Gauthier-Villars, Paris,
1937.
4. B. B. Mandelbrot, J. Business 36, 294 (1963).
5. J. P. Bouchaud and M. Potters, Theorie des Risques Financieres, Alea-Saclay, Eyrolles, 1998.
6. R. N. Mantegna and H. E. Stanley, An Introduction to Econophysics: Correlations
and Complexity in Finance, Cambridge University Press, Cambridge 1999.
7. I. Kondor and J. Kertesz (eds) Econophysics: An Emerging Science, Kluwer, Dordrecht, 1999.
8. R. N. Mantegna (ed), Proceedings of the International Workshop on Econophysics
Statistical Finance, Physica A [special issue] 269, 1 (1999).
242
Downloaded 13 Jan 2009 to [Link]. Redistribution subject to AIP license or copyright; see [Link]
243
Downloaded 13 Jan 2009 to [Link]. Redistribution subject to AIP license or copyright; see [Link]
(1997).
41. Y. Liu, P. Gopikrishnan, P. Cizeau, M. Meyer, C.-K. Peng and H. E. Stanley, Phys.
Rev. 60, xxx (1999), cond-mat/9903369.
42. P. Cizeau, Y. Liu, M. Meyer, C.-K. Peng, and H. E. Stanley, Physica A 245, 441
(1997).
43.
44.
45.
46.
47.
48.
49.
50.
J. Beran, Statistics for Long-Memory Processes (Chapman & Hall, NY, 1994).
R. F. Engle, Econometrica 50, 987 (1982).
R. Bailie, J. Econometrics 17, 223 (1991) .
S. Galluccio, J.-P. Bouchaud, and M. Potters, Physica A 259, 449 (1998); L. Laloux,
P. Cizeau, J.-P. Bouchaud and M. Potters, cond-mat/9810255; Risk 3, 69 (March,
1999).
V. Plerou, P. Gopikrishnan, B. Rosenow, L. A. N. Amaral, and H. E. Stanley, condmat/9903369.
E. P. Wigner, Ann. Math. 53, 36 (1951); E. P. Wigner, in Conference on Neutron
Physics by Time-of-flight (Gatlinburg, Tennessee, 1956).
F. J. Dyson, J. Math. Phys. 3, 140 (1962); F. J. Dyson and M. L. Mehta, J. Math.
Phys. 4, 701 (1963); M. L. Mehta and F. J. Dyson, J. Math. Phys. 4, 713 (1963).
T. Guhr, A. Miiller-Groeling, and H. A. Weidenmiiller, Phys. Rep. 299, 190 (1998);
M. L. Mehta, Random Matrices, Academic Press, Boston, 1991.
244
Downloaded 13 Jan 2009 to [Link]. Redistribution subject to AIP license or copyright; see [Link]
245
Downloaded 13 Jan 2009 to [Link]. Redistribution subject to AIP license or copyright; see [Link]