# The Chaotic Attractor Analysis of DJIA Based on Manifold Embedding and Laplacian Eigenmaps

Authors: Xiaohua Song; Dongxiao Niu; Yulin Zhang
Journal: Mathematical Problems in Engineering
 (2016)

---

## Abstract

By using the techniques of Manifold Embedding and Laplacian Eigenmaps, a novel strategy has been proposed in this paper to detect the chaos of Dow Jones Industrial Average. Firstly, the chaotic attractor of financial time series is assumed to lie on a low-dimensional manifold that is embedded into a high-dimensional Euclidean space. Then, an improved phase space reconstruction method and a nonlinear dimensionality reduction method are introduced to help reveal the structure of the chaotic attractor. Next, the empirical study on the financial time series of Dow Jones Industrial Average shows that there exists an attractor which lies on a manifold constructed by the time sequence of Moving average convergence divergence; finally, Determinism Test, Poincaré section, and translation analysis are used as test approaches to prove both whether it is a chaos and how it works.


---

## Body

## 1. Introduction

Chaos is a concept defined as a deterministic dynamical system that is sensitive to initial conditions and gives rise to an unpredictable behavior in the long term. Although it seems like a paradox, this concept of chaos has been widely accepted in fields such as economics and finance for almost three decades, since the first chaotic phenomenon was reported in economics [1]. Chaos theory promotes the search for a mechanism that allows economists to interpret the observed movements in real economic and financial data and helps to distinguish whether it is exogenous or endogenous. In this sense, it represents a deeper thinking about the essential characteristics of the evolutionary process of an economic and finance system, namely, whether the fluctuation of a system is caused by stochastic phenomena or effected by strong nonlinearities. If the outlier events are caused by some sort of random walk which is unpredictable, we may just neglect them in empirical study. However, if they are caused by internal factors of the system, their characteristics are informative for revealing their generating mechanism, which will be an important issue for forecasting the system.Chaos and its implications for system-forecasting have been enthusiastically and hotly debated in financial markets for a long time [2]. Takens’ theorem has proved that a chaotic dynamical system can be reconstructed from a sequence of observations of the state of a dynamical system [3], thus making financial time series a potentially simple and convenient measurement for chaos detecting. Despite this, it remains difficult to detect an obviously chaotic behavior in financial time series. This is because, from an empirical point of view, it is difficult to distinguish which part of the sequence is provoked by randomness and which part is determined by nonlinearity [4].In view of this, some indirect methods have been introduced to help investigate the natural characteristics of financial chaotic phenomena. For example, the correlation dimension method is a straightforward and fast implement to calculate the fractal dimension of a possibly underlying strange attractor [5, 6]; the Largest Lyapunov Exponents (LLE) proposed by Wolf et al. and Rosenestein et al. is a more appropriate tool to detect chaotic phenomena [7]. Besides, based on the correlation integral, Brock et al. developed a statistical test for independence known as the BDS test to test the serial dependence and nonlinear structures of financial time series. Notably, unlike the correlation dimension method, it assumes the delay time equals one [8].The above described applications for chaos detecting require large, stable, and clean sample data. However, the empirical data from financial markets are often small and noisy, which make it hard to detect the chaos. Thus, to overcome this limit, some topological methods have been introduced to detect the financial chaotic behaviors. These methods are characterized by studying the organization of the strange attractor [9] and exploiting an essential property of chaotic systems, through which the relationship between the empirical data that are not possible to be discovered in the original time sequence can be revealed [10].Although studies on chaotic financial systems are fruitful [11–13], there are still various disputes. Among them, two most controversial views are as follows: Could the financial system be characterized by a low-dimensional chaos [14]? And if there is a chaotic phenomenon, is it able to prove that it is generated from a deterministic system [4]?In this paper, Dow Jones Industrial Average (Dow Jones Industrial Average, US) is analyzed as representative financial time series. Before the studying, a specific topological structure of financial system has been assumed, in which the complicated high-dimensional phase space of financial system can be represented by a few independent degrees of freedom embedded in a low-dimensional nonlinear manifold. Under this assumption, the rest of this paper is organized as follows. In Section2, the theorems of time series reconstruction and dimensionality reduction are probed. Firstly, an improved embedding theorem is adopted to construct a compact but redundant phase space of financial system. Then, Laplacian Eigenmaps is used to denoise and map this space to a low-dimensional manifold, for the sake of extracting the strange attractor hidden behind the complicated chaotic system. In Section 3, empirical data selection and processing are discussed, and the time series of Moving average convergence divergence is applied to substitute the traditional price series to reconstruct the chaotic attractor of financial system. In addition, the wavelet denoising tool is also introduced to smooth the noise of MACD sequence. In Section 4, empirical data of financial market DJIA are studied to find out whether it truly has a chaotic attractor, by virtue of some measures such as Determinism Test and Poincaré section and then translation analysis of chaotic attractor will perform to help understand the internal mechanisms of chaos generation. Finally, some conclusions are drawn in Section 5.


## 2. Phase Space Reconstruction and Manifold Dimensionality Reduction

### 2.1. Phase Space Reconstruction

For a dynamical system, phase space volumes are contracted by the time evolution. The trajectory of such a system typically settles on a subset ofR
n, which is called an attractor [15]. However, in practice, it is usually impossible to measure all the components of the n dimensional vector space. Fortunately, the embedding theorem proved the following.IfM is a d dimensional manifold, φ
:
M
→
M, φ is a smooth diffeomorphism, and y is a twice-differentiable function, let y
:
M
→
R, ∅
(
φ
,
y
)
:
M
→
R
2
d
+
1 where ∅
(
φ
,
y
)
=
(
y
(
x
)
,
y
(
φ
(
x
)
)
,
y
(
φ
2
(
x
)
)
,
…
,
y
(
φ
2
n
(
x
)
)
); then ∅
(
φ
,
y
) is an embedding from M to R
2
d
+
1.Thus, the phase space of dynamical system may be reconstructed by one of the embedding features ofX
⃑ and can be formalized as(1)
X
⃑
t
i
=
X
t
i
,
X
t
i
+
l
Δ
t
,
…
,
X
t
i
+
d
-
1
l
Δ
t
,where d denotes the embedding dimension, Δ
t is the sample time, l is an appropriate integer, and the delay time is τ
d
=
l
Δ
t. A diffeomorphism between the reconstructed and original phase spaces exists if(2)
D
≥
2
d
+
1
,where d is the dimension of the compact manifold containing the attractor. If an appropriate parameter of delay time τ
d and embedding dimension d [16] are selected, the original dynamical system can be characterized by the reconstructed phase space.In this paper, different from traditional methods of phase space reconstruction, the time delayτ
d is fixed to 1 so that a minimum correlation dimension of embedding space with respect to τ
d
=
1 can be decided by the G-P algorithm.The formula follows:(3)
C
m
,
N
,
r
,
τ
d
=
2
M
M
-
1
∑
1
≤
i
≤
j
≤
M
Θ
r
-
x
i
-
x
j
,
r
>
0
,where Θ
α
=
0 if α
<
0, or Θ
α
=
1; N is the size of the data set, M
=
N
-
m
+
1 denotes the number of embedded points in m-dimensional space, and · represents the sup-norm.C
m
,
N
,
r
,
1
≡
C
m
,
N
,
r measures the fraction of the pairs of points x
i, whose sup-norm separation is not bigger than r. If the limit of C
m
,
N
,
r as N
→
∞ exists for each r, then the fraction of all state vector points within r can be denoted by C
m
,
r
=
lim
N
→
∞
⁡
C
m
,
N
,
r. In particular, when r
→
0, the correlation dimension is defined as(4)
C
m
=
lim
r
→
0
⁡
log
⁡
C
m
,
r
log
⁡
r
.Thus, the G-P algorithm is simplified to a single parameter estimation with respect tom. However, the simplification may cause another problem; that is, although the reconstruction of the phase space can preserve the properties of the original dynamical system, there will be a significant redundancy within the reconstructed phase space because of τ
d
=
1. This makes it hard to reveal the true structure of the chaotic attractor embedded in the phase space. To overcome this limit, some ideas of nonlinear dimensionality reduction [17] can be introduced. One approach is to assume the chaotic attractor of the dynamical system lies on an embedded nonlinear manifold within a higher-dimensional Euclidean space.


### 2.2. Manifold Dimensionality Reduction

In this paper, the manifold dimensionality reduction methodology called Laplacian Eigenmaps [18] is used to extract the chaotic attractor from the redundant reconstructed phase space of financial system. Suppose there is a Euclidean space R
D, and the Laplacian Eigenmaps is defined as L
→
Δ
D
-
W, where W is the adjacency matrix of edge weights and D is the diagonal matrix with D
i
i
=
∑
j
ω
i
j. The eigenvalues and eigenvectors of the Laplacian reveal a wealth of information about the phase space R
D. The algorithmic procedure can be illustrated as follows.Firstly, construct the adjacency distance betweenχ
i and χ
j, where χ
i
,
χ
j
∈
R
D by distinguishing whether χ
i
-
χ
j
2
<
ε
(
ε
∈
R
), where the norm is the usual Euclidean norm in R
D.Then, determine the weights of the adjacency matrixW by heat kernel, in which the weight is set as w
i
j
=
e
-
χ
i
-
χ
j
2
/
i if there is an edge between node i and j.Finally, compute Eigenmaps under the assumption that the weighted graphG
=
(
V
,
E
) is connected, and map it into a lower d dimensional space, where y
1
,
y
2
,
…
,
y
k
∈
R
d. This embedding is given by the k
×
d matrix f
=
f
1
,
f
2
,
…
,
f
d where the ith row provides the embedding coordinates of the ith vertex. Thus, the objective function can be formalized as(5)
∑
i
,
j
f
i
-
f
j
ω
i
j
=
t
r
f
T
L
f
=
L
f
,where f
(
i
)
=
[
f
1
(
i
)
,
f
2
(
i
)
,
…
,
f
d
(
i
)
]
T is the d dimensional representation of the ith vertex. In this way, this algorithm can be reduced to finding a
r
g
m
i
n
f
T
L
f
=
1
L
f
.The entire process of this section can be shown in Figure1 (the Lorenz system as an example).Figure 1
Demonstration of the entire process of phase space reconstruction and manifold dimensionality reduction. Note: a one-dimensional time series (a) is firstly mapped to aD-dimensional Euclidean space according to embedding theorem, where we use phase space reconstruction (τ
d
=
1) so that a compact but redundant dynamic space will be built (b). Then, NLDR such as Laplacian Eigenmaps is used to reduce the phase space to ad-dimensional manifold (D
≥
2
d
+
1), which helps to eliminate the redundancy and illustrate the chaotic attractor of the original dynamical system (c).


## 2.1. Phase Space Reconstruction

For a dynamical system, phase space volumes are contracted by the time evolution. The trajectory of such a system typically settles on a subset ofR
n, which is called an attractor [15]. However, in practice, it is usually impossible to measure all the components of the n dimensional vector space. Fortunately, the embedding theorem proved the following.IfM is a d dimensional manifold, φ
:
M
→
M, φ is a smooth diffeomorphism, and y is a twice-differentiable function, let y
:
M
→
R, ∅
(
φ
,
y
)
:
M
→
R
2
d
+
1 where ∅
(
φ
,
y
)
=
(
y
(
x
)
,
y
(
φ
(
x
)
)
,
y
(
φ
2
(
x
)
)
,
…
,
y
(
φ
2
n
(
x
)
)
); then ∅
(
φ
,
y
) is an embedding from M to R
2
d
+
1.Thus, the phase space of dynamical system may be reconstructed by one of the embedding features ofX
⃑ and can be formalized as(1)
X
⃑
t
i
=
X
t
i
,
X
t
i
+
l
Δ
t
,
…
,
X
t
i
+
d
-
1
l
Δ
t
,where d denotes the embedding dimension, Δ
t is the sample time, l is an appropriate integer, and the delay time is τ
d
=
l
Δ
t. A diffeomorphism between the reconstructed and original phase spaces exists if(2)
D
≥
2
d
+
1
,where d is the dimension of the compact manifold containing the attractor. If an appropriate parameter of delay time τ
d and embedding dimension d [16] are selected, the original dynamical system can be characterized by the reconstructed phase space.In this paper, different from traditional methods of phase space reconstruction, the time delayτ
d is fixed to 1 so that a minimum correlation dimension of embedding space with respect to τ
d
=
1 can be decided by the G-P algorithm.The formula follows:(3)
C
m
,
N
,
r
,
τ
d
=
2
M
M
-
1
∑
1
≤
i
≤
j
≤
M
Θ
r
-
x
i
-
x
j
,
r
>
0
,where Θ
α
=
0 if α
<
0, or Θ
α
=
1; N is the size of the data set, M
=
N
-
m
+
1 denotes the number of embedded points in m-dimensional space, and · represents the sup-norm.C
m
,
N
,
r
,
1
≡
C
m
,
N
,
r measures the fraction of the pairs of points x
i, whose sup-norm separation is not bigger than r. If the limit of C
m
,
N
,
r as N
→
∞ exists for each r, then the fraction of all state vector points within r can be denoted by C
m
,
r
=
lim
N
→
∞
⁡
C
m
,
N
,
r. In particular, when r
→
0, the correlation dimension is defined as(4)
C
m
=
lim
r
→
0
⁡
log
⁡
C
m
,
r
log
⁡
r
.Thus, the G-P algorithm is simplified to a single parameter estimation with respect tom. However, the simplification may cause another problem; that is, although the reconstruction of the phase space can preserve the properties of the original dynamical system, there will be a significant redundancy within the reconstructed phase space because of τ
d
=
1. This makes it hard to reveal the true structure of the chaotic attractor embedded in the phase space. To overcome this limit, some ideas of nonlinear dimensionality reduction [17] can be introduced. One approach is to assume the chaotic attractor of the dynamical system lies on an embedded nonlinear manifold within a higher-dimensional Euclidean space.


## 2.2. Manifold Dimensionality Reduction

In this paper, the manifold dimensionality reduction methodology called Laplacian Eigenmaps [18] is used to extract the chaotic attractor from the redundant reconstructed phase space of financial system. Suppose there is a Euclidean space R
D, and the Laplacian Eigenmaps is defined as L
→
Δ
D
-
W, where W is the adjacency matrix of edge weights and D is the diagonal matrix with D
i
i
=
∑
j
ω
i
j. The eigenvalues and eigenvectors of the Laplacian reveal a wealth of information about the phase space R
D. The algorithmic procedure can be illustrated as follows.Firstly, construct the adjacency distance betweenχ
i and χ
j, where χ
i
,
χ
j
∈
R
D by distinguishing whether χ
i
-
χ
j
2
<
ε
(
ε
∈
R
), where the norm is the usual Euclidean norm in R
D.Then, determine the weights of the adjacency matrixW by heat kernel, in which the weight is set as w
i
j
=
e
-
χ
i
-
χ
j
2
/
i if there is an edge between node i and j.Finally, compute Eigenmaps under the assumption that the weighted graphG
=
(
V
,
E
) is connected, and map it into a lower d dimensional space, where y
1
,
y
2
,
…
,
y
k
∈
R
d. This embedding is given by the k
×
d matrix f
=
f
1
,
f
2
,
…
,
f
d where the ith row provides the embedding coordinates of the ith vertex. Thus, the objective function can be formalized as(5)
∑
i
,
j
f
i
-
f
j
ω
i
j
=
t
r
f
T
L
f
=
L
f
,where f
(
i
)
=
[
f
1
(
i
)
,
f
2
(
i
)
,
…
,
f
d
(
i
)
]
T is the d dimensional representation of the ith vertex. In this way, this algorithm can be reduced to finding a
r
g
m
i
n
f
T
L
f
=
1
L
f
.The entire process of this section can be shown in Figure1 (the Lorenz system as an example).Figure 1
Demonstration of the entire process of phase space reconstruction and manifold dimensionality reduction. Note: a one-dimensional time series (a) is firstly mapped to aD-dimensional Euclidean space according to embedding theorem, where we use phase space reconstruction (τ
d
=
1) so that a compact but redundant dynamic space will be built (b). Then, NLDR such as Laplacian Eigenmaps is used to reduce the phase space to ad-dimensional manifold (D
≥
2
d
+
1), which helps to eliminate the redundancy and illustrate the chaotic attractor of the original dynamical system (c).


## 3. Empirical Data Selection and Processing

Although studies on nonlinearities and chaos of financial time series have lasted for almost two decades [19], there is still no convincing support for considering chaos as a consequence of some finite length and high noise level financial sequences. This is because a financial time series, such as stock price, is usually not cyclic [20], and it is difficult to find out a specific oscillation pattern for financial system by observing the movement of financial prices directly.In view of this, a time series that can characterize the periodic or quasi-periodic property of financial system must be discovered. In this paper, the famous trading indicator used in technical analysis of financial market, known as MACD (short for Moving average convergence divergence), is adopted as a new time series response of financial system, since the MACD [21] is supposed to reveal changes in the strength, direction, momentum, and duration of the trend of a financial price.Using the MACD sequence to substitute the traditional price time series will not lose the characteristic of original information of this financial system. This is because, from the mathematical point of view, the transformation is just like a functionf
:
P
r
i
c
e
∈
R
n
→
M
A
C
D
∈
R, and the formula can be given as (6)
M
A
C
D
=
E
M
A
N
,
P
r
i
c
e
_
F
+
E
M
A
N
,
P
r
i
c
e
_
S
,where E
M
A means the Exponential Moving Average which can be represented as a linear transformation E
M
A
=
f
P
r
i
c
e
=
∑
i
=
1
N
ω
i
x
i
=
W
X, where ω
1
=
ω
2
=
⋯
=
ω
N
=
1
/
N are the weight parameters and X
=
{
x
1
,
x
2
,
…
,
x
N
} is the vector of lag prices. As the mapping f
:
P
r
i
c
e
→
M
A
C
D is a bijection, it is not hard to prove there exists an isomorphism between P
r
i
c
e and M
A
C
D, which means the algebraic structure of this function guarantees no information will be lost when transforming it into MACD sequence.Here, 2000 intraday pieces of data of DJIA, from August 2006 to January 2015, are sampled as the original time series. Since the sequence of DJIA continues to be influenced not only by corporate and economic reports but also by domestic and foreign political events such as war and peace, as well as natural disasters that can potentially lead to financial harm, it is a manifest representation of the financial system of the United States. Moreover, the time series of DJIA should be translated into the time series of MACD so that the amplitude of short time trend of financial market can be amplified (seen in Figure2).Figure 2
The price time series of DJIA and its MACD indicator sequence.Yet, there is still a difficulty for financial phase space reconstruction; that is, the MACD series adopted from real-time empirical data of financial prices contains “external” noises which may interfere with the observation of endogenous chaotic phenomenon of financial system. Hence, some denoising measures should be implemented. Here, the wavelet filter is used to smooth the signal series of MACD.A brief introduction of the wavelet denoising [22] is provided as follows.Given a mother waveletψ
(
t
) (which can be considered simply as a basis function of L
2), the continuous wavelet transform (CWT) of a function f
(
t
) is defined as (7)
C
W
T
ψ
f
a
,
b
=
W
f
b
,
a
=
1
a
∫
-
∞
∞
f
t
ψ
t
-
b
a
d
t
,where a
,
b
∈
R, a
≠
0 are dilating and translating coefficients, respectively. The multiplication of 1
/
a is for energy normalization purposes so that the transformed signal has the same energy at every scale. When a
=
2
-
m and b
=
n
2
-
m
,
m
,
n
∈
Z, (7) can be simplified as(8)
D
W
T
ψ
f
m
,
n
=
∫
-
∞
∞
f
t
ψ
m
,
n
t
d
t
,where ψ
m
,
n
t
=
2
-
m
ψ
(
2
m
t
-
n
) is the dilated and translated version of the mother wavelet ψ
(
t
).The relationship among the low-pass filter, the high-pass filter, the scalar function(
t
), and the wavelet ∅
(
t
) can be stated as follows:(9)
A
i
k
=
∑
k
=
-
∞
∞
D
i
+
1
k
·
g
-
n
+
2
k
+
A
i
+
1
k
·
h
-
n
+
2
k
,where A
i
k
=
∑
n
A
i
-
1
t
·
h
[
2
k
-
n
] denotes the approximated coefficient of the signal series and D
i
k
=
∑
n
A
i
-
1
t
·
g
[
2
k
-
n
] denotes the detailed coefficient of the signal series.The wavelet denoise coefficients of MACD time series are given in Figure3.Figure 3
The wavelet denoising on the sampled MACD time series of DJIA. Note: this figure shows the effect of wavelet denoising of MACD series. The Mexican Hat (Mexh) is used as the wavelet base function, whereψ
t
=
(
2
/
3
)
π
-
1
/
4
(
1
-
x
)
e
-
x
2
/
2. The three left plots illustrate the threshold under different scales of details coefficients, the top right plot demonstrates the comparison between original times series and filtered series, and the two right plots below describe the coefficients among different scale levels.As we have seen in Figure3 the Mexh wavelet decomposed the MACD into three levels of details coefficients D
i
k
,
i
=
1,2
,
3. Each D
i
k represents the importance of the original signal under different kinds of scales, and if the value of D
i
k is small enough it may be just considered as noise and filtered. However how to choose suitable threshold to eliminate those noises remains difficult. Here some statistical test indicators S
_
I, N
_
I and S
N
_
I are proposed to measure the performance of noise removing, where S
_
I
=
∑
i
k
(
x
i
-
x
~
)
2 denotes the smoothing of the filtered series, N
_
I
=
∑
i
k
x
i
-
x
i
′ denotes the difference between original time series and noise-filtered sequence, and S
N
_
I
=
S
_
I
+
N
_
I presents the denoising performance. As we know, if we increase the threshold of noise filtering, the value of S
_
I will rise and the value of N
_
I will down; instead if we decrease the threshold, the two indicators will move in the opposite direction. Thus a suitable threshold of each level is balancing these two indicators (S
N
_
I) at a relatively low value.Table1 shows three different common methods of the wavelet denoising. The increasing of threshold at Lev1, Lev2, and Lev3 has caused the rising of S_I and dropping of N_I. But the performance of SN_I had a special change of descending firstly and then ascending. In order to find out the best noise-filtered series, the wavelet filter with a white denoised threshold has been adopted because it has a minimal value of S
N
_
I
=
82.455.Table 1
Comparison between original MACD and noise-filtered series.
Denoise methods
Threshold
Statistical test
Lev1
Lev2
Lev3
S_I
N_I
SN_I
Original signal
—
—
—
85.798
0
85.798
Unscaled denoise
3.865
3.865
3.865
80.871
2.847
83.7180
White denoise
34.039
39.521
54.711
75.311
7.144
82.455
No-white denoise
288.649
121.691
54.039
72.844
15.611
88.455From top right plot of Figure3 the difference between original and denoised signals can be compared. We can see that the wavelet filtered sequence not only reserves the structure of MACD signal but also smoothes the noise of original time series.


## 4. Chaotic Analysis of DJIA’s Financial System

### 4.1. Attractor Reconstruction

In Sections2 and 3 we have discussed the theoretical methodology of chaos detecting. Consequently, in order to confirm the chaotic behavior of financial system, empirical data of financial time series will be studied in this section to confirm the chaotic behavior of the financial system.Firstly, the time series of DJIA should be translated into the time series of MACD so that the amplitude of short time trend of financial market can be amplified, with a designed wavelet filter to smooth the curve of MACD sequence.Secondly, the nonlinearity of MACD sequence is discussed before we do the phase space reconstruction and Laplacian Eigenmaps, and it ensures the assumption about nonlinear system of financial market. In this part, DBS test is applied to discuss this property; the BDS test detects nonlinear serial dependence in time series, after removing the linear structure from the time series by any linear model (such as ARMA(p, q), ARCH(q), GARCH(p, q)), the BDS test is a method to testify whether the time series is linearity. If the null hypothesis is rejected, the time series can be regarded as nonlinearity. BDS statistic can be formalized as follows:(10)
W
m
,
n
ε
=
n
-
m
+
1
C
m
,
n
ε
-
C
1
,
n
-
m
+
1
ε
σ
m
,
n
ε
,where the estimated variance of C
m
,
n
ε
-
C
1
,
n
-
m
+
1
(
ε
) is given by (3), and the variance of σ
m
,
n
(
ε
) can be illustrated as(11)
σ
m
,
n
ε
=
2
k
m
+
2
∑
j
=
1
m
-
1
k
m
-
j
c
2
j
+
m
-
1
2
c
2
m
-
m
2
k
c
2
m
,where parameter c is the first-dimensional correlation integral: c
=
c
1
,
n
(
ε
).Tables2 and 3 show the results of the BDS test of the residual of the MACD sequence under different sup-norm separation r (besides before the BDS test, AR(2) model is used to remove the linear structure of the original time series), and we can see that the probabilities of rejecting the null hypothesis are overwhelming significance which means the MACD time series has a strong nonlinearity.Table 2
BDS test for the RESID of MACD sequence withr
=
1 std.
Dimension
BDS statistic
Std. error
Z-statistic
Prob.
2
0.011587
0.002083
5.562627
0
3
0.0322
0.003313
9.719828
0
4
0.050152
0.003948
12.70234
0
5
0.061936
0.004119
15.03699
0
6
0.068312
0.003976
17.18115
0Table 3
BDS test for the RESID of MACD sequence withr
=
1.5 std.
Dimension
BDS statistic
Std. error
Z-statistic
Prob.
2
0.007182
0.001558
4.611332
0
3
0.013112
0.001428
9.18146
0
4
0.013243
0.000983
13.46789
0
5
0.01034
0.000593
17.42805
0
6
0.007488
0.000332
22.58516
0Thirdly, the redundant high-dimensional Euclidean space of original financial system is built using the technique of phase space reconstruction, with the delay timeτ
d
=
1. According to the correlation integral, the correlation dimension is C
=
4.5, and the embedding dimension is m
=
9 (Figure 4).Figure 4
The plot correlation dimension. Note: when delay time is fixed toτ
d
=
1, the slope of C
m converges gradually after m
=
9, which means the correlation dimension approximately equals C
=
4.5.Finally, the Laplacian Eigenmaps is used to extract the chaotic attractor from the redundant reconstructed phase space of DJIA, so that the attractor of this financial system embedded in a low-dimensional manifold can be decided. Here, two methods of attractor reconstruction (the traditional phase space reconstruction and the Laplacian Eigenmaps) are used to compare the effect of chaos detecting. As can be seen in Figure5, the attractor constructed by Takens’ phase space reconstruction is a spaghetti mess (Figure 5(a)). In contrast, the attractor constructed by Laplacian Eigenmaps embedded in a low-dimensional (3D) manifold looks like a closed ring rotating with the fluctuation of price trends denoted by MACD indicator (Figure 5(b)). Despite coarseness, it can be regarded as a stable quasi-periodic circular orbit which contains price behaviors of the financial market.Figure 5
Attractor reconstruction by normal method and Laplacian Eigenmaps. Note: (a) shows the traditional phase space reconstruction in 3D Euclidean space, in which the embedding dimensionm
=
9 and time delay τ
d
=
9; (b) demonstrates the attractor reconstructed by Laplacian Eigenmaps, where the embedding dimension m
=
9 and time delay τ
d
=
1.
(a)
(b)


### 4.2. Determinism Test

Although we have reconstructed the attractor, it is inappropriate to claim that the chaotic attractor has been found because we have not yet verified the chaotic characteristics of this attractor. A Determinism Test must be performed to assure that the studied attractor is indeed organized by a deterministic system, which enables us to distinguish between deterministic chaos and irregular random behaviors resembling chaos.If a system is deterministic, it can be described by a set of ordinary differential equations [23]. But, for financial system, it is impossible to formalize a group of deterministic differential equations (nonlinear) due to the complexity of its mechanism. Fortunately, the vector field at every point of the attractor is uniquely determined by the differential equations [24], which aids us in constructing the vector field of the system directly from the time series and subsequently testing if the reconstructed vector field is indeed a unique solution in the attractor.The approach is based on the observation that the tangent to the trajectory generated by a deterministic system is a function of position in phase space, and therefore, all the tangents to the trajectory in a given region of phase space have similar orientations. To perform the Determinism Test, the embedding manifold of financial system has to be coarse grained into equal-size boxes, and a vector is assigned to each box that is occupied by the trajectory. Each passi of the trajectory through the jth box generates a unit vector e
i whose direction is determined by the point where the trajectory first enters the box and then passes away. The approximation for the vector field V
j in the jth box of the manifold can be formalized as follows [25]:(12)
V
j
=
1
P
j
∑
i
=
1
P
j
e
i
,where P
j is the number of all passing through the jth box.A determinism factorL
n
d
=
V
j
P
j for the vector field of thed-dimensional financial system can be obtained by combining all occupied boxes containing P
j passes of the trajectory. If the system is deterministic, the obtained vector field should consist solely of vectors that have unit length; namely, the definite measure of the average length L exactly equals 1. In contrast, for a system with a stochastic component, this value is substantially smaller than 1.Figure6 illustrates the result of the Determinism Test where the three-dimensional embedding space is coarse-grained into a 16 × 16 × 16 grid. Although the pertaining determinism factor of the approximated vector field calculated by (12) is L
=
0.81, which is smaller than the unit length, considering the ever-present measurement errors in the observed time series, it can still be declared that there is a deterministic signature within the attractor of DJIA’s financial system.Figure 6
Determinism Test of strange attractor. Note: Determinism Test of the approximated vector field for the attractor of DJIA, which is embedded in a three-dimensional manifold, with the pertaining determinism factorL
=
0.81.


### 4.3. Poincaré Section

Since the deterministic characteristic of the attractor has been confirmed, the chaos phenomenon of the DJIA’s financial system will now be discussed. Instead of traditional methods such as Largest Lyapunov Exponents (LLE) estimation [26], a more convenient topological method called Poincaré Map is applied to detect the chaotic property of DJIA’s financial system [27].In mathematics, particularly in dynamical systems, Poincaré Map is the intersection of a periodic orbit in the state space of a continuous dynamical system with a certain lower-dimensional subspace called Poincaré section, which is transversal to the flow of the system. More precisely, one considers a periodic orbit with initial conditions within a section of the space, which leaves that section afterwards, and observes the point at which this orbit first returns to the section. Then, he/she creates a map to send the first point to the second. The transversality of the Poincaré section means that periodic orbits starting from the subspace flow through it rather than paralleling it. Obviously, it preserves many properties of periodic and quasi-periodic orbits of the original system and has a lower-dimensional state space. Hence, Poincaré section is often used for analyzing the chaos system.Figure7 shows the Poincaré section of the DJIA’s attractor which indicates a discrete dynamical structure of original financial system. As we know, if the fixed points of a Poincaré section are asymptotically stable, or hierarchically distributed, this Poincaré section can be regarded as a chaotic state. Thus, the ribbon scatters projecting on the section imply that there is indeed a chaotic attractor embedded in the original phase space of financial time series.Figure 7
A two-dimensional Poincaré section of the embedded attractor.


### 4.4. The Translation Analyses of Chaotic Attractor

Up to now, the strange attractor of DJIA’s financial time series has been found, but the internal mechanisms of chaotic financial system are still mystery; namely, the generation of chaos cannot be merely understood by detecting of the strange attractor. To solve this issue, the translation process of chaotic attractor must be analyzed. However, due to the complexity of a financial market, it is hard to find out some specific differential equations to represent the behavior of this dynamical system. This means the translation of the chaotic attractor cannot be controlled by adjusting the parameters of the underlying differential equations.Hence, a period span analysis is proposed to study the homotopic process of chaotic attractor. In this paper, the empirical data of DJIA’s time series is divided into five periods (as seen in Figure8):(i)
Oscillator Period (08/2006–11/2007), which means the market is in a stable area without strong volatility.
(ii)
Drop Down (11/2007–09/2008), which denotes the unilateral decline range.
(iii)
Short-Term Rise (09/2008–11/2009), which represents a unilateral up range.
(iv)
Trend Reversal (11/2007–11/2009) that combines both Drop Down and Short-Term Rise periods together illustrates a sharp shock of the financial system.
(v)
Long-Term Rise (11/2009–01/2015), just as its name, implies the long time boost of the financial system.Figure 8
Period span analysis chart of the DJIA’s attractor translation process.The homotopic translation of the attractor is accompanied with the period change of the dynamical financial system. In theOscillator Period, the attractor of DJIA looks like a limit cycle with a periodic trajectory, which implies that the fluctuation of the financial system is stable. However, when there happens to be a trend (Drop Down/Short-Term Rise), with the increasing energy (here, energy can be regarded as the sentiment of the financial market which is produced by the investors of the financial market) of the system, the oscillation of the limit cycle will homotopically translate to a torus with a quasi-periodic trajectory. And this torus will finally turn to be a strange attractor when it absorbs enough energy (as the attractor ofLong-Term Rise shows). Furthermore, the attractor of a financial system will dissipate to a fractal only if a huge volatility occurs in the market such asTrend Reversal. These translations of the attractor reveal that the birth and death process of a chaos phenomenon in the financial system corresponds to the sentiment of the financial market. We can summarize as follows: the chaotic attractor is not always existent; it is generated by the expanding energy of the system and will dissipate until the structure itself cannot withstand this massive energy.


## 4.1. Attractor Reconstruction

In Sections2 and 3 we have discussed the theoretical methodology of chaos detecting. Consequently, in order to confirm the chaotic behavior of financial system, empirical data of financial time series will be studied in this section to confirm the chaotic behavior of the financial system.Firstly, the time series of DJIA should be translated into the time series of MACD so that the amplitude of short time trend of financial market can be amplified, with a designed wavelet filter to smooth the curve of MACD sequence.Secondly, the nonlinearity of MACD sequence is discussed before we do the phase space reconstruction and Laplacian Eigenmaps, and it ensures the assumption about nonlinear system of financial market. In this part, DBS test is applied to discuss this property; the BDS test detects nonlinear serial dependence in time series, after removing the linear structure from the time series by any linear model (such as ARMA(p, q), ARCH(q), GARCH(p, q)), the BDS test is a method to testify whether the time series is linearity. If the null hypothesis is rejected, the time series can be regarded as nonlinearity. BDS statistic can be formalized as follows:(10)
W
m
,
n
ε
=
n
-
m
+
1
C
m
,
n
ε
-
C
1
,
n
-
m
+
1
ε
σ
m
,
n
ε
,where the estimated variance of C
m
,
n
ε
-
C
1
,
n
-
m
+
1
(
ε
) is given by (3), and the variance of σ
m
,
n
(
ε
) can be illustrated as(11)
σ
m
,
n
ε
=
2
k
m
+
2
∑
j
=
1
m
-
1
k
m
-
j
c
2
j
+
m
-
1
2
c
2
m
-
m
2
k
c
2
m
,where parameter c is the first-dimensional correlation integral: c
=
c
1
,
n
(
ε
).Tables2 and 3 show the results of the BDS test of the residual of the MACD sequence under different sup-norm separation r (besides before the BDS test, AR(2) model is used to remove the linear structure of the original time series), and we can see that the probabilities of rejecting the null hypothesis are overwhelming significance which means the MACD time series has a strong nonlinearity.Table 2
BDS test for the RESID of MACD sequence withr
=
1 std.
Dimension
BDS statistic
Std. error
Z-statistic
Prob.
2
0.011587
0.002083
5.562627
0
3
0.0322
0.003313
9.719828
0
4
0.050152
0.003948
12.70234
0
5
0.061936
0.004119
15.03699
0
6
0.068312
0.003976
17.18115
0Table 3
BDS test for the RESID of MACD sequence withr
=
1.5 std.
Dimension
BDS statistic
Std. error
Z-statistic
Prob.
2
0.007182
0.001558
4.611332
0
3
0.013112
0.001428
9.18146
0
4
0.013243
0.000983
13.46789
0
5
0.01034
0.000593
17.42805
0
6
0.007488
0.000332
22.58516
0Thirdly, the redundant high-dimensional Euclidean space of original financial system is built using the technique of phase space reconstruction, with the delay timeτ
d
=
1. According to the correlation integral, the correlation dimension is C
=
4.5, and the embedding dimension is m
=
9 (Figure 4).Figure 4
The plot correlation dimension. Note: when delay time is fixed toτ
d
=
1, the slope of C
m converges gradually after m
=
9, which means the correlation dimension approximately equals C
=
4.5.Finally, the Laplacian Eigenmaps is used to extract the chaotic attractor from the redundant reconstructed phase space of DJIA, so that the attractor of this financial system embedded in a low-dimensional manifold can be decided. Here, two methods of attractor reconstruction (the traditional phase space reconstruction and the Laplacian Eigenmaps) are used to compare the effect of chaos detecting. As can be seen in Figure5, the attractor constructed by Takens’ phase space reconstruction is a spaghetti mess (Figure 5(a)). In contrast, the attractor constructed by Laplacian Eigenmaps embedded in a low-dimensional (3D) manifold looks like a closed ring rotating with the fluctuation of price trends denoted by MACD indicator (Figure 5(b)). Despite coarseness, it can be regarded as a stable quasi-periodic circular orbit which contains price behaviors of the financial market.Figure 5
Attractor reconstruction by normal method and Laplacian Eigenmaps. Note: (a) shows the traditional phase space reconstruction in 3D Euclidean space, in which the embedding dimensionm
=
9 and time delay τ
d
=
9; (b) demonstrates the attractor reconstructed by Laplacian Eigenmaps, where the embedding dimension m
=
9 and time delay τ
d
=
1.
(a)
(b)


## 4.2. Determinism Test

Although we have reconstructed the attractor, it is inappropriate to claim that the chaotic attractor has been found because we have not yet verified the chaotic characteristics of this attractor. A Determinism Test must be performed to assure that the studied attractor is indeed organized by a deterministic system, which enables us to distinguish between deterministic chaos and irregular random behaviors resembling chaos.If a system is deterministic, it can be described by a set of ordinary differential equations [23]. But, for financial system, it is impossible to formalize a group of deterministic differential equations (nonlinear) due to the complexity of its mechanism. Fortunately, the vector field at every point of the attractor is uniquely determined by the differential equations [24], which aids us in constructing the vector field of the system directly from the time series and subsequently testing if the reconstructed vector field is indeed a unique solution in the attractor.The approach is based on the observation that the tangent to the trajectory generated by a deterministic system is a function of position in phase space, and therefore, all the tangents to the trajectory in a given region of phase space have similar orientations. To perform the Determinism Test, the embedding manifold of financial system has to be coarse grained into equal-size boxes, and a vector is assigned to each box that is occupied by the trajectory. Each passi of the trajectory through the jth box generates a unit vector e
i whose direction is determined by the point where the trajectory first enters the box and then passes away. The approximation for the vector field V
j in the jth box of the manifold can be formalized as follows [25]:(12)
V
j
=
1
P
j
∑
i
=
1
P
j
e
i
,where P
j is the number of all passing through the jth box.A determinism factorL
n
d
=
V
j
P
j for the vector field of thed-dimensional financial system can be obtained by combining all occupied boxes containing P
j passes of the trajectory. If the system is deterministic, the obtained vector field should consist solely of vectors that have unit length; namely, the definite measure of the average length L exactly equals 1. In contrast, for a system with a stochastic component, this value is substantially smaller than 1.Figure6 illustrates the result of the Determinism Test where the three-dimensional embedding space is coarse-grained into a 16 × 16 × 16 grid. Although the pertaining determinism factor of the approximated vector field calculated by (12) is L
=
0.81, which is smaller than the unit length, considering the ever-present measurement errors in the observed time series, it can still be declared that there is a deterministic signature within the attractor of DJIA’s financial system.Figure 6
Determinism Test of strange attractor. Note: Determinism Test of the approximated vector field for the attractor of DJIA, which is embedded in a three-dimensional manifold, with the pertaining determinism factorL
=
0.81.


## 4.3. Poincaré Section

Since the deterministic characteristic of the attractor has been confirmed, the chaos phenomenon of the DJIA’s financial system will now be discussed. Instead of traditional methods such as Largest Lyapunov Exponents (LLE) estimation [26], a more convenient topological method called Poincaré Map is applied to detect the chaotic property of DJIA’s financial system [27].In mathematics, particularly in dynamical systems, Poincaré Map is the intersection of a periodic orbit in the state space of a continuous dynamical system with a certain lower-dimensional subspace called Poincaré section, which is transversal to the flow of the system. More precisely, one considers a periodic orbit with initial conditions within a section of the space, which leaves that section afterwards, and observes the point at which this orbit first returns to the section. Then, he/she creates a map to send the first point to the second. The transversality of the Poincaré section means that periodic orbits starting from the subspace flow through it rather than paralleling it. Obviously, it preserves many properties of periodic and quasi-periodic orbits of the original system and has a lower-dimensional state space. Hence, Poincaré section is often used for analyzing the chaos system.Figure7 shows the Poincaré section of the DJIA’s attractor which indicates a discrete dynamical structure of original financial system. As we know, if the fixed points of a Poincaré section are asymptotically stable, or hierarchically distributed, this Poincaré section can be regarded as a chaotic state. Thus, the ribbon scatters projecting on the section imply that there is indeed a chaotic attractor embedded in the original phase space of financial time series.Figure 7
A two-dimensional Poincaré section of the embedded attractor.


## 4.4. The Translation Analyses of Chaotic Attractor

Up to now, the strange attractor of DJIA’s financial time series has been found, but the internal mechanisms of chaotic financial system are still mystery; namely, the generation of chaos cannot be merely understood by detecting of the strange attractor. To solve this issue, the translation process of chaotic attractor must be analyzed. However, due to the complexity of a financial market, it is hard to find out some specific differential equations to represent the behavior of this dynamical system. This means the translation of the chaotic attractor cannot be controlled by adjusting the parameters of the underlying differential equations.Hence, a period span analysis is proposed to study the homotopic process of chaotic attractor. In this paper, the empirical data of DJIA’s time series is divided into five periods (as seen in Figure8):(i)
Oscillator Period (08/2006–11/2007), which means the market is in a stable area without strong volatility.
(ii)
Drop Down (11/2007–09/2008), which denotes the unilateral decline range.
(iii)
Short-Term Rise (09/2008–11/2009), which represents a unilateral up range.
(iv)
Trend Reversal (11/2007–11/2009) that combines both Drop Down and Short-Term Rise periods together illustrates a sharp shock of the financial system.
(v)
Long-Term Rise (11/2009–01/2015), just as its name, implies the long time boost of the financial system.Figure 8
Period span analysis chart of the DJIA’s attractor translation process.The homotopic translation of the attractor is accompanied with the period change of the dynamical financial system. In theOscillator Period, the attractor of DJIA looks like a limit cycle with a periodic trajectory, which implies that the fluctuation of the financial system is stable. However, when there happens to be a trend (Drop Down/Short-Term Rise), with the increasing energy (here, energy can be regarded as the sentiment of the financial market which is produced by the investors of the financial market) of the system, the oscillation of the limit cycle will homotopically translate to a torus with a quasi-periodic trajectory. And this torus will finally turn to be a strange attractor when it absorbs enough energy (as the attractor ofLong-Term Rise shows). Furthermore, the attractor of a financial system will dissipate to a fractal only if a huge volatility occurs in the market such asTrend Reversal. These translations of the attractor reveal that the birth and death process of a chaos phenomenon in the financial system corresponds to the sentiment of the financial market. We can summarize as follows: the chaotic attractor is not always existent; it is generated by the expanding energy of the system and will dissipate until the structure itself cannot withstand this massive energy.


## 5. Conclusions

In summary, this paper discussed how to detect the chaotic attractor in a financial time series and analyzed whether there was a chaos phenomenon in real-time financial time series under the assumption that the chaotic attractor of financial time series lies on a low-dimensional manifold embedded into a high-dimensional Euclidean space.Firstly, a theoretical approach was narrated to demonstrate how to extract the attractor of a complicated dynamical system (Figure1 gives the entire procedure of this approach). Two ways are essential: firstly the embedding theorem was used to reconstruct a high-dimensional phase space with the fixed time delay τ
d
=
1, which constructed a linear but redundant Euclidean space. Then Laplacian Eigenmaps was applied to map the redundant phase space to a low-dimensional manifold so as to find the potential attractor.Secondly, since the traditional time series are not cyclic, the MACD sequence was adopted to substitute price series for the reconstruction of the phase space of DJIA, and the wavelet denoising was used to filter the exogenous interference of MACD series.Finally, based on the empirical study of DJIA, a low-dimensional attractor was captured. By Determinism Test, and Poincaré section, the deterministic and chaotic behaviors of this attractor were identified, respectively, which help us solve the two most controversial views mentioned in Introduction. Moreover, the homotopy study of the chaotic attractor unfolds the internal mechanisms of chaos generation in DJIA’s financial system. Meanwhile, the result of this paper also is self-consistent with the aforementioned assumption about the chaotic attractor of financial time series.



---

*Source: 8087178-2016-06-22.xml*