Scaling in the distribution of intertrade durations of Chinese stocks

The distribution of intertrade durations, defined as the waiting times between two consecutive transactions, is investigated based upon the limit order book data of 23 liquid Chinese stocks listed on the Shenzhen Stock Exchange in the whole year 2003. A scaling pattern is observed in the distributions of intertrade durations, where the empirical density functions of the normalized intertrade durations of all 23 stocks collapse onto a single curve. The scaling pattern is also observed in the intertrade duration distributions for filled and partially filled trades and in the conditional distributions. The ensemble distributions for all stocks are modeled by the Weibull and the Tsallis q-exponential distributions. Maximum likelihood estimation shows that the Weibull distribution outperforms the q-exponential for not-too-large intertrade durations which account for more than 98.5% of the data. Alternatively, nonlinear least-squares estimation selects the q-exponential as a better model, in which the optimization is conducted on the distance between empirical and theoretical values of the logarithmic probability densities. The distribution of intertrade durations is Weibull followed by a power-law tail with an asymptotic tail exponent close to 3.     

Long-term correlations and multifractal nature in the intertrade durations of a liquid Chinese stock and its warrant

The intertrade duration of equities is an important financial measure, characterizing trading activities; it is defined as the waiting time between successive trades of an equity. Using the ultrahigh-frequency data of a liquid Chinese stock and its associated warrant, we perform a comparative investigation of the statistical properties of their intertrade duration time series. The distributions of the two equities can be better described by the shifted power-law form than the Weibull form, and their scaled distributions do not collapse onto a single curve. Although the intertrade durations of the two equities have very different magnitude, their intraday patterns exhibit very similar shapes. Both detrended fluctuation analysis (DFA) and detrending moving average analysis (DMA) show that the 1 min intertrade duration time series of the two equities are strongly correlated. In addition, both multifractal detrended fluctuation analysis (MFDFA) and multifractal detrending moving average analysis (MFDMA) unveil that the 1 min intertrade durations possess multifractal nature. However, the difference between the two singularity spectra of the two equities obtained from the MFDMA is much smaller than that from the MFDFA.     

Message transfer in a communication network

We study message transfer in a 2-d communication network of regular nodes and randomly distributed hubs. We study both single message transfer and multiple message transfer on the lattice. The average travel time for single messages travelling between source and target pairs of fixed separations shows q-exponential behaviour as a function of hub density with a characteristic power-law tail, indicating a rapid drop in the average travel time as a function of hub density. This power-law tail arises as a consequence of the log-normal distribution of travel times seen at high hub densities. When many messages travel on the lattice, a congestion-decongestion transition can be seen. The waiting times of messages in the congested phase show a Gaussian distribution, whereas the decongested phase shows a log-normal distribution. Thus the waiting time distributions carry the signature of congested or decongested behaviour.      

Tail distribution of index fluctuations in World markets

We have investigated the tail distribution of the daily fluctuations in 202 different indices in the stock markets of 59 countries for the time span of the last 20 years. Power law, log-normal, Weibull, exponential and power law with exponential cutoff distributions are considered as possible candidates for the tail distribution of the normalized returns. It is found that the power exponent depends strongly on the choice of the tail threshold and a sizeable number of indices can be better fitted by a distribution function other than the power law at the region that has power law exponent of 3. Also, we have found that the power exponent is not an indicator of the maturity of the market.     

Transition in the waiting-time distribution of price-change events in a global socioeconomic system

The goal of developing a firmer theoretical understanding of inhomogeneous temporal processes–in particular, the waiting times in some collective dynamical system–is attracting significant interest among physicists. Quantifying the deviations between the waiting-time distribution and the distribution generated by a random process may help unravel the feedback mechanisms that drive the underlying dynamics. We analyze the waiting-time distributions of high-frequency foreign exchange data for the best executable bid–ask prices across all major currencies. We find that the lognormal distribution yields a good overall fit for the waiting-time distribution between currency rate changes if both short and long waiting times are included. If we restrict our study to long waiting times, each currency pair’s distribution is consistent with a power-law tail with exponent near to 3.5. However, for short waiting times, the overall distribution resembles one generated by an archetypal complex systems model in which boundedly rational agents compete for limited resources. Our findings suggest that a gradual transition arises in trading behavior between a fast regime in which traders act in a boundedly rational way and a slower one in which traders’ decisions are driven by generic feedback mechanisms across multiple timescales and hence produce similar power-law tails irrespective of currency type.     

Degree distributions of the visibility graphs mapped from fractional Brownian motions and multifractal random walks

The dynamics of a complex system is usually recorded in the form of time series, which can be studied through its visibility graph from a complex network perspective. We investigate the visibility graphs extracted from fractional Brownian motions and multifractal random walks, and find that the degree distributions exhibit power-law behaviors, in which the power-law exponent α is a linear function of the Hurst index H of the time series. We also find that the degree distribution of the visibility graph is mainly determined by the temporal correlation of the original time series with minor influence from the possible multifractal nature. As an example, we study the visibility graphs constructed from three Chinese stock market indexes and unveil that the degree distributions have power-law tails, where the tail exponents of the visibility graphs and the Hurst indexes of the indexes are close to the α∼H linear relationship.     

Detrended fluctuation analysis of intertrade durations

The intraday pattern, long memory, and multifractal nature of the intertrade durations, which are defined as the waiting times between two consecutive transactions, are investigated based upon the limit order book data and order flows of 23 liquid Chinese stocks listed on the Shenzhen Stock Exchange in 2003. An inverse U-shaped intraday pattern in the intertrade durations with an abrupt drop in the first minute of the afternoon trading is observed. Based on a detrended fluctuation analysis, we find a crossover of power-law scaling behaviors for small box sizes (trade numbers in boxes) and large box sizes and strong evidence in favor of long memory in both regimes. In addition, the multifractal nature of intertrade durations in both regimes is confirmed by a multifractal detrended fluctuation analysis for individual stocks with a few exceptions in the small-duration regime. The intraday pattern has little influence on the long memory and multifractality.     

Stretched-Gaussian asymptotics of the truncated Lévy flights for the diffusion in nonhomogeneous media

The Lévy, jumping process, defined in terms of the jumping size distribution and the waiting time distribution, is considered. The jumping rate depends on the process value. The fractional diffusion equation, which contains the variable diffusion coefficient, is solved in the diffusion limit. That solution resolves itself to the stretched Gaussian when the order parameter μ→2. The truncation of the Lévy flights, in the exponential and power-law form, is introduced and the corresponding random walk process is simulated by the Monte Carlo method. The stretched Gaussian tails are found in both cases. The time which is needed to reach the limiting distribution strongly depends on the jumping rate parameter. When the cutoff function falls slowly, the tail of the distribution appears to be algebraic.     

Kinetic market models with single commodity having price fluctuations

We study here numerically the behavior of an ideal gas like model of markets having only one non-consumable commodity. We investigate the behavior of the steady-state distributions of money, commodity and total wealth, as the dynamics of trading or exchange of money and commodity proceeds, with local (in time) fluctuations in the price of the commodity. These distributions are studied in markets with agents having uniform and random saving factors. The self-organizing features in money distribution are similar to the cases without any commodity (or with consumable commodities), while the commodity distribution shows an exponential decay. The wealth distribution shows interesting behavior: gamma like distribution for uniform saving propensity and has the same power-law tail, as that of the money distribution, for a market with agents having random saving propensity.     

供应链型网络中双幂律分布模型

考察了供应链网络的基本特征，提出了节点到达过程是更新过程、新增入边和出边数是具有Bernoulli分布随机变量的供应链型有向网络.研究了这类网络节点的瞬态度分布和稳态平均度分布.利用更新过程理论对这类网络进行了分析，获得了网络节点瞬态度分布和网络稳态平均度分布的解析表达式.分析表明, 虽然这类网络节点的稳态度分布不存在，但是网络的稳态平均度分布具有双向幂律性. 关键词： 复杂网络 入度 出度 度分布     

The Weibull-log Weibull distribution for interoccurrence times of earthquakes

By analyzing the Japan Meteorological Agency (JMA) seismic catalog for different tectonic settings, we have found that the probability distributions of time intervals between successive earthquakes-interoccurrence times-can be described by the superposition of the Weibull distribution and the log-Weibull distribution. In particular, the distribution of large earthquakes obeys the Weibull distribution with the exponent α₁<1, indicating the fact that the sequence of large earthquakes is not a Poisson process. It is found that the ratio of the Weibull distribution to the probability distribution of the interoccurrence time gradually increases with increase in the threshold of magnitude. Our results infer that Weibull statistics and log-Weibull statistics coexist in the interoccurrence time statistics, and that the change of the distribution is considered as the change of the dominant distribution. In this case, the dominant distribution changes from the log-Weibull distribution to the Weibull distribution, allowing us to reinforce the view that the interoccurrence time exhibits the transition from the Weibull regime to the log-Weibull regime.     

On the probability distribution of stock returns in the Mike-Farmer model

Recently, Mike and Farmer have constructed a very powerful and realistic behavioral model to mimick the dynamic process of stock price formation based on the empirical regularities of order placement and cancelation in a purely order-driven market, which can successfully reproduce the whole distribution of returns, not only the well-known power-law tails, together with several other important stylized facts. There are three key ingredients in the Mike-Farmer (MF) model: the long memory of order signs characterized by the Hurst index H_s, the distribution of relative order prices x in reference to the same best price described by a Student distribution (or Tsallis’ q-Gaussian), and the dynamics of order cancelation. They showed that different values of the Hurst index H_s and the freedom degree α_x of the Student distribution can always produce power-law tails in the return distribution f_r(r) with different tail exponent α_r. In this paper, we study the origin of the power-law tails of the return distribution f_r(r) in the MF model, based on extensive simulations with different combinations of the left part L(x) for x < 0 and the right part R(x) for x > 0 of f_x(x). We find that power-law tails appear only when L(x) has a power-law tail, no matter R(x) has a power-law tail or not. In addition, we find that the distributions of returns in the MF model at different timescales can be well modeled by the Student distributions, whose tail exponents are close to the well-known cubic law and increase with the timescale.     

Changes of firm size distribution: The case of Korea

In this paper, the distribution and inequality of firm sizes is evaluated for the Korean firms listed on the stock markets. Using the amount of sales, total assets, capital, and the number of employees, respectively, as a proxy for firm sizes, we find that the upper tail of the Korean firm size distribution can be described by power-law distributions rather than lognormal distributions. Then, we estimate the Zipf parameters of the firm sizes and assess the changes in the magnitude of the exponents. The results show that the calculated Zipf exponents over time increased prior to the financial crisis, but decreased after the crisis. This pattern implies that the degree of inequality in Korean firm sizes had severely deepened prior to the crisis, but lessened after the crisis. Overall, the distribution of Korean firm sizes changes over time, and Zipf’s law is not universal but does hold as a special case.     

Luminescence decays with underlying distributions: General properties and analysis with mathematical functions

In this work, an analysis of the general properties of the luminescence decay law is carried out. The conditions that a luminescence decay law must satisfy in order to correspond to a probability density function of rate constants are established. From an analysis of the general form of the luminescence decay law, it is concluded that the decay must be either exponential or sub-exponential for all times, in order to be represented by a distribution of rate constants H(k). Sub-exponentiality is nevertheless not a sufficient condition. Only decays that are completely monotonic have a probability density function of rate constants. The construction of the decay function from cumulant and moment expansions is studied, as well as the corresponding calculation of H(k) from a cumulant expansion. The asymptotic behavior of the decay laws is considered in detail, and the relation between this behavior and the form of H(k) for small k is explored. Several generalizations of the exponential decay function, namely the Kohlrausch, Becquerel, Mittag-Leffler and Heaviside decay functions, as well as the Weibull and truncated Gaussian rate constant distributions are discussed.     

Statistical properties of daily ensemble variables in the Chinese stock markets

We study dynamical behavior of the Chinese stock markets by investigating the statistical properties of daily ensemble return and variety defined, respectively, as the mean and the standard deviation of the ensemble daily price return of a portfolio of stocks traded in China's stock markets on a given day. The distribution of the daily ensemble return has an exponential form in the center and power-law tails, while the variety distribution is lognormal in the bulk followed by a power-law tail for large variety. Based on detrended fluctuation analysis, R/S analysis and modified R/S analysis, we find evidence of long memory in the ensemble return and strong evidence of long memory in the evolution of variety.     

Runaway events dominate the heavy tail of citation distributions

Statistical distributions with heavy tails are ubiquitous in natural and social phenomena. Since the entries in heavy tail have unproportional significance, the knowledge of its exact shape is very important. Citations of scientific papers form one of the best-known heavy tail distributions. Even in this case there is a considerable debate whether citation distribution follows the log-normal or power-law fit. The goal of our study is to solve this debate by measuring citation distribution for a very large and homogeneous data. We measured citation distribution for 418, 438 Physics papers published in 1980–1989 and cited by 2008. While the log-normal fit deviates too strong from the data, the discrete power-law function with the exponent γ = 3.15 does better and fits 99.955% of the data. However, the extreme tail of the distribution deviates upward even from the power-law fit and exhibits a dramatic “runaway” behavior. The onset of the runaway regime is revealed macroscopically as the paper garners 1000-1500 citations, however the microscopic measurements of autocorrelation in citation rates are able to predict this behavior in advance.     

Predicting the long tail of book sales: Unearthing the power-law exponent

The concept of the long tail has recently been used to explain the phenomenon in e-commerce where the total volume of sales of the items in the tail is comparable to that of the most popular items. In the case of online book sales, the proportion of tail sales has been estimated using regression techniques on the assumption that the data obeys a power-law distribution. Here we propose a different technique for estimation based on a generative model of book sales that results in an asymptotic power-law distribution of sales, but which does not suffer from the problems related to power-law regression techniques. We show that the proportion of tail sales predicted is very sensitive to the estimated power-law exponent. In particular, if we assume that the power-law exponent of the cumulative distribution is closer to 1.1 rather than to 1.2 (estimates published in 2003, calculated using regression by two groups of researchers), then our computations suggest that the tail sales of Amazon.com, rather than being 40% as estimated by Brynjolfsson, Hu and Smith in 2003, are actually closer to 20%, the proportion estimated by its CEO.     

Understanding the Nature of the Long-Range Memory Phenomenon in Socioeconomic Systems

In the face of the upcoming 30th anniversary of econophysics, we review our contributions and other related works on the modeling of the long-range memory phenomenon in physical, economic, and other social complex systems. Our group has shown that the long-range memory phenomenon can be reproduced using various Markov processes, such as point processes, stochastic differential equations, and agent-based models—reproduced well enough to match other statistical properties of the financial markets, such as return and trading activity distributions and first-passage time distributions. Research has lead us to question whether the observed long-range memory is a result of the actual long-range memory process or just a consequence of the non-linearity of Markov processes. As our most recent result, we discuss the long-range memory of the order flow data in the financial markets and other social systems from the perspective of the fractional Lèvy stable motion. We test widely used long-range memory estimators on discrete fractional Lèvy stable motion represented by the auto-regressive fractionally integrated moving average (ARFIMA) sample series. Our newly obtained results seem to indicate that new estimators of self-similarity and long-range memory for analyzing systems with non-Gaussian distributions have to be developed.     

Financial Return Distributions: Past,Present, and COVID-19

We analyze the price return distributions of currency exchange rates, cryptocurrencies, and contracts for differences (CFDs) representing stock indices, stock shares, and commodities. Based on recent data from the years 2017–2020, we model tails of the return distributions at different time scales by using power-law, stretched exponential, and q-Gaussian functions. We focus on the fitted function parameters and how they change over the years by comparing our results with those from earlier studies and find that, on the time horizons of up to a few minutes, the so-called “inverse-cubic power-law” still constitutes an appropriate global reference. However, we no longer observe the hypothesized universal constant acceleration of the market time flow that was manifested before in an ever faster convergence of empirical return distributions towards the normal distribution. Our results do not exclude such a scenario but, rather, suggest that some other short-term processes related to a current market situation alter market dynamics and may mask this scenario. Real market dynamics is associated with a continuous alternation of different regimes with different statistical properties. An example is the COVID-19 pandemic outburst, which had an enormous yet short-time impact on financial markets. We also point out that two factors—speed of the market time flow and the asset cross-correlation magnitude—while related (the larger the speed, the larger the cross-correlations on a given time scale), act in opposite directions with regard to the return distribution tails, which can affect the expected distribution convergence to the normal distribution.     

Strength statistics and the distribution of earthquake interevent times

The Weibull distribution is often used to model the earthquake interevent times distribution (ITD). We propose a link between the earthquake ITD on single faults with the Earth’s crustal shear strength distribution by means of a phenomenological stick–slip model. For single faults or fault systems with homogeneous strength statistics and power-law stress accumulation we obtain the Weibull ITD. We prove that the moduli of the interevent times and crustal shear strength are linearly related, while the time scale is an algebraic function of the scale of crustal shear strength. We also show that logarithmic stress accumulation leads to the log-Weibull ITD. We investigate deviations of the ITD tails from the Weibull model due to sampling bias, magnitude cutoff thresholds, and non-homogeneous strength parameters. Assuming the Gutenberg–Richter law and independence of the Weibull modulus on the magnitude threshold, we deduce that the interevent time scale drops exponentially with the magnitude threshold. We demonstrate that a microearthquake sequence from the island of Crete and a seismic sequence from Southern California conform reasonably well to the Weibull model.