Next Article in Journal
On the Extremality of Harmonic Beltrami Coefficients
Next Article in Special Issue
Analytic Valuation Formula for American Strangle Option in the Mean-Reversion Environment
Previous Article in Journal
Matrix Summability of Walsh–Fourier Series
Previous Article in Special Issue
A Modified Black-Scholes-Merton Model for Option Pricing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fourier Integral Operator Model of Market Liquidity: The Chinese Experience 2009–2010

School of Business Administration, Anglo-American University, 118 00 Prague, Czech Republic
Mathematics 2022, 10(14), 2459; https://doi.org/10.3390/math10142459
Submission received: 8 May 2022 / Revised: 30 June 2022 / Accepted: 5 July 2022 / Published: 14 July 2022
(This article belongs to the Special Issue Application of Mathematical Methods in Financial Economics)

Abstract

:
This paper proposes and motivates a dynamical model of the Chinese stock market based on linear regression in a dual state-space connected to the original state-space of correlations between the volume-at-price buckets by a Fourier transform. We apply our model to the price migration of orders executed by Chinese brokerages in 2009–2010. We use our brokerage tapes to conduct a natural experiment assuming that tapes correspond to randomly assigned, informed, and uninformed traders. Our analysis demonstrates that customers’ orders were tightly correlated—in the highly nonlinear sense of prediction by the neural networks—with Chinese market sentiment, significantly correlated with the returns of the Chinese stock market, and exhibited no correlations with the yield of the bellwether bond of the Bank of China. We did not notice any spike of illiquidity transmitting from the US Flash Crash in May 2010 to trading in China.

Graphical Abstract

1. Introduction

The purpose of this paper is to extend methods of spectral analysis which facilitate visual representation of trading patterns including high-frequency data (see [1]. This paper does not use high-frequency data but analyzes the database of all executed orders by a select anonymous brokerage). The research intends to find new uses of Fourier analysis and deep learning in financial econometrics. Namely, we analyze predictable market dynamics in the state-space dual to the price-to-volume distribution and connected with the state-space by Fourier transform as an alternative of the conventional vector autoregression (VAR). Regression residuals were subjected to analysis using several neural networks. Another novelty is the use of deep learning not to predict the market data but to verify a protocol, which imitates a natural experiment.
The objective of our analysis is to validate the utility of this new method. In principle, it applies to datasets of any size, yet we use a limited dataset of the daily retail trades in the Chinese stock market during the years 2009–2010 so that calculations can be performed on the PC. These datasets were provided by the retail brokerages to the regulator as a matter of compliance.
The application of the state-space approach to the analysis of market microstructure is not new. Hendershott and Mankveld (2014) [2] specifically emphasized this line of research with application to HFT as an alternative to an autoregressive class of models. A standard way to analyze state-space distributions is Kalman filtering, which [3] notably implemented it to distinguish between liquidity-driven and informed trading components of the trading volume for S&P500. The preprocessing stage described below is our new alternative to Kalman filtering.
We employ a correlation measure of the state-space invented by [4], but here apply it to the price bucket in its entirety rather than to individual stocks. Taking correlations as the first stage of data de-noising has been done by [5], in particular, for his studies of the “Flash Crash” of 6 May 2010.
The paper is structured as follows. In Section 2, we run a literature review. In Section 3, we provide summary statistics of databases at our disposal. In Section 4, we describe the state-space of the problem. In Section 5, we compose the model of predictable trade and a description of its inputs and outputs. In Section 6, we provide validation for our model for the predictable variation of trading intensity. The residuals of our prediction model are being analyzed through shallow and deep learning networks simulating the decision process of the traders in response to the new events. We discuss information which can be gleaned by the fictitious traders in the subsequent Section 7. In Section 8, we investigate the prevalence of low-priced stock in the early Chinese stock market. In Section 9, we introduce a dynamic version of the Amihud illiquidity measure which we employ in Section 10 for a single event study. This single event study is a hypothesized influence of the “Flash Crash” in the US, and is analyzed based on the Amihud illiquidity measure.

2. Literature Review

The empirical market microstructure has to deal with many complexities: the latency of execution, and incomplete or deliberately manipulated data. One such complexity is that order flow “lives” in transaction time rather than in physical time [6,7]. Another is that the real trading costs can be hard to estimate and relatively easy to conceal [8]. We partially circumvent limitations of the first kind by using interday correlations of intraday price migrations, using volume buckets as regression panels. Correlations should be free from the absence of transaction time stamps in our database because of the T + 1 rule, unique to the Chinese stock markets described in [9,10].
Many semi-empirical measures have been used to describe market behavior on a microstructure level [11]. The most popular or theoretically well-researched are the ones called VPIN (Volume-Synchronized Probability of Informed Trading [12], volume imbalances [13,14,15], VWAP and its modifications (Volume at Weighted Average Price, [16,17,18,19]), and all the different versions of Amihud measure [20,21]. Note that VPIN or VWAP measures do not distinguish particular stocks, placing them in uniform price buckets. This is the methodology we accept for the current paper; however, a relatively little depth of the 2009-2010 stock market in China, I had to modify it for available data. Consequently, in the low-priced stock segment (below 10 renminbi, further CNY), a bucket can contain a portfolio of similar-priced stocks, while in the high-priced market segment, one or zero stocks are more the norm (see Section 7 for details).
Three types of volume per weighted price distributions were analyzed: the “Buy” volumes, the “Sell” volumes, and the imbalances volumes, which could be called BVWAP, SVWAP, and IVWAP, in deference to extant terminology. The simulations were done with “Buy”, “Sell”, and “Imbalance” separately, but we use IVWAP in the rest of the paper unless explicitly indicated because it is most transparent in terms of interpretation. Our definition of imbalance varies slightly from the standard [16], where it is defined as twice the difference between buy and sell volume divided by their arithmetic mean. We use geometric means for the same purpose (see Section 3).
In the case of imbalances, our measure is similar to the VPIN distribution proposed by Easley, Prado, and O’Hara [12], except that it does not involve the computation of intra-bucket price variance. Instead, we use day-to-day correlations of volumes within a given price bucket. We also employ the measure inspired by [20] to test the contagion between America and China during the days surrounding the Flash Crash. Evidence of contagion had previously been presented by [22].
To test our methodology, we used brokerage tapes of several Chinese brokerages submitted to the mainland Chinese stock exchange as a matter of regulatory compliance during 2009–2010. These tapes contain only completed trades; they do not have timestamps beyond one day, but display most of the trades with “Buy” or “Sell” indicators across the entire price range. Because the tapes divide trades by “Buy” and “Sell” (less than 10% of the records miss this stamp), we do not rely on the algorithmic estimation of this division as in [23]. Whatever incompleteness exists in our data, it lies in the reporting procedures for the brokerages which existed during these years.
To analyze the volumetric data, we combine them into uniform buckets of 0.5 CNY so that a typical number of buckets is around 150–180 during any given day during 2009–2010. We track the migration between the buckets as an indicator of the direction of trading. Then, we build a dual space—connected by the Fourier transform with the original state space—model of the Chinese stock market microstructure, which we further analyze by (deep) learning algorithms.
In the above analysis, we follow in the footsteps of [24] Foster and Viswanathan (1996), who developed a theoretical model of several groups of traders who try to predict the actions of others. Our model allows us to gauge how these predictions could have panned out empirically. We use three metrics of market reaction: the Chinese market sentiment [25,26,27], returns on the Shanghai stock market, and yields on a Bank of China 10-year bond. Using our model, we can directly and relatively parsimoniously explore the conditions prevailing in dark pools, artificially supplying or denying our assumed traders any external information about the activities of their colleagues or the direction of indexes.
The use of neural networks to analyze financial data now seems routine, although only a few substantive papers were published as late as five years ago [28]. Only in 2020 did top journals begin to publish research papers which used neural nets [29]. The difference between the present research and all the papers known to this author (see above-cited papers and [30]) is that the present research does not try to use deep learning to beat forecasts of the market data, which is their conventional purpose.
We certainly cannot match the sophistication of algorithms being used by the modern HFT firms and hedge funds and the computational power available to them [31,32,33], though our “primitive” algorithms could have been closer to the state-of-the-art in 2009–2010. Therefore, the protocol similar to the analysis of natural experiments (see Section 7 for more details) and a small, relatively outdated database were used to compensate for the lack of available resources.

3. Summary Statistics of the Databases

Brokerage tapes provided as Excel files have the following format shown in Figure 1.
The brokerage tapes include date, order price in CNY, type of order “buy” and “sell”, as well as the order’s volume. They did not have a timestamp in the years 2009–2010. Further on, we attribute separate spreadsheets as the “tapes” of fictitious traders from zero to four for all brokerages to conduct a natural experiment. As we can see from the data in Table 1, summary statistics for individual traders are comparable and we treat them as an extra level of randomization of our data. We plot the stock volumes and prices from one of the tapes in Figure 2. Visual inspection confirms an overall growth in trading throughout the two years, which is not surprising for a developing stock market, and includes a few spikes but no other obvious tendencies.
The tapes do not distinguish between individual stocks. Given the comparatively thin volume of trading in the stock markets of Mainland China during 2009–2010, we can only surmise that the trades within a given price bucket belong to one stock or a maximum of two stocks.
Our database included eight brokerages with the symbols rfokp4c, hvw5se4, zbe0rgv0, qwixupca, qguyi05q, gxbmxv0, q1ysmbyz, and 5vuyp3bu. Only the last three of these brokerage tapes contained complete data on the volume, which we further denote by the first symbol as “g”, “q”, and “5”. Most numerical examples in this paper, as well as the data, refer to Brokerage 5.
Even though it is likely that the division of trades between tapes is arbitrary, for later analysis we shall imagine them as belonging to separate “traders”. This corresponds to the intuitive idea that in modern high-frequency trading, the trader is a computer algorithm which arbitrarily parses the state-space.
Our tapes contain (see Table 1) several tens of thousands of trades over the period of two years. For comparison, a characteristic latency of a trading signal is τ ≈ 2–3 ms, which roughly corresponds to the computer messages cycling the circumference of New York City and its vicinity with the speed of light [6]. Inherent latency of trading quotes is even shorter, see [34], Table 1. Thus, if one wants to project this rate on the intensity of modern high-frequency trading, all the tapes of one brokerage would correspond to 7–8 min of wholesale trading. This illustrates the utility of observing emerging markets, where the tendencies requiring very large datasets to analyze can be observed with much less granularity.

4. Formation of the State-Space

We apply three stages of data analysis in our model. In the first stage, we allocate all daily orders to the price buckets. From these price buckets, we construct a state-space from the day-to-day correlations of order volumes, which we use as new vectors of our state-space.
X t ( x ) = ( V a r ( V t ) V a r ( V t + 1 ) ) 1 / 2 i C o v [ V i , t , V i + x , t + 1 ]
In Equation (1), V = { V 0 ,   V 2 ,   V n } is a vector of the trading volumes in the buckets i, numbered from 0 to n. Cov[.,.] is a covariance of the vectors V for the lag = x between the buckets. Most of the covariances for large x are negligible (see below).
The intuition for this definition is that only a reasonably small number of price buckets contribute to this measure because a daily price change for any stock is expected to be small compared to its price. Indeed, only the price changes below 8 CNY were consistently observed on each trading day, although these changes could obviously be larger for some days.
Equation (1) was written with an assumption that day-to-day correlations of order volumes exhibit more stability than the volumes themselves. Heuristically, this assumption is supported by the existence of the unique T + 1 rule in the Chinese stock markets, according to which one has to hold stock one day or more before selling [9] so that intraday noise must be uncorrelated between today and tomorrow.
The predictive power of the correlations (if any) can be used by an informed broker in the following manner. If correlations between buckets ρij were persistent, a broker could predict the migration by x in a given price bucket i, by the formula:
X t i ( x ) = ρ i , i + x V t i i j ρ i , j , t 1
In real life, the brokers can use a next-day predictor ρ ^ i j , to guess their next-day order. Note that in the 100%-efficient markets, or for markets in equilibrium, our state variable is exactly zero.
Covariance matrices could be more consistent from the mathematical point of view, but they are harder to interpret intuitively—particularly because they grow as squares of the volume with more active trading―and visualize. Moreover, covariance matrices, because of their nonlinear growth with average volume, obscure participation of the “high impact trades”, i.e., the trades, which influence price much in excess of their size (Xiaozhou, 2019). In Equation (1) we use three types of variables as “volumes”: a volume of buy trades, a volume of sell trades, and an imbalance volume, which we define as the difference between buy and sell volume at a given price bucket. In the case of imbalance statistics, our measure is reasonably similar to VPIN proposed by [12].
Our state-space is a discrete space of sixteen price buckets, separated by Δ = 0.5 CNY (daily changes of stock prices by more than 16Δ = 8 CNY were seldom observed in the sample). We split the entire trading book into 0.5 CNY buckets by the price change (a few hundred encompassing all the stock price range) but use only the first sixteen buckets. Using a larger number causes a spurious periodicity in our data. A smaller granularity will leave too few events in each bucket to allow confident averaging, while a larger granularity will average over most daily price changes. The division of the trading book into equal buckets allows us to avoid the problem that the trades in our database are not stamped with the name of a particular stock. The trading volume of all stocks experiencing “zero” or “significant” price change goes into the same bucket. Our construction of the phase space potentially allows a two-way analysis: panel analysis which is based on individual buckets and time-series analysis which follows the evolution of buckets through time.
The second stage of our analysis is building a dynamic model of trading. Our only assumption is that the state variables evolve by a linear dynamics, which we estimate from our data:
X i , t + 1 ( x ) = X i , t ( x ) + K ^ i j · X j , t ( x ) + ε j , t ( x ) δ i j
For our analysis, we use a dual state-space obtained by the Fourier transform of the initial state-space of a model. Philosophically, our choice of the state variable is based on the Bochner theorem in functional analysis, which states that covariance of a weakly stationary stochastic process always has a representation as a Fourier integral of a stationary measure [35]. Henceforth, a broad class of stochastic processes can be represented in the form above [36] (Chapters 14 and 15). Here, we only display our model in the form we had used in our analysis.
The Fourier transform of the Equation (4) gives a linear regression in dual state-space:
Δ X t + 1 , ω X t + 1 , ω X t , ω = β ^ ω ω X t , ω + e t + 1 , ω δ ω ω
In Equation (4), because of the Fourier transform, the vectors are assumed to be complex, i.e., with twice the dimensionality of the original state-space. The beta matrix has the dimension of 32 × 32 if we separate real and complex parts. Because our initial state vectors are real, there is covert symmetry in the coefficients and some rows in the beta matrix are identical zeroes. The Kronecker delta in the regression residual assumes that all spurious correlations between volumes disappear in one-two days.
Note that we make no assumptions about the random process governing the price dynamics. The only limitation of Equation (4) is the size of the beta matrix we use to approximate a continuous Fourier integral operator [37,38]. The Inverse Fourier transform of our beta operator is analogous to the Q-operator in [39] Markovian model of the Limit Order Book (LOB).
Original daily state vectors are recovered through the inverse Fourier transform below, where the hat denotes the predicted independent variable. They can have a small complex part because of the finite representation of decimals in the computer, which we ignore.
X t ( x ) , X t ^ ( x ) = F ω x 1 [ X t , ω , X ^ t , ω ]
And
ε t ( x ) = F ω x 1 [ e t , ω ]
In the third stage of our analysis, we employ neural networks to make sense of the regression residuals, i.e., whether they are reflecting real economic surprises or a result of noise trading. We do not know the prediction algorithms being used by the traders and, with time, they might become more complicated than anything we can devise. Therefore, we try an inverted strategy of deep learning. Namely, given an unpredictable part of the day-to-day volume correlations, we try to predict the realized indexes of the Chinese economy. The intuition behind this method is that if there is systematic unexpected buying or selling pressure in the market, it must reflect prevailing market sentiment.

5. First Validation of the Model

We have tested our model’s beta estimator for different traders in our database. Our results are represented by the sets of 484 × 16 matrices (the number of trading days during 2009–2010 times the number of the price buckets). The correlations between columns and rows of the matrix β ^ ω ω for the imbalance volumes in Formula (3) are given in Table 2. (In the series of our tests, we used Buy, Sell, Total Volume, and Imbalance indicators. The results were broadly similar across all selected measures (for instance, see Appendix B). For most of this paper, except Section 7 where we used “Buy” quotes, we selected imbalances to represent our data. In particular, imbalances can be directly compared with the “Cost of Trading” measure [40]. Complex beta matrices have dimensions 32 × 32 because of the real and imaginary parts of Fourier-transformed state vectors. Yet, under inverse Fourier transform, because of the internal symmetry, both the prediction and residual vectors are real.
We display temperature maps of an estimation of a single tape in Figure 3.
Correlations of beta matrices, Equation (4), computed between columns and rows for the temporal correlation of imbalance volumes (see Equation (1)). Tape numbering corresponds to the rows in Table 1. Correlations are symmetric across the diagonal. All the correlations between coefficients are insignificantly different from unity.
In testing regression (3) for the five data tapes, beta matrices are practically identical despite the state vectors being vastly different. This suggests the robustness of our model for the predictable component of the daily correlation of the imbalances. A similar picture was also observed from correlating betas between Buy and Sell tapes.

6. Predictable Component of the Bid-Ask Volume Correlations

As a criterion for the quality of approximation of the Equations (5) and (6), we use the vector error estimator for the predictor and the residuals:
F n = 1 T t = 0 T ε t 2 ( x n ) P n = 1 T t = 0 T X ^ t 2 ( x n ) σ 2 ( 1 F n σ 2 ) σ 2 = ( 1 T t = 0 T X t 2 ( x n ) )
where σ is the empirical volatility of the data and X ^ is the estimate from the regression of Equation (5).
In Equation (8), we retained the multiplier and 1/T, T = 484 (trading days) in the denominator as well as the numerator for clarity. Index n = 1 ÷ 16 numbers a vector of the state-space. A typical plot of variances of the predictor and the residual is given in Figure 4.
We note that the predictor and the residual time series by construction have zero correlation. However, the coincidence of time-weighted variances between the price buckets in Figure 4 is quite impressive and it is typical for a tally of the imbalances.
For quantifying the determination of regression prediction and regression residuals, we used running correlations of panel variances for the 484 trading days in the sample. The matrix of these correlations is provided in Table 3. From this matrix, we observe that 30–40% of the daily variability of the traders’ samples and 50% of the monthly variability is contributed by the prediction variance, and the rest by the regression residuals. The same observation that the variance of the empirical distributions is being split approximately 50:50 between the predictor and the residual can be made from Figure 4.
Table entries are the squares of r i j = C o r r e l a t i o n t [ V a r Δ ( Y ^ i ) ,   V a r Δ ( ε j ) ] where i = 1 ÷ 5 and j = 1 ÷ 5 are individual traders indicating the explanatory power of linear regression.
Correlations between different brokerage tapes are statistically insignificant. This exercise suggests that the trader’s samples are independent in the sense of linear regression. For the trader, it means that processing data from another trader by linear regression does not contribute any valuable information. Individual traders can fairly predict their correlations between today’s and the next day’s volumes, i.e., persistence of their own demand across all price buckets, but not correlations for other traders.

7. Analysis of the Phase Space Regression Residuals

The model of Equations (3) and (4) describes a predictable component in the day-to-day correlations of trading volume of price migrations including the zeroth price bucket (price changes below 0.5 CNY). Residuals contain both microstructure noise and reactions to unpredictable economic events in the market. To analyze the residuals, we employ several methods inspired by neural networks.
We do not know what kind of training algorithms traders might be using, given the quick progress in the algorithmic finance and computation power since 2009–2010 and even as this paper is being written. Henceforth, we employ the following method. Instead of a prediction of the out-of-sample trading data, we attempt backdating market data through a simulated experiment, namely, a neural network trained by our order data attempts to predict the Chinese market sentiment index, returns on the Shanghai stock index, and yields on the bellwether 10-year bond of the Bank of China (for details of the protocol, see below). Because of the monthly periodicity of the sentiment index, for the consistency of our tests, we used monthly stock returns and monthly bond yield as well.
In our case, real-life trading algorithms would have to predict the “unexpected” direction of price changes imprinted in brokerage orders given their information on the markets. Yet, we assign to our imagined traders—represented by the brokerage tapes—a much simpler task of predicting a monthly index given their observation of the day-on-day correlation of orders within a given range of price change. (This reasoning is based on an unproven but intuitive assumption that an economically simpler problem—guessing a “covert” index from proprietary trading data rather than the other way around— is less demanding algorithmically).
Our procedure corresponds to the following stylized situation. We select a randomly chosen “informed” trader who observes orders from her own clients and trains her network by predicting the index. We use her data as a network input and then simulate the behavior of other traders whom we consider uninformed as to the direction of the three chosen indexes but who controllably can observe or be in the dark concerning the actions of their colleagues from the same brokerage (Figure 5).
The situation of “leaky brokers” has been described in [41] in the following terms: “When considering the theoretical soundness of a market equilibrium in which brokers leak order flow information, one may wonder why an informed asset manager is willing to trade with brokers that tend to leak to other market participants… The broker would enforce this cooperative equilibrium across subsequent rounds of trading. In particular, the broker can exclude from the club the managers that never share their private information and reward with more tips the managers that are more willing to share”.
The results were averaged over six or twelve independent runs of the network and were not significantly different.
The best results were obtained by a seven-layer convolutional neural network (CNN), though other options have been explored (Appendix C). Conventionally, CNN is used for image recognition and analysis. Essentially, we used matrices of the residuals of the output regression (depicted as a heat map of the matrix in Figure 3) as if they were digitized information for the visual images to predict the direction of an index. Table 4 displays the trials with randomly selected informed traders in both training and predictive samples, as well as training samples with only “uninformed” traders, i.e., the traders who observe only bid and ask volumes per basket, without access to current or past magnitudes of the index. Unlike the results from Table A1 in Appendix C, the statistically significant results from Table 4 were broadly reproducible on successive runs of the network.
The general conclusion from Table 4 is that CNN can reliably predict Chinese market sentiment from daily imbalances, prediction of the stock market returns is usually significant at 10% but not at 5%, and the yield of BOC bond cannot be inferred from the imbalances.
The implication of insider information does not improve prediction very much for the market sentiment index, somewhat helps to predict the direction of the stock index but well within the assumed 10% statistical dispersion of the results, and is irrelevant for the direction of bond yields. In all cases, there is little difference whether an “uninformed” trader trains her network on the imbalances of her informed colleague or another uninformed trader.

8. Discussion: Possible Sampling Issues

Reporting files do not contain identification of the execution price and a particular stock. Our formation of the price buckets was based on a uniform division of the price range in a trading book into the intervals of 0.5 CNY. We considered this choice optimal because it allows a significant number of price buckets (120 ÷ 400). Furthermore, migration of the price for more than 8 CNY in a given day is rare, and we can use a parsimonious 16×16 matrix approximation for the Fourier evolution operator.
This, or any similar choice―for instance, an arbitrary 0.4317 CNY―assures that there could be one, several, or no trades in a particular price bucket. Yet, the Chinese stock market, which, in 2009–2010 was in its nascence, was dominated by the low-priced stocks (below 10 CNY). At the end of 2009, of the 293 listed equities, only 34 issues (11.6%) had a mean price below 10 CNY, and 259 had a price above that number. Henceforth, the lower price buckets could be systematically different from the upper buckets in that they might contain several stocks, while the price buckets above the average (see Table 1) can contain only one stock or none altogether.
To clarify this problem, we artificially split the trader books into parts comprising the stocks with the price below 10 CNY and above 10 CNY. Of course, there could be some borderline migration of stocks priced slightly above 10 CNY into the first portfolio and stocks priced slightly below into the second. We expect a small influence of this issue on day-to-day trading and we ignored it in our analysis. When we accomplished the procedure described in Section 4, Section 5, Section 6 and Section 7, for our censored trading books, the results were broadly the same as having been observed in Table 4. Namely, if one uses portfolios of low-priced stocks to train the net and then predict the index from the trading books, in which only the high-priced stocks are included, one can confidently infer the sentiment index; the Shanghai stock index is predicted in some runs but with low statistical validity, and there are no correlations between the net trained our traders’ positions and the yields on the 10-year bond of the bank of China. These results do not change much if we train the net on the high-priced stock and leave the prediction to the trader of the low-priced stock. Only the prediction of the sentiment index slightly improves (correlation grows from ~95–97% to ~97–99%), which suggests that the higher-priced stock was more liquid and, henceforth, had a better predictive value. As is the case with all statistical experiments, this observation is only tentative.

9. Empirical Liquidity of the Chinese Stock Market in the Period 2009–2010

A proposed microstructure model of the Chinese stock market allows us to analyze both predictable and unpredictable frictions resulting from two interleaving factors: (1) imperfect balance between buy and sell orders, and (2) securities changing value during trading.
The net cost of trading is computed similarly to [40], though their formula can accept different conventions. Our formula presumes that the brokerage sells an asset in today’s quantity marked to market at yesterday’s buy price and replenishes its inventory sold yesterday at today’s ask price, with the cost to the customer being the same in value and opposite in magnitude. Of course, the signs in Equation (9) are arbitrary.
π t i = p a i , t 1 V b , t p b i , t 1 V s , t
In Equation (9), πt is our definition of the cost of trading, and pa, pb are the ask/bid price buckets. The Vb, Vs are the volumes of buy/sell orders. The index i = 1–16 signifies the price bucket. Note that in market equilibrium, in Equation (7), cost averaged over all buckets is equal to the (constant) bid-ask spread times the daily turnover and is always non-negative. Outside of equilibrium, the sign of π t can be arbitrary because of fluctuating stock prices.
Our analysis by CNN indicates that the net cost of trading is a fair predictor of the market direction in the sense that we have outlined in a previous section (Figure 5). Namely, if we assume that the broker or regulator is “blind” to the order size, she can get a clear idea about the Chinese market sentiment from the trading costs only. Their idea of stock market direction would be imperfect but statistically significant and, finally, there is no connection to the Bank of China bond prices through our model.
While this exercise is purely imaginary as being applied to the Chinese brokerages, we suggest that this conclusion—that, in the observed period, trading costs reflected market sentiment more or less mechanically (see Figure 6)―can help traders and regulators alike in the case of “Dark Pools”. In the latter case, the information about the exchange’s strategy is covert and can be gleaned only indirectly.
The Equation (9) can be recast in the (dynamic) version of [20] liquidity measure. Namely,
λ t i = | π t i | ( V b , t + V a , t 1 ) 2
Further on, we use liquidity lambda to predict the same indexes. The intuitive meaning of our version of the Amihud measure is that it represents the average cost for the agent to make a roundtrip inside the same price bucket with one share. To provide a glimpse of the magnitude and volatility of λ, we display its daily dynamics in Figure 7.

10. Single Event Analysis

Our sample includes a day of the Flash Crash in US stock markets (6 May 2010), which means, dependent on daytime, either trading day 326 or 327 in our sample. There is no visible anomaly in the liquidity of the Chinese stock market during or after that day. It is interesting to analyze this event using our methods.
The spillovers from the established stock markets into the Chinese markets had previousy been studied by [42] using the measure of volatility proposed by [43,44]. They observed that before 2010, the Chinese stock market produced volatility spillovers to Taiwan and Hong Kong, but its influence on the European and North American markets was statistically insignificant. However, the shocks in the US stock market affected every stock market they studied. Beginning in approximately 2007, some pushback from the Chinese stock market could be observed.
The question of whether the mutual influences between the stock markets were macroeconomic or microstructural in nature was investigated by [22]. They also observed asymmetric shocks, i.e., the shocks propagating predominantly from US markets into China but not the other way around. Li and Peng noticed that a structural shock in the US markets usually decreases correlations between the Chinese and American stock markets.
We decided to investigate the influence of the US stock markets on China by our methods of the CNN analysis of the regression residuals. We display the testing strategy in Figure 8. The observation period (years 2009–2010) is split into eight overlapping samples of 60 days each. One sample of two adjacent periods (usually, but not necessarily the first) is used for network training. This constitutes one training and five predictor samples. We then attempt to predict monthly indexes backward from the training data. The null hypothesis is formulated as follows:
Hypothesis 1 (H1).
For i = 1 ÷ 5, Correlation [Index prediction from λi, Index] = Correlation [Average[λ], Index].
Hypothesis 2 (H2).
For at least one i, Correlation [Index prediction from λi, Index] ≠ Correlation [Average[λ], Index].
The intuitive meaning of the null hypothesis above is that predictions obtained from subsamples of illiquidity measures are no different from each other. We, of course, would prefer that the null be rejected for the samples containing the Flash Crash in America.
The statistical significance of the correlations of the predicted monthly indexes for one of the traders (tapes) is shown in Table 5. Note that, in agreement with the results of Section 6, almost no statistically significant correlation of prediction of the return on the stock index and none for the yields on the 10-year bond can be observed. On the contrary, the prediction of the market sentiment from an observed Amihud illiquidity measure is robust. This corresponds to the testing of H1, and of H2 above where the right-hand side is replaced by zero.
As a rule, null cannot be rejected for any of the three indexes. (A) Probabilities for prediction of sentiment index, (B) probabilities for prediction of stock market returns, and (C) probabilities for prediction of bond yields. Highlighted are the subsamples, which include the US Flash Crash of 6 May 2010.
We observe from Table 5 that the null hypothesis―illiquidity in the subsamples that include the date of the American Flash Crash behaves any different from other samples—cannot be rejected. Only one probability in Table 5 is below 10%, and it is not stable with respect to the consecutive runs of the neural network with a different seed.
We cannot discern an influence of a shock from microstructure data, yet the results, e.g., of Li and Peng (2017), suggest that the transmission of shocks was real. This suggests that the transmission between the US and Chinese markets was mainly through macroeconomic fundamentals. As we have seen from Table 5, our deep learning-based method could simply require too long a sample to resolve the spillover.

11. Conclusions

In our paper, we propose a microstructure model of the Chinese stock market. We build it from a state-space of interday correlations of volumes between price buckets. Because of the T + 1 rule in the Chinese stock markets, interday correlations are expected to significantly reduce the microstructure noise. For the 100%-efficient market in equilibrium, our time series would be exactly equal to zero.
This model is analyzed by OLS in a dual state-space, connected to an original state-space by a Fourier transform across the multiple price buckets. This procedure corresponds to an approximation of the predictable fraction of the time evolution of trading volumes by a Fourier integral (pseudodifferential) operator of a general form. In a discrete case, this operator can be represented by a matrix of arbitrary dimension (selected by convenience as 32 × 32) acting on the space of the Fourier coefficients (see Appendix A for details). This method is completely general and can be applied to any time series, which can be grouped into panels, in our case volumes per a given stock price.
The presented model based on the Fourier integral operators, augmented by the deep learning methods, can be calibrated daily by the market makers. It predicts the next day’s volume at a price change x (see Equation (2)) and has no fundamental price inputs. Henceforth, the prediction reveals the dynamics of brokerage demand driven by its clientele without intervening events in the market.
The proposed method (or protocol, in the parlance of the deep learners: (1) preprocessing, filtering using Fourier Integral Operators; (2) processing, recognition of the sequence of residual matrices using neural network; (3) postprocessing, index forecast with a trained network) is universal and can be applied to any trading system as a continually updated model. It can simulate any given market with a level of complexity chosen by the broker and/or the regulator. Namely, each day, a predictable trading forecast is provided by a multivariate filter (Fourier Integral Operator in the current paper, but other methods, e.g., Kalman filter, can be used as well). This stage is “dumb”, in the sense that it does not take into account intervening changes in market fundamentals. Past residuals to the regression can be independently predicted by the neural network and added up to the results of the filtering. In the third stage, different validation strategies can be applied. This protocol can incorporate generic processes such as multivariate ARMA (similar to the one used in the paper) or a specific microstructure model, e.g., Roll’s [4].
The main conclusion of our analysis is that the unpredictable dynamics of trades completed by the Chinese brokers in the period 2009–2010 were tightly correlated with the Chinese market sentiment in a highly nonlinear sense of machine learning. This can be interpreted as investors trading according to the available market information they receive. The returns on the stock index were predictable, but statistically they were barely significant. Because the predictive power of the trades for the actual stock market returns was small, this indicates investors were following a herd mentality. We observed no connection between the trading activity and the yields on the bellwether 10-year bond of the Bank of China. This may signify that the market risk in an emerging market, which describes the Chinese stock market in 2009–2010, has only a small dependence on prevailing borrowing rates.
Finally, we tested whether the Flash Crash of American markets on 6 May 2010, was reflected in the liquidity of the Chinese stock market. For this, we used a dynamic version of the Amihud illiquidity measure λ. The illiquidity measure predicted by the CNN was almost as good a predictor of market sentiment as the VAP (see Figure 6), and we used it to analyze a possible contagion of the American Flash Crash on the Chinese markets.
With our methods, we did not find any statistically significant evidence that liquidity was higher or lower than the average during the period preceding or following the Flash Crash. This indicates a need for a finer measure of contagion between American and Chinese stock markets. In particular, samples of 60 days long can be too crude to resolve the influence of the Flash Crash, which may have lasted only 1–2 trading days.
The limitations of the proposed method are that all regressions at the preprocessing stage must be linear. Non-linear regressions lead to the non-linear integral operators on the right-hand side, which may be more complicated or less computationally stable than the original equations.
The second limitation is that the choice of neural network (discussed in Appendix C) can be rather arbitrary with few hard-wired criteria to guide the selection. Moreover, the main limitation of this study is its restricted array of free data from the Shenzhen Stock Exchange and its limited computer resources. (All computer algorithms were tailored to be run on a conventional PC using not more than 15–30 min of processor time in the case of R algorithms.) However, all of the paper’s methods can be applied to datasets of arbitrary size.

Funding

The paper was completed without an external funding.

Institutional Review Board Statement

Not relevant for this publication.

Informed Consent Statement

Not relevant for this publication.

Data Availability Statement

All datasets and programs used in the current study are available on personal request from this researcher.

Acknowledgments

I thank my former student, Lin Jilei from Wenzhou Kean University (U. Illinois) for his invaluable aid with coding, Oliviero Roggi, chairman of my section, participants of the Risk Society Conference (IRMC, Florence), and an anonymous referee for the helpful suggestion for the Section 8. I also thank Pratiksha Sharma, my discussant at GFC (2021, online), for the literary and computational suggestions for the Section 10 and the organizers of this conference, especially Manuchekhr Sharokhi, who conducted it during such a difficult time as well as the organizers of the Portuguese Finance Society PFN 2021 conference and Vietnam’s Conference on Banking and Finance 2021.

Conflicts of Interest

The author does not recognize any conflict of interest.

Appendix A

Mathematically, one can make the following observation about the temporary correlations as the state-space variables. Let u ˜ t = u t + ϵ t 1 and v ˜ t = v t + ϵ t 2 be our time series, where ut and vt are signals with the correlation coefficient ρ and epsilon terms are microstructure noises with variances σ ϵ 1 2 and σ ϵ 2 2 uncorrelated with the signals and each other: E [ u t , ϵ t 1 ] = E [ u t , ϵ t 2 ] = E [ v t , ϵ t 1 ] = E [ v t , ϵ t 2 ] = E [ ϵ t 1 , ϵ t 2 ] = 0 . Then, their noisy correlation for a small white noise becomes:
ρ ^ = C o r r [ u ˜ t , v ˜ t ] C o v [ u ˜ t , v ˜ t ] V a r [ u ˜ t ]   V a r [ v ˜ t ] = C o v [ u t , v t ] V a r [ u t ] + E [ ( ϵ t 1 ) 2 ]   V a r [ v t ] + E [ ( ϵ t 2 ) 2 ] C o v [ u t , v t ] V a r [ u t ]   V a r [ v t ] ( 1 1 2 ( σ ϵ 1 2 V a r [ u t ] + σ ϵ 2 2 V a r [ v t ] ) ) ρ
This is a downward biased estimator of a true correlation, stationary, as long as the second moments of the noise are small and stationary. Furthermore, the correlation coefficients are concentrated in the [−1, 1] range and, thus, are more amenable to intuitive and graphic interpretation.

Appendix B. Pseudodifferential Operators

Pseudodifferential operators are a particular case of the Integral Fourier Operators [43]. Frequently, financial time series are being described by the AR(n) models, which are, in turn, the discrete analogs of the differential operators with constant coefficients. Pseudodifferential operators can be considered as extensions of ARMA(p,q) models. From that angle, conventional ARMA(p,q) models are discrete pseudodifferential operators with a rational function as a symbol. The formal definition of the pseudodifferential operator goes as follows (e.g., [37,38,45]:
A ( u ) = 1 ( 2 π ) n 2 n e i ( x y ) ξ a ( x , y , ξ ) u ( y ) d y d ξ
One can easily write a solution for a Cauchy problem for the Kolmogorov–Fokker equation, describing Ito diffusion through a pseudodifferential operator. Indeed,
ψ t = ( ψ )
ψ ( x , 0 ) = f ( x )
where is a generator of Ito diffusion of the following form. Here, a is an n-dimensional vector, Σ is the n × n matrix, and i, j = 1 ÷ n—dimensions of the state-space.
= a · i + i Σ ^ j
This problem’s solution expressed in a form of Equation (A1) looks as:
ψ ( x , t ) = 1 ( 2 π ) n 2 n e i k ( x y ) A ( k , t ) · f ( y ) d n y d n k
With
A ( k , t ) = e ( i k · a k T Σ ^ k ) t
In our case, the phase-space regression has a form of Equation (4) of the main text:
X t + 1 , ω X t , ω = β ^ ω ω X t , ω + e t + 1 , ω δ ω ω
To better demonstrate a connection of this regression to the pseudodifferential operators, we replace discrete time steps in Equation (A6) with continuous time. The transfer Equation (A7) becomes:
d X ( t ) ω = β ^ ω ω X ( t ) ω + d ε ( t ) ω δ ω ω
The value of the state vector as a function of a state variable x for an arbitrary time T can be expressed as a moving average-type equation:
X ( x , T ) = ( e β ω ω ( T t ) + i ( ω ω ) t d t ) X ( x , 0 ) + ω , ω e β ω ω ( T t ) + i ( ω ω ) t ε ( t ) ω d t
In Equation (A4), the symbol of the pseudodifferential operator is equal to:
A ( ω , t , T ) = ω , ω e β ω ω ( T t )

Appendix C. Types of Neural Networks Used in Our Estimations

The first test used the prediction of three indexes from the first four average moments using a shallow neural network with one hidden layer. When we feed the network data from the stock imbalances, the shallow learning network predicts a constant answer indicating zero information about the direction of the indexes. This negative result is not all bad because it suggests that the distribution of the residuals is practically indistinguishable from a normal distribution, a nice enough feature for such an artificial model.
Figure A1. Example of estimating of shallow three-layer network with one hidden layer for the prediction of the market sentiment (stock returns, bond yields) from the first four moments of the error distribution. The input data for the network are pre-processed and the data for the first four monthly moments are being fed into the network. Blue arrows from circles marked as ‘1’ signify the training of one sample by the actual market index data.
Figure A1. Example of estimating of shallow three-layer network with one hidden layer for the prediction of the market sentiment (stock returns, bond yields) from the first four moments of the error distribution. The input data for the network are pre-processed and the data for the first four monthly moments are being fed into the network. Blue arrows from circles marked as ‘1’ signify the training of one sample by the actual market index data.
Mathematics 10 02459 g0a1
For our second test, we have used a complete set of residual matrices and the following experimental procedure (Figure 4 in the main text). We trained a relatively deep 10-layer network on a sample of all monthly regression errors from a randomly chosen trader omitting the last day of each month and tested this prediction on the same trader’s data sample of the last days of the month. Then, we tried to predict our indexes (sentiment, stock returns, and bond yields) back from the supposedly unexpected changes—provided by our regression―in other trading tapes.
One of the difficulties in dealing with neural networks is that results frequently represent multidimensional tensors. By their origin, they cannot be listed compactly in two-dimensional tables and their presentation on a sheet of paper or a computer screen is non-intuitive.
We present the results of the 10-layer network in Table A1. Training on buy or sell signals had some limited power for predicting market indexes from a sample of imbalance correlations for the same trader, or imbalances for other traders whom we considered “blind”, i.e., who predict the direction of indexes from their observations of the imbalance quotes. Some outputs from this model are shown in Figure A2.
Figure A2. An example of a 10-layer network. The input of the network indicates a 16-dimensional vector—the size of state-space of the price changes and the scalar output—training on, or prediction of, one of the three monthly indexes.
Figure A2. An example of a 10-layer network. The input of the network indicates a 16-dimensional vector—the size of state-space of the price changes and the scalar output—training on, or prediction of, one of the three monthly indexes.
Mathematics 10 02459 g0a2
The inverse path using the same network, i.e., predicting indexes by the training sets on unexpected changes in imbalances, was added for control. The scalar network results show little dependence on the number of training rounds and functional shapes of individual transmission functions (ReLu—rectified linear unit, hyperbolic tangent, or logit).
Figure A3. Examples of the Chinese sentiment index (blue) and its backward prediction (red) by a 10-layer scalar network from buy quotes (correlation ρ = 0.5841) and sell quotes (correlation ρ = 0.2934) of randomly selected traders. Traders are “uninformed”, i.e., they make predictions based on their observations through a network trained by an informed trader. Results are poorly reproducible.
Figure A3. Examples of the Chinese sentiment index (blue) and its backward prediction (red) by a 10-layer scalar network from buy quotes (correlation ρ = 0.5841) and sell quotes (correlation ρ = 0.2934) of randomly selected traders. Traders are “uninformed”, i.e., they make predictions based on their observations through a network trained by an informed trader. Results are poorly reproducible.
Mathematics 10 02459 g0a3
Table A1. Select runs of a 10-layer neural network for the backward prediction of monthly indexes from traders’ activity. Capital letters B, S, and I mean “Buy”, “Sell”, and “Imbalance” samples of residuals. Letters from a-d refer to a particular trader. An arrow designates training vs. prediction tapes. The symbols r1–3 refer to the correlation coefficients of the neural network predictor and the actual indexes. The explanatory power of the predictions can be inferred from the squares of the correlation coefficient. For instance, if we treat broker “b” as informed, trader “a” could have predicted bond yield for the next month from their imbalances with an explanatory power ρ = r 3 2 = 0.5178 2 26.8 % . Test results were poorly reproducible on successive runs of the network.
Table A1. Select runs of a 10-layer neural network for the backward prediction of monthly indexes from traders’ activity. Capital letters B, S, and I mean “Buy”, “Sell”, and “Imbalance” samples of residuals. Letters from a-d refer to a particular trader. An arrow designates training vs. prediction tapes. The symbols r1–3 refer to the correlation coefficients of the neural network predictor and the actual indexes. The explanatory power of the predictions can be inferred from the squares of the correlation coefficient. For instance, if we treat broker “b” as informed, trader “a” could have predicted bond yield for the next month from their imbalances with an explanatory power ρ = r 3 2 = 0.5178 2 26.8 % . Test results were poorly reproducible on successive runs of the network.
Trader/B, S, ISentiment, r1Market Return, r2BOC Bond, r3 Notes
Bb → Ia 0.03890.4099 (−0.2888)0.5178(Tanh) instead of ReLU
Bb → Ib−0.16630.2607−0.1339
Bb → Ic0.3075−0.3951−0.2713
Bb → Id0.46000.39360.2525
Sb → Bc−0.0431−0.3944−0.2389Control
Ic → Bb−0.2255−0.32820.1444Control
Our third exercise was to use a seven-layer convolutional neural network (CNN) sketched in Figure A1. CNN is conventionally used for image recognition and analysis. Essentially, we used matrices of output regression (Example provided in Figure 2) as if they were digitized information for the visual images to predict the direction of an index. Unlike the results from Table 3, the statistically significant results from Table 4 were broadly reproducible on successive runs of the network.
Figure A4. Sketch of the 7-layer CNN network.
Figure A4. Sketch of the 7-layer CNN network.
Mathematics 10 02459 g0a4
Finally, we tested the Long Short Term Memory (LTSM) network and post-processed the results by several smoothing algorithms (mean, moving average, exponential moving average) for the previous 21 days—the average trading month in China for the years in question—to compare the results with monthly indexes. The box chart of the procedure is portrayed in Figure A5. We plot some of the outputs in Figure A6. The results of this procedure do not depend much on the number of epochs and batch sizes are slightly better than the 10-layer network above but worse than the CNN network, which we use in the main text.
Figure A5. Box chart of the Long Short Term Memory estimation in our paper’s context. Index input signifies one of three monthly indexes (sentiment, stock market returns, and yield on BOC bond). Daily predictions of the indexes are smoothed through post-processing before being compared to the original inputs. The LTSM network was run with default Keras parameters.
Figure A5. Box chart of the Long Short Term Memory estimation in our paper’s context. Index input signifies one of three monthly indexes (sentiment, stock market returns, and yield on BOC bond). Daily predictions of the indexes are smoothed through post-processing before being compared to the original inputs. The LTSM network was run with default Keras parameters.
Mathematics 10 02459 g0a5
Figure A6. Backward LTSM prediction of monthly sentiment index (blue) using data from the broker’s own training tape (orange), and three other brokers (blue, grey, and yellow dash). The general downward tendency of the sentiment for the period but not much else can be predicted from other brokers’ tapes.
Figure A6. Backward LTSM prediction of monthly sentiment index (blue) using data from the broker’s own training tape (orange), and three other brokers (blue, grey, and yellow dash). The general downward tendency of the sentiment for the period but not much else can be predicted from other brokers’ tapes.
Mathematics 10 02459 g0a6

References

  1. Angel, J.J. When Finance Meets Physics: The Impact of The Speed of Light on Financial Markets and Their Regulation. Financ. Rev. 2014, 49, 271–281. [Google Scholar] [CrossRef] [Green Version]
  2. Henderschott, T.; Menkveld, A.J. Price pressures. J. Financ. Econ. 2014, 114, 405–423. [Google Scholar] [CrossRef]
  3. Rzaev, K.; Ibikunle, G. A state-space modeling of the information content of trading volume. J. Financ. Mark. 2019, 46, 100507. [Google Scholar] [CrossRef]
  4. Roll, R. A simple measure of the effective bid-ask spread. J. Financ. 1984, 39, 1127–1139. [Google Scholar] [CrossRef]
  5. Lehalle, C.-A.; Laruelle, S. Market Microstructure in Practice; WSPC: Singapore, 2013. [Google Scholar]
  6. Hasbrouck, J. YouTube FMA Lecture. 2016. Available online: https://www.youtube.com/watch?v=EZCgW1mFRP8&t=1905s+ (accessed on 29 December 2016).
  7. Zhou, X.; Georges, D. Information Environments and High Price Impact Trades: Implications for Volatility and Price Efficiency. In Proceedings of the Paris December 2019 Finance Meeting EUROFIDAI-ESSEC, Paris, France, 19 December 2019. Working Paper HEC 19-3. [Google Scholar]
  8. Keim, D.B.; Madhavan, A. The cost of institutional equity trades: An overview. Financ. Anal. J. 1998, 54, 50–69. [Google Scholar] [CrossRef]
  9. Guo, M.; Li, Z.; Tu, Z. A unique “T + 1 trading rule” in China: Theory and evidence. J. Bank. Financ. 2012, 36, 575–583. [Google Scholar] [CrossRef]
  10. Qiao, K.; Dam, L. The overnight return puzzle and the “T + 1” trading rule in Chinese stock markets. J. Financ. Mark. 2020, 50, 100534. [Google Scholar] [CrossRef]
  11. Fong, K.Y.L.; Holden, C.W.; Trzcinka, C.A. What are the best liquidity proxies for global research? Rev. Financ. 2017, 21, 1355–1401. [Google Scholar] [CrossRef]
  12. Easley, D.; De Prado, M.M.L.; O’Hara, M. Flow toxicity and liquidity in a high-frequency world. Rev. Financ. Stud. 2012, 25, 1457–1493. [Google Scholar] [CrossRef] [Green Version]
  13. Humphery-Jenner, M.L. Optimal VWAP trading under noisy conditions. J. Bank. Financ. 2011, 35, 2319–2329. [Google Scholar] [CrossRef]
  14. Lipton, A.; Pesavento, V.; Sotiropoulos, M.G. Trade arrival dynamics and quote imbalance in a limit order book. arXiv 2013, arXiv:1312.0514. [Google Scholar]
  15. Cartea, A.; Jaimungal, S. A Closed Form Strategy to Target VWAP; IPAM, UCLA: Los Angeles, CA, USA, 2015. [Google Scholar]
  16. Cont, R.; Kukanov, A.; Stoikov, S. The price impact of order book events. J. Financ. Econ. 2013, 12, 47–88. [Google Scholar] [CrossRef] [Green Version]
  17. Cartea, A.; Donnelly, R.F.; Jaimungal, S. Enhancing trading strategies with order book signals. Appl. Math. Financ. 2018, 25, 1–35. [Google Scholar] [CrossRef]
  18. Cartea, A.; Jaimungal, S.; Ricci, J. Buy low, sell high: A high frequency trading perspective. SIAM J. Financ. Math. 2014, 5, 415–444. [Google Scholar] [CrossRef]
  19. Lerner, P. Patience vs. impatience of traders: Formation of the value-at-price distribution through competition for liquidity. Int. J. Financ. Eng. 2015, 2, 1550029. [Google Scholar] [CrossRef]
  20. Amihud, Y.; Mendelson, H. Asset pricing and the bid-ask spread. J. Financ. Econ. 1986, 17, 223–249. [Google Scholar] [CrossRef]
  21. Hasbrouck, J. Empirical Market Microstructure; Oxford University Press: Oxford, UK, 2007. [Google Scholar]
  22. Li, X.-M.; Peng, L. US economic policy uncertainty and co-movements between Chinese and US stock markets. Econ. Model. 2017, 61, 27–39. [Google Scholar] [CrossRef]
  23. Lee, C.; Ready, M. Inferring trade direction from intraday data. J. Financ. 1991, 46, 733–746. [Google Scholar] [CrossRef]
  24. Foster, F.; Viswanathan, S. Strategic trading when agents forecast the forecasts of others. J. Financ. 1996, 51, 1437–1478. [Google Scholar] [CrossRef]
  25. Baker, M.; Wurgler, J. Investor sentiment and the cross-section of stock returns. J. Financ. 2006, 61, 1645–1680. [Google Scholar] [CrossRef] [Green Version]
  26. Hu, C.; Chi, C. Investor sentiment: Rational or irrational. Chin. Rev. Financ. Stud. 2012, 6, 46–62. [Google Scholar]
  27. Liu, C.; An, Y. Investor sentiment and the basis of CSI 300 stock index futures: An empirical study based on QVAR model and quantile regression. Discret. Dyn. Nat. Soc. 2018, 2018, 4783214. [Google Scholar] [CrossRef]
  28. Dixon, M.; Klabjan, D.; Bang, J.H. Classification-based financial markets prediction using deep neural networks. Algorithm. Financ. 2017, 6, 67–77. [Google Scholar] [CrossRef] [Green Version]
  29. Gu, S.; Kelly, B.; Xiu, D. Empirical asset pricing via machine learning. Rev. Financ. Stud. 2020, 33, 2223–2273. [Google Scholar] [CrossRef] [Green Version]
  30. Gu, S.; Kelly, B.; Xiu, D. Autoencoder asset pricing models. J. Econ. 2021, 222, 429–450. [Google Scholar] [CrossRef]
  31. Yang, S.Y.; Qiao, Q.; Beling, P.A.; Scherer, W.T.; Kirilenko, A.A. Gaussian process based trading strategy identification. Quant. Financ. 2015, 10, 1683–1703. [Google Scholar] [CrossRef] [Green Version]
  32. Davey, K. Building Algorithmic Trading Systems: A Trader’s Journey from Data Mining, to Monte Carlo Simulations, to Live Trading; Wiley: New York, NY, USA, 2014. [Google Scholar]
  33. de Prado, L. Advances in Finance Machine Learning; Wiley: New York, NY, USA, 2018. [Google Scholar]
  34. Bartlett, R.P., III; McCrary, J. How rigged are stock markets? Evidence from microsecond timestamps. J. Financ. Mark. 2019, 45, 37–60. [Google Scholar] [CrossRef] [Green Version]
  35. Stein, M.L. Interpolation of Spatial Data; Springer: New York, NY, USA, 1999. [Google Scholar]
  36. Liptser, R.S.; Shiryaev, A.N. Statistics of Random Processes II: Applications; Springer: New York, NY, USA, 1978. [Google Scholar]
  37. Taylor, M.E. Pseudodifferential Operators; Princeton University Press: Princeton, NJ, USA, 1981. [Google Scholar]
  38. Hörmander, L. The Analysis of Linear Partial Differential Operators; Springer: Heidelberg, Germany, 1987. [Google Scholar]
  39. Huang, W.; Lehalle, C.-A.; Rosenbaum, M. Stimulating and analyzing order book data: The queue reactive model. J. Am. Stat. Assoc. 2015, 110, 107–122. [Google Scholar] [CrossRef] [Green Version]
  40. Menkveld, A.J.; Zoican, M. Need for speed? Exchange latency and liquidity. Rev. Financ. Stud. 2017, 30, 1188–1228. [Google Scholar] [CrossRef]
  41. Barbon, A.; Di Maggio, M.; Franzoni, F.; Landier, A. Brokers and flow leakage: Evidence from fire sales. J. Financ. 2019, 74, 2707–2749. [Google Scholar] [CrossRef]
  42. Zhou, X.; Zhang, W.; Zhang, J. Volatility spillovers between the Chinese and world equity markets. Pac. Basin Financ. J. 2012, 20, 247–270. [Google Scholar] [CrossRef]
  43. Diebold, F.X.; Yilmaz, K. Measuring financial asset return and volatility spillovers, with applications to global equity markets. Econ. J. 2009, 199, 158–171. [Google Scholar] [CrossRef] [Green Version]
  44. Yilmaz, K. Return and volatility spillovers among the East Asian equity markets. J. Asian Econ. 2010, 21, 304–313. [Google Scholar] [CrossRef] [Green Version]
  45. Treves, J.-F. Introduction to Pseudodifferential and Fourier Integral Operators; Springer: New York, NY, USA, 1980; Volume 1. [Google Scholar]
Figure 1. Printout from one of the brokers’ tapes.
Figure 1. Printout from one of the brokers’ tapes.
Mathematics 10 02459 g001
Figure 2. Daily turnover and volume from one of the brokerage tapes. Turnover (in CNY) is shown by the blue curve (left axis) and trading volume is shown by an orange curve (right axis).
Figure 2. Daily turnover and volume from one of the brokerage tapes. Turnover (in CNY) is shown by the blue curve (left axis) and trading volume is shown by an orange curve (right axis).
Mathematics 10 02459 g002
Figure 3. Application of the model of state-space regressions (Equations (6) and (7)) to one brokerage tape. The horizontal coordinate is the observation day and the vertical coordinate is the price change bucket from 0 to 8 CNY in units of 0.5 CNY. (A) Temperature map of an original signal. Orange-yellow colors describe a positive day-to-day correlation, while light and dark blue colors describe a negative correlation. (B,D) predictions of the model, (C) model residuals.
Figure 3. Application of the model of state-space regressions (Equations (6) and (7)) to one brokerage tape. The horizontal coordinate is the observation day and the vertical coordinate is the price change bucket from 0 to 8 CNY in units of 0.5 CNY. (A) Temperature map of an original signal. Orange-yellow colors describe a positive day-to-day correlation, while light and dark blue colors describe a negative correlation. (B,D) predictions of the model, (C) model residuals.
Mathematics 10 02459 g003
Figure 4. Respective variances of the state vectors of correlations, n = 1 ÷ 16 for the trade imbalances for one randomly chosen tape in the units of σ2. We display original empirical data (blue), the predictor variance Pn (green), and the variance of the residuals Fn (red dash), all integrated for 484 trading days.
Figure 4. Respective variances of the state vectors of correlations, n = 1 ÷ 16 for the trade imbalances for one randomly chosen tape in the units of σ2. We display original empirical data (blue), the predictor variance Pn (green), and the variance of the residuals Fn (red dash), all integrated for 484 trading days.
Mathematics 10 02459 g004
Figure 5. The protocol of our tests. “Informed” trader observes an index and her own orders and trains the net. “Uninformed” trader uses only his own orders and a trained net to backtest the index.
Figure 5. The protocol of our tests. “Informed” trader observes an index and her own orders and trains the net. “Uninformed” trader uses only his own orders and a trained net to backtest the index.
Mathematics 10 02459 g005
Figure 6. Correlation of CNN predictions of the market sentiment from the unexpected monthly cost-of-trading (TC, brown, solid) with predictions of the unexpected order volume (OV, red, dash). ρε,TCρε,OV ≈ 0.913. For comparability, the scale of market sentiment was increased by 100.
Figure 6. Correlation of CNN predictions of the market sentiment from the unexpected monthly cost-of-trading (TC, brown, solid) with predictions of the unexpected order volume (OV, red, dash). ρε,TCρε,OV ≈ 0.913. For comparability, the scale of market sentiment was increased by 100.
Mathematics 10 02459 g006
Figure 7. Time series of average Amihud λ (Equation (9)) across all price buckets for a single trader.
Figure 7. Time series of average Amihud λ (Equation (9)) across all price buckets for a single trader.
Mathematics 10 02459 g007
Figure 8. Schematic description of the testing procedure. We use the data on λ, which we consider a measure of market friction, from the training sample to predict indexes from λ’s in five other samples. In the drawing above, the training sample comes first in calendar time but its position within the entire sample can be arbitrary.
Figure 8. Schematic description of the testing procedure. We use the data on λ, which we consider a measure of market friction, from the training sample to predict indexes from λ’s in five other samples. In the drawing above, the training sample comes first in calendar time but its position within the entire sample can be arbitrary.
Mathematics 10 02459 g008
Table 1. Summary statistics of Brokerage 5 for all the 484 trading days during 2009–2010. (Table 1 queries the data from original datasets. Traders are marked as 0 through 4. Volume data are rounded to one share. The stock price is expressed in the Chinese Yuan. The index B or S refers to a “Buy” or a “Sell” order on a given tape. We used sample volume variance as an indicator of volume dispersion. The daily standard deviation of the volume can be approximated as √484 times √SVV. There is no systematic difference between the first four tapes, while records in the fifth tape gravitate towards higher-priced stocks.
Table 1. Summary statistics of Brokerage 5 for all the 484 trading days during 2009–2010. (Table 1 queries the data from original datasets. Traders are marked as 0 through 4. Volume data are rounded to one share. The stock price is expressed in the Chinese Yuan. The index B or S refers to a “Buy” or a “Sell” order on a given tape. We used sample volume variance as an indicator of volume dispersion. The daily standard deviation of the volume can be approximated as √484 times √SVV. There is no systematic difference between the first four tapes, while records in the fifth tape gravitate towards higher-priced stocks.
Tape No./Type of OrderNo. TradesMin Price, CNYAvg. Price, CNYMax Price, CNYStd. Price, CNYAvg. Daily VolumeSample Volume Variance (SVV)
05B32,0412.348.77636.1711,04637,535
05S28,5132.368.6361.695.83969133,411
15B31,5502.39.91745.9410,79233,999
15S28,8642.239.5277.955.48996329,251
25B31,9882.579.84966.29950922,946
25S28,5842.539.5895.16.09846818,797
35B31,8242.2912.3669.916.35914429,004
35S26,8672.3111.9667.696.13763226,000
45B19,9753.9516.95123.59.24493938,321
45S15,0953.99.34193.89.34390665,500
Table 2. Correlations in the complex beta matrices of the linear regression.
Table 2. Correlations in the complex beta matrices of the linear regression.
(A) Correlation between columns of the beta matrix of coefficients for different traders.
Tape 0Tape 1Tape 2Tape 3Tape 4
Tape 010.96200.96300.96350.9633
Tape 1 10.96180.95630.9606
Tape 2 10.96010.9678
Tape 3 10.9618
Tape 4 1
(B) Correlation between rows of the beta matrix of coefficients for different traders.
Tape 0Tape 1Tape 2Tape 3Tape 4
Tape 010.94160.94240.95740.9431
Tape 1 10.94400.94210.9479
Tape 2 10.96350.9700
Tape 3 10.9452
Tape 4 1
Table 3. Mutual determination of variances (||ρij||=||rij2||) for the regression predictions and the regression residuals.
Table 3. Mutual determination of variances (||ρij||=||rij2||) for the regression predictions and the regression residuals.
Traders’# Prediction/Residuals12345
10.41260.00500.00030.00020.0006
20.00010.30810.04100.00520.0052
30.00030.02190.36390.03900.0374
40.00070.00390.01630.34080.0364
50.00010.00100.02190.01400.4027
Table 4. Correlation of predictions of the three indexes. Traders “a” and “b” were designated as informed about the index. The first two columns include informed traders in both training and prediction samples. The second pair of columns include data from an informed trader only in the training sample. The third pair of columns does not have an informed trader at all. For comparisons, we show the results obtained by a single 500 round and by changing ReLU into the TANH perceptron function. The lower left rectangle indicates data obtained from six TANH trials.
Table 4. Correlation of predictions of the three indexes. Traders “a” and “b” were designated as informed about the index. The first two columns include informed traders in both training and prediction samples. The second pair of columns include data from an informed trader only in the training sample. The third pair of columns does not have an informed trader at all. For comparisons, we show the results obtained by a single 500 round and by changing ReLU into the TANH perceptron function. The lower left rectangle indicates data obtained from six TANH trials.
(A) Correlations of monthly predictions based on market sentiment with actual data.
Predicted IndexSentiment
Traders Training-Predicta–cb–cb–cc–dc–dd–e
0.97560.93820.95050.95310.9120.9166
0.95270.96430.94740.96690.93460.9014
0.95450.90070.92320.93990.94550.9307
0.96840.95180.94830.96950.95690.966
0.97270.93640.96760.96990.93440.9431
0.96970.93710.94760.96360.95810.9225
Average0.96560.93810.94750.96050.94030.93
t-statistic 10%0.01390.03070.020380.0170.02490.0324
(B) Prediction of the returns on the Shanghai stock market.
Predicted IndexStock Index Return
a–cb–cb–cc–dc–dd–e
0.41690.40560.68820.70650.1820.2083
0.53630.37750.4370.42960.53410.182
0.1017−0.11940.24510.27890.4510.1656
0.05370.32930.27670.42790.25280.0906
0.39470.46360.6770.40510.06930.4107
0.4310.30380.38510.37170.32070.6272
       
Average-60.32230.29330.45150.43660.30170.2807
t-statistic at 10%0.28260.30230.27660.20670.24720.2889
Average-120.33850.271
t-statistic at 10%0.2280.2374
(C) Predictions of the yield of the bellwether Chinese bond.
Predicted IndexThe Yield on a 10-Year BOC Bond
a–cb–cb–cc–dc–dd–e
0.002−0.1425−0.0575−0.0281−0.1188−0.0296
0.0413−0.1958−0.05560.19210.19530.0502
−0.0729−0.1434−0.2029−0.10400.1149−0.1185
0.00190.00060.34160.062−0.05140.1155
0.14270.2798−0.0575−0.170−0.3020−0.4854
−0.5078−0.38180.10610.0081−0.1738−0.0613
       
Average0.17660.10780.05460.12360.10790.1032
t-statistic at 10%0.18190.2990.26230.17930.270.226
500 Rounds0.14150.2063−0.23050.1633−0.1407−0.2493
TANH−0.05170.3209−0.23020.07540.02279−0.0748
t-statistic at 10%, TANH 0.27550.3273
Table 5. Probabilities of null hypothesis H1 from a typical correlation test of index predictions from Amihud λ.
Table 5. Probabilities of null hypothesis H1 from a typical correlation test of index predictions from Amihud λ.
(A)
Trading Days in the SampleP from PearsonP from Spearman
120–2400.3640.932
180–3000.6050.976
240–3600.2270.062
300–4200.8470.342
360–4800.9750.903
(B)
Trading Days in the SampleP from PearsonP from Spearman
120–2400.40.712
180–300NA0.561
240–3600.4520.437
300–4200.8570.806
360–4800.6160.061
(C)
Trading Days in the SampleP from PearsonP from Spearman
120–2400.6220.285
180–3000.7210.565
240–3600.9260.743
300–4200.6880.554
360–4800.1670.143
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lerner, P.B. Fourier Integral Operator Model of Market Liquidity: The Chinese Experience 2009–2010. Mathematics 2022, 10, 2459. https://doi.org/10.3390/math10142459

AMA Style

Lerner PB. Fourier Integral Operator Model of Market Liquidity: The Chinese Experience 2009–2010. Mathematics. 2022; 10(14):2459. https://doi.org/10.3390/math10142459

Chicago/Turabian Style

Lerner, Peter B. 2022. "Fourier Integral Operator Model of Market Liquidity: The Chinese Experience 2009–2010" Mathematics 10, no. 14: 2459. https://doi.org/10.3390/math10142459

APA Style

Lerner, P. B. (2022). Fourier Integral Operator Model of Market Liquidity: The Chinese Experience 2009–2010. Mathematics, 10(14), 2459. https://doi.org/10.3390/math10142459

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop