ORIGINAL RESEARCH article

Front. Phys., 27 March 2014

Sec. Interdisciplinary Physics

Volume 2 - 2014 | https://doi.org/10.3389/fphy.2014.00016

Empirical properties of inter-cancellation durations in the Chinese stock market

  • 1. Department of Finance, School of Business, East China University of Science and Technology Shanghai, China

  • 2. Research Center for Econophysics, East China University of Science and Technology Shanghai, China

  • 3. Institute of Systems Engineering, College of Management and Economics, Tianjin University Tianjin, China

  • 4. China Center for Social Computing and Analytics, Tianjin University Tianjin, China

  • 5. Department of Mathematics, School of Science, East China University of Science and Technology Shanghai, China

Abstract

Order cancellation process plays a crucial role in the dynamics of price formation in order-driven stock markets and is important in the construction and validation of computational finance models. Based on the order flow data of 18 liquid stocks traded on the Shenzhen Stock Exchange in 2003, we investigate the empirical statistical properties of inter-cancellation durations in units of events defined as the waiting times between two consecutive cancellations. The inter-cancellation durations for both buy and sell orders of all the stocks favor a q-exponential distribution when the maximum likelihood estimation method is adopted; In contrast, both cancelled buy orders of 6 stocks and cancelled sell orders of 3 stocks prefer Weibull distribution when the non-linear least-squares estimation is used. Applying detrended fluctuation analysis (DFA), centered detrending moving average (CDMA) and multifractal detrended fluctuation analysis (MF-DFA) methods, we unveil that the inter-cancellation duration time series process long memory and multifractal nature for both buy and sell cancellations of all the stocks. Our findings show that order cancellation processes exhibit long-range correlated bursty behaviors and are thus not Poissonian.

1. Introduction

In an order-driven market, order submission and cancellation play the most important role in the process of price formation. For order submission process, lots of studies have been conducted to investigate the statistical properties of the ingredients of an order including order price [], order size or volume [], order direction [, ], and so on. Special attention has been paid to the probability distribution and memory effect of these ingredients and many stylized facts have been documented.

Order cancellation is a process of removing orders from the limit-order book which is a queue of limit orders waiting to be executed and constructed according to the price-time priority. If all orders placed at the best ask or best bid are cancelled, the mid-price defined as the mean value of the best ask and best bid will change. If cancellation occurs inside the limit order book, it causes changes of the structure of limit order book and has potential impact on price fluctuation.

The motivation of order cancellation is related to the non-execution (NE) risk or free option (FO) risk [, ], and the former is the major reason for cancelling limit orders []. NE risk arises when the current security price moves away from the submitting price. The orders submitted in the front of limit-order book cannot be transacted immediately, which makes the traders suffer opportunity cost. Traders may cancel the stale orders and resubmit more aggressive ones to increase the transaction probability. So in order to reduce NE risk, buy traders potentially drive the security price up, and sell traders drive the price down. FO risk arises when important news arrives. The intrinsic value of asset will be underestimated (for good news) or overestimated (for bad news) according to the current price. In order to prevent to be traded at the unfavorable price, traders will cancel their limit orders and resubmit unaggressive ones. So conversely, in order to reduce FO risk, buy traders potentially drive the price down, and sell traders drives the price up.

Since there are rare cancellation data recorded in the past, only a few literatures investigated the empirical regularities of order cancellation. With the development of information technology and computer science, it is possible to record the order flow data which enables us to analyze the statistical properties of order cancellation and construct cancellation models. Ni et al. [] investigated the empirical regularities of inter-cancellation duration of 22 stocks in the Chinese stock market, and made a conclusion that order cancellation is a non-Poisson process. Liu [] showed a simple model of order revision and cancellation, and found that the frequency of order cancellation is positively related to order submission risk and stock capitalization, but negatively related to bid-ask spread. In an order-driven model, Daniels et al. [] assumed that order cancellation follows a Poisson process, which makes the model having powerful predictions of stylized facts, such as price diffusion, price impact, and so on. In the empirical model proposed by Mike and Farmer [], the order cancellation process is determined by three independent factors, the position in the order book relative to the opposite best price, the imbalance of buy and sell orders in the limit-order book, and the total number of orders stored in the limit-order book. This cancellation model gives an excellent prediction for the life time of cancelled orders.

In financial markets, a widely studied subject is the recurrence interval defined as the waiting time between two consecutive events. Many scholars have analyzed the probability distribution of recurrence intervals of different variables such as returns, volatilities and trading volumes. However, the results are controversial. Power-law distribution [] and stretched exponential distribution [] are mainly selected to fit the probability density function (PDF) in different financial markets. Moreover, other distributions are also proposed for complement [, ]. It is interesting that the recurrence interval time series usually processes long memory [, , , , , , , , ]. Recurrence interval analysis has also been applied to other fields such as the energy dissipation rate in three-dimensional fully developed turbulence [].

Similar to recurrence interval, intertrade duration which is defined as the interval between two consecutive transactions is another important topic. Continuous-time random walk (CTRW) proposed by Montroll and Weiss [] has been widely utilized to deal with the intertrade duration in the financial time series []. Empirical results indicate that intertrade durations might follow a power-law distribution [, ], a stretched exponential or Weibull distribution [62], or a q-exponential distribution [62, 63]. On the other hand, some studies showed that the intertrade durations are neither exponentially distributed [, 61, 64, 65] nor power-law distributed [63].

The goodness of fit for Weibull distribution and q-exponential distribution has been estimated in the distribution of intertrade durations. Jiang et al. [62] studied the limit order data of 18 liquid stocks listed in the Chinese stock market, and showed that Weibull distribution gives better fitting than q-exponential distribution with the maximum likelihood estimation method, while q-exponential distribution outperforms Weibull distribution with the non-linear least-squares estimation method. Politi and Scalas [63] analyzed the tick-by-tick data set of DJIA stocks traded at NYSE in year 1999, and found that q-exponential distribution compares well to the Weibull distribution.

In this paper, we will study the statistical properties of inter-cancellation durations in event time for both cancelled buy and sell orders of 18 stocks listed on the Shenzhen Stock Exchange. The rest of paper is organized as follows. We study the probability distributions of the inter-cancellation durations based on the maximum likelihood estimation and non-linear least-square estimation methods. We further discuss the memory effect and multifractal nature.

2. Materials and methods

2.1. Data sets

Our analysis is based on the order flow data of 18 liquid stocks traded on the Shenzhen Stock Exchange in 2003. There were three periods in a trading day in 2003: opening call auction, cool period and continuous double auction. Opening call auction is held from 9:15 to 9:25 a.m., referring to the process of one-time centralized matching of buy and sell orders accepted during a specified period to generate the opening price at 9:25 in a trading day. Following the opening call auction, cool period is held from 9:25 to 9:30 a.m. when the Exchange is opened to orders routing from investors, but the orders or cancellations are not allowed to be processed. The main trading period is the continuous auction (9:30–11:30 a.m. and 13:00–15:00 p.m.), which refers to the process of continuous matching of both buy and sell orders on a one-by-one basis.

Our database records ultra-high-frequency order flow data whose time stamps are accurate to 0.01s. It contains the details of order placement and order cancellation. For example, the stock Ping An Bank Co., Ltd. (000001) contains 3,925,832 records in the whole year of 2003. So we can rebuild the limit-order book based on the prefect database according the trading rules. On the other hand, it is difficult to obtain this type of prefect database and we only have the data of 23 stocks in the whole year of 2003 among which 5 stocks have wrong records of order cancellation, and we select the rest 18 stocks to study the statistical properties of order cancellation. For the 18 stocks analyzed, they cover 9 CSRC (China Securities Regulatory Commission) Industries, such as finance and insurance, real estate, transportation, machinery, to list a few. On the other hand, in the year of 2003, the Chinese stock index first went up then fell down. Bull market and bear market were both existed in 2003. So the database we studied generally presents the situation of Chinese market.

The paper not only focuses on the cancellation data in the continuous auction, but also includes the cancellation data in the opening call auction and cool period. We count the cancellation numbers NC for both cancelled buy and sell orders of each stock, and then calculate the ratio r of NC to the number of all the orders NA (including both submitted orders and cancelled orders). The results are listed in Table 1. We find that the ratio r fluctuates within a wide range, being [0.097, 0.195] for cancelled buy orders and being [0.087, 0.189] for cancelled sell orders. An interesting feature shows that the ratio of cancelled buy orders approximates to cancelled sell orders for each stock, which implies that a large proportion of buy orders cancelled corresponds to a large number of cancellations for sell orders, and vice versa.

Table 1

StockCancelled buy ordersCancelled sell orders
NCrγdNCrγd
000001320,8720.1570.17412.12277,8780.1480.14813.73
000009185,0180.1760.18011.81188,7250.1640.18311.38
000012115,3350.1940.1879.84107,4430.1890.20010.62
00001660,7240.1700.17911.8959,7230.1540.15212.09
000021158,4700.1880.18210.70154,1570.1740.17310.98
00002443,0130.1770.17912.0346,1300.1560.16011.23
000066110,9790.1910.18810.60107,6410.1730.18210.85
00040694,9980.1640.16011.9292,8410.1570.17312.21
00042937,3320.1610.15812.6336,5850.1400.13812.96
00048833,1830.1480.16713.1834,4780.1450.15012.75
00053927,3220.1270.12414.7627,4870.1270.12814.88
00054119,9380.1490.17813.2920,1120.1290.13513.24
000550123,9480.1880.19310.99131,2700.1810.18210.31
00058127,5700.1520.16913.9929,9070.1310.12413.05
000625124,7360.1780.18711.38133,6270.1790.18710.65
00070966,2910.1520.14513.1664,6970.1360.14513.55
00072016,7670.0970.11617.4414,4410.0870.08920.29
00077844,1130.1510.15014.1748,0410.1310.13813.13

The statistics of inter-cancellation durations for both cancelled buy and sell orders of 18 stocks.

NC is the number of cancellation. r is the ratio of NC to NA. γ is the slope of the fitted line presented in Figure 1. 〈d〉 is the average inter-cancellation duration in units of events.

We further investigate the ratios r in each trading day for all the stocks. Figure 1 presents the linear relations between NC and NA for both cancelled buy and sell orders of two stocks. The slopes γ of the fitted lines are calculated using the least-squares fitting method, and the values of 18 stocks are listed in Table 1. It is evident that the value of γ is close to the value of r for each stock, which implies that the ratios of cancellation in each trading day are almost similar for both cancelled buy and sell orders of each stock.

Figure 1

In the paper, the inter-cancellation duration is defined as the interval between two consecutive cancellations in units of events, which reads

where t(i) is the event time when the i-th cancellation takes place. It is clear that inter-cancellation duration d(i) is the number of orders (including both buy orders and sell orders) submitted between the (i − 1)-th cancellation and the i-th cancellation. We calculate the average values 〈d〉 of inter-cancellation durations for both cancelled buy and sell orders of each stock, and depict the results in Table 1. According to the definition of inter-cancellation duration, we easily obtain the relation r > 1/〈d〉 which is confirmed by the data listed in the table. The reason is that the ratio r is defined based on a certain kind of orders (buy orders or sell orders) while the inter-cancellation duration d is considered as the number of both buy and sell orders submitted between two successive cancellations.

2.2. Multifractal detrended fluctuation analysis

Multifractals are ubiquitous in the nature and society [66], including financial time series. In this section, we investigate the multifractal properties of inter-cancellation durations applying the multifractal detrended fluctuation analysis (MF-DFA) method [67] which is generalized from the DFA method. The MF-DFA algorithm is described as follows.

Step 1. Consider an inter-cancellation duration series d(i), where i = 1, 2,…,N. Construct the cumulative sum sequence y(i) as follows,

Step 2. Divide the sequence y(i) into Ns disjoint segments with the same length s, where Ns = [N/s], and [x] is the largest integer not larger than x. Each segment can be denoted as yv such that yv(j) = y(ℓ + j) for 1 ⩽ js, and ℓ = (v − 1)s. Since the length of the inter-cancellation duration series N might not be a multiple of the segment size s, a remaining part (with the length smaller than s) at the end of sequence y(i) is not covered by the dividing procedure. We will select another Ns disjoint segments from the end of the series for compensating the remaining part, and then consider the 2Ns segments which covers the whole sequence y(i).

Step 3. In each segment yv, a polynomial function is utilized to represent the trend by the least-squares regression. The simplest function could be a line, and in the paper we adopt the linear function with 1 ⩽ js to remove the trend. The residual ϵv(j) in the segment yv can be calculated by

Step 4. The detrended fluctuation function F(v, s) of the segment yv is defined as follows,

Step 5. The q-th order overall fluctuation function Fq(s) is determined through

where q can take any real value except for q = 0. When q = 0, according to L'Hôspital's rule we have

Step 6. Varying the values of segment size s from 10 to [N/6], a power-law relation between the function Fq(s) and the size scale s is determined, which reads

According to the standard multifractal formalism, the multifractal scaling exponent τ(q) characterizes the multifractal nature, which reads

where Df is the fractal dimension of the geometric support of multifractal measure. For one dimensional time series analysis, we have Df = 1. If the scaling exponent τ(q) is a non-linear function of q, the series has multifractal nature. Finally, it is easy to obtain the singularity strength function α(q) and the multifractal spectrum f(α) via the Legendre transform, that is,

3. Results

3.1. Probability distribution

Probability distribution of financial variables has crucial implications on asset pricing and risk management. In this section, we focus on investigating the probability distributions of inter-cancellation durations for both cancelled buy and sell orders of 18 stocks. The probability density functions (PDFs) P(d) of four randomly chosen stocks are presented in Figure 2A.

Figure 2

According to the empirical results shown in the literature and the curve shape presented in Figure 2, we fit the distributions by Weibulls [, 62] and q-exponentials [, 62, 63, 68, 69]. For Weibull distributions, we have

where a is the scale parameter and b is the shape parameter. The q-exponential distributions can be described as follows:

where κ is the scale parameter and q is the shape parameter.

The maximum likelihood estimation (MLE) method is first applied to estimate the parameters of Weibull and q-exponential distributions. As we know, MLE method captures the major part of the fitting data. We find that it accounts for 64.9% (63.7%) in the range d ≤ 10 for cancelled buy (sell) orders. The parameters a and b of the Weibull distribution and the parameters κ and q of the q-exponential distribution estimated with the MLE method are listed in the left panel of Table 2.

Table 2

StockMLENLSE
Weibullq-exponentialWeibullq-exponential
abχWBLκqχqEabχWBLκqχqE
00000111.550.920.00759.471.220.00574.880.540.00927.971.260.0029
00000911.450.940.00649.591.180.00444.730.540.01147.341.280.0013
0000129.570.950.00908.001.180.00664.660.590.00987.071.210.0041
00001611.630.960.00739.991.160.00557.260.680.00629.111.170.0038
00002110.310.930.00858.521.200.00634.960.580.00927.671.220.0042
00002411.560.930.00829.641.200.00656.440.630.00688.461.220.0043
00006610.320.950.00828.671.180.00605.520.620.00887.621.210.0034
00040611.590.950.00699.851.170.00524.350.540.01087.431.260.0016
00042912.380.960.006710.741.150.00517.660.680.00649.291.180.0026
00048812.240.880.00999.671.270.00876.830.600.00519.251.240.0079
00053912.160.760.02167.911.510.02116.310.520.01618.751.360.0222
00054112.690.920.008610.601.200.00756.010.590.00708.601.240.0049
00055010.380.900.01028.321.240.00826.190.620.00678.141.220.0078
00058112.630.840.01139.601.320.01066.050.540.00658.671.310.0094
00062510.630.890.01058.371.270.00866.170.610.00668.471.220.0087
00070912.780.940.006110.881.170.00465.370.560.01048.341.260.0016
00072011.280.630.03944.572.000.03467.820.510.03699.241.470.0429
00077813.690.940.006611.711.170.00567.770.650.006110.251.190.0038
MEAN12.010.890.01049.561.270.00886.500.600.00928.791.250.0076

Parameters of the Weibull and q-exponential distributions based on the MLE and NLSE methods for cancelled buy orders of 18 stocks.

χ is the root mean square of the difference between the best fit and the empirical data.

Since Weibull distribution has the same number of parameters as q-exponential distribution, the root mean square χ of the difference between the best fit and the empirical data is applied to compare the performance of the two distributions, which is presented in Table 2 as well. It is clear that χqE of q-exponential distribution are smaller than χWBL of Weibull distribution (χqE < χWBL). So we conclude that q-exponential distribution outperforms Weibull distribution for the cancelled buy orders of each stock with the MLE method.

In order to capture the tail behavior of the distribution, we then utilize the non-linear least-squares estimation (NLSE) method to fit the distribution of cancelled buy orders, and the parameters of Weibull and q-exponential distributions are listed in the right panel of Table 2. The parameters a and b calculated from the NLSE method are all smaller than those with the MLE method for 18 stocks. For the parameter κ, 3 stocks out of 18 stocks own larger values with NLSE method, while for the parameter q, 12 stocks have larger values with NLSE method. We also select the root mean square χ to compare the performance of the two distributions with the NLSE method. According to the values of χ listed in Table 2, we find that the result is different from the MLE method. There are 6 stocks prefer Weibull distribution, and the rest 12 stocks are better fitted by q-exponential distribution. The mean values of four parameters for cancelled buy orders are also presented in the last row of Table 2. The mean value of the four parameters obtained from the MLE method are larger than those from the NLSE method.

With the same procedure mentioned above, we then analyze the probability distribution of cancelled sell orders with the MLE and NLSE methods, and obtain similar results. The parameters a and b of Weibull distribution and the parameters κ and q of q-exponential distribution are listed in Table 3. For the cancelled sell orders, the relation χqE < χWBL is satisfied for each stock when using the MLE method, which indicates that the distribution prefers q-exponential distribution than Weibull distribution with the MLE method. However, there are 3 stocks out of 18 stocks have smaller values of χWBL and prefer Weibull distribution with the NLSE method.

Table 3

StockMLENLSE
Weibullq-exponentialWeibullq-exponential
abχWBLκqχqEabχWBLκqχqE
00000113.050.900.005810.681.220.00426.040.560.00949.011.270.0015
00000911.110.950.00669.441.170.00464.920.570.01087.971.220.0015
00001210.310.950.00838.641.180.00614.920.590.00947.581.210.0036
00001611.960.970.006310.521.130.00475.120.590.01097.811.220.0021
00002110.720.950.00739.091.170.00534.630.570.01077.491.230.0016
00002411.160.990.00659.831.120.00464.320.560.01226.801.250.0035
00006610.650.950.00729.111.160.00525.800.630.00907.951.190.0025
00040611.880.950.007810.231.160.00655.040.570.00887.881.240.0036
00042912.870.980.006011.451.120.00449.440.770.004510.561.120.0029
00048811.960.890.01259.721.240.01198.670.700.006410.131.170.0124
00053913.630.860.009310.681.280.00887.300.570.00539.241.310.0073
00054112.860.940.007511.111.160.00654.690.520.01027.411.310.0047
00055010.060.940.00788.481.170.00574.220.550.01116.661.250.0016
00058112.430.910.008010.321.210.00705.700.560.00757.881.290.0041
00062510.330.940.00868.701.180.00673.710.530.01086.701.250.0024
00070913.240.950.005611.481.150.00435.220.550.01068.511.250.0023
00072014.370.660.02757.021.820.02599.870.530.024110.931.460.0304
00077812.900.960.006811.361.130.00558.270.710.005010.021.150.0036
MEAN11.980.930.00859.941.200.00706.020.590.00958.291.250.0048

Parameters of Weibull and q-exponential distributions based on the MLE and NLSE methods for cancelled sell orders of 18 stocks.

χ is the root mean square of the difference between the best fit and the empirical data.

Similar to the cancelled buy orders, the parameters a and b of Weibull distribution calculated from the NLSE method are smaller than those from the MLE method for cancelled sell orders of each stock. However, when considering q-exponential distribution, there are 2 stocks (000488 and 000720) having larger values of κ and smaller values of q with the NLSE method. The mean values of four parameters for the cancelled sell orders are presented in the last row of Table 3 as well. It is evident that the mean value of the four parameters obtained from the MLE method are larger than those from the NLSE method, except for the parameter q.

We rescale the inter-cancellation duration d to d/〈d〉 and the probability density function P(d) to P(d)〈d〉, where 〈d〉 is the mean value of inter-cancellation durations d. The rescaled PDFs of inter-cancellation durations for the same four stocks are presented in Figure 2B. We find that four rescaled curves collapse together, showing a perfect scaling behavior. Since the rescaled probability distribution has an excellent scaling, we aggregate all the inter-cancellation durations of 18 stocks together and treat them as an ensemble to obtain a better statistic. The rescaled PDFs of ensemble durations for both cancelled buy and sell orders are shown in Figure 3.

Figure 3

The parameters of Weibull and q-exponential distributions for ensemble inter-cancellation durations are calculated by the MLE and NLSE methods, respectively. When Using the MLE method, we obtain a = 11.21, b = 0.91 and κ = 9.13, q = 1.22 for cancelled buy orders, and a = 11.55, b = 0.93 and κ = 9.64, q = 1.19 for cancelled sell orders. When applying the NLSE method, we have a = 4.79, b = 0.54 and κ = 7.90, q = 1.25 for cancelled buy orders, and a = 4.80, b = 0.54 and κ = 7.91, q = 1.25 for cancelled sell orders. It is clear that the values of the parameters obtained from ensemble inter-cancellation durations are similar to the mean values presented in the last rows of Tables 2, 3, which confirms that the scaling behavior is truly existed. In addition, we find in Figure 3 that Weibull distribution evidently deviates from the empirical data in the tail with the MLE method, which is consistent with the relation χqE < χWBL for both cancelled buy and sell orders.

3.2. Memory effect

Another important issue about financial time series is the memory effect. Many methods have been proposed for quantitatively measuring the memory effort, such as the rescaled range (RS) analysis [70, 71], the fluctuation analysis (FA) [72], the wavelet transform module maxima (WTMM) method [73, 74], the detrended fluctuation analysis (DFA) [75], the detrending moving average (DMA) [76], and so on. Shao et al. [77] compared the performance of the FA, DFA, and DMA methods using different long-range correlated time series, and found that centred detrending moving average (CDMA) has the best performance and DFA is only slightly worse in some situations, while FA performs the worst. In this paper we apply the DFA and CDMA to investigate the memory effect of inter-cancellation duration series for both cancelled buy and sell orders of 18 stocks.

Figure 4 presents the detrended fluctuation functions FDFA(s) with respect to the size scale s using the DFA method for both cancelled buy and sell orders of four stocks, 000001, 000009, 000012, and 000021. Each curve reveals excellent power-law scaling over more than three orders of magnitude.

Figure 4

Applying the DFA method, the Hurst exponents H of 18 stocks are estimated according to the power-law relation FDFA(s) ~ sH, which are the slopes of solid lines shown in the log-log plot of Figure 4. We list the Hurst exponents for both cancelled buy and sell orders of 18 stocks in Table 4. The value of H for cancelled buy orders varies in the range [0.68, 0.82] with the mean value 〈H〉 = 0.76 ± 0.04, and for cancelled sell orders it varies from 0.68 to 0.85 with the mean value 〈H〉 = 0.76 ± 0.04. Since all the Hurst exponents are evidently larger than 0.5, we conclude that the inter-cancellation duration time series of both cancelled buy and sell orders process long memory.

Table 4

StockCancelled buy ordersCancelled sell orders
DFADMADFADMA
HHSFLHHSFLHHSFLHHSFL
0000010.8080.5000.8140.4990.8450.5010.8370.497
0000090.7950.5010.7960.4950.8410.5020.8390.494
0000120.7820.5030.7630.4950.8060.5000.7910.493
0000160.7390.5000.7260.4950.7260.5000.7370.499
0000210.7540.5010.7530.4960.7830.5000.7770.495
0000240.7730.5010.7680.4920.7310.5030.7460.497
0000660.7750.5000.7700.4980.7970.5000.7960.494
0004060.7340.5000.7380.4960.7830.5000.7840.494
0004290.7190.5030.7080.4940.6790.5020.6900.493
0004880.7640.5010.7510.4920.7470.5010.7210.498
0005390.8110.5010.8180.4940.7520.5000.7400.496
0005410.7410.5020.7220.5030.7150.5020.7120.496
0005500.7730.5000.7660.4930.7620.5000.7610.495
0005810.8160.5000.8090.4960.7680.5010.7570.500
0006250.7430.5000.7390.4900.7450.5010.7420.496
0007090.7470.5000.7420.4960.7380.5010.7310.500
0007200.6790.5000.6850.4920.7100.5010.7010.493
0007780.7210.5020.7190.4940.7240.5030.7170.494

Hurst exponents H of inter-cancellation durations for both cancelled buy and sell orders of 18 stocks based on the DFA and DMA methods.

HSFL is the mean Hurst exponent of 100 shuffled inter-cancellation durations.

Comparing with the backward detrending moving average (BDMA) and forward detrending moving average (FDMA), the centered detrending moving average (CDMA) performs better in one dimensional time series [78]. We then choose the CDMA method to estimate the memory effect of inter-cancellation duration series. The detrended fluctuation functions FCDMA(s) calculated from the CDMA method for both cancelled buy and sell orders of four stocks are depicted in Figure 5.

Figure 5

Perfect power-law scalings are observed in the log-log plot which implies that the relation FCDMA(s) ~ sH is well satisfied. The Hurst exponents H are the slopes of solid lines in the log-log plot. Using the least-squares fitting method, we calculate the Hurst exponents for both cancelled buy and sell orders of 18 stocks which are listed in Table 5 as well. It is obvious that the values of H obtained from the CDMA method are close to those calculated from DFA method. With all the Hurst exponents apparently larger than 0.5 with the two methods, we conclude that the inter-cancellation duration time series for both cancelled buy and sell orders of 18 stocks process long memory.

Table 5

StockCancelled buy ordersCancelled sell orders
ΔαΔαSFLRΔαΔαSFLR
0000010.810.29 ± 0.010.520.800.29 ± 0.010.51
0000090.530.32 ± 0.010.210.410.27 ± 0.020.13
0000120.450.30 ± 0.020.150.460.30 ± 0.020.16
0000160.370.26 ± 0.020.100.780.32 ± 0.020.46
0000210.670.28 ± 0.020.390.500.28 ± 0.020.22
0000240.970.32 ± 0.030.650.560.39 ± 0.020.17
0000660.360.28 ± 0.020.080.410.26 ± 0.020.15
0004060.840.37 ± 0.020.470.780.30 ± 0.020.48
0004290.620.29 ± 0.030.330.410.24 ± 0.030.17
0004880.770.36 ± 0.030.410.660.29 ± 0.030.37
0005391.070.53 ± 0.040.540.920.43 ± 0.030.49
0005410.700.44 ± 0.040.260.820.53 ± 0.040.29
0005500.770.28 ± 0.020.490.580.33 ± 0.020.24
0005810.840.47 ± 0.030.370.690.42 ± 0.040.26
0006250.800.30 ± 0.020.500.610.34 ± 0.020.27
0007090.870.34 ± 0.030.540.960.36 ± 0.030.59
0007201.280.70 ± 0.060.581.380.60 ± 0.060.78
0007780.620.31 ± 0.020.310.500.24 ± 0.020.25

The width of the multifractal spectra Δα of inter-cancellation durations for both cancelled buy and sell orders of 18 stocks based on the MF-DFA method.

ΔαSFLis the mean width of 100 shuffled inter-cancellation durations. R is the residual of spectrum width by removing the shuffled width ΔαSFLfrom the original one Δα.

On the other hand, the memory effect might be affected by the distribution of inter-cancellation durations. In order to test this hypothesis, we first shuffle the inter-cancellation duration series of each stock for 100 times, and then calculate the Hurst exponents HSFL for each shuffling series based on both DFA and DMA methods, respectively. The mean Hurst exponents HSFL of 100 shuffling series for both cancelled buy and sell orders of each stock are shown in Table 4. We find that the values of HSFL are all extremely close to 0.5, being significant smaller than the original ones H. So we conclude that the distribution of inter-cancellation duration series has little impact on its memory effect, and confirm that inter-cancellation duration series of both cancelled buy and sell orders truly exhibit significant long memory for all the 18 stocks.

Memory effect presents the time persistence of inter-cancellation durations. It reflects the clustering behavior of order cancellation which is caused by traders' similar reactions to the market. For example, when good news arrives, traders will immediately cancel their limit orders in order to avoid being transacted at the unfavorable price. Many cancellations occur in a short period, which results to the clustering behavior and long memory effect of order cancellation.

3.3. Multifractal nature

We calculate the q-th order fluctuation functions Fq(s) of inter-cancellation durations for both cancelled buy and sell orders of two stocks, 000009 and 000012, and present the fluctuation functions Fq(s) in Figure 6. We find that the function Fq(s) has a excellent power-law scaling with respect to the scale size s. Using the least-squares fitting method, we obtain the slopes h(q) for q = −4, −2, 0, 2, 4, respectively.

Figure 6

Figure 7 presents the scaling exponents τ(q) with respect to the order q and the multifractal spectra f(α) as a function of the singularity strength α for both cancelled buy and sell orders of the two stocks. We observe that the function τ(q) is non-linear with respect to q, which illustrates that the inter-cancellation durations own multifractal nature.

Figure 7

In addition, the strength of multifractal can also be measured by the width of the multifractal spectrum f(α) (Δα = αmax − αmin), and a larger value of Δα corresponds to stronger multifractal. We calculate the singularity width Δα for both cancelled buy and sell orders of 18 stocks and list the results in Table 5. The value Δα of cancelled buy orders varies from 0.37 to 1.28 with the mean value 〈Δα〉 = 0.74 ± 0.24, and for cancelled sell orders the value Δα varies in the range [0.41, 1.38] with the mean value 〈Δα〉 = 0.68 ± 0.25. Since all the values of Δα larger than zero, we consider that the inter-cancellation duration series for both cancelled buy and sell orders of 18 stocks have multifractal nature, which is consistent with the results obtained from the scaling exponent τ(q).

Similar to the memory effect, the probability distribution might have influence upon the multifractal nature of inter-cancellation durations. In order to test the influence of distribution, we shuffle the inter-cancellation durations for 100 times to test this influence. For each shuffling series, the width of the multifractal spectrum ΔαSFL is obtained based on the MF-DFA method. The mean values of 100 shuffled series for 18 stocks are listed in Table 5. The values of ΔαSFL are clearly larger than zero, which indicates that the distribution of inter-cancellation durations reliably generates multifractal. We define the residual of spectrum width R through removing the shuffled width ΔαSFL from the original one Δα, that is, R = Δα − ΔαSFL, and list the values R of 18 stocks in Table 5. Since the values of R are evidently larger than zero, we conclude that inter-cancellation durations process multifractal nature for both cancelled buy and sell orders of all the stocks.

4. Discussion

Order cancellation is an important issue in the dynamics of price formation in financial markets. We have carried out empirical investigations on the statistical properties of inter-cancellation durations (in units of events) using the order flow data of 18 liquid stocks traded on the Shenzhen Stock Exchange in the whole year of 2003.

We first study the probability distributions of inter-cancellation durations for both cancelled buy and sell orders, and find that the rescaled probability density functions have a scaling behavior. When fitting the probability distributions by Weibull and q-exponential distributions, we find that both cancelled buy and sell orders prefer q-exponential distribution with MLE method. However, applying the NLSE method, we find that cancelled buy orders of 6 stocks and cancelled sell orders of 3 stocks prefer Weibull distribution which is different from the result obtained from the MLE method.

We then investigate the memory effect of inter-cancellation durations based on the detrended fluctuation analysis (DFA) and centered detrending moving average (CDMA) methods. Using the DFA method we obtain the average Hurst exponent of 18 stocks 〈H〉 = 0.76 for both cancelled buy and sell orders, and with the CDMA method it is 〈H〉 = 0.75 for both cancelled buy and sell orders. According to the results from these two methods, it is evident that the inter-cancellation duration series processes the same strength of long memory for both cancelled buy and sell orders.

Finally, we investigate the multifractal properties applying the multifractal detrended fluctuation analysis (MF-DFA) method. We find that the average width of multifractal spectrum 〈Δα〉 = 0.74 for cancelled buy orders of 18 stocks and it is 〈Δα〉 = 0.68 for cancelled sell orders. So we conclude that the inter-cancellation duration series has multifractal nature, and inter-cancellation duration series of buy orders has little stronger multifractality than cancelled sell orders. Our findings indicate that the cancellation process has a bursty behavior and possesses long-range correlations. Such non-Poisson behaviors have been unveils in many other human dynamics [79].

Conflict of interest statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Statements

Acknowledgments

This work was partly supported by National Natural Science Foundation of China (Grants No. 71101052, 71131007 and 11075054), Shanghai Rising Star (Follow-up) Program (Grant No. 11QH1400800), and the Fundamental Research Funds for the Central Universities.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  • 1.

    ZovkoIFarmerJD. The power of patience: a behavioural regularity in limit-order placement. Quant Finance (2002) 2:38792. 10.1088/1469-7688/2/5/308

  • 2.

    BouchaudJPMézardMPottersM. Statistical properties of stock order books: empirical results and models. Quant Finance (2002) 2:2516. 10.1088/1469-7688/2/4/301

  • 3.

    PottersMBouchaudJP. More statistical properties of order books and price impact. Physica A (2003) 324:13340. 10.1016/S0378-4371(02)01896-4

  • 4.

    RanaldoA. Order aggressiveness in limit order book markets. J Financial Markets (2004) 7:5374. 10.1016/S1386-4181(02)00069-1

  • 5.

    MaskawaJI. Correlation of coming limit price with order book in stock markets. Physica A (2007) 383:905. 10.1016/j.physa.2007.04.091

  • 6.

    LilloF. Limit order placement as an utility maximization problem and the origin of power law distribution of limit order prices. Eur Phys J B (2007) 55:4539. 10.1140/epjb/e2007-00067-9

  • 7.

    MikeSFarmerJD. An empirical behavioral model of liquidity and volatility. J Econ Dyn Control. (2008) 32:20034. 10.1016/j.jedc.2007.01.025

  • 8.

    GuGFChenWZhouWX. Empirical regularities of order placement in the Chinese stock market. Physica A (2008) 387:317382. 10.1016/j.physa.2008.01.114

  • 9.

    GopikrishnanPPlerouVGabaixXStanleyHE. Statistical properties of share volume traded in financial markets. Phys Rev E (2000) 62:R44936. 10.1103/PhysRevE.62.R4493

  • 10.

    PlerouVGopikrishnanPGabaixXAmaralLANStanleyHE. Price fluctuations, market activity and trading volume. Quant Finance (2001) 1:2629. 10.1088/1469-7688/1/2/308

  • 11.

    FarmerJDLilloF. On the origin of power-law tails in price fluctuations. Quant Finance (2004) 4:C711. 10.1088/1469-7688/4/1/C01

  • 12.

    QueirosSMD. On the emergence of a generalised Gamma distribution: application to traded volume in financial markets. Euro Phys Lett. (2005) 71:33945. 10.1209/epl/i2005-10109-0

  • 13.

    de SouzaJMoyanoLGQueirosSMD. On statistical properties of traded volume in financial markets. Eur Phys J B. (2006) 50:1658. 10.1140/epjb/e2006-00130-1

  • 14.

    PlerouVStanleyHE. Tests of scaling and universality of the distributions of trade size and share volume: evidence from three distinct markets. Phys Rev E (2007) 76:046109. 10.1103/PhysRevE.76.046109

  • 15.

    MuGHChenWKertészJZhouWX. Preferred numbers and the distributions of trade sizes and trading volumes in the Chinese stock market. Eur Phys J B (2009) 68:14552. 10.1140/epjb/e2009-00059-9

  • 16.

    GuGFRenFNiXHChenWZhouWX. Empirical regularities of opening call auction in Chinese stock market. Physica A (2010) 389:27886. 10.1016/j.physa.2009.09.019

  • 17.

    ZhouWX. Determinants of immediate price impacts at the trade level in an emerging order-driven market. New J Phys. (2012) 14:023055. 10.1088/1367-2630/14/2/023055

  • 18.

    ZhouWX. Universal price impact functions of individual trades in an order-driven market. Quant Finance (2012) 12:125363. 10.1080/14697688.2010.504733

  • 19.

    LilloFFarmerJD. The long memory of the efficient market. Stud Nonlin Dyn Econometr. (2004) 8:133. 10.2202/1558-3708.1226

  • 20.

    GuGFZhouWX. Emergence of long memory in stock volatilities from a modified Mike-Farmer model. Euro Phys Lett. (2009) 86:48002. 10.1209/0295-5075/86/48002

  • 21.

    GuGFZhouWX. On the probability distribution of stock returns in the Mike-Farmer model. Eur Phys J B. (2009) 67:58592. 10.1140/epjb/e2009-00052-4

  • 22.

    GriffithsMDSmithBFTurnbullDASWhiteRW. The costs and determinants of order aggressiveness. J Financial Econ. (2000) 56:6588. 10.1016/S0304-405X(99)00059-8

  • 23.

    FongKYLLiuWM. Limit order revisions. J Bank Finance (2010) 34:187385. 10.1016/j.jbankfin.2009.12.010

  • 24.

    NiXHJiangZQGuGFRenFChenWZhouWX. Scaling and memory in the non-Poisson process of limit order cancelation. Physica A (2010) 389:275161. 10.1016/j.physa.2010.02.040

  • 25.

    LiuWM. Monitoring and limit order submission risks. J Financial Markets (2009) 12:10741. 10.1016/j.finmar.2008.04.003

  • 26.

    DanielsMGFarmerJDGillemotLIoriGSmithE. Quantitative model of price diffusion and market friction based on trading as a mechanistic random process. Phys Rev Lett. (2003) 90:108102. 10.1103/PhysRevLett.90.108102

  • 27.

    KaizojiTKaizojiM. Power law for the calm-time interval of price changes. Physica A (2004) 336:56370. 10.1016/j.physa.2003.12.054

  • 28.

    YamasakiKMuchnikLHavlinSBundeAStanleyHE. Scaling and memory in volatility return intervals in financial markets. Proc Natl Acad Sci USA (2005) 102:94248. 10.1073/pnas.0502613102

  • 29.

    LeeJWLeeKERikvoldPA. Waiting-time distribution for Korean stock-market index KOSPI. J Korean Phys Soc. (2006) 48:S1236. 10.3938/jkps.48.123

  • 30.

    RenFZhouWX. Recurrence interval analysis of trading volumes. Phys Rev E (2010) 81:066107. 10.1103/PhysRevE.81.066107

  • 31.

    WangFZYamasakiKHavlinSStanleyHE. Scaling and memory of intraday volatility return intervals in stock markets. Phys Rev E (2006) 73:026117. 10.1103/PhysRevE.73.026117

  • 32.

    WangFWeberPYamasakiKHavlinSStanleyHE. Statistical regularities in the return intervals of volatility. Eur Phys J B (2007) 55:12333. 10.1140/epjb/e2006-00356-9

  • 33.

    Vodenska-ChitkushevIWangFZWeberPYamasakiKHavlinSStanleyHE. Comparison between volatility return intervals of the S&P 500 index and two common models. Eur Phys J B (2008) 61:21723. 10.1140/epjb/e2008-00066-4

  • 34.

    JungWSWangFZHavlinSKaizojiTMoonHTStanleyHE. Volatility return intervals analysis of the Japanese market. Eur Phys J B (2008) 62:1139. 10.1140/epjb/e2008-00123-0

  • 35.

    WangFZYamasakiKHavlinSStanleyHE. Indication of multiscaling in the volatility return intervals of stock markets. Phys Rev E (2008) 77:016109. 10.1103/PhysRevE.77.016109

  • 36.

    QiuTGuoLChenG. Scaling and memory effect in volatility return interval of the Chinese stock market. Physica A (2008) 387:68128. 10.1016/j.physa.2008.09.002

  • 37.

    RenFZhouWX. Multiscaling behavior in the volatility return intervals of Chinese indices. Euro Phys Lett. (2008) 84:68001. 10.1209/0295-5075/84/68001

  • 38.

    RenFGuGFZhouWX. Scaling and memory in the return intervals of realized volatility. Physica A (2009) 388:478796. 10.1016/j.physa.2009.08.009

  • 39.

    RenFGuoLZhouWX. Statistical properties of volatility return intervals of Chinese stocks. Physica A (2009) 388:88190. 10.1016/j.physa.2008.12.005

  • 40.

    WangFWangJ. Statistical analysis and forecasting of return interval for SSE and model by lattice percolation system and neural network. Comput Ind Eng. (2012) 62:198205. 10.1016/j.cie.2011.09.007

  • 41.

    MengHRenFGuGFXiongXZhangYJZhouWXet al. Effects of long memory in the order submission process on the properties of recurrence intervals of large price fluctuations. Euro Phys Lett. (2012) 98:38003. 10.1209/0295-5075/98/38003

  • 42.

    XieWJJiangZQZhouWX. Extreme value statistics and recurrence intervals of NYMEX energy futures volatility. Econ Model. (2014) 36:817. 10.1016/j.econmod.2013.09.011

  • 43.

    ZhangJHWangJShaoJG. Finite-range contact process on the market return intervals distributions. Adv Complex Sys. (2010) 13:64357. 10.1142/S0219525910002797

  • 44.

    JeonWMoonHTOhGYangJSJungWS. Return intervals analysis of the Korean stock market. J Korean Phys Soc. (2010) 56:9225. 10.3938/jkps.56.922

  • 45.

    RenFZhouWX. Recurrence interval analysis of high-frequency financial returns and its application to risk estimation. New J Phys. (2010) 12:075030. 10.1088/1367-2630/12/7/075030

  • 46.

    YamasakiKMuchnikLHavlinSBundeAStanleyHE. Scaling and memory in return loss intervals: application to risk estimation In: Takayasu H, ed Practical Fruits of Econophysics. (Berlin: Springer-Verlag) (2006). p. 4351. 10.1007/4-431-28915-1_7

  • 47.

    LiuCJiangZQRenFZhouWX. Scaling and memory in the return intervals of energy dissipation rate in three-dimensional fully developed turbulence. Phys Rev E (2009) 80:046304. 10.1103/PhysRevE.80.046304

  • 48.

    MontrollEWWeissGH. Random walks on lattices. II. J Math Phys. (1965) 6:16781. 10.1063/1.1704269

  • 49.

    ScalasEGorenfloRMainardiF. Fractional calculus and continuous-time finance. Physica A (2000) 284:37684. 10.1016/S0378-4371(00)00255-7

  • 50.

    MainardiFRabertoMGorenfloRScalasE. Fractional calculus and continuous-time finance II: the waiting-time distribution. Physica A (2000) 287:46881. 10.1016/S0378-4371(00)00386-1

  • 51.

    MasoliverJMonteroMWeissGH. Continuous-time random walk model for financial distribution. Phys Rev E. (2003) 67:021112. 10.1103/PhysRevE.67.021112

  • 52.

    KimKYoonSM. Dynamic behavior of continuous tick data in futures exchange market. Fractals (2003) 11:1316. 10.1142/S0218348X03002142

  • 53.

    ScalasE. The application of continuous-time random walks in finance and economics. Physica A (2006) 362:22539. 10.1016/j.physa.2005.11.024

  • 54.

    MasoliverJMonteroMPerellóJWeissGH. The continunous time random walk formalism in financial markets. J Econ Behav Org. (2006) 61:57798. 10.1016/j.jebo.2004.07.015

  • 55.

    SabatelliLKeatingSDudleyJRichmondP. Waiting time distributions in financial markets. Eur Phys J B (2002) 27:2735. 10.1140/epjb/e20020151

  • 56.

    YoonSMChoiJSLeeCCYumMKKimK. Dynamical volatilities for yen-dollar exchange rates. Physica A (2006) 359:56975. 10.1016/j.physa.2005.05.089

  • 57.

    RabertoMScalasEMainardiF. Waiting-times and returns in high-frequency financial data: an empirical study. Physica A (2002) 314:74955. 10.1016/S0378-4371(02)01048-8

  • 58.

    BartiromoR. Dynamics of stock prices. Phys Rev E (2004) 69:067108. 10.1103/PhysRevE.69.067108

  • 59.

    IvanovPCYuenAPodobnikBLeeYK. Common scaling patterns in intertrade times of U. S. stocks. Phys Rev E (2004) 69:056107. 10.1103/PhysRevE.69.056107

  • 60.

    EislerZKertészJ. Size matters: some stylized facts of the stock market revisited. Eur Phys J B (2006) 51:14554. 10.1140/epjb/e2006-00189-6

  • 61.

    SazukaN. On the gap between an empirical distribution and an exponential distribution of waiting times for price changes in a financial market. Physica A (2007) 376:5006. 10.1016/j.physa.2006.10.094

  • 62.

    JiangZQChenWZhouWX. Scaling in the distribution of intertrade durations of Chinese stocks. Physica A (2008) 387:581825. 10.1016/j.physa.2008.06.039

  • 63.

    PolitiMScalasE. Fitting the empirical distribution of intertrade durations. Physica A (2008) 387:202534. 10.1016/j.physa.2007.11.018

  • 64.

    ScalasEGorenfloRLuckockHMainardiFMantelliMRabertoM. Anomalous waiting times in high-frequency financial data. Quant Finance (2004) 4:695702. 10.1080/14697680500040413

  • 65.

    ScalasEGorenfloRLuckockHMainardiFMantelliMRabertoM. On the intertrade waiting-time distribution. Financ Lett. (2005) 3:695702.

  • 66.

    MandelbrotBB. The Fractal Geometry of Nature. New York, NY: W. H. Freeman; (1983).

  • 67.

    KantelhardtJWZschiegnerSAKoscielny-BundeEHavlinSBundeAStanleyHE. Multifractal detrended fluctuation analysis of nonstationary time series. Physica A (2002) 316:87114. 10.1016/S0378-4371(02)01383-3

  • 68.

    NadarajahSKotzS. q exponential is a Burr distribution. Phys Lett A (2006) 359:5779. 10.1016/j.physleta.2006.07.035

  • 69.

    NadarajahSKotzS. On the q-type distributions. Physica A (2007) 377:4658. 10.1016/j.physa.2006.11.054

  • 70.

    HurstHE. Long-term storage capacity of reservoirs. Trans Amer Soc Civil Eng. (1951) 116:770808.

  • 71.

    MandelbrotBBVan NessJW. Fractional Brownian motions, fractional noises and applications. SIAM Rev. (1968) 10:42237. 10.1137/1010093

  • 72.

    PengCKBuldyrevSVGoldbergerALHavlinSSciortinoFSimonsMet al. Long-range correlations in nucleotide sequences. Nature (1992) 356:168170. 10.1038/356168a0

  • 73.

    HolschneiderM. On the wavelet transformation of fractal objects. J Stat Phys. (1988) 50:96393. 10.1007/BF01019149

  • 74.

    MuzyJFBacryEArnéodoA. Wavelets and multifractal formalism for singular signals: application to turbulence data. Phys Rev Lett. (1991) 67:35158. 10.1103/PhysRevLett.67.3515

  • 75.

    PengCKBuldyrevSVHavlinSSimonsMStanleyHEGoldbergerAL. Mosaic organization of DNA nucleotides. Phys Rev E (1994) 49:16859. 10.1103/PhysRevE.49.1685

  • 76.

    AlessioECarboneACastelliGFrappietroV. Second-order moving average and scaling of stochastic time series. Eur Phys J B (2002) 27:197200. 10.1140/epjb/e20020150

  • 77.

    ShaoYHGuGFJiangZQZhouWXSornetteD. Comparing the performance of FA, DFA and DMA using different synthetic long-range correlated time series. Sci Rep. (2012) 2:835. 10.1038/srep00835

  • 78.

    ArianosSCarboneA. Detrending moving average algorithm: a closed-form approximation of the scaling law. Physica A (2007) 382:915. 10.1016/j.physa.2007.02.074

  • 79.

    JiangZQXieWJLiMXPodobnikBZhouWXStanleyHE. Calling patterns in human communication dynamics. Proc Natl Acad Sci USA. (2013) 110:16005. 10.1073/pnas.1220433110

Summary

Keywords

econophysics, order flow, inter-cancellation duration, probability distribution, memory effect, multifractal nature

Citation

Gu G-F, Xiong X, Zhang W, Zhang Y-J and Zhou W-X (2014) Empirical properties of inter-cancellation durations in the Chinese stock market. Front. Physics 2:16. doi: 10.3389/fphy.2014.00016

Received

19 July 2013

Accepted

05 March 2014

Published

27 March 2014

Volume

2 - 2014

Edited by

Hai-Qing Lin, Beijing Computational Science Research Center, China

Reviewed by

Dror Y. Kenett, Boston University, USA; Shiquan Su, University of Tennessee at Knoxville, USA

Copyright

*Correspondence: Wei-Xing Zhou, School of Business, East China University of Science and Technology, 130 Meilong Road, Shanghai 200237, China e-mail:

This article was submitted to Interdisciplinary Physics, a section of the journal Frontiers in Physics.

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics