Causalities of the Taiwan Stock Market
Julian Juhi-Lian Ting∗

arXiv:cond-mat/0304132v1 [cond-mat.stat-mech] 6 Apr 2003

No.38, Lane 93, Sec.2, Leou-Chuan E. Rd., Taichung, Taiwan 40312, ROC.
(Dated: February 2, 2008)
Volatility, fitting with first order Landau expansion, stationarity, and causality of the Taiwan
stock market (TAIEX) are investigated based on daily records. Instead of consensuses that consider
stock market index change as a random time series we propose the market change as a dual time
series consists of the index and the corresponding volume. Therefore, causalities between these two
time series are investigated.
PACS numbers: 05.45.Tp,89.90.+n,05.40.-a

I.

INTRODUCTION

Physicists are interested in studying stock markets as complex systems. Almost all questions asked can be summarized as searching for price formation theories. Previous researches shown the distribution of price change has
pronounced tail distribution in contrast to Gaussian distribution expected. Furthermore, the auto-correlation of price
change decays exponentially with a characteristic time scale around 5 minutes. Stock crash or rally have also been
identified by physicists as a kind of herd behaviour[1].
Johansen and Sornette considered fitting most stock markets for the bubbles using Landau expansions of the
index[2]. They showed evidence that market crashes as well as large corrections are preceded by speculative bubbles
with two main characteristics: a power law acceleration of the market price decorated with log-periodic oscillations.
For most markets the log-frequency ω/2π is close to unity. However, most data analysis were done for the index of
more mature market like S & P 500. Not all market indexes are defined in the same way. Could different weighting
methods reach the same conclusion? Will government intervention play a role in the conclusion? Johansen and
Sornette further found emergent markets have larger fluctuations. They extended the expansion up to third order
and successfully predicted Nikkei raise in the year 2000[3]. However, why do these fluctuations exist? The logperiodic oscillation appeared in a wide class of out of equilibrium dynamic systems, like ruptures in heterogeneous
media, historic analysis of earthquakes data, and world population. Canessa tried to establish universality for the
exponents from a renormalization group theory[4], and used a stochastic theory to show that the log-periodicities are
a consequence of transient clusters introduced by an entropy-like term. As entropy in thermodynamics corresponds
to the information in an information theory. The possibility to arrive at the log-periodic oscillation therefore suggests
the log-periodic oscillation is a consequence of information exchange between different species of a large system[5].
The effect of volume is less analyzed in the literature. Volume is a measure of market liquidity while index means
the price. Gopikrishnan et al.[6] analyzed the statistical properties of number of shares traded of a particular stock
at a given time interval from an empirical rule saying that it takes volume to push the index. Bonanno et al.[7] also
analyzed the number of shares traded of selected stocks and find a power spectrum of approximately 1/f .
In the present work we considered the volume effects of Taiwan stock market (TAIEX). In particular, the causeeffect relation between the volume time series and the index time series is analyzed. Taiwan stock market is one of
the largest emerging markets. Johansen and Sornette did not studied it because of availability of trading information.
In the present work we tried to make up this missing piece.
When people talk about the volume involved, sometimes they mean the number of shares traded, sometimes the
amount of money involved. We consider money flow to be more important than share number flow if the whole market
is considered, since the person or institute involved should have fixed amount of money.
II.

BASIC PROPERTIES

Taiwan has two stock markets. TAIEX is the major market[8]. TAIEX definition is taken the yearly average of
1966 as 100. All stocks traded are taken into account.

∗ Electronic

address: jlting@chez.com; URL: http://www.chez.com/jlting

2
We analyzed our daily data of TAIEX from January 3, 1991 to December 30, 2000, which include 2814 trading days
for a ten years period. Intra-day time is treated as continuous. Within our data time, several significant events can
be identified. In August 26, 1997, the market reached a local peak 10116.84. The highest index happened at February
17, 2000 with an index 10202.2. The lowest point is at January 7, 1993 with an index 3135.56.
In order to compare Taiwan stock market with other countries we need to know some basic properties studied in
other markets. The first properties are of course daily index and volume histories. They are plotted in Fig.(1).

10000
9000
8000
7000

TAIEX

6000
5000

4000

3000

920930

940630

960402

971223

991102

971223

991102

date

3500

volume (x 10 8NT$)

3000

2500

2000

1500

1000

500

0
920930

940630

960402

date

FIG. 1: Histories of TAIEX daily variations and the corresponding volumes.

One of the most studied properties is volatility. It is found to be a good measure for long time behaviour and risk[9].
It is also the key input parameter for all option pricing models. Many different quantitative definitions of volatility
are used in the literature. Following Liu et al. we define the price change G(t) as the change in the logarithm of the
index,
Z(t + ∆t) − Z(t)
G(t) ≡ ln Z(t + ∆t) − ln Z(t) ∼
,
=
Z(t)

(1)

where ∆t is the sampling time interval. The absolute value of G(t) describes the amplitude of the fluctuation, as
shown in Fig. (2). The large values of |G(t)| correspond to the crashes and big rallies. None of the |G(t)| ever be
larger than 0.07 because government regulation limited maximum daily changes to 7%.
We define the volatility as the average of |G(t)| over a time window T = n · ∆t, i.e.,
VT (t) ≡

t+n−1
1 X
|G(t′ )| ,
n ′
t =t

(2)

3

0.07

0.06

fluctuation |G(t)|

0.05

0.04

0.03

0.02

0.01

0.00
920930

940630

960402

971223

991102

date

FIG. 2: Daily fluctuations amplitude, |G(t)|.

in which n is the moving window size. Fig. (3) shows the calculated volatility VT (t) for a large averaging window
T = 30 days.

0.035

0.030

volatility

0.025

0.020

0.015

0.010

0.005

0.000
921019

940720

960422

980114

991122

date

FIG. 3: Volatility history over a time window T = 30 days. Dates are labeled at the centers of the average windows.

Fig. (4) shows the probability density function P (VT ) of the volatility for several values of T with ∆t = 1 day. The
10-day data is further fitted with Gaussian distribution in Fig. (5). Apparently the data distribution deviates from
the Gaussian curve in many places. The right hand side heel has much higher probability than the Gaussian fit.
An important property need to be investigated for any time series analysis is its stationarity. Most empirical
analysis assumes the time series to be stationary. A time series is stationary if its mean value and its variance do
not vary systematically over time[10]. For TAIEX we choose a window of 1000 trading days, say, and calculate the
corresponding mean value and variance. The result is plotted in Fig. (6). We find the averages of both index and
volume have clear upward trend, while the variances of index oscillate. The volume average and its variance curves
have a turning point around the year 1995, which marks a fundamental change of Taiwan economics. The index curve
is further fitted with a linear curve. The slope is 1.6, which means if an investor holds a index portfolio for over 1000
days her expectation to profit is 1.6 point/day.
PAnother test for stationarity is based on sample autocorrelation function, ρk = γk /γ0 , in which γk =
(Yt − Ȳ )(Yt+k − Ȳ ) is the covariance at lag k. If ρk decays very fast with increasing k, the time series is stationary. The resulting correlogram is plotted in Fig.(7). Apparently, the index time series is more stochastic than the

4

6 days
10 days
21 days

200
180
160
140

probability

120
100
80
60
40
20
0
-20
0.01175

0.02175

0.03175

0.04175

volatility

FIG. 4: Volatility distribution for various time windows with ∆t = 30 min.

180
160
140

frequency count

120
100
80
60
40
20
0
-20
0.01175

0.02175

0.03175

0.04175

volatility

FIG. 5: Gaussian fit for the volatility distribution for T = 10 days and ∆t = 1 days.

volume time series. The results of Fig.(6) and Fig.(7) are consistent.
III.

CRASHES

In this section, we follow Johansen and Sornette to fit the index changes and obtain future trends prediction.
However, our analysis differs from theirs in labeling the date. They convert all dates into decimal numbers. We
labeled only by trading dates. Our labeling method will be more reliable in dating, since trading days spread
irregularly amount years.
Johansen and Sornette wrote down the first, second, and third order Landau expansion of the index starting or
ending at the historical high tc . The first order expansion reads:
p1 (t) ≈ A + Bτ α + Cτ α cos [ω ln (τ ) + φ] .

(3)

in which τ = |t − tc | is the absolute value of time difference from the critical time of the highest index tc . The
parameter φ correspond to the time unit used, while the parameter A, B,and C are units determined by the index as

5

(a)

index average
volume average

3

8x10

3

7x10

3

TAIEX or Volume

6x10

3

5x10

3

4x10

3

3x10

3

2x10

3

1x10

910901

930113

940528

951010

970221

980706

date

2.4
2.2

(b)

index variance
volume variance

2.0
1.8

variance (x 10 6)

1.6
1.4
1.2
1.0
0.8
0.6
0.4
0.2
0.0
910901

930113

940528

951010

970221

980706

date

FIG. 6: Stationarity test for TAIEX with a window size T = 1000 days for (a) average of the index and the corresponding
volume (b) variance of the index and the corresponding volume. The window is labeled by its central date.

well as the historical period. Only tc , α, and ω are the key parameters: α stands for the power law acceleration, while
ω quantifies for the log-periodic oscillation. Numerically, A is roughly the index at tc .
Our fitting start from January 7, 1993 to February 17, 2000. The result is shown in Fig.(8). The fitting does not
appears good, since it is only a first order approximation while the fitting period is about 7 years. There exits a fitting
solution with tc far beyond the highest index point. However, we cannot take that solution because it is obviously
incorrect. The fitted tc has to come before the highest point.
A theory will only be useful if it can also predict for the future. We used the same equation to fit from February
17,2000 in Fig.(9). We find the future is promising. However, in both cases the ω were found to be much higher than
the previous authors’.
IV.

CAUSALITY TEST

It is well known to the practitioner, as an empirical rule, that volume changes come before price change. A rule
of thumb is that higher index will come after huge volume, and if the volume is reaching a minimum the bottom of
the index will not be far. However, how much ahead? What kinds of correlation do they posses? Why is there a lag?
This kind of dual series with one influencing the other has been analyzed in the econometric literature for GNP (gross
national product)-money supply relation, for example[10, 11]. It is called causality as in the literature of physics. We
have already known our time series are non-stationary. In order for the causality test to be meaningful, the two time
series have to be cointegrated, i.e. their wave-lengths of variation have to be of the same order. Fig.(6)(b) and Fig.(7)

6

1.05
1.00
0.95

index

Correlogram

0.90
0.85
0.80
0.75
0.70

volume

0.65
0.60
0

10

20

30

40

50

lag (k)

FIG. 7: Correlogram plot for both index and the corresponding volume time series.

20000

TAIEX

10000
9000
8000
7000
6000
5000
4000
3000

2000
500

1000

1500

2000

2500

date

FIG. 8: TAIEX fit by first order Landau expansion of the index starting from January 7,1993 to February 17,2000. Fitted
parameters are: tc = 2147.999, α = 1.276, ω = 11.013, φ = 7.037, A = 8250.314, B = −0.49, and C = −0.236.

show their time scale differs but not far.
In the econometrics literature a linear functional form is used for the regression of these two time series. Furthermore,
because the future cannot cause the past, we can write down a functional form for the future index, It and the future
trading volume, Vt , as
P
P
(4)
It = i ai Vt−i + j bj It−j + ut
P
P
(5)
Vt = i li Vt−i + j kj It−j + vt .
We fitted for 90-days from the history using the above functional form. The result is plotted in Fig.(10). Both
the index fitting and the volume fitting show stronger dependence on the index history than on the volume history.
However, only the previous trading day shows significance. Another plot is made for the index difference and the
corresponding volume for the future index difference, It and the future trading volume, Vt , as
P
P
(6)
δIt = i ai Vt−i + j bj δIt−j + ut
P
P
(7)
Vt = i li Vt−i + j kj δIt−j + vt .
The result is plotted in Fig.(11). The figure suggests the index difference time series indeed received slightly higher

7

11000

10000

9000

TAIEX

8000

7000

6000

5000

4000
50

100

150

200

250

date

FIG. 9: TAIEX fit by first order Landau expansion of the index starting from February 17, 2000. Fitted parameters are:
tc = −0.491, α = 1.106, ω = 12.650, φ = 5.558, A = 9776.675, B = −9.079, and C = 2.926.

influence from the previous index difference and the previous volume. However, these effects are pretty small. Furthermore, these peaks will be changed with the truncation of the time series.
V.

CONCLUSION

We tried to find the correlation between the volume of transactions involved and stock market fluctuations beyond
simply considering temporal price series, using Taiwan stock market as an example in the present work.
We find TAIEX volatility spread relatively uniformly in the past ten years, which can also be found in the stationarity
test. The distribution of price change shows fat-heel behaviour.
A slightly higher influence coefficient is found in the index difference fitting with previous index and previous
volume. Johansen and Sornette’s result suggests us to try regression with log-periodic functional forms in the future.
However, it does not seem meaningful to go beyond 90-days for the time lag, since companies will have their seasonal
reports by that time.

[1]
[2]
[3]
[4]
[5]
[6]
[7]
[8]
[9]
[10]
[11]

Victor M. Eguiluz and Martin G. Zimmermann, Phys. Rev. Lett., 85, 5659 (2000).
A. Johansen and D. Sornette, Int. J. Theor. and Applied. Finance, cond-mat/9810071.
A. Johansen and D. Sornette, Int. J. Mod. Phys. C, 10, 563 (1999).
E. Canessa, Euro. Phys. J. B, 20, 591 (2001) cond-mat/0008300.
E. Canessa, J. Phys. A, 33, 9131 (2000) cond-mat/0012031.
P. Gopikrishnan, V. Plerou, X. Gabaix, and H.E. Stanley, Phys. Rev. E, 62, R4493 (2000).
G. Bonanno, F. Lillo, and R. N. Mantegna, cond-mat/9912006 Physica A.
Official site for TAIEX located at http://www.tse.com.tw , see also http://www.sfi.org.tw/newsfi/intdb/menu/MainFrame.asp?Mitem=E
Y. Liu, P. Gopikrishnan, P. Cizeau, M. Meyer, C.-K. Peng and H.E. Stanley, Phys. Rev. E, 60, 1390 (1999).
D. N. Gujarati, Basic Econometrics, 3ed, (McGraw-Hill 1995).
G. S. Maddala, Introduction to Econometrics, 2ed, (Macmillan Publishing Co. 1992).

8

bj
ai

1.2
1.0

index coefficients

0.8
0.6
0.4
0.2
0.0
-0.2
0

20

40

60

80

lag/day

0.8

kj
li

0.6

volume coefficients

0.4
0.2
0.0
-0.2
-0.4
-0.6
-0.8
0

20

40

60

80

lag/day

FIG. 10: Causality coefficients for index and for volume.

9

bj
ai

1

standardized index coefficients

0.3

0.2

79
0.1

0.0

-0.1

-0.2

-0.3
0

20

40

60

80

lag/day

kj
li
0.6

standardized volume coefficients

0.5
0.4
0.3
0.2
0.1
0.0
-0.1
0

20

40

60

80

lag/day

FIG. 11: Causality coefficients for index difference and for the corresponding volume.

