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.

Introduction

Sequences in complex systems generally exhibit fluctuations on a wide range of time scales of the data, where the distinct fluctuations follow a scaling relation over several orders of magnitude. Such scaling laws provide a description of the data and the complex system by fractal scaling exponents. Nevertheless, many records do not present a simple monofractal scaling behavior and different parts of the series require different scaling exponents. In even more complicated cases, such different scaling behavior can be observed for many interwoven fractal subsets of the time series, where a multitude of scaling exponents is required for a full characterization of the scaling behavior in the same range of time scales and a multifractal analysis must be applied. In order to observe the fractals and multifractals in time series, algorithms have been developed one by another [1], [2].

Since detrended fluctuation analysis (DFA) has been proposed by Peng et al. [3] to detect the long-range power-law correlations in DNA sequences, it has been successfully applied to diverse fields [4], [5], [6], [7], [8], [9], [10], [11]. Then Podobnik and Stanley [12] generalized DFA and introduced the detrended cross-correlation analysis (DCCA) for two non-stationary time series, which has aroused increasing interest in analysis of long-range cross-correlation and multifractality [1], [13], [14], [15], [16], [17], [18], [19], [20], [21], [22], [23], [24], [25]. Specifically, the analysis is based on the bivariate Hurst exponent h estimation, which is related to an asymptotic power-law decay of the cross-correlation function. A power-law cross-correlated process has the cross-correlation function for and the cross-power spectrum for . is characteristic for the absence of power-law cross-correlation, while processes with h > 0.5 are cross-persistent and h < 0.5 indicates the anti-persistent cross-correlation of the data [2], [12], [25]. Basis of these, the cross-correlation coefficient (σ) was introduced with the objective of quantifying the level of cross-correlation between non-stationary time series [26], [27], [28]. In the DFA–DCCA-σ line, Kristoufek [29] has recently introduced the DFA framework to estimate regression parameters at different scales and under potential non-stationarity and power-law correlations. Alternatively, the detrending moving-average (DMA) method can also be used for fractal or multifractal analysis, found that the backward multifractal detrending moving average (MFDMA) algorithm outperforms the multifractal detrended fluctuation analysis (MFDFA) [30], [31], [32], [33]. As an alternative to σ, based on the detrending moving-average cross-correlation analysis (DMCA), the DMCA coefficient ρ was proposed as well [17], [34]. To investigate the multifractal features of two cross-correlated non-stationary series, several methods have been proposed for this purpose. Zhou [14] further proposed multifractal detrended cross-correlation analysis (MFDCCA) as a generation of DCCA. Jiang and Zhou [35] developed a class of MFDCCA algorithms based on DMA, called MFXDMA whose performances are comparative to MFDCCA’s. Moreover, a new method for the detection of long-range cross-correlations and multifractality based on scaling of qth order covariances, multifractal height cross-correlation analysis (MF-HXA), was proposed by Kristoufek [18]. Meanwhile, statistical moments function or partition function has been widely applied to analyze the multifractal features of financial time series [36], [37], [38], [39], [40]. Based on this approach, Jiang and Zhou [38] investigated the multifractal characteristics in intraday minutely data for four stock market indexes within individual trading days and found the so-called multifractality is merely an illusion. Dissimilarly, when they performed the partition function approach on the 1-min volatility of two indexes and 1139 stocks in the Chinese stock markets, multifractal nature is significant [39]. Extending that from one dimension to two dimensions, Wang et al. [40] introduced the multifractal cross-correlation analysis based on statistical moments (MFSMXA) as an application of the joint multifractal measures [41], finding that MFSMXA has comparative performance with MFXDMA and sometimes perform slightly better than MFDCCA. By far, multifractal analysis has been employed successfully in various provinces, such as human nature [42], [43], [44], financial time series [45], [46], [47], river flow [48] and traffic signals [49], etc.

However, there are limitations in these algorithms and the common one is that they ignored the difference of the measurement in each segment that contributes to the final fluctuation functions unequally theoretically. While in other academic fields, weighting function is selected for improvement, such as Shannon entropy [50], [51], [52], [53], variance or energy [54], and so on [55]. Inspired by their work, basis of MFSMXA, we propose in this paper a weighted method based on Shannon entropy (W-MFSMXA) to analyze the multifractal cross-correlation between two series. Results show that W-MFSMXA is comparable to and in some extent better than MFSMXA.

The remainder of this paper is organized as follows. In Section 2, the methodologies of MFSMXA and W-MFSMXA are introduced. In Section 3, we test the effectiveness of the W-MFSMXA algorithms with three types of artificial time series: two-component ARFIMA stochastic processes, binomial multifractal model and the NBVP time series. Performances are also compared with these in MFSMXA. Application to financial time series is presented in Section 4. At last, Section 5 gives the conclusions and appendix exhibits charts of the finite-size effect for the binomial multifractal model.

Methodologies

MFSMXA method

The MFSMXA method [40] consists of four steps. Consider two time series {X(i)} and {Y(i)} of the same length N, where .

Divide each series into N = int(N/s) non-overlapping segments of equal length s. For a given segment size s, the elements in the vth segment can be described as: where and .

For each segment, we define a quantity u as follows, where is the vth segment. The measure μ in each segment is constructed as follows,

Calculate the statistical moments function F(q, s)

Determine scaling behavior of cross statistical moment functions by analyzing log-log plots of F(q, s) versus s for each value of q. We expect it to scale as where the exponent τ(q) is a scaling exponent function, regarded as a characteristic function of the fractal behavior. Note that . For the method is used to detect the long-range cross-correlations between series, which we call it cross-correlation analysis based on statistical moments (SMXA). When , the above method reduces to the multifractal analysis for separate series (MFSMA) [36], [37], [38], [39].

In light of the relationship of two sets of multifractal scaling exponents, there is an analytic form between τ(q) and the generalized bivariate Hurst exponent h(q): where D is the fractal dimension of the geometric support of the multifractal measure and for time series analysis [1]. If q dependence of the scaling exponent τ(q) is nonlinear, there is multifractality in the signal. The singularity strength λ and its spectrum f(λ) related to τ(q) through a Legendre transform [56] can be defined as, If the singularity width gets close to zero, the series has practically no multifractality. Otherwise, there is a multifratcal nature in the series.

According to Ref. [18], for time series in a standard way, and Further, the left part of Eq. (9) can be rewritten through the definition of covariance: and combining Eqs. (6) and (8), the first part of the right side of Eq. (10) implies i.e, Now consider a new scaling exponent η(q) and relationship If then holds, which is identical with in Ref. [18]. Otherwise, it does not establish unless covariances do not scale with s. Thus, the importance of long-range cross-correlations and multifractality resides in the scaling of covariances between |μ(v, s)| and |μ(v, s)| regarding varying s and q. As a consequence, there are two types of cross-multifractality: (i) cross-multifractality caused by multifractality of the separate series, and (ii) cross-multifractality by scaling of covariances between |μ(v, s)| and |μ(v, s)|.

Weighted MFSMXA method

It can be noted that in the MFSMXA, the product μ of each segment in both sequences contributes to the final function F(q, s) equally, but it is rational to make F(q, s) be proportional to μ in a simplest view. And to a certain extent, the measure μ is kind of frequency or probability. To improve that, we propose a modification of current MFSMXA, called weighted MFSMXA method based on Shannon entropy (W-MFSMXA), to take into account that not all measurements contribute equally when constructing the statistical moments function and to save as much information as possible carried by the series. Specifically, the motivation behind this is to weight proportionally segments with distinct measurements μ but the same final F(q, s). For example, consider a simple case and let and where and . Generally we will get but if we consider the effect of weight on the contrary, and are not equal. Weights here are and according to the procedure described detailedly in the following. The selection of Shannon entropy [57] to be the weight factor makes sure that μ with higher value is assigned more weights. Furthermore, Shannon entropy as a complexity (or regularity or orderliness) measure quantifies the information richness of μ in series and this novel information-based approach has been successfully employed in literary authorship disputes [50], SARS coronavirus [51], biologic signals [52] and financial time series analysis [53].

Therefore, we replace the statistical moments function in step 3 in MFSMXA method with a weighted statistical moments function, which is calculated as follows: where the weighting function ω(v, s) is selected to be the sum of Shannon entropy for the vth segment of series {x(i)} and {y(i)} [50], [51], [52], [53]: and Z is a normalization factor such that . The normalization factor Z in Eq. (16) is given by Then, same as MFSMXA in the step 4, we expect to scale as It reduces to the long-range cross-correlation analysis in a weighted way for (W-SMXA) and the multifractal analysis for individual series when {X(i)}={Y(i)} (W-MFSMA). Note that .

Numerical results for artificial time series

In this section, we use three types of artificial time series, the two-exponent ARFIMA processes [58], binomial multifractal model [56] and the NBVP time series [59], to test the effectiveness of this new procedure and make comparison between two methods. The latter two have been proved the existence of cross-multifractality.

Two-component ARFIMA stochastic processes

The power-law autocorrelations in stochastic variables can be modeled by the fractionally autoregressive integrated moving-average process (ARFIMA) process [60]: where d ∈ (0, 0.5) is a memory parameter related to Hurst exponent as [12], and a(d) are the weights defined by where Γ denotes the Gamma function and n is the time scale.

The two-component ARFIMA stochastic processes can be defined by two variables x and y [58]: where a(d 1) and a(d 2) are the weights defined in Eq. (19) through the scaling parameters d 1 and d 2 (0 < d 1, 2 < 0.5), and W is a free parameter controlling the coupling strength between x and y (0.5 ≤ W ≤ 1). When the process in Eq. (20) reduces to two decoupled ARFIMA processes defined in Eq. (19), and consequently the long-range cross-correlations between x and y vanish while both x and y remain long-range power-law auto-correlated. Additionally, the mixed-correlated ARFIMA process (MC-ARFIMA) has been proposed as a good alternative to the two-component ARFIMA process with analytically solved asymptotic properties [61].

In our simulations, and parameters (d 1, d 2) are set as and separately. Series preprocessed with are then investigated by the proposed W-MFSMXA. Figs. 1 a–c show the numerical results of the process in Eq. (20) with . Fig. 1a illustrates the dependence of the statistical moments function regarding the scale s in log-log plots for q = 4:2:4. Excellent power-law cross-correlations are observed. Figs. 1b and c show the corresponding scaling exponents and singularity spectrum respectively. The linear dependence of with respect to q as well as the singularity widths closer to 0 unveils the fractal property in two-component ARFIMA stochastic processes. Similar results can also be obtained for the process in Eq. (20) with and in Figs. 1d–f. The singularity width Δλ ( under this case is 0.0076.

Fig. 1

Long-range cross-correlation analysis of the two-component ARFIMA processes in Eq. (20) by W-MFSMXA. . (a–c). Processes with . (a) The power-low dependence of the statistical moments function versus the scale s in log-log plots for . The dashed lines are the best power-law fits for . (b) The corresponding scaling exponents versus the moment orders q. The straight line is the linear fit for . (c) The singularity spectrum versus the singularity strength λ. (d–f). Processes with and . Similar to (a–c), (d), (e) and (f) are plotted as well.

Binomial multifractal series

Using the binomial multifractal model (BMF) [56], we construct several binomial multifractal series {x(i)} with the length of . The kth number in series is defined as where 0 < p < 0.5 is a parameter, n(k) is the number of digits equal to 1 in the binary representation of index k. For example, since 14 corresponds to binary 1110. Series generated by this model has known analytic multifractal properties. The scaling exponent τ(q) has an analytic form [62]: satisfying Eq. (14) [35], [40].

In this paper, we generated series with length of 216 with for {x(i)}, for {y(i)} and for {z(i)}, just same with these in Ref. [40]. The two pairs of time series, {x(i)} & {y(i)} and {y(i)} & {z(i)}, are highly correlated with a coefficient of 0.8307 and 0.8316 respectively while {x(i)} & {z(i)} significantly correlated with the coefficient of 0.4709. It derives from the fact that these sequences are constructed according to the same rules despite variations of parameters [35]. Besides, BMF has been successfully implemented in multifractal analysis and the cross-multifractality between two series is significant [14], [18], [35], [40]. Nevertheless, cross-persistence in BMF is mainly due to the persistence of the separate processes and the fact that the series are correlated [18]. For comparison, we applied MFSMXA and W-MFSMXA methods to analyze these series. Results are depicted in the Figs. 2 and 3 .

Fig. 2

For three series generated by BMF in Eq. (21) with and , respectively, statistical moments functions in MFSMXA and W-MFSMXA with respect to s in log-log plots for q from 3 to 3 (from top to bottom) with the step-length of 1. Correspond scaling exponents τ(q) are obtained. (a) F(q, s) and versus s for series {x(i)} and {y(i)}. . (b) F(q, s) and versus s for series {x(i)} and {z(i)}. . (c) F(q, s) and versus s for series {y(i)} and {z(i)}. .

Fig. 3

Multifractal analysis of BMF series when q is from 5 to 5 with the step-length of 0.2. (a) Scaling exponents τ(q) in MFSMXA (filled circles) and W-MFSMXA (open squares). Different colors stand for different types of scaling exponents: τ(q) (blue), τ(q) (black) and τ(q) (red). Theoretical values are calculated by Eq. (22). (b) Scaling exponents τ(q) obtained for separate series in W-MFSMA. To further verify the Eq. (14), and the average of and are also obtained. (c) The singularity spectra f(λ) for the separate series in W-MFSMA. The solid lines are the corresponding theoretical values in Eq. (23). (d) The singularity spectra f(λ) between two series in W-MFSMXA. The solid lines are their average values calculated according to Eq. (23). (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article).

The statistical moments functions F(q, s) of the two series {x(i)} and {y(i)} against s in log-log plot are presented in Fig. 2a for both MFSMXA and W-MFSMXA methods when q is from 3 to 3 with the step-length of 1. Nice power-law properties are observed for both positive and negative q values in three cases. Moreover, what’s intriguing is that the procedure after being weighted did not change the structure of the statistical moments function, which is particularly evident in Fig. 2c. It is noteworthy that in W-MFSMXA is more synchronous with in MFSMXA apparently when q is non-negative and especially when q equals zero in W-MFSMXA while one in MFSMXA.

Analogous phenomena turn up when plotting the scaling exponent function τ(q) versus the moment order q in Fig. 3a. Obviously, are nonlinear and upwards convex, which indicates the multifractal nature between two time series. What’s more, and τ(q) are almost parallel to each other. When shifting the figure of to right with unit 1, and τ(q) coincide pretty well (not shown here), which is identical with closer to . Similar results can be obtained for series {x(i)} and {z(i)} and for series {y(i)} and {z(i)}, which are obvious in Figs. 2 and 3a. Table 1 shows the corresponding singularity widths Δλ captured before and after the weighting procedure. For two algorithms, Δλ are much larger than zero and differences of Δλ deviate insignificantly from zero. Thus, we assume that except for the phase, the multifractal structure of the series will not be changed after being weighted, and it’s reasonable to believe there is a delay in MFSMXA compared to W-MFSMXA.

Table 1

The singularity widths Δλ obtained by MFSMXA and W-MFSMXA algorithms for BMF series.

Δλ		MFSMXA	W-MFSMXA	Diff.
BMF	Δλxy	1.5708	1.5613	0.0095
	Δλxz	1.2401	1.2405	0.0004
	Δλyz	0.8146	0.8134	0.0012

The power-law dependence of and versus s for separate series {x(i)} and {y(i)} based on weighted statistical moments are also displayed in Fig. 3b, which demonstrates that the slope of weighted cross statistical moments function is the average slopes of and satisfying Eq. (14). Same results can also be obtained when computing the slopes of and . As the scaling exponent function for q in W-MFSMXA is synchronous with that in MFSMXA for then we assume and theoretically when Note the difference of the q index is 1 rather than the step-length. The singularity spectrum f(λ) of the theoretical analysis in Eq. (23) is displayed in Fig. 3c. Numerical results are significantly in line with their theoretical values. Fig. 3d shows the singularity spectra between two series and they basically accord with the averages of theoretical values on the basis of Eq. (23).

It is worthwhile to note that short length N of the investigated series or small scale s often results in a spurious detection of multifractal behavior from a monofractal model [19]. To evaluate the finite-size effect [25] of the BMF model, performances of the scaling exponent τ(q) estimators in MFSM(X)A and the proposed W-MFSM(X)A are studied. In all cases, time series with length are analyzed, where . The scales s are given by where and round(X) means rounding the elements of X to the nearest integers. For series with different parameters in Eq. (21), results of the evaluation of (cross-) multifractality are presented in appendix. Significantly, for analyses of separate series in MFSMA and W-MFSMA (see Figs. 9 and 10), two τ(q) estimators, namely τ(q) and have the same performances: (1) For sequences with different length sizes, numerical results of τ(q) and coincide well with the corresponding theoretical values respectively except for ; (2) For in Figs. 9e and 10 e, all series perform well especially in W-MFSMA; (3) Generally, are acceptable while are more reliable, which is accordant with the results of MFSMXA procedure between two time series in Fig. 11. Results of between two time series in W-MFSMXA are plotted in Fig. 12. Yet again, nine pairs of data used here for cross-multifractal analysis are significantly correlated (see Table 2). It is clear to see that deviate significantly from the theoretical values when same as in MFSM(X)A and W-MFSMA. In most of cases, performances of when are comparative. While when and (see Figs. 12b and c), are the best. In general, on the basis of investigations in this section, performances of W-MFSMA are as good as MFSMA’s and are suggested for use in MFSM(X)A or W-MFSM(X)A. Nevertheless, W-MFSMXA outperforms MFSMXA for series with relatively shorter length, i.e., . At this point, W-MFSM(X)A is comparable to MFSM(X)A. In order not to seem rigmarole, we only illustrate the results of one pair of experimental data in following experiments.

Fig. 9

The scaling exponents τ(q) for separate series in MFSMA. Theoretical values are calculated by Eq. (22) for different parameters p.

Fig. 10

The scaling exponents τ(q) for separate series in W-MFSMA. Corresponding theoretical values are calculated according to Eq. 23 for different parameters p.

Fig. 11

The scaling exponents τ(q) between two time series in MFSMXA. Theoretical values are calculated in light of Eqs. (14) and (22).

Fig. 12

The scaling exponents τ(q) between two series in W-MFSMXA. Theoretical values are calculated based on Eqs. (14) and (23).

Table 2

Correlation coefficients of two series generated by the process in Eq. (21) with varied parameters.

Corrcoef		py
		0.1	0.2	0.25	0.3	0.4
px	0.1	1	0.8930	0.7606	0.5980	−
	0.2		1	0.966 2	0.8676	0.5667
	0.25			1	0.9642	0.7257
	0.3				1	0.8689
	0.4					1

In fact, from the log-log plot, the function in Eq. (18) can be defined as follows [44], [63], [64]: so can be the Eq. (5). Combining Eqs. (5), (18) and (24), we can get a new equation: where a(q) is the exponent of the best power law fit for versus F(q, s) in logarithms, which can be regarded as a scaling exponent ratio related to the moment order q and used to explain why there is a delay. As for positive time series {x(i)} and {y(i)}, the values of F(1, s) and F(1, s) are equivalent to the sum of the entire series individually, which is uncorrelated with the scale s. Therefore, according to Eq. (14), and a(1) diverges. Notice .

The a(q) with respect to q for series {x(i)} and {y(i)} is depicted in Fig. 4a. In Fig. 3a, and τ(q) increase with the increasing orders q and for all values of q. Thus, for then 0 < a(q) < 1; for 1 < q < 5, then a(q) > 1; otherwise, and then a(q) < 0. Since the a(q) functions diverges at a(1) is removed for better illustrations as well in the following experiments. Then, the graph of a(q) with respect to q can be approximately viewed as a centrosymmetric hyperbola. For q → ±5, a(q) approximates to 1, demonstrating two methods are practically consistent, while differences are apparent for relatively smaller orders q. Actually, for higher orders q > 5 or a(q) → 1 establishes as well. However, a(q) ≠ 1 holds for all q, consequently delay between W-MFSMXA and MFSMXA is substantial. As we can see in Fig. 4b, the profiles of and τ(q) after multiplying corresponding a(q) are coincident, confirming the Eq. (25).

Fig. 4

Schematic diagram of Eq. (25) for series {x(i)} and {y(i)}. (a) Scaling exponent ratio a(q) against q. For better captions, a(1) is eliminated as it diverges. Insert figure: weighted statistical moment function with respect to F(q, s) in log-log plots for and . The black dash lines are the best fits and the exponents are and . (b) Scaling exponent τ(q) against q. τ(q) and the theoretical value of τ(q) multiplied by a(q) in Eq. (25).

NBVP time series

Deriving from the classical continuous BVP oscillator, a two-dimensional map by using the forward Euler discrete scheme can be written as: where 0 < ρ ≪ 1, 0 < a < 1, 0 < b < 1, μ is a stimulus intensity and 0 < δ < 1 is the step size. Given a set of variables and initial values x 0 and y 0, we can get time series {x(i)} and {y(i)}.

First, we generate two more complex BVP series with the length of which are called NBVP time series, whose initial value μ is a random variable, satisfying μ ∼ Norm(0.2, 1) and other initial values are constant values, with and where the cross-correlation relationship is found to be multifractal by Xue et al. [59]. To our knowledge, no analytic expression has been generated for this model. Series analyzed here are with a given preparation as in Section 3.1, then we apply MFSMXA and W-MFSMXA algorithms to the NBVP series to examine the validity of the assumption about delay in Section 3.2. Results are shown in Fig. 5. Singularity widths in this case are and which are insignificantly deviated. Same as the previous part, except for the position and several discrepant points, the multifractal structure of the series remains unchanged. and τ(q) after multiplying correspond a(q) are nearly coincident and the hyperbolic characteristics of a(q) with respect to q is more obvious. In other words, the τ(q) measurements in both method are apparently different for relatively smaller q but more accordant for higher q while never equal. Proposition about the delay is further verified. It can be seen in Fig. 5d, the scaling of covariances is higher than the average of scaling of separate processes for negative moment orders q, the deviations are nevertheless not significant for higher positive q. In contrast, basically coincides with τ(q), implying that the cross-multifractality of series is mainly caused by the scaling of covariances between |μ(v, s)| and |μ(v, s)| as the construction of the NBVP mixes the long memory of the separate processes.

Fig. 5

Multifractal cross-correlation analysis of the NBVP time series. (a) Statistical moments functions F(q, s) in MFSMXA and W-MFSMXA with respect to s are shown in log-log plots for q from 3 to 3 (from top to bottom) with the step-length of 1. (b) Scaling exponents τ(q) are calculated by the proposed W-MFSMXA and MFSMXA. The theoretical values according to Eq. (25) are also performed. (c) Scaling exponent ratio a(q) as a function of q. (d) Scaling of covariances between |μ(v, s)| and |μ(v, s)| by MFSMXA according to Eq. (13). Scaling of covariances here is plotted as for better comparison. τ(q), τ(q) and their average, and τ(q) are also shown.

Numerical results for financial time series

In order to prove the applicability of the proposed W-MFSMXA method in real world data, we then apply this method to financial time series to analyze their cross-multifractality. Here, we perform two US stock indices: the Dow Jones Industrial Average (DJI) and the National Association of Securities Dealers Automated Quotations (NAS). The original data were obtained from yahoo finance (http://finance.yahoo.com) from Jan. 2, 1997 to Dec. 29, 2014.

To make sure of the synchronization of the time series, we exclude the asynchronous data and then reconnect the remaining parts of the original series. In so doing, the resulting length is 4526. The data we analyze in this paper are defined as the absolute logarithmic returns, this is, where x denote the closing price for a stock index on day t. The daily absolute returns for the two stock indices are shown in Fig. 6. Profiles of two series are very similar to each other.

Fig. 6

The absolute logarithmic returns for DJI and NAS from Jan. 2, 1997 to Dec. 29, 2014.

Let x denotes DJI series and y denotes NAS throughout this section. Fig. 7a depicts the statistical moments function F(q, s) as a function of scale s in log-log plots for q from 3 to 3 with the step-length of 1 by W-MFSMXA and MFSMXA methods. Excellent power laws are also observed for both algorithms. Same as before, in W-MFSMXA is more similar to in MFSMXA, which is conformed when plotting the scaling exponent function τ(q) versus the moment order q in Fig. 7b. The scaling exponent ratio a(q) is plotted in Fig. 7c and clearly approximates a hyperbola. The nonlinearity of τ(q) demonstrates the existence of but weak multifractality in the absolute logarithmic daily returns series. Besides, in Fig. 7d, τ(q) are practically consisitent with the averages of τ(q) and τ(q) while the scaling of covariances is obviously the lowest, indicating that the cross-multifractality in series is mainly dominated by the multifractality of separate series.

Fig. 7

Multifractal cross-correlation analysis for the absolute logarithmic daily returns series. (a) Statistical moments functions in MFSMXA and W-MFSMXA with respect to s. (b) The scaling exponents τ(q) are calculated by the proposed W-MFSMXA and MFSMXA. Theoretical values in light of Eq. (25) are also obtained. (c) Scaling exponent ratio a(q) as a function of q. is removed for its divergence. (d) Scaling of covariances between |μ(v, s)| and |μ(v, s)| by MFSMXA in Eq. (13). τ(q), τ(q) and their average, and τ(q) are also depicted.

To give a more systemic statistical test on these returns series, we adopt the following bootstrapping approach to get the possibility that the empirical multifractality could be spurious [38], [39]. We reshuffle each series 1000 times to remove any potential temporal correlation, then carry out the W-MFSMXA procedure for the surrogate data as for the original data. A very strict null hypothesis is imposed to investigate whether the singularity spectrum f(λ) is wider than those produced by chance. The null hypothesis is the following: The associated probability of false alarm for multifractality is defined by where n is the number of shufflings and counts the number of Δλ whose value is not greater than Δλ. As n → ∞, it is clear that the estimated bootstrap p-value will tend to the ideal bootstrap p-value. Under the conventional significance level of 0.01, the multifractal phenomenon is statistically significant if and only if p ≤ 0.01. When p > 0.01, the null hypothesis cannot be rejected.

Results of the statistical test are presented in Fig. 8. A remarkable deviation of the singularity spectrum f(λ) of the original data from f(λ) of the surrogate data is clear to see in Fig. 8b. Hence, we can infer that cross-multifractality in series is destroyed by the shuffling procedure as a consequence of the eliminated long memory in the separate returns series. Further, we find using the significance level of 0.01, which is a significant evidence for the presence of cross-multifractality in the investigated DJI and NAS data. Comparing the scaling exponent ratio a(q) in Figs. 4a, 5 c and 7 c, it can be found the similar hyperbolic feature of a(q). Moreover, the statistical moments functions as well as the scaling exponents in W-MFSMXA are almost exactly the same as those in MFSMXA if taking no account of their phases or positions.

Fig. 8

Multifractal cross-correlation analysis for the original and shuffled series in W-MFSXMA. τ(q) and f(λ) for the shuffled series are accordingly calculated by the average of results of 1000 surrogate data. (a) Scaling exponents τ(q). (b) Singularity spectra f(λ). . (c) Singularity widths Δλ of the shuffled data. .

Conclusions

In this paper, we have proposed a modified MFSMXA method called weighted MFSMXA method based on Shannon entropy (W-MFSMXA). To test the robustness of this new procedure and compare to the heuristic MFSMXA, we first employed three types of artificial time series: the two-exponent ARFIMA processes with monofractal feature, the binomial multifractal model (BMF) and the NBVP time series, which have been proved the existence of multifractal properties. However, unlike BMF, cross-multifractality in NBVP is mainly caused by the scaling of covariances rather than the correlated separate processes.

Long-range cross-correlation analysis of the two-component ARFIMA stochastic processes by W-MFSMXA uncovered the fractal nature in these power-law cross-correlated series. Results of multifractal analyses show that the statistical moments function versus scale s for q in W-MFSMXA is more similar to that in MFSMXA for . The scaling exponents obtained by two methods are almost parallel to each other. What’s more, differences of the corresponding singularity widths captured by two algorithms do not deviate significantly from zero. Thus, it is rational to believe that the multifractal structure of series has not been changed by the weighting procedure without considering their positions. Basis of these, we generated analytic formulas of the BMF model for W-MFSMXA. Theoretical analysis and finite-size effect test verified the good performances of W-MFSMXA. Specifically, numerical results are well in line with the theoretical values and W-MFSMXA outperforms MFSMXA for analysis of series with relatively shorter length. Meanwhile, we generated the scaling exponent ratio a(q), whose profile is approximately a centrosymmetric hyperbola, to describe the relation between W-MFSMXA and MFSMXA. Scaling exponent estimators in both methods are more consistent for relatively higher orders q while never equal, indicating the delay in MFSMXA compared to W-MFSMXA.

At last, we applied W-MFSMXA to the absolute logarithmic returns of two US stock market indices: the Dow Jones Industrial Average (DJI) and the National Association of Securities Dealers Automated Quotations (NAS). The results obtained are same with those obtained for artificial time series. Besides, weak cross-multifractal nature is found in the returns series and their cross-multifractality is mainly due to the multifractality of separate series. Further, we compared outcomes of the original and shuffled series. Statistical test proved that cross-multifractality in the investigated DJI and NAS data is destroyed by the shuffling procedure as a consequence of the removed long memory in separate returns series.

In general, all these illustrations make we believe that the proposed W-MFSMXA method not only will not change the multifractal structure of series but contains more significant information of series compared to the previous MFSMXA method for the introduced weighting function. In some extent, it’s comparable to and better than MFSMXA. Though there is no criterion for Shannon entropy as the weight factor as the best choice, it did exhibit elegant performances as we can see in the text. This goes nevertheless beyond the scope of current work.