Analysis of clusters formed by the moving average of a long-range correlated time series

Kullback-Leibler cluster entropy to quantify volatility correlation and risk diversity

The Kullback-Leibler cluster entropy D_{C}[P∥Q] is evaluated for the empirical and model probability distributions P and Q of the clusters formed in the realized volatility time series of five assets (S&P500, NASDAQ, DJIA, DAX, and FTSEMIB). The Kullback-Leibler functional D_{C}[P∥Q] provides complementary perspectives about the stochastic volatility process compared to the Shannon functional S_{C}[P]. While D_{C}[P∥Q] is maximum at the short timescales, S_{C}[P] is maximum at the large timescales leading to complementary optimization criteria tracing back respectively to the maximum and minimum relative entropy evolution principles. The realized volatility is modelled as a time-dependent fractional stochastic process characterized by power-law decaying distributions with positive correlation (H>1/2). As a case study, we build a multiperiod portfolio on diversity indexes derived from the Kullback-Leibler entropy measure of the realized volatility. The portfolio is robust and exhibits better performances over the horizon periods. A comparison with the portfolio built either according to the uniform distribution or in the framework of the Markowitz theory is also reported.

Asymmetric multiscale multifractal analysis (AMMA) of heart rate variability

Objective.The physiological activity of the heart is controlled and modulated mostly by the parasympathetic and sympathetic nervous systems. Heart rate variability (HRV) analysis is therefore used to observe fluctuations that reflect changes in the activity in these two branches. Knowing that acceleration and deceleration patterns in heart rate fluctuations are asymmetrically distributed, the ability to analyze HRV asymmetry was introduced into MMA.Approach. The new method is called asymmetric multiscale multifractal analysis (AMMA) and the analysis involved six groups: 36 healthy persons, 103 cases with aortic valve stenosis, 36 with hypertrophic cardiomyopathy, 32 with atrial fibrillation, 59 patients with coronary artery disease (CAD) and 13 with congestive heart failure.Main results. Analyzing the results obtained for the 6 groups of patients based on the AMMA method, i.e. comparing the Hurst surfaces for heart rate decelerations and accelerations, it was noticed that these surfaces differ significantly. And the differences occur in most groups for large fluctuations (multifractal parameterq > 0). In addition, a similarity was found for all groups for the AMMA Hurst surface for decelerations to the MMA Hurst surface-heart rate decelerations (lengthening of the RR intervals) appears to be the main factor determining the shape of the complete Hurst surface and so the multifractal properties of HRV. The differences between the groups, especially for CAD, hypertrophic cardiomyopathy and aortic valve stenosis, are more visible if the Hurst surfaces are analyzed separately for accelerations and decelerations.Significance. The AMMA results presented here may provide additional input for HRV analysis and create a new paradigm for future medical screening. Note that the HRV analysis using MMA (without distinguishing accelerations from decelerations) gave satisfactory screening statistics in our previous studies.

Prediction of Locomotor Activity by Infrared Motion Detector on Sleep-wake State in Mice

Objective

Behavioral assessments that effectively predict sleep-wake states were tried in animal research. This study aimed to examine the prediction power of an infrared locomotion detector on the sleep-wake states in ICR (Institute Cancer Research) mice. We also explored the influence of the durations and ways of data processing on the prediction power.

Methods

The locomotor activities of seven male mice in home cages were recorded by infrared detectors. Their sleep-wake states were assessed by video analysis. Using the receiver operating characteristic curve analysis, the cut-off score was determined, then the area under the curve (AUC) values of the infrared motion detector that predicted sleep-wake states were calculated. In order to improve the prediction power, the four ways of data processing on the prediction power were performed by Matlab 2013b.

Results

In the initial analysis of raw data, the AUC value was 0.785, but it gradually reached to 0.942 after data summation. The simple data averaging and summation among four different methods showed the maximal AUC value. The 10-minute data summation improved sensitivity (0.889) and specificity (0.901) significantly from the baseline value (sensitivity 0.615; specificity 0.936) (p < 0.001).

Conclusion

This study suggests that the locomotor activity measured by an infrared motion detector might be useful to predict the sleep-wake states in ICR mice. It also revealed that only simple data summation may improve the predictive power. Using daily locomotor activities measured by an infrared motion detector is expected to facilitate animal research related to sleep-wake states.

Long-Range Dependence in Financial Markets: A Moving Average Cluster Entropy Approach

A perspective is taken on the intangible complexity of economic and social systems by investigating the dynamical processes producing, storing and transmitting information in financial time series. An extensive analysis based on the moving average cluster entropy approach has evidenced market and horizon dependence in highest-frequency data of real world financial assets. The behavior is scrutinized by applying the moving average cluster entropy approach to long-range correlated stochastic processes as the Autoregressive Fractionally Integrated Moving Average (ARFIMA) and Fractional Brownian motion (FBM). An extensive set of series is generated with a broad range of values of the Hurst exponent H and of the autoregressive, differencing and moving average parameters p,d,q. A systematic relation between moving average cluster entropy and long-range correlation parameters H, d is observed. This study shows that the characteristic behaviour exhibited by the horizon dependence of the cluster entropy is related to long-range positive correlation in financial markets. Specifically, long range positively correlated ARFIMA processes with differencing parameter d≃0.05, d≃0.15 and d≃0.25 are consistent with moving average cluster entropy results obtained in time series of DJIA, S&P500 and NASDAQ. The findings clearly point to a variability of price returns, consistently with a price dynamics involving multiple temporal scales and, thus, short- and long-run volatility components. An important aspect of the proposed approach is the ability to capture detailed horizon dependence over relatively short horizons (one to twelve months) and thus its relevance to define risk analysis indices.

Multifractal analysis of financial markets: a review

Multifractality is ubiquitously observed in complex natural and socioeconomic systems. Multifractal analysis provides powerful tools to understand the complex nonlinear nature of time series in diverse fields. Inspired by its striking analogy with hydrodynamic turbulence, from which the idea of multifractality originated, multifractal analysis of financial markets has bloomed, forming one of the main directions of econophysics. We review the multifractal analysis methods and multifractal models adopted in or invented for financial time series and their subtle properties, which are applicable to time series in other disciplines. We survey the cumulating evidence for the presence of multifractality in financial time series in different markets and at different time periods and discuss the sources of multifractality. The usefulness of multifractal analysis in quantifying market inefficiency, in supporting risk management and in developing other applications is presented. We finally discuss open problems and further directions of multifractal analysis.

Some comments on Bitcoin market (in)efficiency

In this paper, we explore the (in)efficiency of the continuum Bitcoin-USD market in the period ranging from mid 2010 to early 2019. To deal with, we dynamically analyse the evolution of the self-similarity exponent of Bitcoin-USD daily returns via accurate FD4 approach by a 512 day sliding window with overlapping data. Further, we define the memory indicator by the difference between the self-similarity exponent of Bitcoin-USD series and the self-similarity index of its shuffled series. We also carry out additional analyses via FD4 approach by sliding windows of sizes equal to 64, 128, 256, and 1024 days, and also via FD algorithm for values of q equal to 1 and 2 (and sliding windows equal to 512 days). Moreover, we explored the evolution of the self-similarity exponent of actual S&P500 series via FD4 algorithm by sliding windows of sizes equal to 256 and 512 days. In all the cases, the obtained results were found to be similar to our first analysis. We conclude that the self-similarity exponent of the BTC-USD (resp., S&P500) series stands above 0.5. However, this is not due to the presence of significant memory in the series but to its underlying distribution. In fact, it holds that the self-similarity exponent of BTC-USD (resp., S&P500) series is similar or lower than the self-similarity index of a random series with the same distribution. As such, several periods with significant antipersistent memory in BTC-USD (resp., S&P500) series are distinguished.

Theoretical foundation of detrending methods for fluctuation analysis such as detrended fluctuation analysis and detrending moving average

We present a bottom-up derivation of fluctuation analysis with detrending for the detection of long-range correlations in the presence of additive trends or intrinsic nonstationarities. While the well-known detrended fluctuation analysis (DFA) and detrending moving average (DMA) were introduced ad hoc, we claim basic principles for such methods where DFA and DMA are then shown to be specific realizations. The mean-squared displacement of the summed time series contains the same information about long-range correlations as the autocorrelation function but has much better statistical properties for large time lags. However, the scaling exponent of its estimator on a single time series is affected not only by trends on the data but also by intrinsic nonstationarities. We therefore define the fluctuation function as mean-squared displacement with weighting kernel. We require that its estimator be unbiased and exhibit the correct scaling behavior for the random component of a signal, which is only achieved if the weighting kernel implies detrending. We show how DFA and DMA satisfy these requirements and we extract their kernel weights.

Stability of synchrony against local intermittent fluctuations in tree-like power grids

90% of all Renewable Energy Power in Germany is installed in tree-like distribution grids. Intermittent power fluctuations from such sources introduce new dynamics into the lower grid layers. At the same time, distributed resources will have to contribute to stabilize the grid against these fluctuations in the future. In this paper, we model a system of distributed resources as oscillators on a tree-like, lossy power grid and its ability to withstand desynchronization from localized intermittent renewable infeed. We find a remarkable interplay of the network structure and the position of the node at which the fluctuations are fed in. An important precondition for our findings is the presence of losses in distribution grids. Then, the most network central node splits the network into branches with different influence on network stability. Troublemakers, i.e., nodes at which fluctuations are especially exciting the grid, tend to be downstream branches with high net power outflow. For low coupling strength, we also find branches of nodes vulnerable to fluctuations anywhere in the network. These network regions can be predicted at high confidence using an eigenvector based network measure taking the turbulent nature of perturbations into account. While we focus here on tree-like networks, the observed effects also appear, albeit less pronounced, for weakly meshed grids. On the other hand, the observed effects disappear for lossless power grids often studied in the complex system literature.

Temporal and spatial correlation patterns of air pollutants in Chinese cities

As a huge threat to the public health, China’s air pollution has attracted extensive attention and continues to grow in tandem with the economy. Although the real-time air quality report can be utilized to update our knowledge on air quality, questions about how pollutants evolve across time and how pollutants are spatially correlated still remain a puzzle. In view of this point, we adopt the PMFG network method to analyze the six pollutants’ hourly data in 350 Chinese cities in an attempt to find out how these pollutants are correlated temporally and spatially. In terms of time dimension, the results indicate that, except for O₃, the pollutants have a common feature of the strong intraday patterns of which the daily variations are composed of two contraction periods and two expansion periods. Besides, all the time series of the six pollutants possess strong long-term correlations, and this temporal memory effect helps to explain why smoggy days are always followed by one after another. In terms of space dimension, the correlation structure shows that O₃ is characterized by the highest spatial connections. The PMFGs reveal the relationship between this spatial correlation and provincial administrative divisions by filtering the hierarchical structure in the correlation matrix and refining the cliques as the tinny spatial clusters. Finally, we check the stability of the correlation structure and conclude that, except for PM₁₀ and O₃, the other pollutants have an overall stable correlation, and all pollutants have a slight trend to become more divergent in space. These results not only enhance our understanding of the air pollutants’ evolutionary process, but also shed lights on the application of complex network methods into geographic issues.

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.

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.

Fast algorithm for scaling analysis with higher-order detrending moving average method

Among scaling analysis methods based on the root-mean-square deviation from the estimated trend, it has been demonstrated that centered detrending moving average (DMA) analysis with a simple moving average has good performance when characterizing long-range correlation or fractal scaling behavior. Furthermore, higher-order DMA has also been proposed; it is shown to have better detrending capabilities, removing higher-order polynomial trends than original DMA. However, a straightforward implementation of higher-order DMA requires a very high computational cost, which would prevent practical use of this method. To solve this issue, in this study, we introduce a fast algorithm for higher-order DMA, which consists of two techniques: (1) parallel translation of moving averaging windows by a fixed interval; (2) recurrence formulas for the calculation of summations. Our algorithm can significantly reduce computational cost. Monte Carlo experiments show that the computational time of our algorithm is approximately proportional to the data length, although that of the conventional algorithm is proportional to the square of the data length. The efficiency of our algorithm is also shown by a systematic study of the performance of higher-order DMA, such as the range of detectable scaling exponents and detrending capability for removing polynomial trends. In addition, through the analysis of heart-rate variability time series, we discuss possible applications of higher-order DMA.

Establishing a direct connection between detrended fluctuation analysis and Fourier analysis

To understand methodological features of the detrended fluctuation analysis (DFA) using a higher-order polynomial fitting, we establish the direct connection between DFA and Fourier analysis. Based on an exact calculation of the single-frequency response of the DFA, the following facts are shown analytically: (1) in the analysis of stochastic processes exhibiting a power-law scaling of the power spectral density (PSD), S(f)∼f(-β), a higher-order detrending in the DFA has no adverse effect in the estimation of the DFA scaling exponent α, which satisfies the scaling relation α=(β+1)/2; (2) the upper limit of the scaling exponents detectable by the DFA depends on the order of polynomial fit used in the DFA, and is bounded by m+1, where m is the order of the polynomial fit; (3) the relation between the time scale in the DFA and the corresponding frequency in the PSD are distorted depending on both the order of the DFA and the frequency dependence of the PSD. We can improve the scale distortion by introducing the corrected time scale in the DFA corresponding to the inverse of the frequency scale in the PSD. In addition, our analytical approach makes it possible to characterize variants of the DFA using different types of detrending. As an application, properties of the detrending moving average algorithm are discussed.

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.

Information measure for long-range correlated sequences: the case of the 24 human chromosomes

A new approach to estimate the Shannon entropy of a long-range correlated sequence is proposed. The entropy is written as the sum of two terms corresponding respectively to power-law (ordered) and exponentially (disordered) distributed blocks (clusters). The approach is illustrated on the 24 human chromosome sequences by taking the nucleotide composition as the relevant information to be encoded/decoded. Interestingly, the nucleotide composition of the ordered clusters is found, on the average, comparable to the one of the whole analyzed sequence, while that of the disordered clusters fluctuates. From the information theory standpoint, this means that the power-law correlated clusters carry the same information of the whole analysed sequence. Furthermore, the fluctuations of the nucleotide composition of the disordered clusters are linked to relevant biological properties, such as segmental duplications and gene density.

Scaling range of power laws that originate from fluctuation analysis

We extend our previous study of scaling range properties performed for detrended fluctuation analysis (DFA) [Physica A 392, 2384 (2013)] to other techniques of fluctuation analysis (FA). The new technique, called modified detrended moving average analysis (MDMA), is introduced, and its scaling range properties are examined and compared with those of detrended moving average analysis (DMA) and DFA. It is shown that contrary to DFA, DMA and MDMA techniques exhibit power law dependence of the scaling range with respect to the length of the searched signal and with respect to the accuracy R^{2} of the fit to the considered scaling law imposed by DMA or MDMA methods. This power law dependence is satisfied for both uncorrelated and autocorrelated data. We find also a simple generalization of this power law relation for series with a different level of autocorrelations measured in terms of the Hurst exponent. Basic relations between scaling ranges for different techniques are also discussed. Our findings should be particularly useful for local FA in, e.g., econophysics, finances, or physiology, where the huge number of short time series has to be examined at once and wherever the preliminary check of the scaling range regime for each of the series separately is neither effective nor possible.

Evaluation of scale invariance in physiological signals by means of balanced estimation of diffusion entropy

By means of the concept of the balanced estimation of diffusion entropy, we evaluate the reliable scale invariance embedded in different sleep stages and stride records. Segments corresponding to waking, light sleep, rapid eye movement (REM) sleep, and deep sleep stages are extracted from long-term electroencephalogram signals. For each stage the scaling exponent value is distributed over a considerably wide range, which tell us that the scaling behavior is subject and sleep cycle dependent. The average of the scaling exponent values for waking segments is almost the same as that for REM segments (∼0.8). The waking and REM stages have a significantly higher value of the average scaling exponent than that for light sleep stages (∼0.7). For the stride series, the original diffusion entropy (DE) and the balanced estimation of diffusion entropy (BEDE) give almost the same results for detrended series. The evolutions of local scaling invariance show that the physiological states change abruptly, although in the experiments great efforts have been made to keep conditions unchanged. The global behavior of a single physiological signal may lose rich information on physiological states. Methodologically, the BEDE can evaluate with considerable precision the scale invariance in very short time series (∼10^{2}), while the original DE method sometimes may underestimate scale-invariance exponents or even fail in detecting scale-invariant behavior. The BEDE method is sensitive to trends in time series. The existence of trends may lead to an unreasonably high value of the scaling exponent and consequent mistaken conclusions.

Statistical tests for power-law cross-correlated processes

For stationary time series, the cross-covariance and the cross-correlation as functions of time lag n serve to quantify the similarity of two time series. The latter measure is also used to assess whether the cross-correlations are statistically significant. For nonstationary time series, the analogous measures are detrended cross-correlations analysis (DCCA) and the recently proposed detrended cross-correlation coefficient, ρ(DCCA)(T,n), where T is the total length of the time series and n the window size. For ρ(DCCA)(T,n), we numerically calculated the Cauchy inequality -1 ≤ ρ(DCCA)(T,n) ≤ 1. Here we derive -1 ≤ ρ DCCA)(T,n) ≤ 1 for a standard variance-covariance approach and for a detrending approach. For overlapping windows, we find the range of ρ(DCCA) within which the cross-correlations become statistically significant. For overlapping windows we numerically determine-and for nonoverlapping windows we derive--that the standard deviation of ρ(DCCA)(T,n) tends with increasing T to 1/T. Using ρ(DCCA)(T,n) we show that the Chinese financial market's tendency to follow the U.S. market is extremely weak. We also propose an additional statistical test that can be used to quantify the existence of cross-correlations between two power-law correlated time series.

Self-similarity of higher-order moving averages

In this work, higher-order moving average polynomials are defined by straightforward generalization of the standard moving average. The self-similarity of the polynomials is analyzed for fractional Brownian series and quantified in terms of the Hurst exponent H by using the detrending moving average method. We prove that the exponent H of the fractional Brownian series and of the detrending moving average variance asymptotically agree for the first-order polynomial. Such asymptotic values are compared with the results obtained by the simulations. The higher-order polynomials correspond to trend estimates at shorter time scales as the degree of the polynomial increases. Importantly, the increase of polynomial degree does not require to change the moving average window. Thus trends at different time scales can be obtained on data sets with the same size. These polynomials could be interesting for those applications relying on trend estimates over different time horizons (financial markets) or on filtering at different frequencies (image analysis).

Quantifying and modeling long-range cross correlations in multiple time series with applications to world stock indices

We propose a modified time lag random matrix theory in order to study time-lag cross correlations in multiple time series. We apply the method to 48 world indices, one for each of 48 different countries. We find long-range power-law cross correlations in the absolute values of returns that quantify risk, and find that they decay much more slowly than cross correlations between the returns. The magnitude of the cross correlations constitutes "bad news" for international investment managers who may believe that risk is reduced by diversifying across countries. We find that when a market shock is transmitted around the world, the risk decays very slowly. We explain these time-lag cross correlations by introducing a global factor model (GFM) in which all index returns fluctuate in response to a single global factor. For each pair of individual time series of returns, the cross correlations between returns (or magnitudes) can be modeled with the autocorrelations of the global factor returns (or magnitudes). We estimate the global factor using principal component analysis, which minimizes the variance of the residuals after removing the global trend. Using random matrix theory, a significant fraction of the world index cross correlations can be explained by the global factor, which supports the utility of the GFM. We demonstrate applications of the GFM in forecasting risks at the world level, and in finding uncorrelated individual indices. We find ten indices that are practically uncorrelated with the global factor and with the remainder of the world indices, which is relevant information for world managers in reducing their portfolio risk. Finally, we argue that this general method can be applied to a wide range of phenomena in which time series are measured, ranging from seismology and physiology to atmospheric geophysics.

Detrending moving average algorithm for multifractals

The detrending moving average (DMA) algorithm is a widely used technique to quantify the long-term correlations of nonstationary time series and the long-range correlations of fractal surfaces, which contains a parameter θ determining the position of the detrending window. We develop multifractal detrending moving average (MFDMA) algorithms for the analysis of one-dimensional multifractal measures and higher-dimensional multifractals, which is a generalization of the DMA method. The performance of the one-dimensional and two-dimensional MFDMA methods is investigated using synthetic multifractal measures with analytical solutions for backward (θ=0), centered (θ=0.5), and forward (θ=1) detrending windows. We find that the estimated multifractal scaling exponent τ(q) and the singularity spectrum f(α) are in good agreement with the theoretical values. In addition, the backward MFDMA method has the best performance, which provides the most accurate estimates of the scaling exponents with lowest error bars, while the centered MFDMA method has the worse performance. It is found that the backward MFDMA algorithm also outperforms the multifractal detrended fluctuation analysis. The one-dimensional backward MFDMA method is applied to analyzing the time series of Shanghai Stock Exchange Composite Index and its multifractal nature is confirmed.

Nonlinear dynamics of cardiovascular ageing

The application of methods drawn from nonlinear and stochastic dynamics to the analysis of cardiovascular time series is reviewed, with particular reference to the identification of changes associated with ageing. The natural variability of the heart rate (HRV) is considered in detail, including the respiratory sinus arrhythmia (RSA) corresponding to modulation of the instantaneous cardiac frequency by the rhythm of respiration. HRV has been intensively studied using traditional spectral analyses, e.g. by Fourier transform or autoregressive methods, and, because of its complexity, has been used as a paradigm for testing several proposed new methods of complexity analysis. These methods are reviewed. The application of time-frequency methods to HRV is considered, including in particular the wavelet transform which can resolve the time-dependent spectral content of HRV. Attention is focused on the cardio-respiratory interaction by introduction of the respiratory frequency variability signal (RFV), which can be acquired simultaneously with HRV by use of a respiratory effort transducer. Current methods for the analysis of interacting oscillators are reviewed and applied to cardio-respiratory data, including those for the quantification of synchronization and direction of coupling. These reveal the effect of ageing on the cardio-respiratory interaction through changes in the mutual modulation of the instantaneous cardiac and respiratory frequencies. Analyses of blood flow signals recorded with laser Doppler flowmetry are reviewed and related to the current understanding of how endothelial-dependent oscillations evolve with age: the inner lining of the vessels (the endothelium) is shown to be of crucial importance to the emerging picture. It is concluded that analyses of the complex and nonlinear dynamics of the cardiovascular system can illuminate the mechanisms of blood circulation, and that the heart, the lungs and the vascular system function as a single entity in dynamical terms. Clear evidence is found for dynamical ageing.

Fractal heterogeneous media

A method is presented for generating compact fractal disordered media by generalizing the random midpoint displacement algorithm. The obtained structures are invasive stochastic fractals, with the Hurst exponent varying as a continuous parameter, as opposed to lacunar deterministic fractals, such as the Menger sponge. By employing the detrending moving average algorithm [A. Carbone, Phys. Rev. E 76, 056703 (2007)], the Hurst exponent of the generated structure can be subsequently checked. The fractality of such a structure is referred to a property defined over a three-dimensional topology rather than to the topology itself. Consequently, in this framework, the Hurst exponent should be intended as an estimator of compactness rather than of roughness. Applications can be envisaged for simulating and quantifying complex systems characterized by self-similar heterogeneity across space. For example, exploitation areas range from the design and control of multifunctional self-assembled artificial nanostructures and microstructures to the analysis and modeling of complex pattern formation in biology, environmental sciences, geomorphological sciences, etc.

Kinetic equation of linear fractional stable motion and applications to modeling the scaling of intermittent bursts

Lévy flights and fractional Brownian motion have become exemplars of the heavy-tailed jumps and long-ranged memory widely seen in physics. Natural time series frequently combine both effects, and linear fractional stable motion (lfsm) is a model process of this type, combining alpha-stable jumps with a memory kernel. In contrast complex physical spatiotemporal diffusion processes where both the above effects compete have for many years been modeled using the fully fractional kinetic equation for the continuous-time random walk (CTRW), with power laws in the probability density functions of both jump size and waiting time. We derive the analogous kinetic equation for lfsm and show that it has a diffusion coefficient with a power law in time rather than having a fractional time derivative like the CTRW. We discuss some preliminary results on the scaling of burst "sizes" and "durations" in lfsm time series, with applications to modeling existing observations in space physics and elsewhere.

Modeling temporal fluctuations in avalanching systems

We demonstrate how to model the toppling activity in avalanching systems by stochastic differential equations (SDEs). The theory is developed as a generalization of the classical mean-field approach to sandpile dynamics by formulating it as a generalization of Itô's SDE. This equation contains a fractional Gaussian noise term representing the branching of an avalanche into small active clusters and a drift term reflecting the tendency for small avalanches to grow and large avalanches to be constricted by the finite system size. If one defines avalanching to take place when the toppling activity exceeds a certain threshold, the stochastic model allows us to compute the avalanche exponents in the continuum limit as functions of the Hurst exponent of the noise. The results are found to agree well with numerical simulations in the Bak-Tang-Wiesenfeld and Zhang sandpile models. The stochastic model also provides a method for computing the probability density functions of the fluctuations in the toppling activity itself. We show that the sandpiles do not belong to the class of phenomena giving rise to universal non-Gaussian probability density functions for the global activity. Moreover, we demonstrate essential differences between the fluctuations of total kinetic energy in a two-dimensional turbulence simulation and the toppling activity in sandpiles.

Serial correlation of detrended time series

A preliminary essential procedure in time series analysis is the separation of the deterministic component from the random one. If the signal is the result of superposing a noise over a deterministic trend, then the first one must estimate and remove the trend from the signal to obtain an estimation of the stationary random component. The errors accompanying the estimated trend are conveyed as well to the estimated noise, taking the form of detrending errors. Therefore the statistical errors of the estimators of the noise parameters obtained after detrending are larger than the statistical errors characteristic to the noise considered separately. In this paper we study the detrending errors by means of a Monte Carlo method based on automatic numerical algorithms for nonmonotonic trends generation and for construction of estimated polynomial trends alike to those obtained by subjective methods. For a first order autoregressive noise we show that in average the detrending errors of the noise parameters evaluated by means of the autocovariance and autocorrelation function are almost uncorrelated to the statistical errors intrinsic to the noise and they have comparable magnitude. For a real time series with significant trend we discuss a recursive method for computing the errors of the estimated parameters after detrending and we show that the detrending error is larger than the half of the total error.

Algorithm to estimate the Hurst exponent of high-dimensional fractals

We propose an algorithm to estimate the Hurst exponent of high-dimensional fractals, based on a generalized high-dimensional variance around a moving average low-pass filter. As working examples, we consider rough surfaces generated by the random midpoint displacement and by the Cholesky-Levinson factorization algorithms. The surrogate surfaces have Hurst exponents ranging from 0.1 to 0.9 with step 0.1, and different sizes. The computational efficiency and the accuracy of the algorithm are also discussed.

Detrended fluctuation analysis for fractals and multifractals in higher dimensions

One-dimensional detrended fluctuation analysis (DFA) and multifractal detrended fluctuation analysis (MFDFA) are widely used in the scaling analysis of fractal and multifractal time series because they are accurate and easy to implement. In this paper we generalize the one-dimensional DFA and MFDFA to higher-dimensional versions. The generalization works well when tested with synthetic surfaces including fractional Brownian surfaces and multifractal surfaces. The two-dimensional MFDFA is also adopted to analyze two images from nature and experiment, and nice scaling laws are unraveled.

Quantifying signals with power-law correlations: a comparative study of detrended fluctuation analysis and detrended moving average techniques

Detrended fluctuation analysis (DFA) and detrended moving average (DMA) are two scaling analysis methods designed to quantify correlations in noisy nonstationary signals. We systematically study the performance of different variants of the DMA method when applied to artificially generated long-range power-law correlated signals with an a priori known scaling exponent alpha(0) and compare them with the DFA method. We find that the scaling results obtained from different variants of the DMA method strongly depend on the type of the moving average filter. Further, we investigate the optimal scaling regime where the DFA and DMA methods accurately quantify the scaling exponent alpha(0) , and how this regime depends on the correlations in the signal. Finally, we develop a three-dimensional representation to determine how the stability of the scaling curves obtained from the DFA and DMA methods depends on the scale of analysis, the order of detrending, and the order of the moving average we use, as well as on the type of correlations in the signal.

Some properties of the entropy in the natural time

We show that the entropy S , defined as S identical with chi ln chi - chi ln chi [Phys. Rev. E 68, 031106 (2003)] where chi stands for the natural time [Phys. Rev. E 66, 011902 (2002)], exhibits positivity and concavity as well as stability or experimental robustness. Furthermore, the distinction between the seismic electric signal activities and "artificial" noises, based on the classification of their S values, is lost when studying the time-reversed signals. This reveals the profound importance of considering the (true) time arrow.