Multiplicative point process as a model of trading
activity

arXiv:cond-mat/0303089v2 [cond-mat.stat-mech] 27 Dec 2004

V. Gontis and B. Kaulakys1
Institute of Theoretical Physics and Astronomy, Vilnius University
A.Goštauto 12, LT-01108 Vilnius, Lithuania
gontis@ktl.mii.lt
Abstract
Signals consisting of a sequence of pulses show that inherent origin of the 1/f
noise is a Brownian fluctuation of the average interevent time between subsequent
pulses of the pulse sequence. In this paper we generalize the model of interevent
time to reproduce a variety of self-affine time series exhibiting power spectral density
S(f ) scaling as a power of the frequency f . Furthermore, we analyze the relation
between the power-law correlations and the origin of the power-law probability distribution of the signal intensity. We introduce a stochastic multiplicative model for
the time intervals between point events and analyze the statistical properties of the
signal analytically and numerically. Such model system exhibits power-law spectral
density S(f ) ∼ 1/f β for various values of β, including β = 1/2, 1 and 3/2. Explicit
expressions for the power spectra in the low frequency limit and for the distribution
density of the interevent time are obtained. The counting statistics of the events is
analyzed analytically and numerically, as well. The specific interest of our analysis is
related with the financial markets, where long-range correlations of price fluctuations
largely depend on the number of transactions. We analyze the spectral density and
counting statistics of the number of transactions. The model reproduces spectral
properties of the real markets and explains the mechanism of power-law distribution
of trading activity. The study provides evidence that the statistical properties of the
financial markets are enclosed in the statistics of the time interval between trades.
A multiplicative point process serves as a consistent model generating this statistics.
Keywords: Stochastic processes; Econophysics; Financial markets; Point processes; 1/f noise
PACS : 05.40.-a; 02.50.Ey; 89.65.Gh

1

Introduction

Complex collective phenomena usually are responsible for the power-laws which are
universal and independent of the microscopic details of the phenomenon. Examples in
physics are numerous. Power-laws are intrinsic features of the economic and financial
data, as well. The aim of this contribution is to analyze a relationship between the
origin of the power-law distribution and the power-law correlations in financial time
series.
There are numerous studies of power-law probability distributions in various
economic systems [1–15]. The key result in recent findings is that the cumulative
distributions of returns and trading activity can be well described by a power-law
1

Corresponding author. Tel. +370-5-2620859; fax +370-5-2125361.
E-mail address: kaulakys@itpa.lt (B. Kaulakys)

1

asymptotic behavior, characterized by an exponent λ ≈ 3, well outside the stable
Levy regime 0 < λ < 2 [13–15]. Empirical studies confirm the power-law asymptotic behavior with stable exponent for different time scales. Since these empirical
distributions are neither Levy stable nor invariant under addition, new stochastic
models with long range correlations and power-law asymptotic behavior are of great
interest [14].
The random multiplicative process built into the model of wealth distribution
yields Pareto power-law [2]. The generalized Lotka-Volterra dynamics is in the use
for various systems including financial markets [9]. However, these models generically
lead to non-universal exponents and do not explain the power-law autocorrelations
and observable spectral density of financial time series [11]. The time-correlations in
the financial time series are studied extensively as well [13–17]. Recent investigations
[13–16] provide empirical evidence that the long-range correlations for volatility are
due to the trading activity, measured by a number of transactions N . Therefore,
various stochastic models of trading activity or for waiting time between transactions
have to be analysed for reproduction the statistics of volatility in financial markets.
Recently we tried to adapt the model of 1/f noise based on the Brownian motion of time interval between subsequent pulses, proposed in [18–21], to model the
share volume traded in the financial markets [22]. The idea to transfer long time
correlations into the stochastic process of the time interval between trades or time
series of trading activity is in consistence with the detailed studies of the empirical
financial data [13, 16] and can reproduce the spectral properties of the financial time
series [22]. However, the investigation of the model revealed, that the simple additive
Brownian model of time interval between trades failed to reproduce the power-law
probability density function (pdf) of the trading activity [23].
On the other hand, several authors have shown empirically that the fluctuations of
various financial time series possess multifractal statistics [24–28]. Therefore, in this
paper we introduce the stochastic multiplicative model for the time interval between
trades and analyze the statistical properties of the trading activity analytically and
numerically. The model is in consistence with the results of statistical analysis of
the empirical financial time series and reveals the relationship between the power-law
probability distribution and the power-law spectral density of the financial time series.
Moreover, the analytically solvable model [18–21] results only in exact 1/f noise while
in the papers [22, 23] essentially only numerical analysis has been undertaken. Here
we generalize the model defining the time series exhibiting the power spectral density
1/f β with 0.5 . β . 1.5 and the power-law probability density of the trading activity.
Explicit analytical expressions for the power spectra and for pdf of the trading activity
N are obtained and analysed.

2

The model

We consider a signal I(t) as a sequence of the random correlated pulses
X
I(t) =
ak δ(t − tk )

(1)

k

where ak is a contribution to the signal of one pulse at the time moment tk , e.g., a
contribution of one transaction to the financial data. Signal (1) represents a point
process used in a large variety of systems with the flow of point objects or subsequent
actions. When ak = ā is constant, the point process is completely described by

2

the set of times of the events {tk } or equivalently by the set of interevent intervals
{τk = tk+1 − tk }. Various stochastic models of τk can be introduced to define a
stochastic point process. In papers [18–21] it has been shown analytically that the
relatively slow Brownian fluctuations of the interevent time τk yield 1/f fluctuations
of the signal (1).
Power spectral density of the signal I(t) is

S(f ) =

=

Ztf

lim

*

2a2
T

lim

*

kmax
2a2 X
T

T →∞

T →∞

I(t) exp(−i2πf t)dt

ti

kmax
X−k

2+

(2)
+

exp(−i2πf ∆(k; q))

k=kmin q=kmin −k

where the brackets h...i denote the averaging over the realizations of the process,
T = tf − ti is the observation time, ∆(k; q) = tk+q − tk is the time difference between
the pulses occurrence times tk+q and tk , while kmin and kmax are minimal and maximal
values of the index k in the interval of observation T .
We adapt and generalize the model of 1/f noise previously proposed in [18–21]
to model the share volume traded in the financial markets [13, 16]. It is useful to
define a discrete time series with equal time intervals τd . Integration of the signal
I(t) over subsequent time intervals of length τd results in a discrete time series and,
by analogy with financial time series, we will call it the volume Vj ,
Vj =

tZ
j +τd

X

I(t)dt =

ak ,

tj = jτd .

(3)

tj <tk <tj +τd

tj

The number of trades Nj in the time interval [tj , tj + τd ] is defined by the same
equation (3) with ak ≡ 1 and Nj ≡ Vj . Then the power spectral density S(fs ) of the
discrete signal Nj may be calculated by the FFT as
S(fs ) =

*

n
2 X
Nj exp(−i2π(s − 1)(j − 1))
τd n
j=1

2+

.

(4)

Here the discrete frequencies fs = (s − 1)/T , s = 1, 2, ...n and T = τd n. For frequencies fs ≪ 1/τd the power spectral density defined by equation (4) coincides with that
given by equation (2).
In this paper we investigate the statistical properties of the time series Nj , when
the sequence of interevent time τk is generated by a multiplicative stochastic process.
First of all the multiplicativity is an essential feature of processes in economics [1, 7].
Multiplicative stochastic processes yield multifractal intermittency and are able to
produce power-law probability distribution functions. We base our study on the
generic multiplicative process for the interevent times, written as
τk+1 = τk + γτk2µ−1 + τkµ σεk .

(5)

Here the interevent time τk fluctuates due to the external random perturbation by
a sequence of uncorrelated normally distributed random variable {εk } with zero expectation and unit variance, σ denotes the standard deviation of the white noise and
γ ≪ 1 is a damping constant.

3

We will restrict the diffusion of the interevent time according to equation (5) to
the finite interval [τmin , τmax ], i.e., 0 < τmin < τk < τmax . The most simple case
is the pure multiplicativity of τk , i.e., when µ = 1. Other values of µ reproduce
the power-laws as well and explicit expressions can be derived without the loss of
generality.
The iterative relation (5) can be rewritten as Langevine stochastic differential
equation in k-space.
dτk
= γτk2µ−1 + τkµ σξ (k) .
(6)
dk
Here we interpret k as continuous variable while

ξ (k) ξ k′ = δ(k − k′ ).
The steady state solution of the corresponding stationary Fokker-Planck equation
with zero flow gives the probability density function for τk in the k-space (see, e.g.,
[29])
Pk (τk ) = Cτkα , α = 2γ/σ 2 − 2µ
(7)

and C has to be defined from the normalization
Z τmax
Pk (τ )dτ = 1.
τmin

The solution (7) assumes Ito convention involved in the relation between expressions (5), (6) and (7).
As the probability distribution function (7) follows a power-law, we expect that
this would result in the power-law behavior in other statistics, as well. In the next
Section we will obtain the power spectral density defined by equation (2).

3

Power spectral density

The power spectral density is a well-established measure of long time correlations
and is widely used in stochastic systems. As it has been already shown [18–21] the
point process with the Brownian interevent time exhibits 1/f noise.
We will derive the formula for the power spectral density of the multiplicative
stochastic point process model, defined by equations (5) and (6) for the interevent
time.
Let us rewrite equation (2) as [21]
2a2 X
χ∆(k:q) (2πf )
T →∞ T

S(f ) = lim

(8)

k,q

where χ∆(k:q)(2πf ) is the characteristic function for the probability distribution of
∆(k; q). For the normal distribution of ∆(k; q) formula (8) takes the form
2a2 X
2
exp{i2πf h∆(k; q)i − 2π 2 f 2 σ∆
(k; q)}
T →∞ T

S(f ) = lim

(9)

k,q

2 is the variance of ∆(k; q), the time difference between the pulses occurrence
where σ∆
times tk+q and tk . For µ = 1 the time difference ∆(k; q) may be expressed from the
solution of multiplicative stochastic equation (6) as

∆(k; q) =

q
X
l=1

l

X
1
εj }.
τk exp{(γ − σ 2 )l + σ
2
j=1

4

(10)

Averaging over the normal distributions of uncorrelated εj results in explicit ex2 (k; q).
pressions for the mean h∆(k; q)i and variance σ∆
In general, for any µ, the perturbative solution of equation (6 ) yields
h∆(k; q)i = τk q +

γ 2µ−1 2
τ
q + o(γ 2 ),
2 k

(11)

σ2 3
q .
3

2
σ∆
(k; q) = τk2µ

In the low frequency limit f ≪ τk−1 we can replace the summation over k and q
by the integration and take into account only the first order terms of equation (11)
in the expression for the power spectral density (9). This yields
Sµ (f ) =

=

4Ca2
Re
τ
√

τZmax

dτ τ

α

τmin

γ 2µ−1 2
τ
q )}dq,
2



γ
πf



α
3−2µ

Ierf (xmin , xmax ),

(12)



√
α
π
−1
exp −i x −
erfc( −ix)x 3−2µ 2 dx
4

(13)

xZmax

xmin

exp{i2πf (τ q +

0

2Ca2
πτ (3 − 2µ)f

Ierf (xmin , xmax ) = Re

Z∞



where
τ̄ = hτk i =

T
kmax − kmin

is the expectation of τk . Here we introduce the scaled variable x =
xmin =

πf 3−2µ
,
τ
γ min

xmax =

πf 3−2µ
γ τ

and

πf 3−2µ
τ
.
γ max

Expression of the power spectral density (12) is appropriate for the numerical
calculations of the generalized multiplicative point process defined by equations (1)
and (5). In the limit τmin → 0 and τmax → ∞ we obtain an explicit expression for
Sµ (f )
  α
1
α
γ 3−2µ Γ( 2 + 3−2µ )
Ca2
.
(14)
Sµ (f ) = √
πα
)
cos( 2(3−2µ)
πτ (3 − 2µ)f πf
Equation (14) reveals that the multiplicative point process (5) results in the power
spectral density S(f ) ∼ 1/f β with the scaling exponent
β =1+

2γ/σ 2 − 2µ
.
3 − 2µ

(15)

Let us compare our analytical results (12), (13) and (14) with the numerical
calculations of the power spectral density according to equations (2) and (5). In
Fig. 1 we present the numerically calculated power spectral density S(f ) of the signal
I(t) for µ = 1 and α = 2γ/σ 2 − 2 = 0, -0.4 and +0.5. Numerical results confirm
that the multiplicative point process exhibits the power spectral density scaled as
S(f ) ∼ 1/f β . Equation (12) describes the model power spectral density very well in
a wide range of parameters. The explicit formula (14) gives a good approximation of
power spectral density for the parameters when β ≃ 1. In Fig. 2 we present numerical

5

SHfL

SHfL
103

a

b

3

10

102

102

Μ=1.
Σ =0.1
Γ =0.008
Β =0.6

1
10-3

f
-3

10
SHfL

-2

10

Μ=1.
Σ =0.1
Γ =0.01
Β =1.

101

-1

10

1

f
10-2

10-1

1

c
102
101
1

Μ=1.
Σ =0.1
Γ =0.0125
Β =1.5

10-3

10-2

f
10-1

1

Figure 1: The power spectral density S(f ) of the signal I(t) versus frequency f calculated
from the model described by equations (1), (2) and (5) with µ = 1. The sinous curves
represent the results of the numerical simulations averaged over 5 realizations of 100000
impulse sequences with σ = 0.1, the straight lines represent the analytical approximation
given by integral formula (12) with τmin = 10−6 ; τmax = 1 and the dashed lines show results
of explicit formula (14). (a) γ = 0.008, (b) γ = 0.01 and (c) γ = 0.0125.
results of the power spectral density calculated for the multiplicative point process
with µ = 0.5. In this case analytical expressions (12), (13) and (14) describe power
spectral density very well in a wide range of the parameter γ. These results confirm
the earlier finding [18–21] that the power spectral density is related to the probability
distribution of the interevent time τk and 1/f noise occurs when this distribution is
flat, i.e., when α = 0.
It is likely that such a stochastic model with parameters in the region 0.5 ≤ β ≤
1.5 may be adaptable for a wide variety of different systems. In this paper we will
investigate applicability of the model for the financial market.

4

Power-law probability distribution function

The statistics of the trading activity N , i.e., the number of trades per time interval
τd is the major task of this paper. The consistency of any model of financial markets
depends on its ability to reproduce empirically measured probability distribution of
the number of trades N , defined by equation (3). Intuitively, the pdf of N is related
with pdf of the rate ν = 1/τ and with pdf of τ . The pdf of τk in k-space is given by
expression (7). It is obvious that in the actual time t the pdf of τ may be written as
′

Pt (τ ) = Pk (τ )τ /τ̄ = C τ α+1

6

(16)

SHfL

SHfL

104

103

a

103

b

102
Μ=0.5
Σ =0.1
Γ =0.001
Β =0.6

102

10-3
SHfL

10-2

Μ=0.5
Σ =0.1
Γ =0.0025
Β =0.75

101
f
10-1

10-3
SHfL

1

10-2

f
10-1

102

c

1

d

102
101
101
1

Μ=0.5
Σ =0.1
Γ =0.005
Β =1.

1
10-1

Μ=0.5
Σ =0.1
Γ =0.0075
Β =1.25

f
10-3

10-2

10-1

f
10-3

1

10-2

10-1

1

Figure 2: The same as in Fig. 1 but for µ = 0.5. (a) γ = 0.001, (b) γ = 0.0025, (c)
γ = 0.005 and (d) γ = 0.0075.
′

where C = C/τ̄ is a new normalization constant.
For the pure multiplicative model with µ = 1 equations (10) and (11) define a
relationship between N and τ after the substitution h∆(k; q)i → τd , τk → τ and
q → N , i.e.,
γ
(17)
τd = τ N + τ N 2 .
2
This relationship may be used for definition of pdf P (N ) from the relation
P (N )dN = Pt (τ )dτ . We have
P (N ) =

C ′ τd2+α (1 + γN )
.
N 3+α (1 + γ2 N )3+α

(18)

This equation yields limiting cumulative distributions for N ,
+∞

Z
P (N )dN ∼
P> (N ) =
N

1
N 2+α ,
1
N 4+2α ,

N ≪ γ −1 ,
N ≫ γ −1 .

(19)

Let us compare probability distribution function for N obtained from the numerical simulation of the model with the analytical results (18) and (19). It is convenient
to use the amplitude versus rank plot (see, e.g., [30] ) for definition of the exponent
λ of the cumulative distribution P> (N ) ∼ 1/N λ . We plot the value of the n-th rank
variable N (n) as a function of its rank n. If P> (N ) behaves asymptotically when
n → 0 as a power-law, one gets a straight line in log-log coordinates of N and n,
with the slope equal to −1/λ. In Fig. 3 we present the plot for N (n) calculated
numerically for pure multiplicative model (µ = 1) with the parameters yielding the
best fit to the empirical data.

7

N
2000
1000
500
200
n
1

10

100

1000

Figure 3: Amplitude versus rank, N(n), plot defining cumulative distribution P> (N) ∼
1/N λ with the parameters of the model µ = 1, σ = 0.1 and γ = 0.0085. Dashed line
represents asymptotic plot N(n) ∼ n−1/λ with λ = 3.4.

5

Discussion and conclusions

We have introduced a multiplicative stochastic model for the time intervals between
events of point process. Such a model of time series has only a few parameters
defining the statistical properties of the system, i.e., the power-law behavior of the
distribution function and the scaled power spectral density of the signal. The ability
of the model to simulate 1/f noise as well as to reproduce signals with the values
of power spectral density slope β between 0.5 and 1.5 promises a wide variety of
applications of the model.
Let us present shortly the possible interpretations of the empirical data of the
trading activity in the financial markets. With a very natural assumption of transactions in the financial markets as point events we can model the number of transactions
Nj in equal time intervals τd as the outcome of the described multiplicative point
process. We already know from available studies [13] that the empirical data exhibit
power spectral density in the low frequency limit with the slope β ≃ 0.7. For the pure
multiplicative model µ = 1 this corresponds to the case when α = 2γ/σ 2 −2µ ≃ −0.3.
The corresponding cumulative distribution of N in the tail of high values (see equation (19)) has the exponent λ = 4 + 2α = 3.4. This is in an excellent agreement with
the empirical cumulative distribution exponent 3.4 defined in [13] for 1000 stocks of
the three major US stock markets.
The numerical results confirm that the multiplicative stochastic model of the
time interval between trades in the financial market is able to reproduce the main
statistical properties of trading activity N and its power spectral density. The powerlaw exponents of the pdf of the interevent time, α, and the cumulative distribution
of the trading activity, λ, as well as the slope of power spectral density, β, are
defined just by one parameter of the model 2γ/σ 2 . The model suggests a simple
mechanism of the power-law statistics of trading activity in the financial markets.
We expect that multiplicative model of the time interval between the trades with
more specific restrictions for the diffusion and more precisely adjusted parameters
lies in the background of the financial markets statistics and may be useful in the
financial time series analysis.
Acknowledgements
We acknowledge support by the Lithuanian State and Studies Foundation.

8

References
[1] R.N. Mantegna, H.E. Stanley, Nature 376 (1995) 46.
[2] M. Levy, S. Solomon, Int.J.Mod.Phys. C 7 (1996) 595.
[3] M. Levy, S. Solomon, Physica A 242 (1997) 90.
[4] D. Sornette, R. Cont, J. Phys. I (France) 7 (1997) 431.
[5] R.N. Mantegna, H.E. Stanley, J. Stat. Phys. 89 (1997) 469.
[6] B.B. Mandelbrot, Fractals and Scaling in Finance, Springer, Berlin, 1997.
[7] P. Gopikrishnan, M. Meyer, L.A.N. Amaral, H.E. Stanley, Eur. Phys. J. B 3
(1998) 139.
[8] B.B. Mandelbrot, Physica A 263 (1999) 477.
[9] A. Blank, S. Solomon, Physica A 287 (2000) 279.
[10] P. Richmond, Eur. Phys. J. B 20 (2001) 523.
[11] J.P. Bouchaud, e-print: cond-mat/0008103, 2000.
[12] G. Bonanno, F. Lillo, R. N. Mantegna, Physica A 280 (2000) 136.
[13] P. Gopikrishnan, V. Plerou, Y. Liu, L.A.N. Amaral, X. Gabaix, H.E. Stanley,
Physica A 287 (2000) 362.
[14] R.N. Mantegna, H.E. Stanley, An Introduction to Econophysics: Correlations
and Complexity in Finance, Cambridge University Press, Cambridge, 2000.
[15] X. Gabaix, P. Gopikrishnan, V. Plerou, H.E. Stanley, Nature 423 (2003) 267.
[16] P. Gopikrishnan, V. Plerou, X. Gabaix, H.E. Stanley, Phys. Rev. E 62 (2000)
R4493.
[17] P. Cizeau, M. Potters, J.-P. Bouchaud, Quant. Finance 1 (2001) 217.
[18] B. Kaulakys, T. Meškauskas, Phys. Rev. E 58 (1998) 7013.
[19] B. Kaulakys, Phys. Lett. A 257 (1999) 37.
[20] B. Kaulakys, T. Meškauskas, Nonlin. Anal.: Mod. Contr., Vilnius 4 (1999) 87.
[21] B. Kaulakys, Microel. Reliab. 40 (2000) 1787; cond-mat/0305067, 2003.
[22] V. Gontis, Lithuanian J. Phys. 41 (2001) 551; cond-mat/0201514, 2002.
[23] V. Gontis, Nonlin. Anal.: Mod. Contr., Vilnius 7 (2002) 43; cond-mat/0211317,
2002.
[24] F. Schmitt, D. Schertzer, S. Lovejoy, Appl. Stoch. Mod. Data Anal. 15 (1999)
29.
[25] F. Schmitt, D. Schertzer, S. Lovejoy, Int. J. Theor. Appl. Fin. 3 (2000) 361.
[26] N. Vandewalle, M. Ausloos, Eur. Phys. J. B 4 (1998) 257.
[27] P. Richmond, S. Solomon, e-print: cond-mat/0010222, 2000.
[28] B.D. Malamud, D.L. Turcotte, Stat. Planing and Inference, 80 (1999) 173.
[29] C.W. Gardiner, Handbook of Stochastic Methods, Springer, Berlin, 1986.
[30] J.-P. Bouchaud, M. Potters, Theory of Financial Risks, Cambridge University
Press, Cambridge, 2000.

9

