Download 4028-0-Introduction

Survey
yes no Was this document useful for you?
   Thank you for your participation!

* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project

Document related concepts
no text concepts found
Transcript
Get Folder
 Network
Neighbourhood
 Tcd.ie
•Ntserver-usr
•Get
• richmond
Econophysics
Physics and Finance
(IOP UK)
Socio-physics
(GPS)
Molecules > people
Physics World October 2003
http://www.helbing.org/
Complexity
Arises from interaction
Disorder & order
Cooperation & competition
Stochastic Processes
Random movements
Statistical Physics
cooperative phenomena
Describes complex, random behaviour in terms of basic elements and
interactions
Physics and Finance-history

Bankers



Gamblers


Pascal, Bernoulli
Actuaries


Newton
Gauss
Halley
Speculators, investors


Bachelier
Black Scholes >Nobel prize for economics
Books –
Econophysics
• Statistical Mechanics of Financial Markets
• J Voit Springer
• Patterns of Speculation; A study in
Observational Econophysics
• BM Roehner Cambridge
• Introduction to Econophysics
• HE Stanley and R Mantegna Cambridge
• Theory of Financial Risk: From Statistical
Physics to Risk Management
• JP Bouchaud & M Potters Cambridge
• Financial Market Complexity
• Johnson, Jefferies & Minh Hui Oxford
Books – Financial math

Options, Futures & Other
Derivatives
• JC Hull

Mainly concerned with solution of
Black Scholes equation
• Applied math (HPC, DCU, UCD)
Books –
Statistical Physics

Stochastic Processes
•
•
•
•
•
•

Quantum Field Theory (Chapter 3) Zimm Justin
Langevin equations
Fokker Planck equations
Chapman Kolmogorov Schmulochowski
Weiner processes; diffusion
Gaussian & Levy distributions
Random Walks & Transport
• Statistical Dynamics, chapter 12, R Balescu

Topics also discussed in Voit
Read the business press
Financial Times
 Investors Chronicle
 General Business pages

Fundamental & technical analysis
 Web sites

• http://www.digitallook.com/
• http://www.fool.co.uk/
Motivation
Perhaps you want to become an
actuary.
 Or perhaps you want to learn about
investing?

What
happened
next?
5000000
0
20/09/04
20/09/01
20/09/98
20/09/95
20/09/92
20/09/89
20/09/86
20/09/83
23/09/04
23/09/01
20/09/80
23/09/98
23/09/95
20/09/77
23/09/92
23/09/89
20/09/74
23/09/86
23/09/83
20/09/71
23/09/80
10000000
23/09/77
20/09/68
20000000
23/09/74
30000000
23/09/71
23/09/68
DJ Closing Price
12000
10000
8000
6000
Volume of stock traded
4000
25000000
2000
15000000
0
FTSE Closing Price
8000
7000
6000
5000
4000
3000
2000
1000
0
1990-05- 1993-01- 1995-10- 1998-07- 2001-04- 2004-01- 2006-1007
31
28
24
19
14
10
Date
Questions

Can we earn money during both upward and
downward moves?
• Speculators

What statistical laws do changes obey? What is
frequency, smoothness of jumps?
• Investors & physicists/mathematicians


What is risk associated with investment?
What factors determine moves in a market?
• Economists, politicians

Can price changes (booms or crashes) be
predicted?
• Almost everyone….but tough problem!
Why physics?

Statistical physics
• Describes complex behaviour in terms of
basic elements and interaction laws

Complexity
• Arises from interaction
• Disorder & order
• Cooperation & competition
Financial Markets
Elements = agents (investors)
 Interaction laws = forces governing
investment decisions

• (buy sell do nothing)

Trading is increasingly automated using
computers
Social Imitation
Theory of Social Imitation Callen & Shapiro Physics Today July 1974
Profiting from Chaos Tonis Vaga McGraw Hill 1994
buy
Hold
Sell
Are there parallels with statistical physics?
E.g. The Ising model of a magnet
Focus on spin I:
Sees local force field,
Yi, due to other spins
{sj}
plus external field, h
I
yi   J ij s j  h
j
si  sgn[ yi ]
V ( si )   yi si
h
Mean Field theory

Gibbsian statistical mechanics
si  si 
 s p(s )
s  , 
i
 yi / kT
i
e
/ kT
 yi / kT
e
j J ij s j  h
 tanh
kT
yi / kT
e
 y
e i
yi   J ij s j  h
j
sgn[ x]  tanh[ x]
Jij=J>0 Total alignment
(Ferromagnet)

Look for solutions <σi>= σ
σ = tanh[(J σ + h)/kT]
+1
-h/J
y= tanh[(J σ+h)/kT]
σ*>0
y= σ
-1
Orientation as function of h
y= tanh[(J σ+h)/kT] ~sgn [J σ+h]
+1
Increasing h
-1
Spontaneous orientation (h=0)
below T=Tc
Suppose h ~ 0; K  J / kT
1
tanh x  x  x3 / 3! ...
6(1  K )

K
 ~ [Tc  T ]1/ 2 T  Tc
T<Tc
+1
T>Tc
 0 T  Tc
σ*
Increasing T
Social imitation

Herding – large number of agents
coordinate their action
Direct influence between traders
through exchange of information
 Feedback of price changes onto
themselves

Cooperative phenomena
Non-linear complexity
Opinion changes
K Dahmen and J P Sethna Phys Rev B53 1996 14872
J-P Bouchaud Quantitative Finance 1 2001 105


magnets si
field h
 trader’s position φi (+ -?)
 time dependent random
a priori opinion hi(t)
• h>0 – propensity to buy
• h<0 – propensity to sell
• J – connectivity matrix
Confidence?

hi is random variable
 <hi>=h(t); <[hi-h(t)]2>=Δ2
 h(t) represents confidence
• Economy strong: h(t)>0
• expect recession: h(t)<0
• Leads to non zero average for pessimism or
optimism
Need mechanism for
changing mind

Need a dynamics
• Eg G Iori
N
i (t  1)  sgn[hi (t )   J ij j (t )]?
i 1


Basic concepts of stocks and investors
Stochastic dynamics
Topics
• Langevin equations; Fokker Planck equations; Chapman,
Kolmogorov, Schmulochowski; Weiner processes;
diffusion




Bachelier’s model of stock price dynamics
Options
Risk
Empirical and ‘stylised’ facts about stock
data
• Non Gaussian
• Levy distributions

The Minority Game
• or how economists discovered the scientific method!

Some simple agent models
• Booms and crashes

Stock portfolios
• Correlations; taxonomy
Basic material

What is a stock?
• Fundamentals; prices and value;
• Nature of stock data
• Price, returns & volatility
Empirical indicators used by
‘professionals’
 How do investors behave?

Normal v Log-normal
distributions
Probability distribution density
functions p(x)
characterises occurrence of random variable, X
For all values of x:
p(x) is positive

p(x) is normalised, ie: -/0  p(x)dx =1
p(x)x is probability that x < X < x+x
a

b
p(x)dx is probability that x lies between a
and b
Cumulative probability function
C(x) = Probability that x<X
x
= -  P(x)dx
= P<(x)
P>(x) = 1- P<(x)
C() = 1; C(-) = 0
Average and expected
values
For string of values x1, x2…xN
average or expected value of any function
f(x) is

N
1
f  x   Lt  f ( xn)   f ( x) p( x)dx
N  N n 1




f ( x)dC ( x)

In statistics & economics literature, often
find E[ f x ] instead of f x


Moments and the ‘volatility’
m
n > =  p(x) x

<
x
n
Mean: m 1 = m
n
dx
Standard deviation, Root mean square (RMS)
variance or ‘volatility’ :
 2 = < (x-m)2 > =  p(x) (x-m)2dx
= m2 – m 2
NB For mn and hence  to be meaningful, integrals
have to converge and p(x) must decrease
sufficiently rapidly for large values of x.
Gaussian (Normal) distributions
PG(x) ≡ (1/ (2π)½σ) exp(-(x-m)2/22)
All moments exist
For symmetric distribution
m=0; m2n+1= 0
and
m2n = (2n-1)(2n-3)…. 2n
Note for Gaussian: m4=34 =3m22
m4 is ‘kurtosis’
Some other Distributions
Log normal
PLN(x) ≡ (1/(2π)½ xσ) exp(-log2(x/x0)/22)
mn = x0nexp(n22/2)
Cauchy
PC (x) ≡ /{1/(2 +x2)}
Power law tail
(Variance diverges)
Levy distributions
NB Bouchaud uses 
instead of 
Curves that have
narrower peaks
and fatter tails
than Gaussians
are said to exhibit
‘Leptokurtosis’
Simple example

Suppose orders arrive sequentially at random with mean waiting
time of 3 minutes and standard deviation of 2 minutes. Consider
the waiting time for 100 orders to arrive. What is the approximate
probability that this will be greater than 400 minutes?



Assume events are independent.
For large number of events, use central limit theorem to obtain m and .
Thus
• Mean waiting time, m, for 100 events is ~ 100*3 = 300 minutes
• Average standard deviation for 100 events is  ~ 2/100 = 0.2 minutes


Model distribution by Gaussian, p(x) = 1/[(2)½] exp(-[x-m]2/22)
Answer required is

•
•
•
•


P(x>400) = 400 dx p(x) ~ 400 dx 1/((2)½) exp(-x2/22)
= 1/()½ z dy exp(-y2)
where z = 400/0.04*2 ~ 7*10+3
=1/2{ Erfc (7.103)} = ½ {1 – Erf (7.103)}
Information given: 2/ * z dy exp(-y2) = 1-Erf (x)
and tables of functions containing values for Erf(x) and or Erfc(x)