Detrended cross-correlation analysis: a new method for analyzing two nonstationary time series

Investigation on law and economics of listed companies' financing preference based on complex network theory

In this paper, complex network theory is used to make time-series analysis of key indicators of governance structure and financing data. We analyze scientific listed companies' governance data from 2010 to 2014 and divide them into groups in accordance with the similarity they share. Then we select sample companies to analyze their financing data and explore the influence of governance structure on financing decision and the financing preference they display. This paper reviews relevant laws and regulations of financing from the perspective of law and economics, then proposes reasonable suggestions to consummate the law for the purpose of regulating listed companies' financing. The research provides a reference for making qualitative analysis on companies' financing.

Evolution of scaling behaviors embedded in sentence series from A Story of the Stone

The novel entitled A Story of the Stone provides us precise details of life and social structure of the 18th century China. Its writing lasted a long duration of about 10 years, in which the author’s habit may change significantly. It had been published anonymously up to the beginning of the 20th century, which left a mystery of the author’s attribution. In the present work we focus our attention on scaling behavior embedded in the sentence series from this novel, hope to find how the ideas are organized from single sentences to the whole text. Especially we are interested in the evolution of scale invariance to monitor the changes of the author’s language habit and to find some clues on the author’s attribution. The sentence series are separated into a total of 69 non-overlapping segments with a length of 500 sentences each. The correlation dependent balanced estimation of diffusion entropy (cBEDE) is employed to evaluate the scaling behaviors embedded in the short segments. It is found that the total, the part attributed currently to Xueqin Cao (X-part), and the other part attributed to E Gao (E-part), display scale invariance in a large scale up to 10³ sentences, while their scaling exponents are almost identical. All the segments behave scale invariant in considerable wide scales, most of which reach one third of the length. In the curve of scaling exponent versus segment number, the X-part has rich patterns with averagely larger values, while the E-part has a U-shape with a significant low bottom. This finding is a new clue to support the attribution of the E-part to E Gao.

A generalization of random matrix theory and its application to statistical physics

To study the statistical structure of crosscorrelations in empirical data, we generalize random matrix theory and propose a new method of cross-correlation analysis, known as autoregressive random matrix theory (ARRMT). ARRMT takes into account the influence of auto-correlations in the study of cross-correlations in multiple time series. We first analytically and numerically determine how auto-correlations affect the eigenvalue distribution of the correlation matrix. Then we introduce ARRMT with a detailed procedure of how to implement the method. Finally, we illustrate the method using two examples taken from inflation rates for air pressure data for 95 US cities.

Dynamic Portfolio Strategy Using Clustering Approach

The problem of portfolio optimization is one of the most important issues in asset management. We here propose a new dynamic portfolio strategy based on the time-varying structures of MST networks in Chinese stock markets, where the market condition is further considered when using the optimal portfolios for investment. A portfolio strategy comprises two stages: First, select the portfolios by choosing central and peripheral stocks in the selection horizon using five topological parameters, namely degree, betweenness centrality, distance on degree criterion, distance on correlation criterion and distance on distance criterion. Second, use the portfolios for investment in the investment horizon. The optimal portfolio is chosen by comparing central and peripheral portfolios under different combinations of market conditions in the selection and investment horizons. Market conditions in our paper are identified by the ratios of the number of trading days with rising index to the total number of trading days, or the sum of the amplitudes of the trading days with rising index to the sum of the amplitudes of the total trading days. We find that central portfolios outperform peripheral portfolios when the market is under a drawup condition, or when the market is stable or drawup in the selection horizon and is under a stable condition in the investment horizon. We also find that peripheral portfolios gain more than central portfolios when the market is stable in the selection horizon and is drawdown in the investment horizon. Empirical tests are carried out based on the optimal portfolio strategy. Among all possible optimal portfolio strategies based on different parameters to select portfolios and different criteria to identify market conditions, 65% of our optimal portfolio strategies outperform the random strategy for the Shanghai A-Share market while the proportion is 70% for the Shenzhen A-Share market.

Multifractality and Network Analysis of Phase Transition

Many models and real complex systems possess critical thresholds at which the systems shift dramatically from one sate to another. The discovery of early-warnings in the vicinity of critical points are of great importance to estimate how far the systems are away from the critical states. Multifractal Detrended Fluctuation analysis (MF-DFA) and visibility graph method have been employed to investigate the multifractal and geometrical properties of the magnetization time series of the two-dimensional Ising model. Multifractality of the time series near the critical point has been uncovered from the generalized Hurst exponents and singularity spectrum. Both long-term correlation and broad probability density function are identified to be the sources of multifractality. Heterogeneous nature of the networks constructed from magnetization time series have validated the fractal properties. Evolution of the topological quantities of the visibility graph, along with the variation of multifractality, serve as new early-warnings of phase transition. Those methods and results may provide new insights about the analysis of phase transition problems and can be used as early-warnings for a variety of complex systems.

Multifractals embedded in short time series: An unbiased estimation of probability moment

An exact estimation of probability moments is the base for several essential concepts, such as the multifractals, the Tsallis entropy, and the transfer entropy. By means of approximation theory we propose a new method called factorial-moment-based estimation of probability moments. Theoretical prediction and computational results show that it can provide us an unbiased estimation of the probability moments of continuous order. Calculations on probability redistribution model verify that it can extract exactly multifractal behaviors from several hundred recordings. Its powerfulness in monitoring evolution of scaling behaviors is exemplified by two empirical cases, i.e., the gait time series for fast, normal, and slow trials of a healthy volunteer, and the closing price series for Shanghai stock market. By using short time series with several hundred lengths, a comparison with the well-established tools displays significant advantages of its performance over the other methods. The factorial-moment-based estimation can evaluate correctly the scaling behaviors in a scale range about three generations wider than the multifractal detrended fluctuation analysis and the basic estimation. The estimation of partition function given by the wavelet transform modulus maxima has unacceptable fluctuations. Besides the scaling invariance focused in the present paper, the proposed factorial moment of continuous order can find its various uses, such as finding nonextensive behaviors of a complex system and reconstructing the causality relationship network between elements of a complex system.

Quantifying distinct associations on different temporal scales: comparison of DCCA and Pearson methods

Cross-correlation between pairs of variables takes multi-time scale characteristic, and it can be totally different on different time scales (changing from positive correlation to negative one), e.g., the associations between mean air temperature and relative humidity over regions to the east of Taihang mountain in China. Therefore, how to correctly unveil these correlations on different time scales is really of great importance since we actually do not know if the correlation varies with scales in advance. Here, we compare two methods, i.e. Detrended Cross-Correlation Analysis (DCCA for short) and Pearson correlation, in quantifying scale-dependent correlations directly to raw observed records and artificially generated sequences with known cross-correlation features. Studies show that 1) DCCA related methods can indeed quantify scale-dependent correlations, but not Pearson method; 2) the correlation features from DCCA related methods are robust to contaminated noises, however, the results from Pearson method are sensitive to noise; 3) the scale-dependent correlation results from DCCA related methods are robust to the amplitude ratio between slow and fast components, while Pearson method may be sensitive to the amplitude ratio. All these features indicate that DCCA related methods take some advantages in correctly quantifying scale-dependent correlations, which results from different physical processes.

Effects of Tactile Sensitivity on Structural Variability of Digit Forces during Stable Precision Grip

This study investigated the effects of fingertip tactile sensitivity on the structural variability of thumb and index finger forces during stable precision grip. Thirty right-handed healthy subjects participated in the experiment. Transient perturbation of tactile afferents was achieved by wrapping up the distal pads of the thumb or index finger with transparent polyethylene films. The time-dependent structure of each digit force and the variability of interdigit force correlation were examined by detrended fluctuation analysis (DFA) and detrended cross-correlation analysis (DCCA), respectively. Results showed that the tactile sensitivity affected α_DFA of the vertical shear force Fx (F_3,239 = 6.814, p < 0.001) and α_DCCA of Fx (χ² = 16.440, p < 0.001). No significant difference was observed in α_DFA or α_DCCA of the normal forces produced by the thumb or index finger. These results suggested that with blurred tactile sensory inputs the central nervous system might decrease the vertical shear force flexibility and increase the interdigit shear force coupling in order to guarantee a stable grip control of an object against gravity. This study shed light on the feedback and feed-forward strategies involved in digit force control and the role of SA-II afferent fibers in regulation of vertical shear force variability for precision grip.

Multiscale limited penetrable horizontal visibility graph for analyzing nonlinear time series

Visibility graph has established itself as a powerful tool for analyzing time series. We in this paper develop a novel multiscale limited penetrable horizontal visibility graph (MLPHVG). We use nonlinear time series from two typical complex systems, i.e., EEG signals and two-phase flow signals, to demonstrate the effectiveness of our method. Combining MLPHVG and support vector machine, we detect epileptic seizures from the EEG signals recorded from healthy subjects and epilepsy patients and the classification accuracy is 100%. In addition, we derive MLPHVGs from oil-water two-phase flow signals and find that the average clustering coefficient at different scales allows faithfully identifying and characterizing three typical oil-water flow patterns. These findings render our MLPHVG method particularly useful for analyzing nonlinear time series from the perspective of multiscale network analysis.

Detrended cross-correlation coefficient: Application to predict apoptosis protein subcellular localization

Apoptosis, or programed cell death, plays a central role in the development and homeostasis of an organism. Obtaining information on subcellular location of apoptosis proteins is very helpful for understanding the apoptosis mechanism. The prediction of subcellular localization of an apoptosis protein is still a challenging task, and existing methods mainly based on protein primary sequences. In this paper, we introduce a new position-specific scoring matrix (PSSM)-based method by using detrended cross-correlation (DCCA) coefficient of non-overlapping windows. Then a 190-dimensional (190D) feature vector is constructed on two widely used datasets: CL317 and ZD98, and support vector machine is adopted as classifier. To evaluate the proposed method, objective and rigorous jackknife cross-validation tests are performed on the two datasets. The results show that our approach offers a novel and reliable PSSM-based tool for prediction of apoptosis protein subcellular localization.

Patterns for the waiting time in the context of discrete-time stochastic processes

The aim of this study is to extend the scope and applicability of the level-crossing method to discrete-time stochastic processes and generalize it to enable us to study multiple discrete-time stochastic processes. In previous versions of the level-crossing method, problems with it correspond to the fact that this method had been developed for analyzing a continuous-time process or at most a multiple continuous-time process in an individual manner. However, since all empirical processes are discrete in time, the already-established level-crossing method may not prove adequate for studying empirical processes. Beyond this, due to the fact that most empirical processes are coupled; their individual study could lead to vague results. To achieve the objectives of this study, we first find an analytical expression for the average frequency of crossing a level in a discrete-time process, giving the measure of the time experienced for two consecutive crossings named as the "waiting time." We then introduce the generalized level-crossing method by which the consideration of coupling between the components of a multiple process becomes possible. Finally, we provide an analytic solution when the components of a multiple stochastic process are independent Gaussian white noises. The comparison of the results obtained for coupled and uncoupled processes measures the strength and efficiency of the coupling, justifying our model and analysis. The advantage of the proposed method is its sensitivity to the slightest coupling and shortest correlation length.

Time and frequency domain characteristics of detrending-operation-based scaling analysis: Exact DFA and DMA frequency responses

We develop a general framework to study the time and frequency domain characteristics of detrending-operation-based scaling analysis methods, such as detrended fluctuation analysis (DFA) and detrending moving average (DMA) analysis. In this framework, using either the time or frequency domain approach, the frequency responses of detrending operations are calculated analytically. Although the frequency domain approach based on conventional linear analysis techniques is only applicable to linear detrending operations, the time domain approach presented here is applicable to both linear and nonlinear detrending operations. Furthermore, using the relationship between the time and frequency domain representations of the frequency responses, the frequency domain characteristics of nonlinear detrending operations can be obtained. Based on the calculated frequency responses, it is possible to establish a direct connection between the root-mean-square deviation of the detrending-operation-based scaling analysis and the power spectrum for linear stochastic processes. Here, by applying our methods to DFA and DMA, including higher-order cases, exact frequency responses are calculated. In addition, we analytically investigate the cutoff frequencies of DFA and DMA detrending operations and show that these frequencies are not optimally adjusted to coincide with the corresponding time scale.

Market Confidence Predicts Stock Price: Beyond Supply and Demand

Stock price prediction is an important and challenging problem in stock market analysis. Existing prediction methods either exploit autocorrelation of stock price and its correlation with the supply and demand of stock, or explore predictive indictors exogenous to stock market. In this paper, using transaction record of stocks with identifier of traders, we introduce an index to characterize market confidence, i.e., the ratio of the number of traders who is active in two successive trading days to the number of active traders in a certain trading day. Strong Granger causality is found between the index of market confidence and stock price. We further predict stock price by incorporating the index of market confidence into a neural network based on time series of stock price. Experimental results on 50 stocks in two Chinese Stock Exchanges demonstrate that the accuracy of stock price prediction is significantly improved by the inclusion of the market confidence index. This study sheds light on using cross-day trading behavior to characterize market confidence and to predict stock price.

Suppressing disease spreading by using information diffusion on multiplex networks

Although there is always an interplay between the dynamics of information diffusion and disease spreading, the empirical research on the systemic coevolution mechanisms connecting these two spreading dynamics is still lacking. Here we investigate the coevolution mechanisms and dynamics between information and disease spreading by utilizing real data and a proposed spreading model on multiplex network. Our empirical analysis finds asymmetrical interactions between the information and disease spreading dynamics. Our results obtained from both the theoretical framework and extensive stochastic numerical simulations suggest that an information outbreak can be triggered in a communication network by its own spreading dynamics or by a disease outbreak on a contact network, but that the disease threshold is not affected by information spreading. Our key finding is that there is an optimal information transmission rate that markedly suppresses the disease spreading. We find that the time evolution of the dynamics in the proposed model qualitatively agrees with the real-world spreading processes at the optimal information transmission rate.

A novel way to detect correlations on multi-time scales, with temporal evolution and for multi-variables

In this paper, two new methods, Temporal evolution of Detrended Cross-Correlation Analysis (TDCCA) and Temporal evolution of Detrended Partial-Cross-Correlation Analysis (TDPCCA), are proposed by generalizing DCCA and DPCCA. Applying TDCCA/TDPCCA, it is possible to study correlations on multi-time scales and over different periods. To illustrate their properties, we used two climatological examples: i) Global Sea Level (GSL) versus North Atlantic Oscillation (NAO); and ii) Summer Rainfall over Yangtze River (SRYR) versus previous winter Pacific Decadal Oscillation (PDO). We find significant correlations between GSL and NAO on time scales of 60 to 140 years, but the correlations are non-significant between 1865–1875. As for SRYR and PDO, significant correlations are found on time scales of 30 to 35 years, but the correlations are more pronounced during the recent 30 years. By combining TDCCA/TDPCCA and DCCA/DPCCA, we proposed a new correlation-detection system, which compared to traditional methods, can objectively show how two time series are related (on which time scale, during which time period). These are important not only for diagnosis of complex system, but also for better designs of prediction models. Therefore, the new methods offer new opportunities for applications in natural sciences, such as ecology, economy, sociology and other research fields.

Robust Statistical Detection of Power-Law Cross-Correlation

We show that widely used approaches in statistical physics incorrectly indicate the existence of power-law cross-correlations between financial stock market fluctuations measured over several years and the neuronal activity of the human brain lasting for only a few minutes. While such cross-correlations are nonsensical, no current methodology allows them to be reliably discarded, leaving researchers at greater risk when the spurious nature of cross-correlations is not clear from the unrelated origin of the time series and rather requires careful statistical estimation. Here we propose a theory and method (PLCC-test) which allows us to rigorously and robustly test for power-law cross-correlations, correctly detecting genuine and discarding spurious cross-correlations, thus establishing meaningful relationships between processes in complex physical systems. Our method reveals for the first time the presence of power-law cross-correlations between amplitudes of the alpha and beta frequency ranges of the human electroencephalogram.

Detrending moving average algorithm: Frequency response and scaling performances

The Detrending Moving Average (DMA) algorithm has been widely used in its several variants for characterizing long-range correlations of random signals and sets (one-dimensional sequences or high-dimensional arrays) over either time or space. In this paper, mainly based on analytical arguments, the scaling performances of the centered DMA, including higher-order ones, are investigated by means of a continuous time approximation and a frequency response approach. Our results are also confirmed by numerical tests. The study is carried out for higher-order DMA operating with moving average polynomials of different degree. In particular, detrending power degree, frequency response, asymptotic scaling, upper limit of the detectable scaling exponent, and finite scale range behavior will be discussed.

A novel coefficient for detecting and quantifying asymmetry of California electricity market based on asymmetric detrended cross-correlation analysis

In order to detect and quantify asymmetry of two time series, a novel cross-correlation coefficient is proposed based on recent asymmetric detrended cross-correlation analysis (A-DXA), which we called A-DXA coefficient. The A-DXA coefficient, as an important extension of DXA coefficient ρDXA, contains two directional asymmetric cross-correlated indexes, describing upwards and downwards asymmetric cross-correlations, respectively. By using the information of directional covariance function of two time series and directional variance function of each series itself instead of power-law between the covariance function and time scale, the proposed A-DXA coefficient can well detect asymmetry between the two series no matter whether the cross-correlation is significant or not. By means of the proposed A-DXA coefficient conducted over the asymmetry for California electricity market, we found that the asymmetry between the prices and loads is not significant for daily average data in 1999 yr market (before electricity crisis) but extremely significant for those in 2000 yr market (during the crisis). To further uncover the difference of asymmetry between the years 1999 and 2000, a modified H statistic (MH) and ΔMH statistic are proposed. One of the present contributions is that the high MH values calculated for hourly data exist in majority months in 2000 market. Another important conclusion is that the cross-correlation with downwards dominates over the whole 1999 yr in contrast to the cross-correlation with upwards dominates over the 2000 yr.

"Intrinsic" correlations and their temporal evolutions between winter-time PNA/EPW and winter drought in the west United States

In this study, relations between winter-time Pacific-Northern America pattern (PNA)/East Pacific wave-train (EPW) and winter-time drought in the west United States over the period of 1951–2010 are analyzed. Considering traditional Pearson’s Correlation Coefficient can be influenced by non-stationarity and nonlinearity, a recently proposed method, Detrended Partial-Cross-Correlation Analysis (DPCCA) is applied. With DPCCA, we analyzed the “intrinsic” correlations between PNA/EPW and the winter drought with possible effects of ENSO and PDO removed. We found, i) significant negative correlations between PNA/EPW and drought on time scales of 5–6 years after removing the effects of ENSO, ii) and significant negative correlations between PNA/EPW and drought on time scales of 15–25 years after removing the effects of PDO. By further studying the temporal evolutions of the “intrinsic” correlations, we found on time scales of 5–6 years, the “intrinsic” correlations between PNA/EPW and drought can vary severely with time, but for most time, the correlations are negative. While on interdecadal (15–25 years) time scales, after the effects of PDO removed, unlike the relations between PNA and drought, the “intrinsic” correlations between EPW and drought takes nearly homogeneous-sign over the whole period, indicating a better model can be designed by using EPW.

Influence of Sub-Daily Variation on Multi-Fractal Detrended Fluctuation Analysis of Wind Speed Time Series

Using multi-fractal detrended fluctuation analysis (MF-DFA), the scaling features of wind speed time series (WSTS) could be explored. In this paper, we discuss the influence of sub-daily variation, which is a natural feature of wind, in MF-DFA of WSTS. First, the choice of the lower bound of the segment length, a significant parameter of MF-DFA, was studied. The results of expanding the lower bound into sub-daily scope shows that an abrupt declination and discrepancy of scaling exponents is caused by the inability to keep the whole diel process of wind in one single segment. Additionally, the specific value, which is effected by the sub-daily feature of local meteo-climatic, might be different. Second, the intra-day temporal order of wind was shuffled to determine the impact of diel variation on scaling exponents of MF-DFA. The results illustrate that disregarding diel variation leads to errors in scaling. We propose that during the MF-DFA of WSTS, the segment length should be longer than 1 day and the diel variation of wind should be maintained to avoid abnormal phenomena and discrepancy in scaling exponents.

Weighted multifractal cross-correlation analysis based on Shannon entropy

In this paper, we propose a modification of multifractal cross-correlation analysis based on statistical moments (MFSMXA) method, called weighted MFSMXA method based on Shannon entropy (W-MFSMXA), to investigate cross-correlations and cross-multifractality between time series. Robustness of this method is verified by numerical experiments with both artificial and stock returns series. Results show that the proposed W-MFSMXA method not only keep the multifractal structure unchanged, but contains more significant information of series compared to the previous MFSMXA method. Furthermore, analytic formulas of the binomial multifractal model are generated for W-MFSMXA. Theoretical analysis and finite-size effect test demonstrate that W-MFSMXA slightly outperforms MFSMXA for relatively shorter series. We further generate the scaling exponent ratio to describe the relation of two methods, whose profile is found approximating a centrosymmetric hyperbola. Cross-multifractality is found in returns series but then destroyed after being shuffled as a consequence of the removed long memory in separate series.

Gene Expression Switching of Receptor Subunits in Human Brain Development

Synaptic receptors in the human brain consist of multiple protein subunits, many of which have multiple variants, coded by different genes, and are differentially expressed across brain regions and developmental stages. The brain can tune the electrophysiological properties of synapses to regulate plasticity and information processing by switching from one protein variant to another. Such condition-dependent variant switch during development has been demonstrated in several neurotransmitter systems including NMDA and GABA. Here we systematically detect pairs of receptor-subunit variants that switch during the lifetime of the human brain by analyzing postmortem expression data collected in a population of donors at various ages and brain regions measured using microarray and RNA-seq. To further detect variant pairs that co-vary across subjects, we present a method to quantify age-corrected expression correlation in face of strong temporal trends. This is achieved by computing the correlations in the residual expression beyond a cubic-spline model of the population temporal trend, and can be seen as a nonlinear version of partial correlations. Using these methods, we detect multiple new pairs of context dependent variants. For instance, we find a switch from GLRA2 to GLRA3 that differs from the known switch in the rat. We also detect an early switch from HTR1A to HTR5A whose trends are negatively correlated and find that their age-corrected expression is strongly positively correlated. Finally, we observe that GRIN2B switch to GRIN2A occurs mostly during embryonic development, presumably earlier than observed in rodents. These results provide a systematic map of developmental switching in the neurotransmitter systems of the human brain.

Synapses change their properties during development affecting information processing and learning. Most synaptic receptors consist of several proteins, each having several variants coded by closely related genes. These protein variants are similar in structure, yet often differ slightly in their biophysical attributes. Switching a synapse from using one variant to another provides the brain with a way to fine-tune electrophysiological properties of synapses and has been described in NMDA and GABA receptors. Here we describe a systematic approach to detect pairs of context-dependent variants at a genome-wide scale based on a set of post-mortem expression measurements taken from brains at multiple ages. We take into account both the profile of expression as it changes along life and also the detrended age-corrected correlation among genes. This method characterizes the landscape of developmental switches in brain transcriptome, putting forward new candidates pairs for deeper analysis. The abundance of switching between context-dependent variants through life suggests that it is a major mechanism by which the brain tunes its plasticity and information processing.

Detrended fluctuation analysis made flexible to detect range of cross-correlated fluctuations

The detrended cross-correlation coefficient ρ(DCCA) has recently been proposed to quantify the strength of cross-correlations on different temporal scales in bivariate, nonstationary time series. It is based on the detrended cross-correlation and detrended fluctuation analyses (DCCA and DFA, respectively) and can be viewed as an analog of the Pearson coefficient in the case of the fluctuation analysis. The coefficient ρ(DCCA) works well in many practical situations but by construction its applicability is limited to detection of whether two signals are generally cross-correlated, without the possibility to obtain information on the amplitude of fluctuations that are responsible for those cross-correlations. In order to introduce some related flexibility, here we propose an extension of ρ(DCCA) that exploits the multifractal versions of DFA and DCCA: multifractal detrended fluctuation analysis and multifractal detrended cross-correlation analysis, respectively. The resulting new coefficient ρ(q) not only is able to quantify the strength of correlations but also allows one to identify the range of detrended fluctuation amplitudes that are correlated in two signals under study. We show how the coefficient ρ(q) works in practical situations by applying it to stochastic time series representing processes with long memory: autoregressive and multiplicative ones. Such processes are often used to model signals recorded from complex systems and complex physical phenomena like turbulence, so we are convinced that this new measure can successfully be applied in time-series analysis. In particular, we present an example of such application to highly complex empirical data from financial markets. The present formulation can straightforwardly be extended to multivariate data in terms of the q-dependent counterpart of the correlation matrices and then to the network representation.

Detrended partial cross-correlation analysis of two nonstationary time series influenced by common external forces

When common factors strongly influence two power-law cross-correlated time series recorded in complex natural or social systems, using detrended cross-correlation analysis (DCCA) without considering these common factors will bias the results. We use detrended partial cross-correlation analysis (DPXA) to uncover the intrinsic power-law cross correlations between two simultaneously recorded time series in the presence of nonstationarity after removing the effects of other time series acting as common forces. The DPXA method is a generalization of the detrended cross-correlation analysis that takes into account partial correlation analysis. We demonstrate the method by using bivariate fractional Brownian motions contaminated with a fractional Brownian motion. We find that the DPXA is able to recover the analytical cross Hurst indices, and thus the multiscale DPXA coefficients are a viable alternative to the conventional cross-correlation coefficient. We demonstrate the advantage of the DPXA coefficients over the DCCA coefficients by analyzing contaminated bivariate fractional Brownian motions. We calculate the DPXA coefficients and use them to extract the intrinsic cross correlation between crude oil and gold futures by taking into consideration the impact of the U.S. dollar index. We develop the multifractal DPXA (MF-DPXA) method in order to generalize the DPXA method and investigate multifractal time series. We analyze multifractal binomial measures masked with strong white noises and find that the MF-DPXA method quantifies the hidden multifractal nature while the multifractal DCCA method fails.

The effect of the underlying distribution in Hurst exponent estimation

In this paper, a heavy-tailed distribution approach is considered in order to explore the behavior of actual financial time series. We show that this kind of distribution allows to properly fit the empirical distribution of the stocks from S&P500 index. In addition to that, we explain in detail why the underlying distribution of the random process under study should be taken into account before using its self-similarity exponent as a reliable tool to state whether that financial series displays long-range dependence or not. Finally, we show that, under this model, no stocks from S&P500 index show persistent memory, whereas some of them do present anti-persistent memory and most of them present no memory at all.

A new methodology of spatial cross-correlation analysis

Spatial correlation modeling comprises both spatial autocorrelation and spatial cross-correlation processes. The spatial autocorrelation theory has been well-developed. It is necessary to advance the method of spatial cross-correlation analysis to supplement the autocorrelation analysis. This paper presents a set of models and analytical procedures for spatial cross-correlation analysis. By analogy with Moran's index newly expressed in a spatial quadratic form, a theoretical framework is derived for geographical cross-correlation modeling. First, two sets of spatial cross-correlation coefficients are defined, including a global spatial cross-correlation coefficient and local spatial cross-correlation coefficients. Second, a pair of scatterplots of spatial cross-correlation is proposed, and the plots can be used to visually reveal the causality behind spatial systems. Based on the global cross-correlation coefficient, Pearson's correlation coefficient can be decomposed into two parts: direct correlation (partial correlation) and indirect correlation (spatial cross-correlation). As an example, the methodology is applied to the relationships between China's urbanization and economic development to illustrate how to model spatial cross-correlation phenomena. This study is an introduction to developing the theory of spatial cross-correlation, and future geographical spatial analysis might benefit from these models and indexes.

Generalized permutation entropy analysis based on the two-index entropic form Sq,δ

Permutation entropy (PE) is a novel measure to quantify the complexity of nonlinear time series. In this paper, we propose a generalized permutation entropy ( PEq,δ) based on the recently postulated entropic form, Sq,δ, which was proposed as an unification of the well-known Sq of nonextensive-statistical mechanics and Sδ, a possibly appropriate candidate for the black-hole entropy. We find that PEq,δ with appropriate parameters can amplify minor changes and trends of complexities in comparison to PE. Experiments with this generalized permutation entropy method are performed with both synthetic and stock data showing its power. Results show that PEq,δ is an exponential function of q and the power ( k(δ)) is a constant if δ is determined. Some discussions about k(δ) are provided. Besides, we also find some interesting results about power law.

Coordination of digit force variability during dominant and non-dominant sustained precision pinch

This study examined the effects of handedness on the inter-digit coordination of force variability with and without concurrent visual feedback during sustained precision pinch. Twenty-four right-handed subjects were instructed to pinch an instrumented apparatus with their dominant and non-dominant hands, separately. During the pinch, the subjects were required to maintain a stable force output at 5 N for 1 min. Visual feedback was given for the first 30 s and removed for the second 30 s. Coefficient of variation and detrended fluctuation analysis were employed to examine the amount and structural variability of the thumb and index finger forces. Similarly, correlation coefficient and detrended cross-correlation analysis were applied to quantify the inter-digit correlation of force amount and structural variability. Results showed that, compared to the non-dominant hand, the dominant hand had higher inter-digit difference in the amount of digit force variability. Without visual feedback, the dominant hand exhibited lower digit force structural variability but higher inter-digit force structural correlation than the non-dominant hand. These results implied that the dominant hand would be more independent, less flexible and with lower dynamic degrees of freedom than the non-dominant hand in coordination of the thumb and index finger forces during sustained precision pinch. The effects of handedness on inter-digit force coordination were dependent on sensory condition, which shed light on higher-level sensorimotor mechanisms that may be responsible for the asymmetries in coordination of digit force variability.

Detecting and interpreting distortions in hierarchical organization of complex time series

Hierarchical organization is a cornerstone of complexity and multifractality constitutes its central quantifying concept. For model uniform cascades the corresponding singularity spectra are symmetric while those extracted from empirical data are often asymmetric. Using selected time series representing such diverse phenomena as price changes and intertransaction times in financial markets, sentence length variability in narrative texts, Missouri River discharge, and sunspot number variability as examples, we show that the resulting singularity spectra appear strongly asymmetric, more often left sided but in some cases also right sided. We present a unified view on the origin of such effects and indicate that they may be crucially informative for identifying the composition of the time series. One particularly intriguing case of this latter kind of asymmetry is detected in the daily reported sunspot number variability. This signals that either the commonly used famous Wolf formula distorts the real dynamics in expressing the largest sunspot numbers or, if not, that their dynamics is governed by a somewhat different mechanism.

Detrended fluctuation analysis as a regression framework: estimating dependence at different scales

We propose a framework combining detrended fluctuation analysis with standard regression methodology. The method is built on detrended variances and covariances and it is designed to estimate regression parameters at different scales and under potential nonstationarity and power-law correlations. The former feature allows for distinguishing between effects for a pair of variables from different temporal perspectives. The latter ones make the method a significant improvement over the standard least squares estimation. Theoretical claims are supported by Monte Carlo simulations. The method is then applied on selected examples from physics, finance, environmental science, and epidemiology. For most of the studied cases, the relationship between variables of interest varies strongly across scales.

Detrended partial-cross-correlation analysis: a new method for analyzing correlations in complex system

In this paper, a new method, detrended partial-cross-correlation analysis (DPCCA), is proposed. Based on detrended cross-correlation analysis (DCCA), this method is improved by including partial-correlation technique, which can be applied to quantify the relations of two non-stationary signals (with influences of other signals removed) on different time scales. We illustrate the advantages of this method by performing two numerical tests. Test I shows the advantages of DPCCA in handling non-stationary signals, while Test II reveals the "intrinsic" relations between two considered time series with potential influences of other unconsidered signals removed. To further show the utility of DPCCA in natural complex systems, we provide new evidence on the winter-time Pacific Decadal Oscillation (PDO) and the winter-time Nino3 Sea Surface Temperature Anomaly (Nino3-SSTA) affecting the Summer Rainfall over the middle-lower reaches of the Yangtze River (SRYR). By applying DPCCA, better significant correlations between SRYR and Nino3-SSTA on time scales of 6 ~ 8 years are found over the period 1951 ~ 2012, while significant correlations between SRYR and PDO on time scales of 35 years arise. With these physically explainable results, we have confidence that DPCCA is an useful method in addressing complex systems.

Existence of anticorrelations for local field potentials recorded from mice reared in standard condition and environmental enrichment

In the present paper, we analyze local field potentials (LFPs) recorded from the secondary motor cortex (M2) and primary visual cortex (V1) of freely moving mice reared in environmental enrichment (EE) and standard condition (SC). We focus on the scaling properties of the signals by using an integrated approach combining three different techniques: the Higuchi method, detrended fluctuation analysis, and power spectrum. Each technique provides direct or indirect estimations of the Hurst exponent H and this prevents spurious identification of scaling properties in time-series analysis. It is well known that the power spectrum of an LFP signal scales as 1/f(β) with β>0. Our results indicate the existence of a particular power spectrum scaling law 1/f(β) with β<0 for low frequencies (f<4 Hz) for both SC and EE rearing conditions. This type of scaling behavior is associated to the presence of anticorrelation in the corresponding LFP signals. Moreover, since EE is an experimental protocol based on the enhancement of sensorimotor stimulation, we study the possible effects of EE on the scaling properties of secondary motor cortex (M2) and primary visual cortex (V1). Notably, the difference between Hurst's exponents in EE and SC for individual cortical regions (M2) and (V1) is not statistically significant. On the other hand, using the detrended cross-correlation coefficient, we find that EE significantly reduces the functional coupling between secondary motor cortex (M2) and visual cortex (V1).

Evaluation of scaling invariance embedded in short time series

Scaling invariance of time series has been making great contributions in diverse research fields. But how to evaluate scaling exponent from a real-world series is still an open problem. Finite length of time series may induce unacceptable fluctuation and bias to statistical quantities and consequent invalidation of currently used standard methods. In this paper a new concept called correlation-dependent balanced estimation of diffusion entropy is developed to evaluate scale-invariance in very short time series with length ~10(2). Calculations with specified Hurst exponent values of 0.2,0.3,...,0.9 show that by using the standard central moving average de-trending procedure this method can evaluate the scaling exponents for short time series with ignorable bias (≤0.03) and sharp confidential interval (standard deviation ≤0.05). Considering the stride series from ten volunteers along an approximate oval path of a specified length, we observe that though the averages and deviations of scaling exponents are close, their evolutionary behaviors display rich patterns. It has potential use in analyzing physiological signals, detecting early warning signals, and so on. As an emphasis, the our core contribution is that by means of the proposed method one can estimate precisely shannon entropy from limited records.

Chaotic behavior of light-assisted physical aging in arsenoselenide glasses

The theory of strange attractors is shown to be adequately applicable for analyzing the kinetics of light-assisted physical aging revealed in structural relaxation of Se-rich As-Se glasses below glass transition. Kinetics of enthalpy losses is used to determine the phase space reconstruction parameters. Observed chaotic behaviour (involving chaos and fractal consideration such as detrended fluctuation analysis, attractor identification using phase space representation, delay coordinates, mutual information, false nearest neighbours, etc.) reconstructed via the TISEAN program package is treated within a microstructure model describing multistage aging behaviour in arsenoselenide glasses. This simulation testifies that photoexposure acts as an initiating factor only at the beginning stage of physical aging, thus facilitating further atomic shrinkage of a glassy backbone.

Spectrum-based estimators of the bivariate Hurst exponent

We discuss two alternate spectrum-based estimators of the bivariate Hurst exponent in the power-law cross-correlations setting, the cross-periodogram and local X-Whittle estimators, as generalizations of their univariate counterparts. As the spectrum-based estimators are dependent on a part of the spectrum taken into consideration during estimation, a simulation study showing performance of the estimators under varying bandwidth parameter as well as correlation between processes and their specification is provided as well. These estimators are less biased than the already existent averaged periodogram estimator, which, however, has slightly lower variance. The spectrum-based estimators can serve as a good complement to the popular time domain estimators.

Spectral domain of large nonsymmetric correlated Wishart matrices

We study complex eigenvalues of the Wishart model for nonsymmetric correlation matrices. The model is defined for two statistically equivalent but different Gaussian real matrices, as C=AB(t)/T, where B(t) is the transpose of B and both matrices A and B are of dimensions N×T. If A and B are uncorrelated, or equivalently if C vanishes on average, it is known that at large matrix dimension the domain of the eigenvalues of C is a circle centered-at-origin and the eigenvalue density depends only on the radial distances. We consider actual correlation in A and B and derive a result for the contour describing the domain of the bulk of the eigenvalues of C in the limit of large N and T where the ratio N/T is finite. In particular, we show that the eigenvalue domain is sensitive to the correlations. For example, when C is diagonal on average with the same element c≠0, the contour is no longer a circle centered at origin but a shifted ellipse. In this case we explicitly derive a result for the spectral density which again depends only on the radial distances. For more general cases, we show that the contour depends on the symmetric and antisymmetric parts of the correlation matrix resulting from the ensemble-averaged C. If the correlation matrix is normal, then the contour depends only on its spectrum. We also provide numerics to justify our analytics.

Asymmetric multiscale detrended cross-correlation analysis of financial time series

We propose the asymmetric multiscale detrended cross-correlation analysis (MS-ADCCA) method and apply MS-ADCCA method to explore the existence of asymmetric cross-correlation for daily price returns in US and Chinese stock markets and to assess the properties of these asymmetric cross-correlations. The results all show the existences of asymmetric cross-correlations, while small asymmetries at small scales and larger asymmetries at larger scales are also displayed. There is a strong similarity between S&P500 and DJI, and we reveal that the asymmetries depend more on the cross-correlations of S&P500 vs. DJI, S&P500 vs. NQCI, DJI vs. NQCI, and ShangZheng vs. ShenCheng when the market is falling than rising, respectively. By comparing the spectra of S&P500 vs. NQCI and DJI vs. NQCI with uptrends and downtrends, we detect some new characteristics which lead to some new conclusions. Likewise, some new conclusions also can be drawn by the new characteristics displayed through the comparison between the spectra of ShangZheng vs. HSI and ShenCheng vs. HSI. Obviously, we conclude that although the overall spectra are similar and one market has the same effect when it is rising and falling in the study of asymmetric cross-correlations between it and different markets, the cross-correlations and asymmetries on the trends of the different markets are all different. MS-ADCCA method can detect the differences on the asymmetric cross-correlations by different trends of markets. Moreover, the uniqueness of cross-correlation between NQCI and HSI can be detected in the study of the asymmetric cross-correlations, which confirms that HSI is unique in the Chinese stock markets and NQCI is unique in the US stock markets further.

Cross-correlation asymmetries and causal relationships between stock and market risk

We study historical correlations and lead-lag relationships between individual stock risk (volatility of daily stock returns) and market risk (volatility of daily returns of a market-representative portfolio) in the US stock market. We consider the cross-correlation functions averaged over all stocks, using 71 stock prices from the Standard & Poor's 500 index for 1994-2013. We focus on the behavior of the cross-correlations at the times of financial crises with significant jumps of market volatility. The observed historical dynamics showed that the dependence between the risks was almost linear during the US stock market downturn of 2002 and after the US housing bubble in 2007, remaining at that level until 2013. Moreover, the averaged cross-correlation function often had an asymmetric shape with respect to zero lag in the periods of high correlation. We develop the analysis by the application of the linear response formalism to study underlying causal relations. The calculated response functions suggest the presence of characteristic regimes near financial crashes, when the volatility of an individual stock follows the market volatility and vice versa.

Synchronization in human musical rhythms and mutually interacting complex systems

Though the music produced by an ensemble is influenced by multiple factors, including musical genre, musician skill, and individual interpretation, rhythmic synchronization is at the foundation of musical interaction. Here, we study the statistical nature of the mutual interaction between two humans synchronizing rhythms. We find that the interbeat intervals of both laypeople and professional musicians exhibit scale-free (power law) cross-correlations. Surprisingly, the next beat to be played by one person is dependent on the entire history of the other person's interbeat intervals on timescales up to several minutes. To understand this finding, we propose a general stochastic model for mutually interacting complex systems, which suggests a physiologically motivated explanation for the occurrence of scale-free cross-correlations. We show that the observed long-term memory phenomenon in rhythmic synchronization can be imitated by fractal coupling of separately recorded or synthesized audio tracks and thus applied in electronic music. Though this study provides an understanding of fundamental characteristics of timing and synchronization at the interbrain level, the mutually interacting complex systems model may also be applied to study the dynamics of other complex systems where scale-free cross-correlations have been observed, including econophysics, physiological time series, and collective behavior of animal flocks.

Modified multidimensional scaling approach to analyze financial markets

Detrended cross-correlation coefficient (σDCCA) and dynamic time warping (DTW) are introduced as the dissimilarity measures, respectively, while multidimensional scaling (MDS) is employed to translate the dissimilarities between daily price returns of 24 stock markets. We first propose MDS based on σDCCA dissimilarity and MDS based on DTW dissimilarity creatively, while MDS based on Euclidean dissimilarity is also employed to provide a reference for comparisons. We apply these methods in order to further visualize the clustering between stock markets. Moreover, we decide to confront MDS with an alternative visualization method, "Unweighed Average" clustering method, for comparison. The MDS analysis and "Unweighed Average" clustering method are employed based on the same dissimilarity. Through the results, we find that MDS gives us a more intuitive mapping for observing stable or emerging clusters of stock markets with similar behavior, while the MDS analysis based on σDCCA dissimilarity can provide more clear, detailed, and accurate information on the classification of the stock markets than the MDS analysis based on Euclidean dissimilarity. The MDS analysis based on DTW dissimilarity indicates more knowledge about the correlations between stock markets particularly and interestingly. Meanwhile, it reflects more abundant results on the clustering of stock markets and is much more intensive than the MDS analysis based on Euclidean dissimilarity. In addition, the graphs, originated from applying MDS methods based on σDCCA dissimilarity and DTW dissimilarity, may also guide the construction of multivariate econometric models.

Scaling analysis of stock markets

In this paper, we apply the detrended fluctuation analysis (DFA), local scaling detrended fluctuation analysis (LSDFA), and detrended cross-correlation analysis (DCCA) to investigate correlations of several stock markets. DFA method is for the detection of long-range correlations used in time series. LSDFA method is to show more local properties by using local scale exponents. DCCA method is a developed method to quantify the cross-correlation of two non-stationary time series. We report the results of auto-correlation and cross-correlation behaviors in three western countries and three Chinese stock markets in periods 2004-2006 (before the global financial crisis), 2007-2009 (during the global financial crisis), and 2010-2012 (after the global financial crisis) by using DFA, LSDFA, and DCCA method. The findings are that correlations of stocks are influenced by the economic systems of different countries and the financial crisis. The results indicate that there are stronger auto-correlations in Chinese stocks than western stocks in any period and stronger auto-correlations after the global financial crisis for every stock except Shen Cheng; The LSDFA shows more comprehensive and detailed features than traditional DFA method and the integration of China and the world in economy after the global financial crisis; When it turns to cross-correlations, it shows different properties for six stock markets, while for three Chinese stocks, it reaches the weakest cross-correlations during the global financial crisis.

Multiscale multifractal analysis of traffic signals to uncover richer structures

Multifractal detrended fluctuation analysis (MF-DFA) is the most popular method to detect multifractal characteristics of considerable signals such as traffic signals. When fractal properties vary from point to point along the series, it leads to multifractality. In this study, we concentrate not only on the fact that traffic signals have multifractal properties, but also that such properties depend on the time scale in which the multifractality is computed. Via the multiscale multifractal analysis (MMA), traffic signals appear to be far more complex and contain more information which MF-DFA cannot explore by using a fixed time scale. More importantly, we do not have to avoid data sets with crossovers or narrow the investigated time scales, which may lead to biased results. Instead, the Hurst surface provides a spectrum of local scaling exponents at different scale ranges, which helps us to easily position these crossovers. Through comparing Hurst surfaces for signals before and after removing periodical trends, we find periodicities of traffic signals are the main source of the crossovers. Besides, the Hurst surface of the weekday series behaves differently from that of the weekend series. Results also show that multifractality of traffic signals is mainly due to both broad probability density function and correlations. The effects of data loss are also discussed, which suggests that we should carefully handle MMA results when the percentage of data loss is larger than 40%.

Detrended cross-correlation analysis consistently extended to multifractality

We propose an algorithm, multifractal cross-correlation analysis (MFCCA), which constitutes a consistent extension of the detrended cross-correlation analysis and is able to properly identify and quantify subtle characteristics of multifractal cross-correlations between two time series. Our motivation for introducing this algorithm is that the already existing methods, like multifractal extension, have at best serious limitations for most of the signals describing complex natural processes and often indicate multifractal cross-correlations when there are none. The principal component of the present extension is proper incorporation of the sign of fluctuations to their generalized moments. Furthermore, we present a broad analysis of the model fractal stochastic processes as well as of the real-world signals and show that MFCCA is a robust and selective tool at the same time and therefore allows for a reliable quantification of the cross-correlative structure of analyzed processes. In particular, it allows one to identify the boundaries of the multifractal scaling and to analyze a relation between the generalized Hurst exponent and the multifractal scaling parameter λ(q). This relation provides information about the character of potential multifractality in cross-correlations and thus enables a deeper insight into dynamics of the analyzed processes than allowed by any other related method available so far. By using examples of time series from the stock market, we show that financial fluctuations typically cross-correlate multifractally only for relatively large fluctuations, whereas small fluctuations remain mutually independent even at maximum of such cross-correlations. Finally, we indicate possible utility of MFCCA to study effects of the time-lagged cross-correlations.

Trading network predicts stock price

Stock price prediction is an important and challenging problem for studying financial markets. Existing studies are mainly based on the time series of stock price or the operation performance of listed company. In this paper, we propose to predict stock price based on investors' trading behavior. For each stock, we characterize the daily trading relationship among its investors using a trading network. We then classify the nodes of trading network into three roles according to their connectivity pattern. Strong Granger causality is found between stock price and trading relationship indices, i.e., the fraction of trading relationship among nodes with different roles. We further predict stock price by incorporating these trading relationship indices into a neural network based on time series of stock price. Experimental results on 51 stocks in two Chinese Stock Exchanges demonstrate the accuracy of stock price prediction is significantly improved by the inclusion of trading relationship indices.