Complexity
Volume 2018 (2018), Article ID 7015721, 13 pages
https://doi.org/10.1155/2018/7015721

Research Article

Dynamical Variety of Shapes in Financial Multifractality

Stanisław Drożdż,1,2 Rafał Kowalski,¹ Paweł Oświȩcimka,¹ Rafał Rak,1,3 and Robert Gȩbarowski²

¹Complex Systems Theory Department, Institute of Nuclear Physics, Polish Academy of Sciences, ul. Radzikowskiego 152, 31-342 Kraków, Poland
²Faculty of Physics, Mathematics and Computer Science, Cracow University of Technology, ul. Warszawska 24, 31-155 Kraków, Poland
³Faculty of Mathematics and Natural Sciences, University of Rzeszów, ul. Pigonia 1, 35-310 Rzeszów, Poland

Correspondence should be addressed to Stanisław Drożdż; stanislaw.drozdz@ifj.edu.pl

Received 16 March 2018; Revised 5 July 2018; Accepted 17 July 2018; Published 16 September 2018

Academic Editor: Lingzhong Guo

Copyright © 2018 Stanisław Drożdż et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The concept of multifractality offers a powerful formal tool to filter out a multitude of the most relevant characteristics of complex time series. The related studies thus far presented in the scientific literature typically limit themselves to evaluation of whether a time series is multifractal, and width of the resulting singularity spectrum is considered a measure of the degree of complexity involved. However, the character of the complexity of time series generated by the natural processes usually appears much more intricate than such a bare statement can reflect. As an example, based on the long-term records of the S&P500 and NASDAQ—the two world-leading stock market indices—the present study shows that they indeed develop the multifractal features, but these features evolve through a variety of shapes, most often strongly asymmetric, whose changes typically are correlated with the historically most significant events experienced by the world economy. Relating at the same time the index multifractal singularity spectra to those of the component stocks that form this index reflects the varying degree of correlations involved among the stocks.

1. Introduction

Multifractality is a concept that is central to the science of complexity. The related multiscale approach [1–3] aims at bridging the wide range of time and length scales that are inherent in a number of complex natural phenomena, and as such, it pervades essentially all scientific disciplines [4]. By now, it finds applications in essentially all areas of the scientific activity, including physics [5, 6], biology [7–9], chemistry [10, 11], geophysics [12, 13], hydrology [14], atmospheric physics [15], quantitative linguistics [16, 17], behavioural sciences [18], cognitive structures [19], music [20, 21], songbird rhythms [22], physiology [23, 25], human behaviour [24, 26, 27], social psychology [28], and even ecological sciences [29], but especially frequently in economic and in financial contexts [30–41] as stimulated by practical aspects and by needs to develop models of the financial dynamics based on multifractality [31, 42–45] such that they help in making predictions. Indeed, the multifractal analyses of the financial time series have provided so far most of the quantitative evidence for the factors that induce the genuine multifractality, such as the temporal long-range nonlinear correlations and, only when such correlations are present, the fat tails in the distribution of fluctuations [46]. In order to unambiguously identify action of such factors and to suppress potential spurious multifractality, the time series under study have to be, however, sufficiently long [47]. In addition, the realistic time series, as generated by the natural phenomena, even if of multifractal character, are typically more involved in composition than the model mathematical uniform multifractals, and they may contain several components of different multifractality characteristics. In such frequent cases, the global hierarchical organization of the series gets distorted, and the multifractal spectrum becomes asymmetric, either left- or right-sided, as recently demonstrated in ref. [48]. Detecting such effects may provide even more valuable information about the mechanism that governs dynamics of a particular time series than just a bare statement that it is multifractal. Such effects of asymmetry are, for instance, already found to constitute a very helpful formal tool in identifying a specific organization of complex networks [49]. Furthermore, directions of the relevant distortions may vary in time parallel to changes of weight of the constituent components in a series. The most straightforward candidate to experience this kind of impact is the stock market index which, by construction, is already a sum, most often weighted, of prices of the constituent companies, and those companies themselves may react differently for the same external news depending of the sector they belong to. It is primarily for this reason that below, the world largest stock market indices are studied. Of course another, more specific, market-oriented reason for this study is to broaden our historical perspective on evolution of the stock market multiscale characteristics over periods comprising the global crashes or transitions due to the technological revolution in trading.

2. Multifractal Formalism

At present, there exist two distinct, commonly accepted, and complementary computational methods that serve quantification of the multifractal characteristics of the time series. One of them—the wavelet transform modulus maxima (WTMM) [3]—makes use of the wavelet expansion of the time series under consideration, and the other one—the multifractal detrended fluctuation analysis (MFDFA) [50]—is based on inspecting the scaling properties of the varying order moments of fluctuations evaluated after an appropriate trend removal. While the former of those techniques allows a better visualization of the underlying patterns in the time series, the latter one often appears more accurate and more stable numerically, and it will therefore be used here. Furthermore, at present, there exists a consistent generalization of MFDFA such that it even allows to properly identify and quantify the multifractal aspects of cross-correlations between two time series [51–53]. This novel method, termed multifractal cross-correlation analysis (MFCCA), consists of several steps that at the beginning are common to all the methods based on detrending.

One thus considers two time series and , where . The signal profile is then calculated for each of them:

where denotes averaging over the entire time series. Next, both these signal profiles are split up into () disjoint segments of length starting both from the beginning and the end of the profile, and in each , the assumed trend is estimated by fitting a polynomial of order ( for and for ). In typical cases, an optimal choice corresponds to [54]. This trend is subtracted from the series, and the detrended cross-covariance within each segment is calculated:

Since can assume both positive and negative values, the th-order covariance function is defined by the following equation:

where denotes the sign of . The parameter in (3) can take any real number except zero. However, for , the logarithmic version of this equation can be employed [50]:

Fractal cross-dependencies between the time series and manifest themselves in the scaling relations:

or for , where is the corresponding scaling exponent whose range of dependence on quantifies the degree of the complexity involved. Scaling with the -dependent exponents reflects a richer, multifractal character of correlations in the time series as compared to the monofractal case when is -independent.

The conventional MFDFA procedure of calculating the singularity spectra for single time series can be considered a special case of the above MFCCA procedure and corresponds to taking and as identical. Equation (3) then reduces to

and to a corresponding counterpart of (4) for . The signatures of multifractality (monofractality) are then reflected, analogously to (5), by

where denotes the generalized Hurst exponent. The singularity spectrum (also referred to as multifractal spectrum) is then calculated from the following relations:

where denotes the Hölder exponent characterizing the singularity strength and reflects the fractal dimension of support of the set of data points whose Hölder exponent equals . In the case of multifractals, the shape of the singularity spectrum typically resembles an inverted parabola and the degree of their complexity is straightforwardly quantified by the width of :

where and correspond to the opposite ends of the values as projected out by different -moments (see (6)). For monofractal signals, the spectrum converges to a single point, though in practice, this often turns out to be a subtle matter [47]. Another important feature of the multifractal spectrum is its asymmetry (skewness), which can be quantified by the asymmetry coefficient [48]:

where and , and for , the spectrum assumes maximum. The positive value of reflects the left-sided asymmetry of ; that is, its left arm is stretched with respect to the right one and, thus, more developed multifractality on the level of large fluctuations in the time series. Negative , on the other hand, reflects the right-sided asymmetry of the spectrum and indicates temporal organization of the smaller fluctuations as the main source of multifractality.

A family of the fluctuation functions as defined by (3) can also be used to define a -dependent detrended cross-correlation (qDCCA) [55] coefficient

which allows to quantify the degree of cross-correlations between two time series and after detrending and at varying time scales . Furthermore, by varying the parameter , one is able to identify the range of detrended fluctuation amplitudes that are correlated the most in the two signals under study [55]. This filtering ability of constitutes an important advantage as cross-correlations among time series typically are not uniformly distributed over their fluctuations of different magnitude [56].

3. Data Specification

In the present study, two sets of data are used:

(i)Daily prices of the S&P500 and NASDAQ indices covering the period January 03, 1950–December 29, 2016 (16,496 data points). The values of the NASDAQ before 1971 (official launching date of the index is February 05, 1971) were reconstructed from the historical data [57].(ii)Daily prices of 9 stocks listed on the NYSE over the period from January 1, 1962, to July 07, 2017 (13,812 points). The analysed companies are GE (General Electric), AA (Alcoa), IBM (International Business Machines), KO (Coca-Cola), BA (Boeing), CAT (Caterpillar), DIS (Walt Disney), HPQ (Hewlett-Packard), and DD (DuPont). These in fact are the only stocks that participate in the Dow Jones Industrial Average (DJIA) over such a long period of time and thus also in the S&P500. They, however, represent a large spectrum of the economy sectors and may thus be considered as a reasonable representation for the larger American indices.

For each time series, the logarithmic returns are calculated according to the equation:

where denotes the stock price or index value and stands for time interval (). All time series are normalized to have unit variance and zero mean.

4. Results

4.1. S&P500 and NASDAQ

The MFDFA multifractal spectra for the S&P500 and NASDAQ indices are shown in Figure 1. For both these indices, the fluctuation functions reveal a convincing power-law behaviour over almost two decades, which is shown in the corresponding lower-right insets; thus, is determined unambiguously. The parameters are taken within the interval , which is common in financial applications because it allows to safely avoid the danger of divergent moments when the fluctuation functions are computed. Cumulative distributions of the return fluctuations for the two indices considered here are shown in the corresponding upper-left panels of Figure 1 and can be seen not to develop thicker tails than the inverse cubic power law [58, 59], and there is thus no danger of the divergent moments. The width of the resulting spectra for the S&P500 and 0.32 for the NASDAQ, correspondingly. The significance of this result is also tested against the two null hypotheses of calculated from (i) the series obtained from the original ones by a random shuffling, thus destroying all the temporal correlations (green triangles), and (ii) Fourier-phase randomized counterparts of the original series which destroys the nonlinear correlations (blue squares). Clearly, the spectra in these two tests get shrank to a form characteristic to monofractals. An additional form of surrogates tested here are time series with the Gaussianized pdfs. In the latter case, the original pdf is replaced by a Gaussian distribution while the amplitude ranks of fluctuations remain preserved. The resulting multifractal spectra appear only slightly narrower than the original ones, and therefore, they are not shown in Figure 1. All these tests thus provide a convincing evidence for quite a rich multifractality of the original time series and, moreover, corroborate the fact that this multifractality is, as expected [47], due to the nonlinear temporal correlations. The obtained multifractal spectra are at the same time visibly left-sided asymmetric [48]. The asymmetry coefficient for the S&P500 and 0.31 for the NASDAQ. The left side of is determined by the positive values which filter out larger events, and the opposite applies to the right side of this spectrum. In the present context, this thus means that it is the dynamics of the large returns which develop more pronounced multifractal organization than those of the small returns.

Figure 1: Main panels: multifractal spectra calculated for the S&P500 and NASDAQ returns (black dots) covering the period January 03, 1950–December 29, 2016. Average spectra obtained for the Fourier phase-randomized surrogates and for the randomly shuffled time series are denoted by blue squares and green triangles, respectively. Upper-left insets display cumulative distributions of return fluctuations, and lower-right insets display the fluctuation functions calculated for the original S&P500 and NASDAQ series.

Figure 1 shows the result of calculations over the entire time span where the time series are taken. It appears that probing this period with a shorter window rolling in time reveals a nontrivial and a very interesting time dependence of the corresponding multifractal spectra. Here, the window size is taken over 5000 data points (equivalent to about 20 years) which in the presence of temporal correlations is sufficiently long to guaranty stability of the result [47] (absence of such correlations demands significantly longer series [47]), and the window is moved with the step of 20 points (approximately one calendar month). The results of such a procedure are highlighted in Figure 2 for the S&P500 and in Figure 3 for the NASDAQ. Panels (a) in these figures show sequences of the singularity spectra , calculated within such windows consecutively, and the calendar date assigned to each corresponds to the end point within a window. Thus, for the time series which begin, as here, in January 1950, the first date appearing in Figures 2 and 3 corresponds to January 1969. In order to better visualize evolution of and of , the panels (b) in these figures show projections of onto the time plane. The three historically most recognized events that influenced the world financial markets are indicated by the vertical dashed lines. These are the Black Monday of October 19, 1987; burst of the Dot-com bubble in March 10, 2000; and bankruptcy of Lehman Brothers in September 15, 2008. Clearly, evolution of in such a 20-year time window reveals sizeable changes in the width of and in its asymmetry, both going somewhat differently in the two indices, however. For S&P500 until about 1985, the spectrum is comparatively broad and then starts quick narrowing, but this narrowing primarily results from shrinkage of the right arm in . For the time window ending in around 1993, this arm almost disappears and starts recovering only in recent years. Interestingly, the left side of got broadened even a few years earlier. The NASDAQ spectrum also experiences sizeable changes in time but differently than the one for the S&P500. On average, this spectrum is broader and strongly asymmetric for time windows ending between about Black Monday and the burst of the Dot-com bubble in 2000, but here, this asymmetry results from a sudden stretching of the left side of while the right side does not experience much changes.

(a)

(b)

Figure 2: (a) For the S&P500 from January 03, 1950, to December 29, 2016, the sequence of singularity spectra calculated within a rolling 20-year window. The calendar date assigned to each corresponds to the end point within a window. This window is moved with the step of 20 points which corresponds approximately to one calendar month. The black sold line corresponds to the left and the blue line to the right side of . (b) Projections of of (a) onto the time plane. The red line illustrates displacement of the maxima of in the consecutive windows. Vertical dashed lines indicate the Black Monday of October 19, 1987; burst of the Dot-com bubble in March 10, 2000; and bankruptcy of Lehman Brothers in September 15, 2008.

(a)

(b)

Figure 3: The same as in Figure 2, but analysis is carried out for the NASDAQ data.

In Figures 2 and 3, one also sees changes in location of the maxima of which are related to a degree of persistence in time series. A parameter that directly quantifies this property is the Hurst exponent . Figure 4 shows the Hurst exponents H, the widths , the asymmetry coefficients , and the widths , for the time sequence of the multifractal spectra already presented in Figure 2 (S&P500), whereas Figure 5 shows these characteristics corresponding to Figure 3 (NASDAQ). The two dates seen to be related to almost discontinuous changes in some of these quantities are the Black Monday of October 19, 1987, and the bankruptcy of Lehman Brothers in September 15, 2008, and these two dates are indicated by the vertical dashed lines. While Black Monday affected the NASDAQ much more spectacularly than the S&P500, though the latter started assuming similar trends already some 2 years earlier, the effect of the bankruptcy of Lehman Brothers was just the opposite. This time, it is the S&P500 which reveals a sudden increase of by a factor of about 2, but remarkably, this increase is entirely due to stretching of the left arm of . A partial identification of the origin of these S&P500 versus NASDAQ differences comes from Figure 6, which displays fluctuations of the daily returns of these two indices and, as the most informative, the time dependence of the local (in the rolling window of trading days) detrended variance. In around the Black Monday, this variance is much larger for the NASDAQ than for the S&P500, and this goes in parallel with a sharp stretching of the left arm in for the NASDAQ. On the other hand, in around the bankruptcy of Lehman Brothers, even though the NASDAQ detrended variance still is somewhat larger than the one of the S&P500, it is much smaller than around the period of the Dot-com burst. In the S&P500 case, the corresponding development is just reversed, and larger variance accompanies the bankruptcy of Lehman Brothers. Thus, in this latter period, the detrended variance of the NASDAQ decreases while the one of the S&P500 increases, and it is in this period when the left arm of for the S&P500 experiences a sudden stretching. Worth noticing is also the fact that the Hurst exponents of these two indices on average decrease when going from past to present and in recent years assume values even lower than 0.5, which indicates antipersistence [60]. Especially monotonic in this respect is the S&P500—one of the most significant global indicators of the world economy—whose Hurst exponent on average systematically decreases over the whole time span considered, and in the last couple of years, it even steadily dropped down below 0.4. There are presumptions [60, 61] that such values of indicate proximity to a crash zone. In view of this result, the log-periodic scenario [62] indicating danger of a much larger world economic decline in around 2025 than anything the World has experienced so far needs to be taken into consideration more and more seriously.

Figure 4: For the S&P500, the Hurst exponents , the widths , the asymmetry coefficients , and the widths for the time sequence of multifractal spectra in the 20-year windows of Figure 2.

Figure 5: For the NASDAQ, the Hurst exponents , the widths , the asymmetry coefficients , and the widths for the time sequence of multifractal spectra in the 20-year windows of Figure 3.

Figure 6: Two upper panels: daily returns for the S&P500 and for the NASDAQ over the period January 03, 1950–December 29, 2016. The bottom panel: the corresponding detrended variance for the S&P500 (red line) and for the NASDAQ (black line).

The window-probed multifractal spectra of Figures 2 and 3 for the S&P500 and for the NASDAQ resemble each other more in the first half of the entire considered interval, until about the mid-1980s, than in the following second half. This similarity or dissimilarity appears to occur even on the deeper level of their multifractal synchrony as reflected by the appropriate cross-correlation measures expressed by (5). The two approximately 20-year-long time periods taken from inside of these halves are selected as September 25, 1957–August 26, 1977, and May 19, 1989–March 20, 2008; the cross-correlation fluctuation functions between the S&P500 and the NASDAQ calculated according to (5), and the result is shown in Figure 7. It is very interesting to see that in the first of these periods, the fluctuation functions display a clear tendency to scaling, which indicates cross-correlations between the two indices even on the level of their multifractal organization. This holds down to the level of their small fluctuations as measured by the negative values. In the second of these time intervals, while for the positive values one may still see some remnants of scaling, for the negative values there is none; thus, the indices are systematically losing their multifractal synchrony, and on the level of the small fluctuations, this synchrony is lost completely.

Figure 7: Cross-correlation fluctuation functions between the S&P500 and NASDAQ calculated according to (5) for in two periods: September 25, 1957–August 26, 1977, and May 19, 1989–March 20, 2009.

4.2. Index versus Companies

It is natural to expect that significant changes in time of the multifractal features of the two indices seen in the previous subsection reflect different market phases, and such phases vary in a degree of coupling among the component shares [63]. These are the individual stocks which are traded, and only a superposition of their multifractal characteristics, not necessarily identical, determines of an index. It is clear that in an uncorrelated sum of many multifractal time series, the multifractality gradually disappears when the number of component series increases, and in addition, this limiting case is typically approached asymmetrically [48]. One may thus anticipate that stronger coupling among the companies that form a basket of an index favours multifractality of that index as well. In the present context, in order to study such effects in more detail, by summing up prices of the 9 companies listed in Section 3, a proxy of the DJIA is formed. It, however, amazingly accurately follows changes in the full DJIA and even all the significant moves in the S&P500, as can be seen from Figure 8. This is likely to reflect the fact that the 9 companies are dispersed over different market sectors, and in total, they well represent the global DJIA market.

Figure 8: Daily prices of the S&P500, of the Dow Jones, and of the sum of 9 DJIA stocks listed on the NYSE over the period from January 1, 1962, to July 07, 2017 (13,812 points). The companies included are GE (General Electric), AA (Alcoa), IBM (International Business Machines), KO (Coca-Cola), BA (Boeing), CAT (Caterpillar), DIS (Walt Disney), HPQ (Hewlett-Packard), and DD (DuPont).

The results of calculations relating to the multifractal spectra , projected onto the time plane, of these companies labelled by i (thus, here ), for illustrative clarity represented by one average and of the index constructed from these 9 companies, in the same rolling window as before, are displayed in panels (a) and (b) of Figure 9, correspondingly. Several interesting observations based on these results can be made. One main finding is that the width of is never smaller than that of the global 9-company index, which is understandable because equality is expected in the case of perfect correlation among prices of all the participating companies. Some decorrelation, which is always the case in real markets, should result in narrowing of the global, here 9 companies, index. A significantly larger difference between the widths of multifractal spectra in the two cases considered is observed for the time period between the Black Monday and the bankruptcy of Lehman Brothers, and the transition is nearly sharp. This difference originates, however, from a sudden stretching of the left side in within that period, which indicates that multifractality of the price changes of individual companies is much more pronounced on the level of larger fluctuations than on the level of small ones. When prices of these companies are summed up to form a global 9-company index, this huge left-side stretching is significantly reduced, which indicates that the large fluctuations of individual stocks are not fully correlated among themselves. Still, within this most volatile period (Figure 6) in the market, even the global index preserves the left-sided asymmetry in indicating dominance of nonlinear correlations on the level of large fluctuations.

Figure 9: Projections onto the time plane of the sequence of singularity spectra calculated within a rolling 20-year window for both the average spectrum (a) and (b) the artificial index of the 9 companies of Figure 8. The calendar date assigned to each corresponds to the ending point within a window. This window is moved with the step of 20 points which corresponds to approximately one calendar month. The red line illustrates displacement of the maxima of in the consecutive windows. Vertical dashed lines indicate the Black Monday of October 19, 1987, and bankruptcy of the Lehman Brothers in September 15, 2008, while the dotted ones indicate October 1990 and April 1994. The bottom three panels display the corresponding Hurst exponents (c1), widths (c2), and the asymmetry coefficients (c3).

An especially interesting related case occurs in the period between October 1990 and April 1994 indicated in Figure 9 by the two vertical dotted lines. In this period, the multifractal spectra of the individual companies on average develop broad multifractal spectra while of the corresponding global 9-company index is so narrow that it can be considered as monofractal. One possible reason for such a result is a substantial suppression of cross-correlations among price changes of the component stocks [48]. Such a possibility is verified using the correlation matrix

where denotes a rectangular matrix formed from time series of length . Entries of the matrix thus correspond to the conventional Pearson correlation coefficients. By diagonalizing (), one obtains the eigenvalues and the corresponding eigenvectors . In the limiting case of entirely random signals, the density of eigenvalues is known analytically [64, 65] as

where the lower and upper bounds of this distribution are given by

In this expression, and is equal to the variance of the time series. The degree of departure of the largest eigenvalue above is a measure of the strength of correlations among the time series participating [4, 66].

Changes of the magnitude of the largest eigenvalue in the rolling time window of length trading days for the present versus the noise regime as set by and for these particular values of and are shown in Figure 10. Furthermore, in the same figure, changes of the largest eigenvalue of an analogous matrix composed of the coefficients as defined by (11) taking for are also shown. Clearly, in both these measures, the largest eigenvalues assume the lowest values in the period of interest, just between October 1990 and April 1994. At one point, the value even touches the border of purely random series. Thus, the scenario of the least correlated 9 companies here studied in this time period applies, indeed, which explains a narrow of the global 9-company index.

Figure 10: The blue line displays the time dependence of the largest eigenvalue of the correlation matrix constructed from the time series representing daily returns of the 9 companies of Figure 9 in a rolling window of size trading days. The red line displays the largest eigenvalues of an analogous rolling window matrix composed of the coefficients as defined by (11) taking .

5. Conclusions

Quantification of the complex time series in terms of multifractality nowadays finds a multitude of applications in diverse areas. Thus far, however, majority of the related studies presented in the scientific literature limit themselves to a sole estimation of the singularity spectrum, and if found multifractal, it usually is treated as evidence of the hierarchical organization of such series, and the width of such a spectrum is considered a measure of the degree of complexity involved. While this indicates some kind of a cascade-like, hierarchical organization indeed, in realistic cases, such an organization is rarely uniform. The time series generated by natural processes may include many convoluted components with different hierarchy generators each, which results in asymmetry of the singularity spectra. Even more, contribution of such components may vary in time, and this thus may introduce further dynamical variability. Definitely, the financial markets constantly functioning in evolving external conditions represent a natural candidate to become a subject of such effects. This can be anticipated to apply almost straightforwardly to the stock market indices as they by construction constitute an average (typically weighted but not always) of the prices of selected stocks representing different economy sectors, thus not necessarily obeying the same multiscaling characteristics. The degree of correlations among such stocks is also known to depend on the global market phases. In the present paper, based on over half a century daily recordings of the S&P500 and NASDAQ, the two world-leading stock market indices, it is shown that they reveal the multiscaling features which expressed in terms of the multifractal spectrum evolve through a variety of shapes whose changes typically appear correlated with the historically most significant events experienced by the world economy. From a more general perspective, these results indicate that the form of the multifractal spectrum, and especially its departures from the model mathematical cases of the uniform cascades, contains richness of information that, if properly interpreted and potentially disentangled, may provide very valuable insight into the underlying dynamics which may be of crucial value for a more accurate modelling of the financial markets. Taking into consideration the effects exposed here may also be very helpful for market regulators and policymakers in stabilizing markets as well as for a flexible portfolio optimization.

Finally, the methodology introduced in Section 4.2 of relating the global (here, index) multifractal spectrum to the corresponding multifractal spectra of subsystems (here, companies) provides an appropriate quantitative tool with potential applications extending far beyond the financial context when various questions related to the so-called complexity matching [67] are addressed and studied empirically as, for instance, those in a psychological/cognitive domain [24, 68–71]. Differences between widths—as an example in Figure 9 shows—of such spectra reflect strength of the underlying complexity matching between subsystems, and this strength may vary in time. The weakest matching, for instance, corresponds to the period between October 1990 and April 1994. Furthermore, appreciating the relative changes in asymmetry of may allow to selectively scan the varying strengths of such a matching for different ranges of fluctuations. Of course, as far as the world financial markets are concerned, one may rely on observations only since, by their very nature, there exists no realistic possibility to set up the world financial experiments. Since phenomena belonging to the domain of social psychology definitely constitute a significant factor driving the markets, a properly coordinated joint multidisciplinary effort may crucially help in understanding the cross-scale dependencies and information flows in the financial markets and in other complex systems as well.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was supported in part by PLGrid Infrastructure.

References

T. C. Halsey, M. H. Jensen, L. P. Kadanoff, I. Procaccia, and B. I. Shraiman, “Fractal measures and their singularities—the characterization of strange sets,” Physical Review A, vol. 33, no. 2, pp. 1141–1151, 1986.

B. B. Mandelbrot, “Multifractal measures, especially for the geophysicist,” Pure and Applied Geophysics, vol. 131, no. 1-2, pp. 5–42, 1989.

J.-F. Muzy, E. Bacry, and A. Arneodo, “The multifractal formalism revisited with wavelets,” International Journal of Bifurcation and Chaos, vol. 4, no. 2, pp. 245–302, 1994.

J. Kwapień and S. Drożdż, “Physical approach to complex systems,” Physics Reports, vol. 515, no. 3-4, pp. 115–226, 2012.

J. F. Muzy, E. Bacry, R. Baile, and P. Poggi, “Uncovering latent singularities from multifractal scaling laws in mixed asymptotic regime. Application to turbulence,” EPL (Europhysics Letters), vol. 82, no. 6, article 60007, 2008.

A. R. Subramaniam, I. A. Gruzberg, and A. W. W. Ludwig, “Boundary criticality and multifractality at the two-dimensional spin quantum Hall transition,” Physical Review B, vol. 78, no. 24, article 245105, 2008.

P. C. Ivanov, L. A. N. Amaral, A. L. Goldberger et al., “Multifractality in human heartbeat dynamics,” Nature, vol. 399, no. 6735, pp. 461–465, 1999.

D. Makowiec, A. Dudkowska, R. Gała̧ska, and A. Rynkiewicz, “Multifractal estimates of monofractality in RR-heart series in power spectrum ranges,” Physica A, vol. 388, no. 17, pp. 3486–3502, 2009.

A. Rosas, E. Nogueira Jr, and J. F. Fontanari, “Multifractal analysis of DNA walks and trails,” Physical Review E, vol. 66, no. 6, article 061906, 2002.

H. E. Stanley and P. Meakin, “Multifractal phenomena in physics and chemistry,” Nature, vol. 335, no. 6189, pp. 405–409, 1988.

V. V. Udovichenko and P. E. Strizhak, “Multifractal properties of copper sulfide film formed in self-organizing chemical system,” Theoretical and Experimental Chemistry, vol. 38, no. 4, pp. 259–262, 2002.

A. Witt and B. D. Malamud, “Quantification of long-range persistence in geophysical time series: conventional and benchmark-based improvement techniques,” Surveys in Geophysics, vol. 34, no. 5, pp. 541–651, 2013.

L. Telesca, V. Lapenna, and M. Macchiato, “Multifractal fluctuations in earthquake-related geoelectrical signals,” New Journal of Physics, vol. 7, p. 214, 2005.

E. Koscielny-Bunde, J. W. Kantelhardt, P. Braun, A. Bunde, and S. Havlin, “Long-term persistence and multifractality of river runoff records: detrended fluctuation studies,” Journal of Hydrology, vol. 322, no. 1–4, pp. 120–137, 2006.

J. W. Kantelhardt, E. Koscielny-Bunde, D. Rybski, P. Braun, A. Bunde, and S. Havlin, “Long-term persistence and multifractality of precipitation and river runoff records,” Journal of Geophysical Research, vol. 111, no. D1, article D01106, 2006.

M. Ausloos, “Generalized Hurst exponent and multifractal function of original and translated texts mapped into frequency and length time series,” Physical Review E, vol. 86, no. 3, article 031108, 2012.

S. Drożdż, P. Oświȩcimka, A. Kulig et al., “Quantifying origin and character of long-range correlations in narrative texts,” Information Sciences, vol. 331, pp. 32–44, 2016.

E. A. F. Ihlen and B. Vereijken, “Multifractal formalisms of human behavior,” Human Movement Science, vol. 32, no. 4, pp. 633–651, 2013.

J. A. Dixon, J. G. Holden, D. Mirman, and D. G. Stephen, “Multifractal dynamics in the emergence of cognitive structure,” Topics in Cognitive Science, vol. 4, no. 1, pp. 51–62, 2012.

G. R. Jafari, P. Pedram, and L. Hedayatifar, “Long-range correlation and multifractality in Bach’s inventions pitches,” Journal of Statistical Mechanics: Theory and Experiment, vol. 2007, no. 4, article P04012, 2007.

P. Oświȩcimka, J. Kwapień, I. Celińska, S. Drożdż, and R. Rak, “Computational approach to multifractal music,” 2011, http://arxiv.org/abs/1106.2902.

T. C. Roeske, D. Kelty-Stephen, and S. Wallot, “Multifractal analysis reveals music-like dynamic structure in songbird rhythms,” Scientific Reports, vol. 8, no. 1, article 4570, 2018.

Z. Nagy, P. Mukli, P. Herman, and A. Eke, “Decomposing multifractal crossovers,” Frontiers in Physiology, vol. 8, article 533, 2017.

D. G. Stephen, W.-H. Hsu, D. Young et al., “Multifractal fluctuations in joint angles during infant spontaneous kicking reveal multiplicativity-driven coordination,” Chaos, Solitons and Fractals, vol. 45, no. 9-10, pp. 1201–1219, 2012.

D. Ghosh, S. Dutta, and S. Chakraborty, “Multifractal detrended cross-correlation analysis for epileptic patient in seizure and seizure free status,” Chaos, Solitons and Fractals, vol. 67, pp. 1–10, 2014.

C.-H. You, D.-I. Lee, and K. Kim, “Analysis of multifractal strengths in game behaviors,” Journal of the Korean Physical Society, vol. 66, no. 10, pp. 1617–1622, 2015.

P. D. Domański, “Multifractal properties of process control variables,” International Journal of Bifurcation and Chaos, vol. 27, no. 06, article 1750094, 2017.

D. G. Kelty-Stephen, “Threading a multifractal social psychology through within-organism coordination to within-group interactions: a tale of coordination in three acts,” Chaos, Solitons and Fractals, vol. 104, pp. 363–370, 2017.

D. G. Kelty-Stephen, K. Palatinus, E. Saltzman, and J. A. Dixon, “A tutorial on multifractality, cascades, and interactivity for empirical time series in ecological science,” Ecological Psychology, vol. 25, no. 1, pp. 1–62, 2013.

M. Ausloos and K. Ivanova, “Multifractal nature of stock exchange prices,” Computer Physics Communications, vol. 147, no. 1-2, pp. 582–585, 2002.

L. Calvet and A. Fisher, “Multifractality in asset returns: theory and evidence,” The Review of Economics and Statistics, vol. 84, no. 3, pp. 381–406, 2002.

A. Turiel and C. J. Perez-Vicente, “Role of multifractal sources in the analysis of stock market time series,” Physica A, vol. 355, no. 2–4, pp. 475–496, 2005.

P. Oświȩcimka, J. Kwapień, and S. Drożdż, “Multifractality in the stock market: price increments versus waiting times,” Physica A, vol. 347, pp. 626–638, 2005.

W.-X. Zhou, “The components of empirical multifractality in financial returns,” EPL (Europhysics Letters), vol. 88, no. 2, article 28004, 2009.

M. I. Bogachev and A. Bunde, “Improved risk estimation in multifractal records: application to the value at risk in finance,” Physical Review E, vol. 80, no. 2, article 026131, 2009.

Z.-Y. Su, Y.-T. Wang, and H.-Y. Huang, “A multifractal detrended fluctuation analysis of Taiwan’s stock exchange,” Journal of the Korean Physical Society, vol. 54, no. 4, pp. 1395–1402, 2009.

S. Drożdż, J. Kwapień, P. Oświȩcimka, and R. Rak, “The foreign exchange market: return distributions, multifractality, anomalous multifractality and the Epps effect,” New Journal of Physics, vol. 12, no. 10, article 105003, 2010.

P. Oświȩcimka, S. Drożdż, J. Kwapień, and A. Z. Górski, “Effect of detrending on multifractal characteristics,” Acta Physica Polonica A, vol. 123, no. 3, pp. 597–603, 2013.

S. Dutta, D. Ghosh, and S. Chatterjee, “Multifractal detrended cross correlation analysis of foreign exchange and SENSEX fluctuation in Indian perspective,” Physica A, vol. 463, pp. 188–201, 2016.

D. Grech, “Alternative measure of multifractal content and its application in finance,” Chaos, Solitons and Fractals, vol. 88, pp. 183–195, 2016.

S. Bayraci, “Testing for multi-fractality and efficiency in selected sovereign bond markets: a multi-fractal detrended moving average (MF-DMA) analysis,” International Journal of Computational Economics and Econometrics, vol. 8, no. 1, pp. 95–120, 2018.

E. Bacry, J. Delour, and J. F. Muzy, “Modelling financial time series using multifractal random walks,” Physica A, vol. 299, no. 1-2, pp. 84–92, 2001.

P. Oświȩcimka, J. Kwapień, S. Drożdż, A. Z. Górski, and R. Rak, “Multifractal model of asset returns versus real stock market dynamics,” Acta Physica Polonica B, vol. 37, pp. 3083–3092, 2006.

T. Lux, “The Markov-switching multifractal model of asset returns,” Journal of Business & Economic Statistics, vol. 26, no. 2, pp. 194–210, 2008.

J. Perelló, J. Masoliver, A. Kasprzak, and R. Kutner, “Model for interevent times with long tails and multifractality in human communications: an application to financial trading,” Physical Review E, vol. 78, no. 3, article 036108, 2008.

W.-X. Zhou, “Finite-size effect and the components of multifractality in financial volatility,” Chaos, Solitons and Fractals, vol. 45, no. 2, pp. 147–155, 2012.

S. Drożdż, J. Kwapień, P. Oświȩcimka, and R. Rak, “Quantitative features of multifractal subtleties in time series,” EPL (Europhysics Letters), vol. 88, no. 6, article 60003, 2009.

S. Drożdż and P. Oświȩcimka, “Detecting and interpreting distortions in hierarchical organization of complex time series,” Physical Review E, vol. 91, no. 3, article 030902(R), 2015.

P. Oświȩcimka, L. Livi, and S. Drożdż, “Right-side-stretched multifractal spectra indicate small-worldness in networks,” Communications in Nonlinear Science and Numerical Simulation, vol. 57, pp. 231–245, 2018.

J. W. Kantelhardt, S. A. Zschiegner, E. Koscielny-Bunde, S. Havlin, A. Bunde, and H. E. Stanley, “Multifractal detrended fluctuation analysis of nonstationary time series,” Physica A, vol. 316, no. 1–4, pp. 87–114, 2002.

B. Podobnik and H. E. Stanley, “Detrended cross-correlation analysis: a new method for analyzing two nonstationary time series,” Physical Review Letters, vol. 100, no. 8, article 084102, 2008.

W.-X. Zhou, “Multifractal detrended cross-correlation analysis for two nonstationary signals,” Physical Review E, vol. 77, no. 6, article 066211, 2008.

P. Oświȩcimka, S. Drożdż, M. Forczek, S. Jadach, and J. Kwapień, “Detrended cross-correlation analysis consistently extended to multifractality,” Physical Review E, vol. 89, no. 2, article 023305, 2014.

P. Oświȩcimka, J. Kwapień, and S. Drożdż, “Wavelet versus detrended fluctuation analysis of multifractal structures,” Physical Review E, vol. 74, no. 1, article 016103, 2006.

J. Kwapień, P. Oświęcimka, and S. Drożdż, “Detrended fluctuation analysis made flexible to detect range of cross-correlated fluctuations,” Physical Review E, vol. 92, no. 5, article 052815, 2015.

J. Kwapień, P. Oświȩcimka, M. Forczek, and S. Drożdż, “Minimum spanning tree filtering of correlations for varying time scales and size of fluctuations,” Physical Review E, vol. 95, no. 5, article 052313, 2017.

Data source: https://stooq.pl/.

P. Gopikrishnan, M. Meyer, L. A. N. Amaral, and H. E. Stanley, “Inverse cubic law for the distribution of stock price variations,” European Physical Journal B, vol. 3, no. 2, pp. 139-140, 1998.

S. Drożdż, M. Forczek, J. Kwapień, P. Oświȩcimka, and R. Rak, “Stock market return distributions: from past to present,” Physica A, vol. 383, no. 1, pp. 59–64, 2007.

D. Grech and Z. Mazur, “Can one make any crash prediction in finance using the local Hurst exponent idea?” Physica A, vol. 336, no. 1-2, pp. 133–145, 2004.

L. Czarnecki, D. Grech, and G. Pamuła, “Comparison study of global and local approaches describing critical phenomena on the Polish stock exchange market,” Physica A, vol. 387, no. 27, pp. 6801–6811, 2008.

S. Drożdż, F. Grümmer, F. Ruf, and J. Speth, “Log-periodic self-similarity: an emerging financial law?” Physica A, vol. 324, no. 1-2, pp. 174–182, 2003.

S. Drożdż, F. Grümmer, A. Z. Górski, F. Ruf, and J. Speth, “Dynamics of competition between collectivity and noise in the stock market,” Physica A, vol. 287, no. 3-4, pp. 440–449, 2000.

V. A. Marčenko and L. A. Pastur, “Distribution of eigenvalues for some sets of random matrices,” Mathematics of the USSR-Sbornik, vol. 1, no. 4, pp. 457–483, 1967.

A. Edelman, “Eigenvalues and condition numbers of random matrices,” SIAM Journal on Matrix Analysis and Applications, vol. 9, no. 4, pp. 543–560, 1988.

S. Drożdż, J. Kwapień, F. Grümmer, F. Ruf, and J. Speth, “Quantifying the dynamics of financial correlations,” Physica A, vol. 299, no. 1-2, pp. 144–153, 2001.

B. J. West, E. L. Geneston, and P. Grigolini, “Maximizing information exchange between complex networks,” Physics Reports, vol. 468, no. 1–3, pp. 1–99, 2008.

D. G. Stephen, N. Stepp, J. A. Dixon, and M. T. Turvey, “Strong anticipation: sensitivity to long-range correlations in synchronization behavior,” Physica A, vol. 387, no. 21, pp. 5271–5278, 2008.

D. G. Stephen and J. A. Dixon, “Strong anticipation: multifractal cascade dynamics modulate scaling in synchronization behaviors,” Chaos, Solitons and Fractals, vol. 44, no. 1–3, pp. 160–168, 2011.

D. H. Abney, C. T. Kello, and A. S. Warlaumont, “Production and convergence of multiscale clustering in speech,” Ecological Psychology, vol. 27, no. 3, pp. 222–235, 2015.

D. Delignières, Z. M. H. Almurad, C. Roume, and V. Marmelat, “Multifractal signatures of complexity matching,” Experimental Brain Research, vol. 234, no. 10, pp. 2773–2785, 2016.