Abstract

In this paper, we develop a new method to measure the nonlinear interactions between nonstationary time series based on the detrended cross-correlation coefficient analysis. We describe how a nonlinear interaction may be obtained by eliminating the influence of other variables on two simultaneous time series. By applying two artificially generated signals, we show that the new method is working reliably for determining the cross-correlation behavior of two signals. We also illustrate the application of this method in finance and aeroengine systems. These analyses suggest that the proposed measure, derived from the detrended cross-correlation coefficient analysis, may be used to remove the influence of other variables on the cross-correlation between two simultaneous time series.

1. Introduction

There are numerous real-world systems where the output signals are nonstationary and exhibit complex self-correlation or cross-correlation over a broad range of time scales. The output signals can be characterized by power-law correlations. One method, which has proved to be quite useful to detect the degree of interrelation between two stationary variables, is Pearson’s correlation coefficient [1]:where is the arithmetic average of and is its standard deviation and likewise for . Proposition of Pearson's correlation coefficient (PCC) has achieved great success in multivariate analysis, such as the principal component analysis [2], random matrix theory [3], and singular value decomposition [4].

Nevertheless, in real-world systems, nonlinear and nonstationary characteristics are present. Therefore, PCC may not be suitable to describe the interrelation between two variables that are nonlinear and nonstationary. For dealing with the drawbacks of PCC, the detrended cross-correlation analysis (DCCA) method and the DCCA coefficient are proposed by Stanley and Podobnik [5, 6]. The advantage of the DCCA method is that it allows the detection of cross-correlations between noisy signals with embedded polynomial trends, which can mask the true cross-correlations in the fluctuations of signals. The DCCA method is widely applied to measure the cross-correlations in different fields, such as social sciences [7], biology [8], climatology [9], geophysics [10, 11], transportation [12, 13], seismic signals [11, 14], economics [1520], and aeroengine dynamics [2124].

Recently, multifractal analysis is one of the major interests for researchers from interdisciplinary domains to uncover the scaling properties and understand the hidden information. Among these researchers, many of them applied the multifractal analysis to meteorology [2527], electroencephalography [28], and economics [2931]. Later, as some researchers thought of extending the research of multifractal analysis to the detrended cross-correlations between time series, the multifractal detrended cross-correlation analysis (MFDXA) was proposed [3234].

The cross-correlation between two variables may be influenced by other variables. Hence, we have to be alert to the possibilities of spurious correlation while investigating the cross-correlation. Then, the methods of partial correlation and partial correlation coefficient are therefore proposed to measure the degree of association between two random variables [35, 36]. The linear effect may be removed using the partial correlation coefficient (partial CC):where and to minimize the mean and likewise for . If n additional variables are to be accounted for, say , the nth-order partial CC can be computed by [36]

Lately, the detrended partial cross-correlation analysis and multifractal detrended partial cross-correlation analysis (MFDPXA) which can measure cross-correlations between nonlinear time series influenced by common external forces is proposed [37, 38].

In order to remove the spurious correlation and improve the estimation performance for quantifying the intrinsic interactions between two nonstationary time series, this paper proposes the method of nth-order multifractal detrended partial cross-correlation analysis by incorporating the partial correlation coefficient with the multifractal detrended cross-correlation analysis.

The rest of the paper is organized as follows. In the next section, we introduce the multifractal DCCA coefficient method and propose the method of nth-order multifractal detrended partial cross-correlation analysis. In Section 3, we show the data results for the randomly generated dataset and stock and engine dataset by the proposed methods. Finally, we draw some conclusions in Section 4.

2. Methodologies

2.1. Multifractal Detrended Partial Cross-Correlation Analysis

For the sake of clarity, we begin with a summary of the multifractal DCCA coefficient algorithm. For two series and with equal length , where , the computational procedure of the multifractal DCCA coefficient is as follows:Step 1: construct the profile of each series by eliminating the mean value:where and are the average values of and , respectively.Step 2: divide the profiles and into nonoverlapping units of equal length . Considering that is usually not a multiple of the time scale , we repeat the same procedure by starting from the opposite end of the sequence in order to take the whole series into account. Thus, we obtain segments of equal length . In this paper, we follow the previous literature practice and set .Step 3: for each segment , the local trends and are estimated on the basis of a least-squares fit of the sequences and , respectively. The corresponding detrended covariance for isand for iswhere and are the fitting polynomials in the segment .Step 4: calculate the average of multifractal detrended covariance fluctuation function over all segments:Generally, can take any real value, except zero. For , the equation becomesFor , is equal to the detrended cross-correlation fluctuation function .Step 5: estimate the multifractal DCCA coefficient:For , the standard DCCA coefficient is retrieved.Step 6: compute the multifractal detrended partial cross-correlation coefficient between and by eliminating the influence of the controlling variable on and analogous to the generalization of the correlation coefficient to partial correlation coefficient:named the first-order multifractal detrended partial cross-correlation coefficient (first-order MFDPCC coefficient), where are random variables, is the controlling variable, and represent the mean of MFDCCA coefficients for and , and , and and , respectively.

For , the first-order detrended partial cross-correlation coefficient (first-order DPCC coefficient) is retrieved.

2.2. The nth-Order Multifractal Detrended Partial Cross-Correlation Analysis and n-Controlling-Variables Detrended Partial Cross-Correlation Coefficient

Considering the cross-correlation between and affected by more than one variable in complex systems, we define the second-order multifractal detrended partial cross-correlation coefficient (second-order MFDPCC coefficient) by using the partial correlation method [36]:where are random variables, controlling variables are not related to each other, and , and are first-order MFDPCC coefficients.

Generally, the nth-order multifractal detrended partial cross-correlation coefficient (nth-order MFDPCC coefficient) is as follows:where , , and are (n − 1) th-order MFDPCC coefficients and controlling variables are not related to each other. For , the nth-order detrended partial cross-correlation coefficient (nth-order DPCC coefficient) is retrieved.

In general, the nth-order partial cross-correlation is necessary when these controlling variables are not related to each other. Nevertheless, in real-world systems, the variables generated by large number of interacting units are cross-correlated. Therefore, we define the n-controlling-variables multifractal detrended partial cross-correlation coefficient (n-variables MFDPCC) by equation (12) for related controlling variables . Note that when the controlling variables are not related to each other, the n-variables MFDPCC is equivalent to the nth-order MFDPCC.

3. Data and Analysis

3.1. Two-Component ARFIMA Process

In order to test the robustness of the proposed n-controlling-variables MFDPCC coefficient method, power-law cross-correlated time series and are generated by using the two-component ARFIMA stochastic process in this section [18, 39, 40]. In this model, the series is defined bywhere is weight, is a free parameter to control the coupling strength between and (0.5 ≤ ≤ 1), and and are independent and identically distributed (i.i.d.) Gaussian variables with and [18, 39]. For different values of , the different coupling strength between the variables and is . In this section, the two-component ARFIMA series and with parameter and  = 0.5, denoted by and , are employed to detect the interactions between two time series. Then, the effect of white noise sequence on the cross-correlation of the two series and is tested to investigate the validity of the n-controlling-variables MFDPCC coefficient analysis mentioned in this paper. For this purpose, we study the difference between the mean of the MFDCCA coefficient and the n-controlling-variables MFDPCC coefficient for any parameter q by using the influence degree function . The influence degree function is defined as

We calculate the influence degree function of the synthetical signals using the proposed first-order MFDPCC coefficient and present the influence degree function vs. parameter q in Figure 1. The results of the influence degree values of different q are just about nil, which indicates that there is hardly any effect of white noise sequence on cross-correlation of the two series and .

3.2. Stock Market

To further exemplify the potential utility of the n-controlling-variables MFDPCC coefficient method for analyzing real-world data, we study daily closing prices of fifteen stock markets including the São Paulo Index (IBOV), the Dow Jones Index (DJI), the NASDAQ Index (IXIC), the Standard & Poor 500 Composite Stock Price Index (SPX), the FTSE Global Equity Index Series (FISE), the French CAC 40 (FCHI), German DAX Index (GDAXI), Nikkei 255 Index (N255), Korea Composite Index (KS11), Hang Seng Index (HSI), Australian Standard & Poor’s 200 (AS51), Mumbai Index (SENSEX), Russian Index (RTS), Shanghai Composite Index (SSEC), and Shenzhen Composite Index (SZI). Datasets are from January 04, 1993, to January 03, 2019.

Figure 2 shows the mean of DCCA coefficients for the stock series. The mean of DCCA coefficient between DJI and SPX is 0.97, which performs relevantly different from other DCCA coefficients. It indicates the close cross-correlation between the American stock markets. The next largest DCCA coefficient is obtained by SSEC and SZI, which indicates the close cross-correlation in Chinese mainland stock markets.

The mean of DCCA coefficients between SZI and stock markets in developed countries (GDAXI, N225, KS11, and AS51) is less than 0.3. It shows that SZI has a weak relationship with stock markets in developed countries. The mean of DCCA coefficients between SZI and HSI is in an intermediate state, which indicates the existence of cross-correlation in Chinese stock markets.

Next, we analyze the effect of the other thirteen stock markets on cross-correlation characteristics between SSEC and SZI, by applying the influence degree of first-order DPCC coefficient. For the effect on cross-correlation characteristics between the SSEC and SZI, the largest influence degree I = 0.05 is obtained by HSI, which shows the information exchange between the Chinese stock market, as seen in Figure 3. The next largest I = 0.04 is acquired by SENSEX, which indicates the association between the stock markets in developing countries (Indian and Chinese stock markets). The I values of other stock time series are less than 0.1, which indicates little information exchange between the Chinese mainland stock market and other stock markets. The influence degree values of 13 stock markets for first-order MFDPCC coefficient with are also demonstrated in the upper left of Figure 3.

During the analysis, we observe the effect of HSI on cross-correlation characteristics between SSEC and SZI from influence degree function that decreases as the scale q increases. And this infers the change of multifractal cross-correlation.

In order to capture the change of multifractal cross-correlation between two nonstationary time series influenced by common external forces, multifractal detrended partial cross-correlation analysis (MFDPXA) is employed [38]. We also investigate the multifractal behavior between the bivariate time series through MFDPXA method for comparison. The result shows that both the corresponding spectra and are wide, but the latter is narrower than the former, which is presented in Figure 4.

Here, we perform cross-correlation analysis using MFDPXA method and give the multifractal spectrum for SZI and SSEC time series in which HSI shows significant influence on multifractal spectrum, as seen in Figure 4. We compare the obtained influence degrees with the aforementioned method and infer that the HSI has significant influence on SZI and SSEC time series. These similar results imply that the partial cross-correlation method is quite efficient in eliminating external common influence factor.

Applied to scalar variables, the first-order MFDPCC will detect the intrinsic interactions by removing the correlations of controlling variables. When variables are time series, this application is equivalent to removal of zero delay correlations, whereas delayed correlations are not considered [36, 37, 40, 41]. Therefore, we investigate the delayed effect of variable on the correlation between variables and . Because the two variables and in question may themselves be correlated at nonzero delays, we write the multifractal detrended partial cross-correlation between and , given , as a function of two time delays:where is the delay between variables and and is the delay between variables and .

In this section, we estimate the delayed effect of HSI on the correlation between SSEC and SZI by using the time delay influence degree . Figure 5 shows the time delay influence degree for . The effect of on influence degree is weaker than that of on influence degree.

We now analyze the 2-controlling-variables effect of the other thirteen stock markets on cross-correlation characteristics between SSEC and SZI, by giving a set of two controlling variables. In Figure 6, we illustrate the comparative relation of the influence degree of 2-controlling-variables DPCC coefficients for elements by the matrix diagram.

We note that the structure of the matrix is symmetrical and that element at the intersection of row i and column j represents the influence of controlling variables on the cross-correlation of SSEC and SZI, where the 2-controlling-variables are the stock time series from IBOV, DJI, IXIC, SPX, FISE, FCHI, GDAXI, N255, KS11, HSI, AS51, SENSEX, and RTS. Therefore, we analyze the top left corner of the matrix. It can be seen that the largest element is the intersection of row 2 and column 10, i.e., SENSEX and HSI, which indicates the association between the Indian and Chinese stock markets. This is consistent with our result of first-order MFDPCC coefficient method.

Concerning the influence degree of the 2-controlling-variables MFDPCC, we demonstrate 5 cases (HSI and SENSEX, HSI and RTS, SENSEX and KS11, FCHI and N255, and IXIC and FISE) for in Figure 7. The largest influence degree is the case SENSEX and HSI, which is consistent with the 2-controlling-variables DPCC method.

3.3. Aeroengine Time Series

Previous research studies show that the aeroengine gas path parameters such as low-pressure rotor speed (N1), high-pressure rotor speed (N2), and fuel flow (WF) play an important role in understanding the aeroengine system [21, 42]. The mean of DCCA coefficients for the aeroengine time series is shown in Figure 8, where the average DCCA coefficient between N1 and N2 is 0.85, which shows the close cross-correlation between N1 and N2.

We here investigate the partial correlation between N1 and N2 given a set of eight controlling variables, including WF, exhaust gas temperature (EGT), N2 tracked vibration channel B (N2TB), inlet air pressure (P2), outlet temperature of high-pressure compressor (T3), outlet temperature of low-pressure compressor (T2.5), and other temperatures (T2 and T2.95).

In Figure 9, we plot the influence degree of first-order DPCC coefficient, investigating the effect of the other eight controlling variables on cross-correlation characteristics between N1 and N2. The largest influence degree I = 0.51, obtained by T3, shows the information exchange between the outlet temperature of high-pressure compressor and the rotor speed system. The next largest I = 0.22 is acquired by WF, which indicates the association between the fuel flow system and rotor speed system.

The result of the influence degree for eight aeroengine parameters applying by first-order MFDPCC coefficient with is also demonstrated in the upper left of Figure 9. The effect of T3 on cross-correlation characteristics, observed from influence degree function , decreases as the scale q increases. It indicates that the multifractal cross-correlation differs across values of q.

Further, we apply the MFDPXA method on aforementioned N1 and N2 time series considering the T3 as common influencing factor. It is observed from Figure 10 that the corresponding spectra and are wide which shows the strength of multifractal behavior in analyzed time series. We observe that the width of singularity spectrum is narrower, and this implies the strength of multifractal nature is weak in analyzed bivariate time series.

Here, we estimate the delayed effect of T3 on the correlation between N1 and N2 by using the time delay influence degree . Figure 11 shows the time delay influence degree for . It is obvious that the time delay influence degree gradually increases and then declines as a single-peak curve when remains constant. As increases, the peak value of time delay influence degree shifts rightward.

The next observation concerns the influence degree of 2-controlling-variables DPCC coefficient in the aeroengine system. We now analyze the influence of two controlling parameters on the cross-correlation between N1 and N2. In Figure 12, we illustrate the comparative relation of the influence degree of 2-controlling-variables DPCC coefficient for aeroengine system. It can be seen that the larger elements in the symmetrical matrix are located at row 3 or column 6, which denote T3 has a greater impact on the correlation between N1 and N2.

Concerning the influence degree of the 2-controlling-variables MFDPCC, we demonstrate 7 cases (T3 and WF, T3 and N2TB, T3 and T2, T3 and T2.5, T2.95 and N2TB, T2.95 and P2, and T2.95 and T2) for in Figure 13. Larger influence degrees exist in the cases with the presence of T3 (T3 and WF, T3 and N2TB, T3 and T2, and T3 and T2.5), which is consistent with the 2-controlling-variables DPCC method, as seen in Figure 12.

For the aeroengine, the parameters N1 and N2 are chosen to indicate the engine thrust which depends on the throttle lever angle. Hence, the cross-correlation between them is strong. The temperature and pressure parameters are linked with many factors, including the compressor power, combustion efficiency, throttle lever angle, etc. Therefore, the dynamic interaction of these three groups makes the aeroengine function. These results estimate the influence of temperature and pressure parameters on the cross-correlation between N1 and N2.

4. Conclusion

In this paper, we propose the nth-order multifractal detrended partial cross-correlation analysis method and the n-controlling-variables multifractal detrended partial cross-correlation analysis method for understanding the interactions between two nonstationary time series. For comparing these new methods with classical measures, we introduce the influence degree function. We then apply the n-controlling-variables multifractal detrended partial cross-correlation analysis of stock markets and aeroengine performance parameters and measure the influence degree function of the partial cross-correlation in a dynamic system.

To understand the numerous real-world systems where the output signals exhibit complex cross-correlation, both cross-correlation and partial correlation are subjects of investigation. The information of n-variables MFDPCC helps people to research information exchange in complex systems. This paper gives two examples, stock markets and aeroengine systems. For stock time series, our results indicate that, concerning closing index values, there is little information exchange between the Chinese stock markets and the American-European stock markets, whereas the SSEC, SZI, and HSI, by first-order MFDPCC method and 2-controlling-variables MFDPCC, show frequent and abundant information exchange in Chinese stock markets. For aeroengine performance parameters, our results show that there is some information exchange between the engine rotor system and the aeroengine parameters, such as the outlet temperature of the high-pressure compressor and the fuel flow.

We believe that the MFDPCC method can be used to detect the intrinsic interactions among multiple dynamical systems, and therefore it can be widely applied to many research fields such as the aeroengine health monitoring systems and the investment portfolio where the covariance is employed to explore the interaction of assets income.

The multifractal detrended partial cross-correlation analysis is used to delete the possible indirect correlation, but it may also delete valuable information. This problem required further investigation, both experimental and theoretical. Hence, the results of this paper should be considered as preliminary results on the multifractal detrended partial cross-correlation analysis. Therefore, we hope that this study will be extended to analyze the filtered information.

Data Availability

The stock market data used to support the findings of this study are available from the corresponding author upon request. The aeroengine data used to support the findings of this study have not been made available because of commercial secrets.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The financial support from the funds of the MOE (Ministry of Education in China) Project of Humanities and Social Sciences under grant no. 19YJC910001 and the Fundamental Research Funds for the Central Universities under grant no. 3122014K013 is gratefully acknowledged.