Fractional Interaction of Financial Agents in a Stock Market Network

Mehmet Ali Balcı

Journal & Issue Details

Applied Mathematics and Nonlinear Sciences

Journal Details

Format: Journal
eISSN: 2444-8656
First Published: 01 Jan 2016
Publication timeframe: 2 times per year
Languages: English

Complex systems are mathematical structures involving interacting agents at different levels. These interactions emerge from the financial, chemical, social, and computer system entities. In the realm of computational finance, a financial market can be viewed as interacting group of boundedly-rational agents and its fluctuation represent strong nonlinearity and persistent memory. The mathematical tools such as network and graph theories can be used to understand and analyze these systems [1, 2]. There are several models expressed in the terms of differential equations in biological complex systems. For instance, the virus models that classify individuals and hosts can be used to analyze spread of a contagion [3, 4, 5]. Besides, bursting electrical activity in the pancreatic β-cell, population models, and unilingual-bilingual interactions [6], and the interaction of biological species living together [7] can also be modelled by differential equations. However, these models are not only restricted to biological systems. Recent studies show that complex systems involving financial agents have similar structures as systems involving biological agents [8, 9]. Therefore, it is reasonable to model the interaction of financial agents as we model the interaction of biological agents.

In a financial market, heterogeneous agents interact through simple investment strategies driven by the investors. In a perfect rational market, information is transmitted continuously and agents adopt their behavior accordingly. Besides, asset prices reflect economic fundamentals. Agents are considered as they only interact though price system. Hence, a complex network where agents are expressed as vertices and edges are formed by the correlation of price fluctuations emerges as a powerful mathematical tool to model such a financial system in the traditional way. In contrast to Keynesian approach, such traditional way takes account of prices of assets as they are only driven by market fundamentals and the role of market psychology is neglected. Even though we use the traditional way to express our model in this study, we need to point out two important classes of investors which are called chartists and fundamentalists in the traditional way of interaction of agents [10]. Chartists tend to look for simple patterns such as trends, past prices, and base and make their investment upon those patterns. Conversely, fundamentalists make their decisions upon the expectation of asset price as moving towards its fundamental value. The fundamentalist investors buy or sell assets that are under or overvalued. The market tends to be dominated by one of fundamentalists or chartists. However, since the behavior of the agents is persistent, the majority of agents switches to the other view at certain point [11, 12].

Our approach in this study aims to model interaction of agents in a stock market network in the traditional way. We first use a threshold method to construct a network model where vertices are the companies operating in a stock market and edges are formed by the correlation distance of daily logarithmic returns of stock prices. The dimensionality of the resulting network model would be really high and the patterns that yield power law of degree distributions would be disappeared, however it would involve optimally many edges to characterize community structures. By maximization of the modularity of edges in the network, we can cluster agents into densely connected vertex sets. Then, each cluster has its subdominant ultra-metric structure that is a hierarchical structure with at least one leading actor. We set the number of cluster to two, then assume the investors, even chartists or fundamentalists, start to invest one cluster regarding to factors such that merging, capital augmentation, public flotation, etc. Then, the investors in the other cluster start to sell assets to get the capital to invest increasing valued assets. Therefore, the price fluctuation spread within each cluster by conducting leading actors. However, at certain time, the profit realizations start within the asset price increasing cluster, and then the capital emergent by the profit realization is used to invest assets in price decreasing cluster. Eventually, the two clusters find an equilibrium state.

In the complex network model of financial agents the interactions are modelled by the correlations of long time series [13, 14, 15, 16, 17, 18, 19]. Beside the other types of complex systems, financial systems have the strong memory and heredity properties. Therefore, while using differential equations in financial models it is much more useful to get fractional calculus involved. Fractional calculus is the extension of the integer order differential and integral operators to fractional orders [20,21]. The dynamic memory in a financial process can be defined as the averaged characteristic that describes the dependence of a process in the past. Such memory assumes withitness of financial agents about the history of the process. In formal way, the information on the state of the process {t, χ(t)} does not only affect the behavior of financial agents, but also the information about the process state {τ, χ(τ)} also has effect at τ ∈ [0,t]. This effect is related with the fact that the change of the factors can leads to different amount of change in indicators that is there exist multivalent dependencies among variables. One type of such memory of financial agents is called the fading memory and have range application area in physical sciences [22, 23, 24, 25, 26, 27, 28, 29, 30]. In this study, we assume that financial agents can remember the previous changes of investments and the impact of these changes on the output by following fading memory by using Caputo’s definition of fractional derivative.

In this study, we present our model with fractional derivative as similar as the model describing biological species living together. This model of biological species is first given in [7] as the following usual integro-differential equations: (1) $\frac{d μ_{1}}{dt} = μ_{1} (t) [k_{1} - γ_{1} μ_{2} (t) - \int_{t - T_{0}}^{t} f_{1} (t - τ) μ_{2} (τ) d τ], k_{1} > 0$ {{d{\mu _1}} \over {dt}} = {\mu _1}(t)\left[ {{k_1} - {\gamma _1}{\mu _2}(t) - \int_{t - {T_0}}^t {f_1}(t - \tau ){\mu _2}(\tau )d\tau } \right],\;\;{k_1} > 0(2) $\frac{d μ_{2}}{dt} = μ_{2} (t) [- k_{2} + γ_{2} μ_{1} (t) + \int_{t - T_{0}}^{t} f_{2} (t - τ) μ_{1} (τ) d τ], k_{2} > 0 .$ {{d{\mu _2}} \over {dt}} = {\mu _2}(t)\left[ { - {k_2} + {\gamma _2}{\mu _1}(t) + \int_{t - {T_0}}^t {f_2}(t - \tau ){\mu _1}(\tau )d\tau } \right],\;\;{k_2} > 0.

Several solution methods are also presented to study this model [31, 32, 33]. The characterization of the fractional order of the model is also studied in [34].

The rest of the paper is organized as follow: In Section 2, we present the preliminaries about to fractional calculus and graph theoretical concepts that we use throughout the paper. We start our analysis by first determining the financial agents in Section 3. The stock market we choose to study is Borsa Istanbul Stock Exchange Market (XU100). The agents are the companies operating in XU100 and expressed with the time series of the time span of working days from 2013 to 2015. Afterwards, we determine the clusters of financial agents that have strong correlations by using high modularity method. Afterwards we introduce the Adomian decomposition method for solving the system numerically. In Section 5, we present the results by solving the model with Adomian decomposition method. And finally, in Section 6, we deeply discuss the computational results.

Preliminaries

Fractional calculus is an efficient mathematical tool to express complex system phenomenon which involve memory effect. Hence, we use fractional derivatives and integrals to study the model we present in this study. In this section, we give some basics about the fractional calculus in Caputo sense. We also introduce some basics about the graph theory which is the fundamental tool for network modelling. Throughout the paper we let Γ to represent Gamma function that is an extension of the factorial function.

2.1

Fractional Calculus

The generalization of the integer order differentiation and integration to the fractional order is called fractional calculus [20, 21]. The basic definitions and properties of fractional calculus theory is given as follows:

Definition 1

[20] For f (x) ∈ 𝒞 (a,b) and n − 1 < α ≤ n, the Caputo fractional derivative operator of order α is given as $_{a}^{C} D_{t}^{α} f (t) = \frac{1}{Γ (n - 1)} \int_{a}^{t} \frac{f^{(n)} (τ)}{{(t - τ)}^{α + 1 - n}} d τ .$ \ _a^CD_t^\alpha f(t) = {1 \over {\Gamma (n - 1)}}\int_a^t {{{f^{(n)}}(\tau )} \over {{{(t - \tau )}^{\alpha + 1 - n}}}}\;d\tau .

Throughout this paper, we denote the Caputo fractional derivative operator as $_{a}^{C} D_{t}^{α} = D_{a}^{α}$ \ _a^CD_t^\alpha = D_a^\alpha . We also let a = 0 since our formulation only involves the initial conditions as t = 0.

Definition 2

[20] The Riemann-Liouville fractional integral operator of order α ≥ 0 of a function f is defined as $\begin{array}{l} J^{α} f (t) = \frac{1}{Γ (α)} \int_{0}^{t} {(x - τ)}^{α - 1} f (τ) d τ, α > 0, t > 0 \\ J^{0} f (t) = f (t) . \end{array}$ \matrix{ {{J^\alpha }f(t) = {1 \over {\Gamma (\alpha )}}\int_0^t {{(x - \tau )}^{\alpha - 1}}f(\tau )\;d\tau ,\;\;\alpha > 0,t > 0} \hfill \cr {{J^0}f(t) = f(t).} \hfill}

The several properties of the Riemann-Liouville fractional integral operator can be found in [35, 36, 37]. Since the Caputo fractional derivative allows traditional initial and boundary conditions to be included in the formulation of the problem [38], we present our model in the sense of Caputo fractional derivative. By the introduction of the J^α, the $D_{a}^{α}$ D_a^\alpha can also be expressed as $D_{0}^{α} f (t) = J^{m - α} D^{m} f (t),$ D_0^\alpha f(t) = {J^{m - \alpha }}{D^m}f(t), where m − 1 < α ≤ m, m ∈ ℕ,t > 0.

Also, the following two basic properties of the entwined relations among Caputo and Riemann-Liouville fractional operators are needed to present the solution of the fractional differential equations.

Lemma 1

[35] If m − 1 < α ≤ m, m ∈ ℕ, then $D_{0}^{α} J^{α} f (t) = f (t)$ D_0^\alpha {J^\alpha }f(t) = f(t)and $J^{α} D_{0}^{α} f (t) = f (t) - \sum_{k = 0}^{m - 1} f^{(k)} (0^{+}) \frac{t^{k}}{k!}, t > 0 .$ {J^\alpha }D_0^\alpha f(t) = f(t) - \sum\limits_{k = 0}^{m - 1} {f^{(k)}}({0^ + }){{{t^k}} \over {k!}},\;\;t > 0.

2.2

Graph Theory

The real–world problems are often expressed with the relations of interacting individuals. One of the efficient mathematical tools to represent such relations is the simple graphs. In the Stock Market Networks, interactions of financial agents can be modelled by simple graphs. Let V be the set of the interacting individuals and E be the set of relations, then a simple graph is a tuple G = (V,E). Here we call V as the vertex set and E as the set of edges. The each element of E is the unordered pair of vertices as e_k = (v_i,v_j), where v_i,v_j ∈ V for all i, j,k ∈ ℕ. The number of edges incident to a vertex v is called as the degree of v and we denote the degree as d_v.

A sequence of edges between the vertices v_i and v_j is called a path, and if there is a path between any vertices of the graph G, then G is called as connected. If there is an edge between all elements of V, then G is called as a complete graph. The k-clique of the graph G is the complete subgraph of G which involves k-many vertices of G.

For the simple graph G = (V,E) with the unordered edges, a binary matrix which has the entities as $A_{G} (i, j) = {\begin{array}{l} 1, & if (v_{i}, v_{j}) \in E \\ 0, & otherwise \end{array}$ {A_G}(i,j) = \left\{ {\matrix{ {1,} \hfill & {if\;({v_i},{v_j}) \in E} \hfill \cr {0,} \hfill & {otherwise} \hfill \cr } } \right. is used to represent the relations and called as adjacency matrix. A_G is symmetric by definition.

Now let |V| = n and D_G be the diagonal degree matrix of G defined as D_G = diag[d_v₁,...,d_vn]. The matrix L_G = A_G − D_G is called the Laplacian matrix of G. The spectrum of the L_G encodes structural properties of G. The one that we use in this study is helpful to construct a threshold network of the financial agents. All eigenvalues of L_G are positive semi-definite with the least one 0. The multiplicity of the 0 eigenvalue equals the connected components number of G [39].

Several types of subgraphs also involve the information about the network which is expressed as a simple graph G = (V,E). One of them is the tree structure that has minimum weight. Such subgraphs are called as Minimum Spanning Tree (MST) and involve the junction vertices which are dominant in the flow of information and comes up with subdominant ultra-metric structure [40, 41]. In the case of financial agents are the vertices of the network, MST gives the hierarchical structure of the financial network [42, 43, 44]. A planar graph is a simple graph that can be embedded on the plane, which is none of the graph edges intersect. Trees like planar graph that involve cliques are also useful to extract information about the network. Such tree like planar graphs have the same hierarchical structure as MST but they contain larger amount of information about the relation among the interacting agents [45, 46, 47]. In [45], authors present a method to obtain a planar graph with maximum non-crossing edges among the agents of a network and called it Planar Maximally Filtered Graph (PMFG).

Model

This study involves 93 companies that have been operating in Borsa Istanbul 100 Stock Exchange Index (XU100) from January 2013 to January 2015. The Pearson correlation coefficient of time series assumes the equality of the length of time series. Hence, even though XU100 has 100 operating companies, we only consider 93 of them which have the equal time length. Trading hours for the stocks are held by two sessions on business days with mid-day break, and one session in some official holidays [44, 48]. The tickers of the companies operating in XU100 and that are considered in this study is given in Table 1. The more details on the data can be found in [44].

Table 1

Sectors of each considered stock

FinancialsAKBNK, SKBNK, SNGYO, TSKB, TEKST, TRGYO, VKGYO, ALGYO, ISGYO, GARAN, ALBRK, GLYHO, ISCTR, YKBNK, SAHOL, GOZDE, HALKB, VAKBN, ECZYT, SAFGY, EKGYO, SAHOL, GSDHO

IndustrialsASELS, TAVHL, TKFEN, TTRAK, CLEBI

Consumer DiscretionaryASUZU, TKNSA, TOASO, YAZIC, AKSA, ARCLK, GSRAY, KARSN, THYAO, BRISA, DOAS, FENER, MNDRS, METRO, VESBE, ADEL, BJKAS, NTTUR, GOODY, OTKAR, TMSN, EGEEN, FROTO, IHLAS

EnergyAYGAZ, TUPRS, IPEKE, KCHOL

TechnologyNETAS, VESTL

MaterialsSASA, AFYON, ANACM, BAGFS, CIMSA, KONYA, KOZAA, ERBOS, KRDMD, PRKME, SISE, ALKIM, TRKCM, GUBRF, KOZAL, BRSAN, KARTN, PETKM, GOLTS, EREGL

CommunicationsTTKOM, TCELL, DOHOL, HURGZ

Consumer StaplesAEFES, CCOLA, BIZIM, ECILC, BIMAS, MGROS, SODA, ULKER

UtilitiesAKSEN, ALARK, TRCAS, ZOREN, ENKAI

The data we use is available with sessional closure price, therefore we calculate sessional closure price logarithmic return as $C l_{i} = log P_{i} (t + 1) - log P_{i} (t),$ C{l_i} = \log {P_i}(t + 1) - \log {P_i}(t), where P_i(t) is the closure price of the i-th stock at the session t. To represent the relation between stock pairs, we use the Pearson correlation coefficient of stocks as $ρ_{ij} = \frac{< C l_{i} C l_{j} > - < C l_{i} > < C l_{j} >}{\sqrt{(< C l_{i}^{2} > - < C l_{i} >^{2}) (< C l_{j}^{2} > - < C l_{j} >^{2})}},$ {\rho _{ij}} = {{ < C{l_i}C{l_j} > - < C{l_i} > < C{l_j} > } \over {\sqrt {( < Cl_i^2 > - < C{l_i}{ > ^2})( < Cl_j^2 > - < C{l_j}{ > ^2})} }}, where < ··· > is the temporal average performed on the trading days.

Pearson correlation coefficient varies between −1 and 1. ρ_ij = 1 indicates that the stocks i and j are completely correlated whilst ρ_ij indicates that the stocks i and j are completely uncorrelated [42]. Hence, it is also possible to introduce a new distance $d_{Corr} : = \sqrt{2 (1 - ρ_{ij})} / 2$ {d_{Corr}}: = \sqrt {2(1 - {\rho _{ij}})} /2 as in [13,44]. We can conclude that if d_Corr(i, j) = 0, then the stocks i and j are completely correlated, and if d_Corr(i, j) = 1, then the stocks i and j are completely uncorrelated.

This distances based on Pearson correlation is helpful to us for edge determination on the network. By using an empirical threshold value among the stocks, it is possible to determine edges representing strong relations as $A_{G} (i, j) = 1 iff d_{Corr} (i, j) \leq ThV,$ {A_G}(i,j) = 1\;\;iff\;\;{d_{Corr}}(i,j) \le ThV, where T hV is the threshold value. The threshold value can be determined by the subdivision of the interval [0,1], where the boundaries are the extremal values of d_Corr, into h many subintervals. The details on the algorithm of network construction and computational complexities can be found in [13, 44].

The model we proposed in this study first deal with the two clusters of financial agents of the network. In the literature, the cluster of densely connected vertices of a network is called graph communities [49, 50]. This densely connection is internal and can be used to analyze the relations that are represented by edges on the network. There are several methods to detect communities in a network [51, 52, 53, 54]. To find the graph communities in the network; we use the Modularity Maximization Method which is based on the maximizing the Newman modularity index [51] defined as $Q = \sum_{i = 1}^{N_{C}} [\frac{E_{k}}{m} - \frac{1}{4 m^{2}} {(\sum_{j \in V_{k}} d_{j})}^{2}],$ Q = \sum\limits_{i = 1}^{{N_C}} \left[ {{{{E_k}} \over m} - {1 \over {4{m^2}}}{{\left( {\sum\limits_{j \in {V_k}} {d_j}} \right)}^2}} \right], where E_k is the number of edges in the k-th module, N_C is the total number of modules, m is the total number of edges and d_j is the vertex degree. Since the resulting communities are non-overlapping and this method let us to determine final number of the communities, we choose it as an efficient tool.

Now let us consider the two communities of financial agents with the total number of investments µ₁(t) and µ₂(t), respectively, at time t. Let us assume the investment in the first community is increasing with the coefficient of increase k₁ and in the second community is decreasing with the coefficient of decrease k₂. Both coefficients k₁ and k₂ are positive reals. If the two communities are left separate; i.e. they are non-overlapping, then the fractional growth of the first can be represented by (3) $D_{0}^{α} μ_{1} (t) = k_{1} μ_{1} (t)$ D_0^\alpha {\mu _1}(t) = {k_1}{\mu _1}(t) and the decline of the second community can be represented by (4) $D_{0}^{α} μ_{2} (t) = - k_{2} μ_{2} (t) .$ D_0^\alpha {\mu _2}(t) = - {k_2}{\mu _2}(t).

The neoclassical liberal economy states that markets always look for the equilibrium state. Hence, if we put these two communities together in the corresponding stock market environment, the decrease of the rate of the increase of first community is proportional to µ₂(t) and vice versa. Therefore, it is reasonable to assume the increase and decrease coefficient as (5) ${k^{'}}_{1} = k_{1} - γ_{1} μ_{2} (t)$ {k'_1} = {k_1} - {\gamma _1}{\mu _2}(t) and (6) ${k^{'}}_{2} = k_{2} + γ_{2} μ_{1} (t),$ {k'_2} = {k_2} + {\gamma _2}{\mu _1}(t), where γ₁ and γ₂ are the proportionality constants which depend on other investor behavior, respectively. The actual decrease and increase of the investments in the communities are due not only to the presence of the other community but also to all previous presences for the whole time interval t − T₀ < τ < t, where T₀ is the finite heredity duration of both communities. In addition to the present γ₁ and γ₂ factors, we may have the record of decrease as f₁(τ) and increase as f₂(τ). Therefore, by considering the heredity duration of both communities, the total decrease in k₁ in the time interval T₀ is (7) $δ k_{1} = - \int_{t - T_{0}}^{t} f_{1} (t - s) μ_{2} (s) ds$ \delta {k_1} = - \int_{t - {T_0}}^t {f_1}(t - s){\mu _2}(s)\;ds the total fractional increase in k₂ is (8) $δ k_{2} = \int_{t - T_{0}}^{t} f_{2} (t - s) μ_{1} (s) ds .$ \delta {k_2} = \int_{t - {T_0}}^t {f_2}(t - s){\mu _1}(s)\;ds.

Now, by considering effective values of the k₁ and k₂ and the equations 3–8, the fractional model with fading memory of the equilibrium state of the two communities of financial agents in same stock market can be given as the system of fractional integro-differential equations as follows: (9) $D_{0}^{α} μ_{1} (t) = μ_{1} (t) [k_{1} - γ_{1} μ_{2} (t) - \int_{t - T_{0}}^{t} f_{1} (t - s) μ_{2} (s) ds],$ D_0^\alpha {\mu _1}(t) = {\mu _1}(t)\left[ {{k_1} - {\gamma _1}{\mu _2}(t) - \int_{t - {T_0}}^t {f_1}(t - s){\mu _2}(s)\;ds} \right],(10) $D_{0}^{α} μ_{2} (t) = μ_{2} (t) [- k_{2} + γ_{2} μ_{1} (t) + \int_{t - T_{0}}^{t} f_{2} (t - s) μ_{1} (s) ds],$ D_0^\alpha {\mu _2}(t) = {\mu _2}(t)\left[ { - {k_2} + {\gamma _2}{\mu _1}(t) + \int_{t - {T_0}}^t {f_2}(t - s){\mu _1}(s)\;ds} \right],(11) $μ_{1} (0) = N_{1}, μ_{2} (0) = N_{2},$ {\mu _1}(0) = {N_1},\,\,\,\,{\mu _2}(0) = {N_2}, where N₁ and N₂ are the initial conditions.

Method

In this section, we present the graph theoretical methods to determine parameters k₁, k₂, γ₁, and γ₂ of the fractional integro-differential equation model in the initial value problem 9–11 and the numerical solution that is based on Adomian Decomposition Method.

4.1

Parameter Estimation

The parameters in network models can be estimated by using graph theoretical concepts. For the system 9–11, to determine the coefficients of increase and decrease, we use an interpretation of the displacement of a vertex in a network from its equilibrium state while the network is under a thermal bath. This thermal bath can be seen as the change of investment strategies on given network. This procedure is called vibrational potential and first presented in [55]. The later studies consider vibrational potential as an efficient measure to vertex centrality [56, 57, 58]. The main idea to compute the vibrational potential of a network is embedding vertices to n-dimensional Euclidean space by using the Moore-Penrose pseudo-inverse of the Laplacian, where n = |V|. Within the hierarchical structure of each community, each stock market tends to be adjacent to junction vertices. Therefore, the change of investment on the junction vertices directly affect the corresponding leaves. Therefore, we correspond the increase/decrease coefficients with vibrational potential of the network. However, instead of direct computing vibrational potential of the network, we compute vibrational potential of each vertex respect to its neighborhood graph.

For this purpose we present the vertex displacement in vibrational potential of a vertex respect to its neighboring vertices with (12) $V ({\vec{x}}_{v}) = \frac{k}{2} {\vec{x}}_{v}^{T} L_{N} {\vec{x}}_{v},$ V({\vec x_v}) = {k \over 2}{\vec x_v}^T{L_N}{\vec x_v}, where k is the spring constant, L_N is the Laplacian of the neighboring graph G_N of the vertex v in G, and ${\vec{x}}_{v}$ {\vec x_v} is the vector whose i-th entry is the displacement of v. The mean displacement of the vertex v can be computed with the reverse temperature β as (13) $Δ x_{i} = \sqrt{\int x_{i}^{2} P ({\vec{x}}_{v}) d {\vec{x}}_{v}},$ \Delta {x_i} = \sqrt {\int x_i^2P({{\vec x}_v})d{{\vec x}_v}} , where the probability distribution $P ({\vec{x}}_{v})$ P({\vec x_v}) is (14) $P ({\vec{x}}_{v}) = exp (- \frac{β k}{2}) / \int exp (- \frac{β k}{2}) d {\vec{x}}_{v} .$ P({\vec x_v}) = \exp \left( { - {{\beta k} \over 2}} \right)/\int \exp \left( { - {{\beta k} \over 2}} \right)d{\vec x_v}.

Similarly the displacement correlation of the vertices in same neighborhood can be defined as (15) $< x_{i}, x_{j} > = \int x_{i} x_{j} P ({\vec{x}}_{v}) d {\vec{x}}_{v},$ < {x_i},{x_j} > = \int {x_i}{x_j}P({\vec x_v})d{\vec x_v}, where < ··· > is the thermal average.

Let $0 = λ_{1}^{N} < λ_{2}^{N} \leq \dots \leq λ_{n}^{N}$ 0 = \lambda _1^N < \lambda _2^N \le \ldots \le \lambda _n^N be the spectrum of L_N respect to eigenvalues $λ_{μ}^{N}$ \lambda _\mu ^N . Since the quantity respect to 0 eigenvalue is the center of mass, the 0 eigenvalue does not affect the vertex displacement. Then the integral measure can be transformed by (16) $d {\vec{x}}_{v} = \prod_{i = 1}^{n} d x_{i} = | det U_{N} | \prod_{i = 1}^{n} d ξ_{i} = d {\vec{ξ}}_{v}$ d{\vec x_v} = \prod\limits_{i = 1}^n d{x_i} = |\det {U_N}|\prod\limits_{i = 1}^n d{\xi _i} = d{\vec \xi _v} where U_N is the matrix formed by the orthogonal eigenvectors of L_N. By the introduction of this transform the new probability distribution can be obtained as (17) $\begin{array}{l} P ({\vec{ξ}}_{v}) & = exp (- \frac{β k}{2} {\vec{ξ}}_{v}^{T} Λ_{N} {\vec{ξ}}_{v}) / \int exp (- \frac{β k}{2} {\vec{ξ}}_{v}^{T} Λ_{N} {\vec{ξ}}_{v}) d {\vec{ξ}}_{v} \\ = exp (- \frac{β k}{2} {\vec{ξ}}_{v}^{T} Λ_{N} {\vec{ξ}}_{v}) / \prod_{μ = 1}^{n} \int_{- \infty}^{\infty} exp (- \frac{β k}{2} λ_{μ}^{N} ξ_{μ}^{2}), \end{array}$ \matrix{ {P({{\vec \xi }_v})} \hfill & { = \exp \left( { - {{\beta k} \over 2}{{\vec \xi }_v}^T{\Lambda _N}{{\vec \xi }_v}} \right)/\int \exp \left( { - {{\beta k} \over 2}{{\vec \xi }_v}^T{\Lambda _N}{{\vec \xi }_v}} \right)d{{\vec \xi }_v}} \hfill \cr {} \hfill & { = \exp \left( { - {{\beta k} \over 2}{{\vec \xi }_v}^T{\Lambda _N}{{\vec \xi }_v}} \right)/\prod\limits_{\mu = 1}^n \int_{ - \infty }^\infty \exp \left( { - {{\beta k} \over 2}\lambda _\mu ^N\xi _\mu ^2} \right),} \hfill} where the diagonal matrix Λ_N involves the eigenvalues $λ_{μ}^{N}$ \lambda _\mu ^N .

Since 0 eigenvalue does not effect the vertex displacement, we can remove the component µ = 1 from the Equation 17, and the probability distribution can be computed as (18) $\begin{array}{l} P ({\vec{ξ}}_{v}) & = exp (- \frac{β k}{2} {\vec{ξ}}_{v}^{T} Λ_{N} {\vec{ξ}}_{v}) / \prod_{μ = 2}^{n} \int_{- \infty}^{\infty} exp (- \frac{β k}{2} λ_{μ}^{N} ξ_{μ}^{2}) \\ = exp (- \frac{β k}{2} {\vec{ξ}}_{v}^{T} Λ_{N} {\vec{ξ}}_{v}) / \prod_{μ = 2}^{n} \sqrt{\frac{2 π}{β k λ_{μ}}} . \end{array}$ \matrix{ {P({{\vec \xi }_v})} \hfill & { = \exp \left( { - {{\beta k} \over 2}{{\vec \xi }_v}^T{\Lambda _N}{{\vec \xi }_v}} \right)/\prod\limits_{\mu = 2}^n \int_{ - \infty }^\infty \exp \left( { - {{\beta k} \over 2}\lambda _\mu ^N\xi _\mu ^2} \right)} \hfill \cr {} \hfill & { = \exp \left( { - {{\beta k} \over 2}{{\vec \xi }_v}^T{\Lambda _N}{{\vec \xi }_v}} \right)/\prod\limits_{\mu = 2}^n \sqrt {{{2\pi } \over {\beta k{\lambda _\mu }}}} .} \hfill}

Hence, by using the probability distribution obtained in Equation 18, it is possible to compute the Equation 13 as (19) $Δ x_{i} = \sqrt{< x_{i}^{2} >} = \sqrt{ℒ_{i} / \prod_{μ = 2}^{n} \sqrt{\frac{2 π}{β k λ_{μ}}}},$ \Delta {x_i} = \sqrt { < x_i^2 > } = \sqrt {{{\cal L}_i}/\prod\limits_{\mu = 2}^n \sqrt {{{2\pi } \over {\beta k{\lambda _\mu }}}} } , where $\begin{array}{l} ℒ_{i} & = \sum_{j = 2}^{n} \int_{- \infty}^{\infty} {(U_{ij} ξ_{j})}^{2} exp (- \frac{β k}{2} λ_{j} ξ_{j}^{2}) d ξ_{j} \times \prod_{μ = 2, μ \neq j}^{n} \int_{- \infty}^{\infty} exp (- \frac{β k}{2} λ_{μ} ξ_{μ}^{2}) d ξ_{μ} \\ = \prod_{μ = 2}^{n} \sqrt{\frac{2 π}{β k λ_{μ}}} \times \sum_{j = 2}^{n} \frac{U_{ij}^{2}}{β k λ_{j}} . \end{array}$ \matrix{ {{{\cal L}_i}} \hfill & { = \sum\limits_{j = 2}^n \int_{ - \infty }^\infty {{({U_{ij}}{\xi _j})}^2}\exp \left( { - {{\beta k} \over 2}{\lambda _j}\xi _j^2} \right)d{\xi _j} \times \prod\limits_{\mu = 2,\mu \ne j}^n \int_{ - \infty }^\infty \exp \left( { - {{\beta k} \over 2}{\lambda _\mu }\xi _\mu ^2} \right)d{\xi _\mu }} \hfill \cr {} \hfill & { = \prod\limits_{\mu = 2}^n \sqrt {{{2\pi } \over {\beta k{\lambda _\mu }}}} \times \sum\limits_{j = 2}^n {{U_{ij}^2} \over {\beta k{\lambda _j}}}.} \hfill \cr }

Therefore, the mean displacement of a vertex from its neighborhood can be computed as $Δ x_{i} = \sqrt{\sum_{j = 2}^{n} \frac{U_{ij}^{2}}{β k λ_{j}}} .$ \Delta {x_i} = \sqrt {\sum\limits_{j = 2}^n {{U_{ij}^2} \over {\beta k{\lambda _j}}}} .

By the introduction of the Moore-Penrose pseudo inverse $L_{ij}^{N +}$ L_{ij}^{N + } of L_N as in [59, 60], it is also possible to compute the mean displacement of a vertex from its neighborhood as (20) $Δ x_{i} = \frac{1}{β k} (L_{ii}^{N +}) .$ \Delta {x_i} = {1 \over {\beta k}}(L_{ii}^{N + }).

We also note that the displacement correlation of the vertices in the same neighborhood that is given in the Equation 15 can be computed in the terms of Moore-Penrose pseudo inverse as (21) $< x_{i}, x_{j} > = \frac{1}{β k} (L_{ij}^{N +}) .$ < {x_i},{x_j} > = {1 \over {\beta k}}(L_{ij}^{N + }).

Another parameters we need to estimate in the system 9–11 are γ₁ and γ₂ which are the proportionality values. The proportionality values control how much financial agent in the same community affect each other. Hence, they can be measured as how strong each agents are connected internally. This measurement is naturally arise from the PMFG of each communities. PMFG structure allows cliques which are the topological subgraph structure representing strong relations. Since PMFG also has the information about the hierarchical structure, it is reasonable to measure internal connectedness of communities by using PMFG. For this measurement we follow the way presented in [45]. The mean disparity measurement < y > of a PMFG can be defined as the mean of (22) $y (i) = \sum_{i \neq j, j \in Clique} {(\frac{d_{Corr} (i, j)}{s_{i}})}^{2},$ y(i) = \sum\limits_{i \ne j,j \in Clique} {\left( {{{{d_{Corr}}(i,j)} \over {{s_i}}}} \right)^2}, where i is the generic element of the clique and (23) $s_{i} = \sum_{i \neq j, j \in Clique} d_{Corr} (i, j) .$ {s_i} = \sum\limits_{i \ne j,j \in Clique} {d_{Corr}}(i,j).

4.2

Adomian Decomposition Method

It is well known that many nonlinear differential equations exhibit strange attractors and their solutions have been discovered to move toward strange attractors. If these strange attractors are examined deeply, it can be seen that these are fractals. Therefore, we aim to deal with fractal nonlinear differential equations rather than classical forms of them. Hence we shall extend the Adomian decomposition method to be used for solving fractional nonlinear equations. For the solution of the system 9–11, we use an efficient decomposition method for approximating the solution of systems of fractional integro-differential equation that are given in Caputo sense. The approximate solutions are calculated in the terms of a convergent series as in [34].

Now let us consider the system 9–11 with 0 < α 1. By following the decomposition idea we may state that (24) $D_{0}^{α} μ_{1} (t) = m_{1} (t),$ D_0^\alpha {\mu _1}(t) = {m_1}(t),(25) $D_{0}^{α} μ_{2} (t) = m_{2} (t) .$ D_0^\alpha {\mu _2}(t) = {m_2}(t).

This equations lead us to the integral equations (26) $μ_{1} (t) = μ_{1} (0) + J^{α} (m_{1} (t))$ {\mu _1}(t) = {\mu _1}(0) + {J^\alpha }({m_1}(t))(27) $= N_{1} + J_{0}^{α} m_{1} (s) ds,$ = {N_1} + J_0^\alpha {m_1}(s)ds,(28) $μ_{2} (t) = μ_{2} (0) + J^{α} (m_{2} (t))$ {\mu _2}(t) = {\mu _2}(0) + {J^\alpha }({m_2}(t))(29) $= N_{2} + J_{0}^{α} m_{2} (s) ds,$ = {N_2} + J_0^\alpha {m_2}(s)ds,(30) $\begin{array}{l} m_{1} (t) & = μ_{1} (t) (k_{1} - γ_{1} μ_{2} (t) - \int_{t - T_{0}}^{t} f_{1} (t - s) μ_{2} (s) ds) \\ = k_{1} (N_{1} + \int_{0}^{t} m_{1} (s) ds) - μ_{1} (t) (γ_{1} μ_{2} (t) - \int_{t - T_{0}}^{t} f_{1} (t - s) μ_{2} (s) ds), \end{array}$ \matrix{ {{m_1}(t)} \hfill & { = {\mu _1}(t)\left( {{k_1} - {\gamma _1}{\mu _2}(t) - \int_{t - {T_0}}^t {f_1}(t - s){\mu _2}(s)ds} \right)} \hfill \cr {} \hfill & { = {k_1}\left( {{N_1} + \int_0^t {m_1}(s)ds} \right) - {\mu _1}(t)\left( {{\gamma _1}{\mu _2}(t) - \int_{t - {T_0}}^t {f_1}(t - s){\mu _2}(s)ds} \right),} \hfill}(31) $\begin{array}{l} m_{2} (t) & = μ_{2} (t) (- k_{2} + γ_{2} μ_{1} (t) - \int_{t - T_{0}}^{t} f_{2} (t - s) μ_{1} (s) ds) \\ = k_{2} (N_{2} + \int_{0}^{t} m_{2} (s) ds) + μ_{2} (t) (γ_{2} μ_{1} (t) + \int_{t - T_{0}}^{t} f_{2} (t - s) μ_{1} (s) ds) . \end{array}$ \matrix{ {{m_2}(t)} \hfill & { = {\mu _2}(t)\left( { - {k_2} + {\gamma _2}{\mu _1}(t) - \int_{t - {T_0}}^t {f_2}(t - s){\mu _1}(s)ds} \right)} \hfill \cr {} \hfill & { = {k_2}\left( {{N_2} + \int_0^t {m_2}(s)ds} \right) + {\mu _2}(t)\left( {{\gamma _2}{\mu _1}(t) + \int_{t - {T_0}}^t {f_2}(t - s){\mu _1}(s)ds} \right).} \hfill}

Afterwards, the Adomian process will be as follows: (32) $μ_{1, 0} = N_{1}, μ_{2, 0} = N_{2},$ {\mu _{1,0}} = {N_1},\;{\mu _{2,0}} = {N_2},(33) $m_{1, 0} = k_{1} N_{1}, m_{2, 0} = k_{2} N_{2},$ {m_{1,0}} = {k_1}{N_1},\;{m_{2,0}} = {k_2}{N_2},(34) $μ_{1, j + 1} = J_{0}^{α} m_{1, j} (s) ds, μ_{2, j + 1} = J_{0}^{α} m_{2, j} (s) ds,$ {\mu _{1,j + 1}} = J_0^\alpha {m_{1,j}}(s)ds,\;{\mu _{2,j + 1}} = J_0^\alpha {m_{2,j}}(s)ds,(35) $\begin{array}{l} m_{1, j + 1} & = k_{1} \int_{0}^{t} m_{1, j} (s) ds - γ_{1} \sum_{k = 0}^{j} μ_{1, k} (t) μ_{2, j - k} (t) \\ - \int_{t - T_{0}}^{t} f_{1} (t - s) (\sum_{k = 0}^{j} μ_{1, k} (t) μ_{2, j - k} (t)) ds \end{array}$ \matrix{ {{m_{1,j + 1}}} \hfill & { = {k_1}\int_0^t {m_{1,j}}(s)ds - {\gamma _1}\sum\limits_{k = 0}^j {\mu _{1,k}}(t){\mu _{2,j - k}}(t)} \hfill \cr {} \hfill & { - \int_{t - {T_0}}^t {f_1}(t - s)\left( {\sum\limits_{k = 0}^j {\mu _{1,k}}(t){\mu _{2,j - k}}(t)} \right)ds} \hfill}(36) $\begin{array}{l} m_{2, j + 1} & = k_{2} \int_{0}^{t} m_{2, j} (s) ds - γ_{2} \sum_{k = 0}^{j} μ_{1, k} (t) μ_{2, j - k} (t) \\ + \int_{t - T_{0}}^{t} f_{2} (t - s) (\sum_{k = 0}^{j} μ_{1, k} (t) μ_{2, j - k} (t)) ds . \end{array}$ \matrix{ {{m_{2,j + 1}}} \hfill & { = {k_2}\int_0^t {m_{2,j}}(s)ds - {\gamma _2}\sum\limits_{k = 0}^j {\mu _{1,k}}(t){\mu _{2,j - k}}(t)} \hfill \cr {} \hfill & { + \int_{t - {T_0}}^t {f_2}(t - s)\left( {\sum\limits_{k = 0}^j {\mu _{1,k}}(t){\mu _{2,j - k}}(t)} \right)ds.} \hfill}

Results

In order to study the proposed model in the Borsa Istanbul Stock Exchange, we first apply our algorithm to the data set to obtain stock market network. For the fraction size h = 10000, the algorithm determines the control parameter as 0.6854. The vertices are sorted from 1 to 93 by the alphabetical order in Table 1. The formed network is presented in Figure 1. The vertices with maximum vertex degree are ADEL, AKBNK, AKSEN, ALBRK, ALGYO, ALKIM, ARCLK, ASELS, BRISA, DOAS, GARAN, GOLTS, HALKB, ISCTR, KARTN, KCHOL, KONYA, KRDMD, MGROS, OTKAR, PRKME, SAHOL, SISE, SNGYO, TKFEN, TKNSA, TMSN, TOASO, TRCAS, VAKBN, and YKBNK with the degree number 92, and the vertex with minimum vertex degree is VKGYO. The maximum and mean correlation distances among the companies are 0.7211 and 0.5798, respectively. Hence it can also be concluded this network has strong internal connectedness. The correlation distance matrix of the agents are presented in Figure 2 with temperature mapping.

The network of XU 100 with the threshold value 0.6854.

The matrix of correlation distance among the financial agents of XU 100. The darker points are closer to 1 while the lighter ones are closer to 0.

To determine the two non-overlapping clusters of financial actors we use the high modularity method. The resulting communities are presented in Figure 3. The agent numbers of each community give us the initial conditions as N₁ = 66 and N = 27.

The community plot of the resulting network.

As aforementioned, parameters of the model described by the system 9–11 are obtained by vibrational potential respect to neighborhood graphs and mean disparity measures of each communities. We need to note that vibrational potentials of each vertices are tend to form internal clusters; i.e., some of them have higher values and some of them have lower values. Therefore, while determining k₁ and k₂ values, we choose the mean value of each vibrational potentials. Afterwards forming the PMFG of each community, it becomes possible to obtain disparity measures respect to 4-cliques, which are the representation of the strongest connections. The resulting parameters γ₁ = 0.3342 and γ₂ = 0.3388. The parameters are close to the value 1/3 which also states that the clustering method we choose is reasonable [45]. To interpret the results, we present MSTs and PMFGs of both two clusters in Figures 4–5

MST (above) and PMFG (below) filtering of Community 1.

MST (above) and PMFG (below) filtering of Community 2.

In the light of these computed parameters we may now state the system 9–11 with T₀=100, f₁(t) = f₂(t) = e^t as (37) $D_{0}^{α} μ_{1} (t) = μ_{1} (t) [5.95 - 0.3342 μ_{2} (t) - \int_{t - 100}^{t} e^{(t - s)} μ_{2} (s) ds],$ D_0^\alpha {\mu _1}(t) = {\mu _1}(t)\left[ {5.95 - 0.3342{\mu _2}(t) - \int_{t - 100}^t {e^{(t - s)}}{\mu _2}(s)\;ds} \right],(38) $D_{0}^{α} μ_{2} (t) = μ_{2} (t) [- 4.39 + 0.3388 μ_{1} (t) + \int_{t - 100}^{t} e^{(t - s)} μ_{1} (s) ds],$ D_0^\alpha {\mu _2}(t) = {\mu _2}(t)\left[ { - 4.39 + 0.3388{\mu _1}(t) + \int_{t - 100}^t {e^{(t - s)}}{\mu _1}(s)\;ds} \right],(39) $μ_{1} (0) = 66, μ_{2} (0) = 27 .$ {\mu _1}(0) = 66,\,\,\,\,{\mu _2}(0) = 27.

By applying the Adomian process obtained in Section 4.2, the solution of the initial value problem 37–39 can be obtained as (40) $μ_{1} (t) = N_{1} + \frac{k_{1} N_{1} t^{α}}{Γ (1 + α)} + \frac{N_{1} t^{α} ((- 1 + e^{- T_{0}}) N_{2} + 4^{- α} k_{1}^{2} \sqrt{π} t^{α} / Γ (0.5 + α) - N_{2} γ_{1})}{Γ (1 + α)},$ {\mu _1}(t) = {N_1} + {{{k_1}{N_1}{t^\alpha }} \over {\Gamma (1 + \alpha )}} + {{{N_1}{t^\alpha }\left( {\left( { - 1 + {e^{ - {T_0}}}} \right){N_2} + {4^{ - \alpha }}k_1^2\sqrt \pi {t^\alpha }/\Gamma (0.5 + \alpha ) - {N_2}{\gamma _1}} \right)} \over {\Gamma (1 + \alpha )}},(41) $μ_{2} (t) = N_{2} + \frac{k_{2} N_{2} t^{α}}{Γ (1 + α)} + \frac{N_{2} t^{α} ((1 - e^{- T_{0}}) N_{1} + 4^{- α} k_{2}^{2} \sqrt{π} t^{α} / Γ (0.5 + α) - N_{1} γ_{2})}{Γ (1 + α)} .$ {\mu _2}(t) = {N_2} + {{{k_2}{N_2}{t^\alpha }} \over {\Gamma (1 + \alpha )}} + {{{N_2}{t^\alpha }\left( {\left( {1 - {e^{ - {T_0}}}} \right){N_1} + {4^{ - \alpha }}k_2^2\sqrt \pi {t^\alpha }/\Gamma (0.5 + \alpha ) - {N_1}{\gamma _2}} \right)} \over {\Gamma (1 + \alpha )}}. with a three-term approximation.

The plots of the solution functions (40–41) are presented in Figures 6–15 for the different α values.

The solutions of the initial value problem 40–41 for α = 0.1

The solutions of the initial value problem 40–41 for α = 0.2

The solutions of the initial value problem 40–41 for α = 0.3

The solutions of the initial value problem 40–41 for α = 0.4

The solutions of the initial value problem 40–41 for α = 0.5

The solutions of the initial value problem 40–41 for α = 0.6

The solutions of the initial value problem 40–41 for α = 0.7

The solutions of the initial value problem 40–41 for α = 0.8

The solutions of the initial value problem 40–41 for α = 0.9

The solutions of the initial value problem 40–41 for α = 0.10

Conclusions

Ordinary differential equations are the most common mathematical tool to represent real world problems. But ordinary differential equations become less effective whenever the problem involves memory effect. Complex systems that representing financial agents have the memory effect, hence it is reasonable to model such systems by using the idea of fractional differential.

In this study, we propose a model which is represents the fractional interaction of financial agents. The interaction of the agents is determined within a complex network of a stock market. We express the model as a system of fractional integro-differential equations in Caputo sense. Hence, we keep the fading memory of the financial interaction. Our model considers two clusters of agents where one cluster tends to get investment flow. To determine the clusters we use maximization of the edge modularity in stock market network. The resulting clusters are consistent with the structure of Borsa Istanbul. Both MST and PMFG filtration of the clusters involve agents of Financials sector as the leading elements. To estimate the parameters of the model, we use graph theoretical concepts such as vibrational potentials and disparity measure of respected PMFGs.

By the computed parameters, we use Adomian decomposition method to obtain a solution of the model. This solution show us that for different fractional dimensions α, the model always reaches to an equilibrium state. For the lesser values of fraction rate α, agents reach to an equilibrium state relatively slower. Besides, the flows of investments tend to be in same characteristics. For the greater α values, agents reach to an equilibrium state faster and similarly the flows of investments tend to be in same characteristics. The model keeps the memory of the investment in best for 0.4 ≤ α ≤ 0.6. This results shows us that the fractional interaction of financial agents is consistent with reality when autocorrelations are discarded.

As neoclassical liberal theory of economics states, markets always seek for an equilibrium state. Hence, the model we present with fractional derivative is consistent with the real data of Borsa Istanbul Stock Exchange Market. We also believe that these kind of models can provide useful information for understanding and prediction of the global economic crisis.

Figures & Tables

Sectors of each considered stock

FinancialsAKBNK, SKBNK, SNGYO, TSKB, TEKST, TRGYO, VKGYO, ALGYO, ISGYO, GARAN, ALBRK, GLYHO, ISCTR, YKBNK, SAHOL, GOZDE, HALKB, VAKBN, ECZYT, SAFGY, EKGYO, SAHOL, GSDHO

IndustrialsASELS, TAVHL, TKFEN, TTRAK, CLEBI

Consumer DiscretionaryASUZU, TKNSA, TOASO, YAZIC, AKSA, ARCLK, GSRAY, KARSN, THYAO, BRISA, DOAS, FENER, MNDRS, METRO, VESBE, ADEL, BJKAS, NTTUR, GOODY, OTKAR, TMSN, EGEEN, FROTO, IHLAS

EnergyAYGAZ, TUPRS, IPEKE, KCHOL

TechnologyNETAS, VESTL

MaterialsSASA, AFYON, ANACM, BAGFS, CIMSA, KONYA, KOZAA, ERBOS, KRDMD, PRKME, SISE, ALKIM, TRKCM, GUBRF, KOZAL, BRSAN, KARTN, PETKM, GOLTS, EREGL

CommunicationsTTKOM, TCELL, DOHOL, HURGZ

Consumer StaplesAEFES, CCOLA, BIZIM, ECILC, BIMAS, MGROS, SODA, ULKER

UtilitiesAKSEN, ALARK, TRCAS, ZOREN, ENKAI

References

[1]Haken H., Jumarie G. (2006) A macroscopic approach to complex system, Springer-Verlag, Berlin, Heidelberg, New York, 3th EditionHakenH.JumarieG.2006A macroscopic approach to complex systemSpringer-VerlagBerlin, Heidelberg, New York3th EditionSearch in Google Scholar

[2]Milo R., Shen-Orr S. (2002) Itzkovitz S., Kashtan N., Chklovskii D., Alon U., Network motifs: simple building blocks of complex networks. Science, 298(5594): 824–827MiloR.Shen-OrrS.ItzkovitzS.KashtanN.ChklovskiiD.AlonU.2002Network motifs: simple building blocks of complex networksScience298559482482710.1515/9781400841356.217Search in Google Scholar

[3]Castellano C., Pastor-Storras R. (2010) Thresholds for epidemic spreading in networks, Phys. Rev. Kett., 105(21): 218701CastellanoC.Pastor-StorrasR.2010Thresholds for epidemic spreading in networksPhys. Rev. Kett.1052121870110.1103/PhysRevLett.105.21870121231361Search in Google Scholar

[4]Read J. M., Eames K. T., Edmunds W. J. (2008) Dynamic social networks and the implications for the spread of infectious disease, J. R. Soc. Interface, 6–6(26): 1001–1007ReadJ. M.EamesK. T.EdmundsW. J.2008Dynamic social networks and the implications for the spread of infectious diseaseJ. R. Soc. Interface6–6261001100710.1098/rsif.2008.0013260743318319209Search in Google Scholar

[5]Keeling M. J. (2008) Rohani P, Modelling Infectious Diaseases in Human and Animals, Princeton University PressKeelingM. J.RohaniP2008Modelling Infectious Diaseases in Human and AnimalsPrinceton University Press10.1515/9781400841035Search in Google Scholar

[6]Busenberg S., Martelli M. (Eds.) (1990) Differential Equations Models in Biology, Epidemiology and Ecology, Proceedings of a Conference Held in Claremont California, January 13–16 (Vol. 92). Springer Science and Business Media, 2002BusenbergS.MartelliM.(Eds.)1990Differential Equations Models in Biology, Epidemiology and EcologyProceedings of a Conference Held in Claremont CaliforniaJanuary 13–1692Springer Science and Business Media, 2002Search in Google Scholar

[7]Jerri A. J. (1999) Introduction to integral equations with applicationsJerriA. J.1999Introduction to integral equations with applicationsSearch in Google Scholar

[8]Rosenthal J., Gilliam D. S. (Eds.) (2012) Mathematical systems theory in biology, communications, computation and finance (Vol. 134). Springer Science and Business MediaRosenthalJ.GilliamD. S.(Eds.)2012Mathematical systems theory in biology, communications, computation and finance134Springer Science and Business MediaSearch in Google Scholar

[9]Capasso V., Bakstein D. (2015) An introduction to continuous-time stochastic processes: theory, models, and applications to finance, biology, and medicine. BirkhäuserCapassoV.BaksteinD.2015An introduction to continuous-time stochastic processes: theory, models, and applications to finance, biology, and medicineBirkhäuser10.1007/978-1-4939-2757-9Search in Google Scholar

[10]Hens T., Schenk-Hoppé K. R. (Eds.) (2009) Handbook of financial markets: dynamics and evolution. ElsevierHensT.Schenk-HoppéK. R.(Eds.)2009Handbook of financial markets: dynamics and evolutionElsevierSearch in Google Scholar

[11]Deffuant G., Neau D., Amblard F. (2000) Weisbuch G, Mixing beliefs among interacting agents. Advances in Complex Systems, 3(01n04): 87–98DeffuantG.NeauD.AmblardF.WeisbuchG2000Mixing beliefs among interacting agentsAdvances in Complex Systems3(01n04):879810.1142/S0219525900000078Search in Google Scholar

[12]Samanidou E., Zschischang E., Stauffer D., Lux T. (2007) Agent-based models of financial markets. Reports on Progress in Physics, 70(3): 409.SamanidouE.ZschischangE.StaufferD.LuxT.2007Agent-based models of financial marketsReports on Progress in Physics70340910.1088/0034-4885/70/3/R03Search in Google Scholar

[13]Balcı M. A. (2016) Fractional virus epidemic model on financial networks. Open Mathematics, 4(1): 1074–1086BalcıM. A.2016Fractional virus epidemic model on financial networksOpen Mathematics411074108610.1515/math-2016-0098Search in Google Scholar

[14]Onnela J. P., Kaski K., Kertész J. (2004) Clustering and information in correlation based financial networks. The European Physical Journal B, 38(2): 353–362OnnelaJ. P.KaskiK.KertészJ.2004Clustering and information in correlation based financial networksThe European Physical Journal B38235336210.1140/epjb/e2004-00128-7Search in Google Scholar

[15]Kumar S., Deo N. (2012) Correlation and network analysis of global financial indices. Physical Review E, 86(2): 026101KumarS.DeoN.2012Correlation and network analysis of global financial indicesPhysical Review E86202610110.1103/PhysRevE.86.02610123005819Search in Google Scholar

[16]Mizuno T., Takayasu H., Takayasu M. (2006) Correlation networks among currencies. Physica A: Statistical Mechanics and its Applications, 364:336–342MizunoT.TakayasuH.TakayasuM.2006Correlation networks among currenciesPhysica A: Statistical Mechanics and its Applications36433634210.1016/j.physa.2005.08.079Search in Google Scholar

[17]Balcı M. A., Akgüller, Ö. (2016) Soft Vibrational Force on Stock Market Networks. Library Journal, 3, e3050.BalcıM. A.AkgüllerÖ.2016Soft Vibrational Force on Stock Market NetworksLibrary Journal3e3050Search in Google Scholar

[18]Akgüller Ö., Balcı, M.A. (2018) Geodetic convex boundary curvatures of the communities in stock market networks. Physica A: Statistical Mechanics and its Applications, 505, pp.569–581.AkgüllerÖ.BalcıM.A.2018Geodetic convex boundary curvatures of the communities in stock market networksPhysica A: Statistical Mechanics and its Applications50556958110.1016/j.physa.2018.03.087Search in Google Scholar

[19]Akgüller, Ö (2017) Geometric Soft Sets. Hittite Journal of Science & Engineering, 4(2).AkgüllerÖ2017Geometric Soft SetsHittite Journal of Science & Engineering4210.17350/HJSE19030000063Search in Google Scholar

[20]Podlubny I. (1999) Fractional Differential Equation, Academic Press, New YorkPodlubnyI.1999Fractional Differential EquationAcademic PressNew YorkSearch in Google Scholar

[21]Gorenflo R., Mainardi F. (1997) Fractional calculus: integral and differential equations of fractional order, Springer, New YorkGorenfloR.MainardiF.1997Fractional calculus: integral and differential equations of fractional orderSpringerNew York10.1007/978-3-7091-2664-6_5Search in Google Scholar

[22]Wang C. C. (1965) The principle of fading memory. Archive for Rational Mechanics and Analysis, 18(5): 343–366WangC. C.1965The principle of fading memoryArchive for Rational Mechanics and Analysis18534336610.1007/BF00281325Search in Google Scholar

[23]Tarasova V. V., Tarasov V. E. (2018) Concept of dynamic memory in economics. Communications in Nonlinear Science and Numerical Simulations, 55:127–145TarasovaV. V.TarasovV. E.2018Concept of dynamic memory in economicsCommunications in Nonlinear Science and Numerical Simulations5512714510.1016/j.cnsns.2017.06.032Search in Google Scholar

[24]David S. A., Fischer C., Machado J. T. (2018) Fractional electronic circuit simulation of a nonlinear macroeconomic model. AEU-International Journal of Electronics and Communications, 84: 210–220.DavidS. A.FischerC.MachadoJ. T.2018Fractional electronic circuit simulation of a nonlinear macroeconomic modelAEU-International Journal of Electronics and Communications8421022010.1016/j.aeue.2017.11.019Search in Google Scholar

[25]Babenkov M. B., Vitokhin E. Y. (2017) Thermoelastic Waves in a Medium with Heat-Flux Relaxation. Encyclopedia of Continuum Mechanics, 1–10BabenkovM. B.VitokhinE. Y.2017Thermoelastic Waves in a Medium with Heat-Flux RelaxationEncyclopedia of Continuum Mechanics110Search in Google Scholar

[26]Bas E., Acay B., Ozarslan R. (2019). The price adjustment equation with different types of conformable derivatives in market equilibrium. AIMS Mathematics, 4(3): 805–820.BasE.AcayB.OzarslanR.2019The price adjustment equation with different types of conformable derivatives in market equilibriumAIMS Mathematics4380582010.3934/math.2019.3.805Search in Google Scholar

[27]Ozarslan R., Ercan A., Bas E. (2019). Novel Fractional Models Compatible with Real World Problems. Fractal and Fractional, 3(2), 15.OzarslanR.ErcanA.BasE.2019Novel Fractional Models Compatible with Real World ProblemsFractal and Fractional321510.3390/fractalfract3020015Search in Google Scholar

[28]Bas E., Acay B., Ozarslan R. (2019). Fractional models with singular and non-singular kernels for energy efficient buildings. Chaos: An Interdisciplinary Journal of Nonlinear Science, 29(2), 023110.BasE.AcayB.OzarslanR.2019Fractional models with singular and non-singular kernels for energy efficient buildingsChaos: An Interdisciplinary Journal of Nonlinear Science29202311010.1063/1.508239030823708Search in Google Scholar

[29]Bas E., Ozarslan R., Baleanu D., Ercan A. (2018). Comparative simulations for solutions of fractional Sturm–Liouville problems with non-singular operators. Advances in Difference Equations, 2018(1), 350.BasE.OzarslanR.BaleanuD.ErcanA.2018Comparative simulations for solutions of fractional Sturm–Liouville problems with non-singular operatorsAdvances in Difference Equations2018135010.1186/s13662-018-1803-8Search in Google Scholar

[30]Bas E., Ozarslan R. (2018). Real world applications of fractional models by Atangana aleanu fractional derivative. Chaos, Solitons & Fractals, 116, 121–125.BasE.OzarslanR.2018Real world applications of fractional models by Atangana aleanu fractional derivativeChaos, Solitons & Fractals11612112510.1016/j.chaos.2018.09.019Search in Google Scholar

[31]Shakeri F., Dehghan M. (2008) Solution of a model describing biological species living together using the variational iteration method. Mathematical and Computer Modelling, 48(5): 685–699ShakeriF.DehghanM.2008Solution of a model describing biological species living together using the variational iteration methodMathematical and Computer Modelling48568569910.1016/j.mcm.2007.11.012Search in Google Scholar

[32]Babolian E., Biazar J. (2002) Solving the problem of biological species living together by Adomian decomposition method. Applied Mathematics and Computation, 129(2): 339–343BabolianE.BiazarJ.2002Solving the problem of biological species living together by Adomian decomposition methodApplied Mathematics and Computation129233934310.1016/S0096-3003(01)00043-1Search in Google Scholar

[33]Yousefi S. A. (2011) Numerical solution of a model describing biological species living together by using Legendre multiwavelet method. International Journal of Nonlinear Science, 11(1): 109–113YousefiS. A.2011Numerical solution of a model describing biological species living together by using Legendre multiwavelet methodInternational Journal of Nonlinear Science111109113Search in Google Scholar

[34]Momani S., Qaralleh R. (2006) An Efficient Method for Solving Sysyems of Fractional Integro-Differential Equations, Computers and Mathematics with Applications, 52: 459–470MomaniS.QarallehR.2006An Efficient Method for Solving Sysyems of Fractional Integro-Differential EquationsComputers and Mathematics with Applications5245947010.1016/j.camwa.2006.02.011Search in Google Scholar

[35]Luchko A., Groneflo R. (1998) The initial value problem for some fractional differential equations with the Capoto derivative, Preprint Series A08-98, FAchbreich Mathematik und Informatik, Freie Universitat BerlinLuchkoA.GronefloR.1998The initial value problem for some fractional differential equations with the Capoto derivativePreprint Series A08-98,FAchbreich Mathematik und Informatik, Freie Universitat BerlinSearch in Google Scholar

[36]Miller K. S., Ross B. (1993) An Introduction to the Fractional Calculus and Fractional Differential Equations, John Wiley and Sons, New YorkMillerK. S.RossB.1993An Introduction to the Fractional Calculus and Fractional Differential EquationsJohn Wiley and SonsNew YorkSearch in Google Scholar

[37]Oldhami K. B., Spanier J. (1974) The Fractional Calculus, Academic Press, New YorkOldhamiK. B.SpanierJ.1974The Fractional CalculusAcademic PressNew YorkSearch in Google Scholar

[38]Ren R. F., Li H. B., Jiang W., Song M. Y. (2013) An efficient Chebyshev-tau method for solving the space fractional diffusion equations. Applied Mathematics and Computation, 224: 259–267RenR. F.LiH. B.JiangW.SongM. Y.2013An efficient Chebyshev-tau method for solving the space fractional diffusion equationsApplied Mathematics and Computation22425926710.1016/j.amc.2013.08.073Search in Google Scholar

[39]Van Dam E. R., Haemers W. H. (2003) Which graphs are determined by their spectrum?, Linear Algebra Appl., 373: 241–272Van DamE. R.HaemersW. H.2003Which graphs are determined by their spectrum?Linear Algebra Appl.37324127210.1016/S0024-3795(03)00483-XSearch in Google Scholar

[40]Cyman J., Lemańska M., Raczek J. (2006) On the doubly connected domination number of a graph, Open Mathematics, 4(1): 34–45CymanJ.LemańskaM.RaczekJ.2006On the doubly connected domination number of a graphOpen Mathematics41344510.1007/s11533-005-0003-4Search in Google Scholar

[41]Graham R. L., Hell P. (1985) On the history of the minimum spanning tree problem, Annals of the History of Computing, 7(1): 43–57GrahamR. L.HellP.1985On the history of the minimum spanning tree problemAnnals of the History of Computing71435710.1109/MAHC.1985.10011Search in Google Scholar

[42]Mantegna R. N. (1999) Hierarchical structure in financial markets. The European Physical Journal B-Condensed Matter and Complex Systems, 11(1):193–197MantegnaR. N.1999Hierarchical structure in financial marketsThe European Physical Journal B-Condensed Matter and Complex Systems11119319710.1007/s100510050929Search in Google Scholar

[43]Naylor M. J., Rose L. C., Moyle B. J. (2007) Topology of foreign exchange markets using hierarchical structure methods. Physica A: Statistical Mechanics and its Applications, 382(1): 199–208NaylorM. J.RoseL. C.MoyleB. J.2007Topology of foreign exchange markets using hierarchical structure methodsPhysica A: Statistical Mechanics and its Applications382119920810.1016/j.physa.2007.02.019Search in Google Scholar

[44]Balcı, M. A. (2018). Hierarchies in communities of Borsa Istanbul Stock Exchange. Hacettepe Journal of Mathematics and Statistics, 47(4), 921–936.BalcıM. A.2018Hierarchies in communities of Borsa Istanbul Stock ExchangeHacettepe Journal of Mathematics and Statistics474921936Search in Google Scholar

[45]Tumminello M., Aste T., Di Matteo T., Mantegna R. N. (2005) A tool for filtering information in complex systems. Proceedings of the National Academy of Sciences of the United States of America, 102(30): 10421–10426TumminelloM.AsteT.Di MatteoT.MantegnaR. N.2005A tool for filtering information in complex systemsProceedings of the National Academy of Sciences of the United States of America10230104211042610.1073/pnas.0500298102118075416027373Search in Google Scholar

[46]Barfuss W., Massara G. P., Di Matteo T., Aste T. (2016) Parsimonious modeling with information filtering networks. Physical Review E, 94(6): 062306BarfussW.MassaraG. P.Di MatteoT.AsteT.2016Parsimonious modeling with information filtering networksPhysical Review E94606230610.1103/PhysRevE.94.06230628085404Search in Google Scholar

[47]Wang G. J., Xie C., He K., Stanley H. E. (2017) Extreme risk spillover network: application to financial institutions. Quantitative Finance, 1–17WangG. J.XieC.HeK.StanleyH. E.2017Extreme risk spillover network: application to financial institutionsQuantitative Finance11710.1080/14697688.2016.1272762Search in Google Scholar

[48]Aksu M., Kosedag A. (2006) Transparency and disclosure scores and their determinants in the Istanbul Stock Exchange. Corporate Governance: An International Review, 14(4): 277–296.AksuM.KosedagA.2006Transparency and disclosure scores and their determinants in the Istanbul Stock ExchangeCorporate Governance: An International Review14427729610.1111/j.1467-8683.2006.00507.xSearch in Google Scholar

[49]Clauset A., Moore C., Newman M. E. J. (2008) Hierarchical structure and the prediction of missing links in networks. Nature, 453(7191): 98–101ClausetA.MooreC.NewmanM. E. J.2008Hierarchical structure and the prediction of missing links in networksNature45371919810110.1038/nature0683018451861Search in Google Scholar

[50]Lü L., Zhou T. (2011) Link prediction in complex networks: A survey. Physica A: statistical mechanics and its applications, 390(6): 1150–1170LüL.ZhouT.2011Link prediction in complex networks: A surveyPhysica A: statistical mechanics and its applications39061150117010.1016/j.physa.2010.11.027Search in Google Scholar

[51]Newman M. E. J. (2004) Detecting community structure in network, The European Physical Journal B-Condensed Matter and Complex Systems, 38(2): 321–330NewmanM. E. J.2004Detecting community structure in networkThe European Physical Journal B-Condensed Matter and Complex Systems38232133010.1140/epjb/e2004-00124-ySearch in Google Scholar

[52]Lancichinetti A., Fortunato S., Kertész J. (2009) Detecting the overlapping and hierarchical community structure in complex networks. New Journal of Physics, 11(3): 033015LancichinettiA.FortunatoS.KertészJ.2009Detecting the overlapping and hierarchical community structure in complex networksNew Journal of Physics11303301510.1088/1367-2630/11/3/033015Search in Google Scholar

[53]Newman M. E. J., Girvan M. (2004) Finding and evaluating community structure in networks. Physical review E, 69(2): 026113NewmanM. E. J.GirvanM.2004Finding and evaluating community structure in networksPhysical review E69202611310.1103/PhysRevE.69.02611314995526Search in Google Scholar

[54]Agarwal G., Kempe D. (2008) Modularity-maximizing graph communities via mathematical programming. The European Physical Journal B, 66(3): 409–418AgarwalG.KempeD.2008Modularity-maximizing graph communities via mathematical programmingThe European Physical Journal B66340941810.1140/epjb/e2008-00425-1Search in Google Scholar

[55]Estrada E., Hatano N. (2010) A vibrational approach to node centrality and vulnerability in complex networks. Physica A: Statistical Mechanics and its Applications, 389(17): 3648–3660EstradaE.HatanoN.2010A vibrational approach to node centrality and vulnerability in complex networksPhysica A: Statistical Mechanics and its Applications389173648366010.1016/j.physa.2010.03.030Search in Google Scholar

[56]Estrada E., Hatano N., Benzi M. (2012) The physics of communicability in complex networks. Physics reports, 514(3): 89–119EstradaE.HatanoN.BenziM.2012The physics of communicability in complex networksPhysics reports51438911910.1016/j.physrep.2012.01.006Search in Google Scholar

[57]Ranjan G., Zhang Z. L. (2013) Geometry of complex networks and topological centrality. Physica A: Statistical Mechanics and its Applications, 392(17): 3833–3845RanjanG.ZhangZ. L.2013Geometry of complex networks and topological centralityPhysica A: Statistical Mechanics and its Applications392173833384510.1016/j.physa.2013.04.013Search in Google Scholar

[58]Feng L., Bhan B. (2015) Understanding dynamic social grouping behaviors of pedestrians. IEEE Journal of Selected Topics in Signal Processing, 9(2): 317–329FengL.BhanB.2015Understanding dynamic social grouping behaviors of pedestriansIEEE Journal of Selected Topics in Signal Processing9231732910.1109/JSTSP.2014.2365765Search in Google Scholar

[59]Davidson J., Savaliya S., Shah J. J. (2013) Least-squares fit of measured points for square line-profiles, Procedia CIRP, 10: 203–210DavidsonJ.SavaliyaS.ShahJ. J.2013Least-squares fit of measured points for square line-profilesProcedia CIRP1020321010.1016/j.procir.2013.08.032Search in Google Scholar

[60]Grant W. S., Voorhies R. C., Itti L. (2013) Finding planes in LiDAR point clouds for real-time registration, In Intelligent Robots and Systems (IROS), 2013 IEEE/RSJ International Conference, 4347–4354GrantW. S.VoorhiesR. C.IttiL.2013Finding planes in LiDAR point clouds for real-time registrationIn Intelligent Robots and Systems (IROS), 2013 IEEE/RSJ International Conference4347435410.1109/IROS.2013.6696980Search in Google Scholar

Applied Mathematics and Nonlinear Sciences

Fractional Interaction of Financial Agents in a Stock Market Network

Mehmet Ali BalcıMehmet Ali Balcı

Published Online: 31 Mar 2020

Volume & Issue: Volume 5 (2020) - Issue 1 (January 2020)

Page range: 317 - 336

Received: 14 Jun 2019

Accepted: 20 Aug 2019

Journal & Issue Details

Applied Mathematics and Nonlinear Sciences

PDF Preview

Article

Fig. 1

Fig. 2

Fig. 3

Fig. 4

Fig. 5

Fig. 6

Fig. 7

Fig. 8

Fig. 9

Fig. 10

Fig. 11

Fig. 12

Fig. 13

Fig. 14

Fig. 15

Figures & Tables

Fig. 1

Fig. 2

Fig. 3

Fig. 4

Fig. 5

Fig. 6

Fig. 7

Fig. 8

Fig. 9

Fig. 10

Fig. 11

Fig. 12

Fig. 13

Fig. 14

Fig. 15

Sectors of each considered stock

References

Recent Articles

Plan your remote conference with Sciendo

Mehmet Ali Balcı
Mehmet Ali Balcı