Next Article in Journal
Hermite Finite Element Method for Time Fractional Order Damping Beam Vibration Problem
Next Article in Special Issue
Inner Multifractal Dynamics in the Jumps of Cryptocurrency and Forex Markets
Previous Article in Journal
Effects of Competition Intensities and R&D Spillovers on a Cournot Duopoly Game of Digital Economies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Gold and Sustainable Stocks in the US and EU: Nonlinear Analysis Based on Multifractal Detrended Cross-Correlation Analysis and Granger Causality

Institute of Economic Sciences, 11000 Belgrade, Serbia
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(10), 738; https://doi.org/10.3390/fractalfract7100738
Submission received: 11 September 2023 / Revised: 25 September 2023 / Accepted: 29 September 2023 / Published: 7 October 2023
(This article belongs to the Special Issue Complex Dynamics and Multifractal Analysis of Financial Markets)

Abstract

:
Geopolitical risks and conflicts wield substantial influence on the global economy and financial markets, fostering uncertainty and volatility. This study investigates the intricate relationship between gold and representatives of green and sustainable stocks in the US and EU during the Russia-Ukraine conflict, employing multifractal detrended cross-correlation analysis (MF-DCCA) and nonlinear Granger causality. MF-DCCA reveals significant multifractal properties and nonlinear cross-correlations across all time series pairs. Notably, conflict weakened the multifractal cross-correlations between US stocks and gold, except for the TESLA/gold pair. A similar significant change in the EU market’s multifractal properties was observed during the conflict. In conjunction with MF-DCCA, Granger causality tests indicate bidirectional nonlinear relationships between gold and green/sustainable stock markets in the USA and EU. Gold’s past movements significantly influence changes in all the considered green and sustainable stocks, enabling predictions of their behavior. These findings shed light on multifractal dynamics during geopolitical crises and emphasize the bidirectional relationships between gold and green and sustainable stock markets. This comprehensive analysis informs investors and policy makers, enhancing their understanding of financial market behavior amid geopolitical instability.

1. Introduction

The conflict between Russia and Ukraine, which began in February 2014 over a dispute about the official status of Crimea and Donbas, escalated on 24 February 2022, and eventually became the most significant war in Europe since World War II [1]. This conflict had severe implications for global financial markets. Investor reactions to news about geopolitical risks and conflicts are often exaggerated [2], which does not help the overall market stability. As an example, Abbassi et al. [3] demonstrated that stock prices are vulnerable to geopolitical risks and trade dependence, underscoring the vulnerability of the stock markets in response to such events.
The economic consequences of the Russia-Ukraine (RU-UA) conflict range from supply chain shortages (especially in base metals and food) to higher inflationary pressures, weaker growth prospects, and other impacts as the battle has escalated [4]. Since the outbreak of the RU-UA conflict, global equity markets have experienced negative cumulative abnormal returns with significant heterogeneous effects [5], accompanied by substantial volatility [6]. As inflationary pressures continue to mount and the conflict shows no signs of abating, investor sentiment appears to have come under significant pressure. Bossman and Gubareva’s [7] analysis of the financial impact of geopolitical risks triggered by the Russia-Ukraine conflict on the seven major stock markets in emerging and developed economies shows that the effect is asymmetric and market-specific, with Brazil, China, Russia, and Turkey being resilient to geopolitical risks during bear markets. The RU-UA conflict and the lingering consequences of COVID-19 are forcing central banks to raise interest rates to support currencies and social cohesion, fueling global economic concerns of a possible sharp slowdown or recession, which, combined with inflation, would lead to a highly adverse situation [8].
Stocks and gold prices are two important indicators of the health of the economy. However, the relationship between stocks and gold prices is complex and can be influenced by a number of factors. When the stock market rises, investors are generally optimistic about the future and willing to invest in companies. On the other hand, when the stock market falls, investors are more cautious and pull their money out of the market. Conversely, when there is economic uncertainty, gold is often seen as a safe haven for investors. When the price of gold rises, it indicates that investors are concerned about inflation or other economic factors that could negatively impact the value of traditional investments like stocks.
However, green and sustainable stocks are relatively less studied in the academic literature than traditional stocks, primarily due to the evolving nature of sustainable finance and environmental, social, and governance (ESG) investing. Traditional finance research has long focused on conventional metrics and analysis of standard stocks, for which well-established data sources and historical performance evidence exist. In contrast, green and sustainable stocks operate at the intersection of finance, environmental impact, and social considerations. As interest in green and sustainable investing continues to grow, academic research in these areas gradually expands; however, this growth has yet to fully bridge the relative scarcity of scholarly work. As sustainability and ESG factors gain prominence in investment decisions, we expect more research to be shedding light on the performance, risks, and impacts of green and sustainable stocks.
The main objective of this study is to examine how the Russia-Ukraine conflict has affected the relationship between the stock prices of selected green and sustainable companies in the US and the EU and the price of gold. To achieve that goal, we employ the multifractal detrended cross-correlation analysis (MF-DCCA) method to assess the cross-correlation between gold and sustainable stocks in the US and EU. Subsequently, we apply the nonlinear Granger causality test method to examine the causal relationship between these financial instruments. Furthermore, this study contributes to closing the existing research gap by exploring the impact of the Russia-Ukraine conflict on green and sustainable stocks and their relationship with gold prices. This analysis enhances our understanding of sustainable finance under geopolitical stress.
The MF-DCCA technique was first proposed by Zhou [9] and has since been successfully applied in a variety of fields, including financial markets [10,11,12,13,14], cryptocurrencies [15], Internet of Things [16], and various environmental phenomena [17,18,19]. Moreover, gold has also been an interesting research object when speaking about MF-DDCA. For instance, gold was analyzed in world markets [20], Chinese markets [21], and with Bitcoin and crude oil [22]. Multifractal analysis is commonly used to study crises’ impact on markets. Bentes [23] investigated the safe haven attributes of gold in CIVETS countries and observed that, prior to the pandemic, gold did indeed exhibit safe haven characteristics; however, as the crises of COVID-19 and the RU-UA conflict unfolded, an increase in multifractality was detected, signifying the diminishing safe haven property of gold during these crises times. Aslam [24] applied multifractal analysis to evaluate the impact of the COVID-19 pandemic on agricultural futures markets and found significant changes in intraday multifractal properties and market efficiency. Mensi [25] used multifractal analysis to study Islamic stock markets and found different levels of efficiency across sectors and changes after the global financial crisis. In addition, multifractal analysis was used to evaluate the impact of the Russia-Ukraine conflict on the intraday efficiency of energy markets [26] and to investigate the multifractal characteristics of the cross-correlation between geopolitical risk and energy markets [27].
In addition to multifractal techniques, Granger causality analysis is widely used to explore nonlinear causal relationships among financial variables that provide insight into how different factors interact with and influence market dynamics. Linear models often cannot capture the complex dynamics during financial crises or significant geopolitical events such as the Russia-Ukraine conflict. Researchers can gain valuable insights into how financial markets nonlinearly respond to such events, providing important information for investors and decision makers. For example, Fernandez [28] used the nonlinear Granger causality test to study the spillover effects of the US subprime crisis on Asian and European economies and its impact on currency and stock markets. Similarly, Alzahrani et al. [29] examined the nonlinear Granger causality between wavelet-transformed spots and future oil prices. This method extends across various domains, encompassing stock markets [30] and their interactions with oil markets [31,32], investor sentiment [33], geopolitical risk assessments [34], responses to events like COVID-19 [35], and considerations of mental health [36].
The following sections of this paper are organized as follows: Section 2 explains the materials and methods used, Section 3 presents the results, and Section 4 presents our conclusions.

2. Methodology and Data

This section presents a concise overview of the MF-DCCA and nonlinear Granger causality test methodologies, outlining their key objectives and primary steps.

2.1. MF-DCCA

The MF-DCCA is one of the most predominant methods for examining multifractal patterns in cross-correlations of time series. Multifractal cross-correlation implies that shared underlying processes or factors influence both time series, but their relationship is more intricate and nuanced than a straightforward linear correlation can depict. MF-DCCA offers distinct advantages over comparable methods when it comes to examining the relationship between two time series. MF-DCCA excels at discovering hidden nonlinear relationships that are often missing in linear methods. Its ability to handle scale invariance makes it suitable for capturing both short-term fluctuations and long-term trends. This is particularly valuable in finance, where multifractal scaling is common due to complex interactions. MF-DCCA is also robust against noise, provides statistical significance testing, and can be used outside of finance, increasing its versatility. However, it requires computational resources that are sensitive to parameter settings and may not always provide direct causal insights.
To investigate the multifractal cross-correlation between the stock prices of selected green and sustainable companies in the US and EU and the price of gold, the following steps are taken.
  • Step 1. Let x t and y t be two time series of the same length N . The cumulative deviation series, X ( t ) and Y ( t ) , are calculated for both time series using (1), where x ¯ and y ¯ represent the averages of the time series, and t = 1 ,   2 ,   ,   N :
    X   t = k = 1 t   x k x ¯     a n d   Y   t = k = 1 t   y k y ¯  
  • Step 2. The X ( t ) and Y ( t ) are divided into   N s non-overlapping segments, where s is the time scale. This process is also performed from the end to the start of the series to ensure all information is used. As a result, two sets of N s non-overlapping segments are obtained.
  • Step 3. Compute local trends with a k th-order polynomial fit via the least squares method for each sub-segment ν
    x ν i = a 1 i k + + a k i + a k + 1 y ν i = b 1 i k + + b k i + b k + 1
    where i = 1,2 , , s ; k = 1,2 , ; ν = 1 , ,   N s .
  • Step 4. For each of the 2 N s segments, we determine the local variance
    F 2 v , s = 1 s i = 1 s X v 1 s + i x v i Y v 1 s + i y v i
    for segments v = 1,2 , , N s and
    F 2 v , s = 1 s i = 1 s | X N v N s s + i x v i | | Y N v N s s + i y v i |
    for segments v = N s + 1 , N s + 2 , , 2 N s .
  • Step 5. Then, we compute the q th order wave function
    F q s = 1 2 N s v = 1 2 N s   F 2 v , s q 2 1 q ,   q   0 F 0 ( s ) = exp { 1 4 N s v = 1 2 N s   log F 2 v , s } ,   q = 0 .
If the series are long-range power-law correlated, F q s increases when s increases. In this case, the power-law correlations satisfy F q s ~ s h x y q , where h x y q denotes the generalized Hurst exponent. If there is only short-range or no correlation in the time series, then h x y 2 = 0.5 . If there is a long-range power-law correlation, then h x y 2 0.5 . Furthermore, if h x y 2 > 0.5 , the long-range auto-correlations are persistent, while if h x y 2 < 0.5 , long-range auto-correlations are anti-persistent. To measure the strength of multifractal features, we compute Δ h x y = h x y   q m i n h x y q m a x . The larger Δ H x y , the stronger the multifractal feature.
The generalized Hurst exponent h x y is directly related to the multifractal scaling exponent τ ( q ) via the following relation: τ ( q ) = q h x y ( q ) 1 . A linear form for the multifractal scaling exponent characterizes the monofractal time series.
The singularity strength α x y = h x y q + q h x y ( q ) describes the singular degree of each segment in a complex system. The singularity spectrum f x y ( α ) = q a x y q h x y + 1 describes the fractal dimension. The spectrum gives information about the relative dominance of various fractal exponents in the series. The wider the range of singularity strength Δ α x y = α x y m a x α x y m i n , the more intense data fluctuation exists.

2.2. Nonlinear Granger Causality

The Granger causality test method proposed by Granger [37] is widely used to analyze the risk conduction effect. Existing research proposes the nonlinear Granger Causality test method [38] based on the linear Granger causality test method, which compensates for the limitation of the linear Granger causality test method. Hence, the key advantage of the nonlinear Granger causality test is its ability to uncover complex, nonlinear causal relationships in data.
The process of conducting a nonlinear Granger causality test is outlined as demonstrated by [30,36]. Consider two time series, { X t } and { Y t } , which are assumed to exhibit strict stationarity and weak dependence. The m -length leading vector of { X t } is represented as { X t m } , while the lag vectors of α -length and β -length for { X t } and { Y t } are denoted as X t α α and Y t β β , respectively.
Formally, we define that
X t m     X t ,   X t + 1 ,   X t + 2 , ,   X t + m 1 Y t m     Y t ,   Y t + 1 ,   Y t + 2 , ,   Y t + m 1 X t α α     X t ,   X t + 1 ,   X t + 2 , ,   X t α Y t β β     Y t ,   Y t + 1 ,   Y t + 2 , ,   Y t β .
For given values of m ,   α , and β (each greater than or equal to 1 ) and a positive constant e , the nonlinear Granger causality test indicates that { Y t } fails to Granger-cause { X t } if the following condition is met [30,36]:
P r ( | | X t m     X s m | | < e | | X t α α     X s α α | | < e ,   | | Y t β β     Y s β β   | | < e ) = P r ( X t m     X s m | | < e | |   X t α α X s α α | < e .
Here, P r ( · ) denotes the probability and · denotes the maximum norm.
The left-hand side probability signifies the likelihood that two arbitrary m -length leading vectors of { X t } are within a distance of e from each other, provided their corresponding α -length lag vectors of { X t } are at a distance shorter than e from each other.
Testing the null hypothesis H 0 :   { Y t } does not exhibit nonlinear Granger causality with respect to { X t } but becomes more insightful when expressing the conditional probabilities through ratios of their corresponding joint probabilities.
We utilize the ratios of joint probabilities, denoted as C 1 ( m + α , β , e )   C 2 ( α , β , e ) and C 3 ( m + α , e ) C 4 α , e , to assess the left-hand and right-hand sides. We define joint probabilities as follows [30,36]:
C 1 ( m + α ,   β ,   e )     P r ( | | X t α m + α X   s α m + α | | < e ,   | | Y t β β     Y s β β   | | < e ) C 2 ( α ,   β ,   e )     P r ( | | X t α α     X s α α   | | < e ,   | | Y t β β     Y s β β   | < e ) C 3 ( m + α ,   e )     P r ( | | X t α m + α X   s α m + α | | < e ) C 4 ( α ,   e )     P r ( | | X t α α     X s α α   | | < e )
Given the provided values of m ,   α , and β (all of which are greater than or equal to 1) and a positive constant e , the null hypothesis H 0 can be formulated as follows [30,36]:
H 0   :   C 1 ( m + α ,   β ,   e ) C 2 ( α ,   β ,   e ) = C 3 ( m + α ,   e ) C 4 ( α ,   e )
Let us define I as a kernel, denoted as I = ( Z 1 ,   Z 2 ,   e ) , which takes on the value of 1 when two compatible vectors Z 1 and Z 2 are within a maximum norm distance of e from each other, and 0 otherwise. This kernel function plays a crucial role in determining proximity based on the specified conditions, [30,36].
I = ( Z 1 ,   Z 2 ,   e ) = 1 ,   Z 1 Z 2     e 0 ,   Z 1 Z 2 > e  
The correlation-integral estimators of the joint probabilities can subsequently be expressed as follows [30,36]:
C 1 ( m + α ,   β ,   e ,   λ )     2 λ ( λ 1 ) t < s I ( x t α m + α   ,   x s α m + α   ,   e )   ·   I ( y t β β   ,   y s β β   ,   e ) C 2 ( α ,   β ,   e ,   λ )   2 λ ( λ 1 ) t < s I ( x t α α   ,   x s α α   ,   e )   ·   I ( y t β β   ,   y s β β   ,   e ) C 3 ( α ,   β ,   e ,   λ )   2 λ ( λ 1 ) t < s I ( x t α m + α   ,   x s β m + β   ,   e ) C 4 ( α ,   e ,   n )   2 λ ( λ 1 ) t < s I ( x t α α   ,   x s β β   ,   e )
where t ,   s = m a x ( α ,   β ) + 1 ,   .   .   .   ,   T m + 1 ,   n = T m a x ( α ,   β ) m + 1 .
When considering the provided values of m     1 ,   L x     1 ,   L y     1 , and e > 0 , the test statistics are defined utilizing the joint probability estimators as follows [30,36]:
n C 1 m + α ,   β ,   e C 2 α ,   β ,   e C 3 m + α ,   e C 4 α ,   e α   N 0 ,   σ 2 m ,   α ,   β ,   e .
In this context, the variance is consistently estimated [39]. Assuming the values of m ,   α , and β     1 and e > 0 , if { X t } does not satisfy the nonlinear Granger cause with respect to { Y t } , then the test statistic follows an asymptotic normal distribution with a zero mean and a constant variance. The conventional critical values are applicable when used to test the null hypothesis that the stock price X t does not nonlinearly Granger-cause { Y t } , due to the asymptotic normality of the test statistics. A similar procedure is employed to test the hypothesis that Y t does not nonlinearly Granger-cause X t .

2.3. Data

We selected 8 representatives of stocks of the greenest and sustainable companies from the United States and Europe. In selecting companies for this study, we were guided by several key criteria. First, we sought diversity across industries, encompassing renewable energy, waste management, utilities, and real estate, thus addressing various environmental challenges. Second, we focused on industry leaders renowned for their innovation, providing valuable insights into the intersection of innovation and sustainability. Third, we prioritized companies with significant global impact, actively shaping sustainability trends and supply chains. Fourth, to facilitate comparisons, we considered companies operating in both the United States and Europe, regions that are at the forefront of sustainability. Finally, and most crucially, the availability of robust and reliable data, critical to scientific research, strongly influenced our selection, leading us to the eight chosen companies. The data were obtained from the Yahoo Finance website. The selected time series spanned from August 4, 2014, to February 21, 2023, and a description of the gold and stocks series can be found in Table 1.
The evolution and development of prices are depicted in Figure 1, while the descriptive statistics and stationarity tests are presented in Table 2.
Descriptive statistics for eight stocks and gold were analyzed for three periods: the whole period, before the RU-UA conflict, and after the conflict. The analysis revealed some interesting findings about the behavior of these assets during the different periods.
First, the average price of most stocks and gold increased after the RU-UA conflict compared with the period before. This can be attributed to the market’s response to geopolitical events, where investors tend to favor safe assets such as gold and utility stocks during periods of uncertainty. In contrast, the average price of Vonovia SE (VNA.DE) declined after the RU-UA conflict, which could be attributed to factors such as increased competition, supply chain disruptions, and/or regulatory issues. The median values are not significantly different from the mean values.
Second, the standard deviation (SD) for all stocks and gold declined after the RU-UA conflict compared with the pre-conflict period. This can be explained by the significantly lower number of observations after the conflict and the fact that they are green and sustainable companies since, from a theoretical point of view, we could expect increased volatility in the market due to geopolitical tensions.
The results also suggest that investors should carefully consider the risks associated with investing in individual stocks, especially those with high volatility. These conclusions are important for both investors and policy makers, as they highlight the need for careful analysis and monitoring of market behavior in times of uncertainty.
Unit root tests determine whether a time series data set is stationary or not. A stationary series has a constant mean, variance, and autocorrelation, facilitating its modelling and analysis. The Augmented Dickey-Fuller (ADF) and Phillips-Perron (PP) tests are commonly used unit root tests. From the results, it can be seen that all time series are stationary over the entire period, except Duke Energy (DUK) for the whole period and before the conflict and Waste Management (WM) after the conflict.

3. Results

3.1. MF-DCCA

In this section, we present results of MF-DCCA used to examine the cross-correlation between selected US (TSLA, FLSR, NEE, WM, DUK) and EU stock prices (VNA.DE, VWS.CO, SU.PA) and the gold price. The analysis was conducted for the entire period and separately for the series before and after the RU-UA conflict in order to examine the impact of the conflict. The reason is that previous crises, such as COVID-19, have significantly affected the gold price and may have affected the multifractality and persistence of both the gold market and green and sustainable stocks.
We first present results from MF-DCCA for the entire period to determine whether the degree of multifractal cross-correlation between the price of gold and stocks increases with the time scale. Figure 2 and Figure 3 illustrate the fluctuation functions for the US and EU companies, respectively. We only show the results for q-orders −10, −6, −2, 2, 6, and 10. All fluctuation functions exhibit an upward slope, suggesting that the cross-correlated multifractal behavior of assets increases with the time scale. This is the first evidence that the gold price amplifies the multifractal behavior of both US and EU stocks.
The values of the cross-correlation generalized Hurst exponents for both the US and EU markets decline as the values of q increase, validating a strong multifractal behavior, as shown in Figure 4a,d. Moreover, we can observe persistence, since h x y ( 2 ) is greater than 0.5. Another indication of multifractality is provided by the scaling exponent properties, which exhibit a nonlinear relationship with q (Figure 4b,f). Lastly, we examine the strength and spectrum of multifractality to analyze the time series pair. The width of the multifractal spectra is considerably larger than zero, which demonstrates that the series are multifractal (Figure 4c,f).
Now, we compare time series pairs in the pre-conflict and conflict periods. First, we analyze the results for US stocks (Figure 5). Focusing on green and sustainable stocks in the US market during the pre-conflict period, we found that the generalized Hurst exponent decreases with increasing q , indicating that each series pair possesses a multifractal property (Figure 5a,b). Furthermore, when q = 2 , the exponents are larger than 0.5, indicating that all time series pairs have persistence. During the conflict period, we observed a lower decay rate of the exponent for all stocks except for DUK. The exponent for DUK did not follow the same pattern as in the pre-conflict period. The values of the generalized Hurst exponents during the conflict period were larger than those in the pre-conflict period, indicating that the cross-correlation was more persistent during the conflict period. The scaling exponents for US markets are nonlinearly dependent on q , showing further evidence of multifractality for both periods (Figure 5c,d). Lastly, we use multifractal strength and spectra to examine time series pairs in the pre- and post-conflict periods. The widths of multifractal spectra are significantly nonzero, indicating that all the series are multifractal (Figure 5e,f).
Overall, we can conclude that multifractal cross-correlation exists between US stocks and gold prices in both periods. These results mean that the time series are correlated in a complex, nonlinear way that a single correlation coefficient cannot fully describe.
Regarding the EU market, we found that during the pre-conflict period, there was evidence of multifractal cross-correlation between gold and the selected green and sustainable stocks (Figure 6). Specifically, we observed that the scaling exponents decrease with increasing q , indicating multifractality in the cross-relations (Figure 6a,b). The scaling exponents for q < 0 are larger than those for q > 0 , although they are all larger than 0.5. This suggests that the cross-correlated behavior of small fluctuations is more persistent than that of large fluctuations. The nonlinear dependency and multifractality of the analyzed relationships are further supported by the fact that the multifractal exponents are nonlinearly dependent on q . To better understand the nonlinear relationship between gold and the selected green and sustainable stocks in the EU market, we conducted a further analysis using multifractal spectra. Our results showed that during the pre-conflict period, there were clear departures from a random walk process for all cases, supporting the presence of multifractality in the cross-correlations. However, during the conflict period, we observed significant changes in the results and cross-correlation properties, suggesting a possible impact on the underlying dynamics of the system—specifically, the generalized Hurst exponent of VNA.DE/Gold and VWS. CO/Gold was no longer a decaying function of q for q > 0 . This suggests that the conflict may have affected the strength or nature of the interactions between the components of the system. Furthermore, our analysis showed that the functions were monotonically increasing for 0 < q < 2 , indicating an increase in the persistence of the cross-correlations, with the “normal” fluctuations enhancing. The multifractal exponents did not show an obvious nonlinear dependence with q (Figure 6c,d), while the multifractal spectra no longer had a reversed U shape (Figure 6e,f). Regarding the pair SU. PA/Gold, the results still exhibited multifractal properties as in the pre-conflict period, but they were much weaker. Overall, our findings suggest that the conflict significantly impacted the multifractal nature and cross-correlations of the relationships between gold and the selected green and sustainable stocks in the EU market.
A consistent fall in the generalized Hurst exponents further strengthens the multifractality of time fluctuations of the cross-correlation between the gold price and the US green and sustainable stock markets. Turning to the specific case of q-order 2, it is clearly seen that the generalized Hurst exponents are higher in the conflict than in the pre-conflict period. The implication is that the conflict intensifies the impact of gold prices on the persistence of all the stock markets, making them less efficient than before the conflict. All time series pairs have a larger Δ h x y and a larger Δ α x y in the pre-conflict period, suggesting stronger multifractality and greater cross-correlations, except TSLA, where all the values are larger in the conflict period. The implication is that the conflict intensifies the impact of gold prices on the persistence of all the stock markets, making them less efficient than before the conflict. All the stock markets’ generalized Hurst exponents rise significantly above 0.5. Thus, the null hypothesis of random walk is rejected in favor of persistence and market inefficiency.
The prices of selected EU green and sustainable companies show the same dynamics concerning gold prices. The multifractal features of cross-correlation are weaker in the period after the conflict. This can be viewed in terms of lower values of Δ h x y and Δ α x y .
The results in Table 3 and Table 4 could be interpreted [in several ways. First, the time series pairs are still correlated after the conflict, but the nature of the correlation has changed. The larger value of Δ h x y and Δ α x y in the pre-conflict period could indicate a stronger, more complex correlation between the time series pairs during that time, which has since weakened or become more linear. Second, the conflict itself may have affected the dynamics of one or both of the time series, leading to changes in their multifractal properties. Finally, it is important to consider the possibility of spurious correlations or other confounding factors that could affect the MF-DCCA analysis. For example, if other significant events or trends are happening around the same time as the conflict, these could influence the analysis results. In any case, further analysis and contextual information would be needed to fully interpret the results of the MF-DCCA analysis and understand the implications of the results before and after the RU-UA conflict.

3.2. Nonlinear Granger

While the MF-DCCA methodology is proficient in detecting cross-correlations, it does not possess the capacity to discern the direction of a relationship. Hence, we employ the nonlinear Granger causality test method to establish and characterize causal relationships, subsequently comparing the obtained results using the multifractal approach on the whole considered period. It is imperative to maintain stationarity in the time series utilized for the Granger causality test in order to mitigate the influence of any autoregressive phenomena.
Table 5 and Table 6 show that bidirectional nonlinear Granger causality relationships exist between gold and two sustainable and green stock markets (EU and USA), where α = β = 1, 2, 3, 4.
The past values of TSLA, FSLR, NEE, WM, and DUK do not significantly Granger-cause changes in Gold’s value for most lag orders. The p-values for these cases are above 0.01, indicating a lack of significant causal influence. On the other hand, gold’s past values significantly Granger-cause fluctuations in TSLA, FSLR, NEE, WM, and DUK for most lag orders (p < 0.01). This suggests that gold’s past movements can predict changes in these stocks. As the lag order decreases, the p-values generally remain below 0.01 for gold’s influence on the stocks, confirming the consistent predictive relationships.
The nonlinear Granger causality test reveals significant relationships between gold and EU stocks, investigated in both directions. Our analysis indicates that gold’s past values exhibit a robust and nonlinear influence on changes in the value of all three stocks, VNA.DE, VWS.CO, and SU.PA. The past values of VNA.DE significantly and nonlinearly influence changes in gold’s value, suggesting a predictive relationship. On the other hand, p-values for the relationships VWS.CO→Gold and SU.PA→Gold are consistently greater than 0.01 for all lag orders. This indicates that the past values of VWS.CO and SU.PA do not significantly Granger-cause changes in gold’s value, suggesting an absence of predictive influence.

4. Conclusions

In this paper, we investigate the impact of gold prices on the persistence and efficiency of green and sustainable stocks in the US and EU markets. Our results are of great importance to investors, financial institutions, governments, and policy makers as they shed light on the complex relationship between these variables and their implications for financial markets and economic policies.
First, we use MF-DCCA to uncover multifractal cross-correlations between the prices of selected green and sustainable stocks and the price of gold. Our analysis reveals significant multifractal properties and nonlinear cross-correlations between all pairs of time series examined. In particular, we observe a distinct shift in multifractal properties within the EU market, especially during the period of the RU-UA conflict. This highlights the differential resilience of green and sustainable stock markets during geopolitical events and underscores the need to assess market dynamics and their potential impact on asset prices during periods of uncertainty.
To further explore these cross-correlations, we use the nonlinear Granger causality test method, which uncovered bidirectional nonlinear Granger causality relationships between gold and both US and EU green and sustainable stocks. This finding deepens our understanding of the complex linkages within financial markets and underscores the importance of considering gold as a relevant risk factor in portfolio construction. Consequently, this finding has direct implications for the risk management and diversification strategies of investors and portfolio managers.
Financial institutions should take these insights into account when designing and managing financial products and services, as they play a central role in facilitating investment in these markets. Recognizing the potential substitution effect between the two markets can inform the development of financial instruments that are aligned with investors’ risk exposures and sustainability preferences. In addition, institutions need to actively monitor and manage risks arising from market spillover effects to ensure portfolio stability.
For governments and policy makers, especially those committed to sustainable finance and climate-friendly investing, this study offers valuable insights. Policy makers should consider potential spillover effects and contagion risks when creating regulatory frameworks and incentives to encourage investment in green and sustainable stocks. Recognizing that these markets are interconnected is essential to a holistic and adaptive approach to sustainable finance. This allows investors and policy makers to navigate evolving market dynamics while aligning financial goals with global sustainability goals, ensuring the stability of financial markets.
It is vital to explore the causes of these dynamic changes for a more comprehensive understanding. Therefore, future research should examine these series using additional methods such as wavelet analysis or empirical mode decomposition to gain further insight into the changes in time series properties before and after the conflict.

Author Contributions

Conceptualization, M.K. and P.M.; methodology, M.K. and P.M.; software, M.K; validation, M.K., P.M. and J.M.; formal analysis, M.K. and P.M.; investigation, M.K. and P.M.; data curation, P.M.; writing—original draft preparation, M.K., P.M. and J.M.; writing—review and editing, M.K. and P.M.; visualization, M.K.; supervision, P.M.; funding acquisition, M.K., P.M. and J.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Science, Technological Development and Innovation of the Republic of Serbia under contract number 451-03-47/2023-01/200005.

Data Availability Statement

Publicly available datasets were analyzed in this study. This data can be found here: https://finance.yahoo.com (accessed on 29 September 2023).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Umar, M.; Riaz, Y.; Yousaf, I. Impact of Russian-Ukraine war on clean energy, conventional energy, and metal markets: Evidence from event study approach. Resour. Policy 2022, 79, 102966. [Google Scholar] [CrossRef]
  2. Zaremba, A.; Cakici, N.; Demir, E.; Long, H. When bad news is good news: Geopolitical risk and the cross-section of emerging market stock returns. J. Financ. Stab. 2022, 58, 100964. [Google Scholar] [CrossRef]
  3. Abbassi, W.; Kumari, V.; Pandey, D.K. What makes firms vulnerable to the Russia-Ukraine crisis? J. Risk Financ. 2023, 24, 24–39. [Google Scholar] [CrossRef]
  4. KPMG. Russia and Ukraine Conflict: Economic Implications. Available online: https://kpmg.com/xx/en/home/insights/2022/03/russia-and-ukraine-conflict-chief-economist-team.html (accessed on 11 July 2023).
  5. Boubaker, S.; Goodell, J.W.; Pandey, D.K.; Kumari, V. Heterogeneous impacts of wars on global equity markets: Evidence from the invasion of Ukraine. Financ. Res. Lett. 2022, 48, 102934. [Google Scholar] [CrossRef]
  6. Sun, M.; Zhang, C. Comprehensive analysis of global stock market reactions to the Russia-Ukraine war. Appl. Econ. Lett. 2022, 30, 2673–2680. [Google Scholar] [CrossRef]
  7. Bossman, A.; Gubareva, M. Asymmetric impacts of geopolitical risk on stock markets: A comparative analysis of the E7 and G7 equities during the Russian-Ukrainian conflict. Heliyon 2023, 9, e13626. [Google Scholar] [CrossRef] [PubMed]
  8. Basdekis, C.; Christopoulos, A.; Katsampoxakis, I.; Nastas, V. The Impact of the Ukrainian War on Stock and Energy Markets: A Wavelet Coherence Analysis. Energies 2022, 15, 8174. [Google Scholar] [CrossRef]
  9. Zhou, W.X. Multifractal detrended cross-correlation analysis for two nonstationary signals. Phys. Rev. E 2008, 77, 066211. [Google Scholar] [CrossRef]
  10. Ruan, Q.; Yang, H.; Lv, D.; Zhang, S. Cross-correlations between individual investor sentiment and Chinese stock market return: New perspective based on MF-DCCA. Phys. A Stat. Mech. 2018, 503, 243–256. [Google Scholar] [CrossRef]
  11. Lahmiri, S.; Bekiros, S. Nonlinear analysis of Casablanca stock exchange, Dow Jones, and S&P500 industrial sectors with a comparison. Phys. A Stat. Mech. 2020, 539, 122923. [Google Scholar] [CrossRef]
  12. Li, S.; Lu, X.; Liu, X. Dynamic relationship between Chinese RMB exchange rate index and market anxiety: A new perspective based on MF-DCCA. Phys. D Nonlinear Phenom. 2020, 541, 123405. [Google Scholar] [CrossRef]
  13. Yuan, X.; Sun, Y.; Lu, X. SHIBOR Fluctuations and Stock Market Liquidity: An MF-DCCA Approach. Emerg. Mark. Financ. Trade 2022, 58, 2050–2065. [Google Scholar] [CrossRef]
  14. Fernandes, L.H.; Silva, J.W.L.; De Araujo, F.H.A.; Tabak, B.M. Multifractal cross-correlations between green bonds and financial assets. Financ. Res. Lett. 2023, 53, 103603. [Google Scholar] [CrossRef]
  15. Lahmiri, S.; Bekiros, S. Big data analytics using multifractal wavelet leaders in high-frequency Bitcoin markets. Chaos Solit. Fractals 2020, 131, 109472. [Google Scholar] [CrossRef]
  16. Miao, R.; Zhang, B. Analysis on Time-Series Data from Movie Using MF-DCCA Method and Recurrent Neural Network Model under the Internet of Things. Comput. Intell. 2022, 2022, 7400833. [Google Scholar] [CrossRef] [PubMed]
  17. Zhang, C.; Ni, Z.; Ni, L. Multifractal detrended cross-correlation analysis between PM2.5 and meteorological factors. Phys. A Stat. Mech. 2015, 438, 114–123. [Google Scholar] [CrossRef]
  18. Kojić, M.; Schlüter, S.; Mitić, P.; Hanić, A. Economy-environment nexus in developed European countries: Evidence from multifractal and wavelet analysis. Chaos Solit. Fractals 2022, 160, 112189. [Google Scholar] [CrossRef]
  19. Wang, J.; Shao, W.; Kim, J. Cross-correlations between bacterial foodborne diseases and meteorological factors based on MF-DCCA: A case in South Korea. Fractals 2020, 28, 2050012. [Google Scholar] [CrossRef]
  20. Lin, M.; Wang, G.; Xie, C.; Stanley, H.E. Cross-correlations and influence in world gold markets. Phys. A Stat. Mech. 2018, 490, 504–512. [Google Scholar] [CrossRef]
  21. Ruan, Q.; Wei, J.; Guo, M. Cross-correlations between price and volume in Chinese gold markets. Phys. A Stat. Mech. 2016, 451, 10–22. [Google Scholar] [CrossRef]
  22. Jin, J.; Jiang, Y.; Hu, Y.; Shang, Y. Which one is more informative in determining price movements of hedging assets? Evidence from Bitcoin, gold and crude oil markets. Phys. A Stat. Mech. 2019, 527, 121121. [Google Scholar] [CrossRef]
  23. Bentes, S.R. Is gold a safe haven for the CIVETS countries under extremely adverse market conditions? Some new evidence from the MF-DCCA analysis. Phys. A Stat. Mech. 2023, 623, 128898. [Google Scholar] [CrossRef]
  24. Aslam, F.; Ferreira, P.; Ali, H. Analysis of the Impact of COVID-19 Pandemic on the Intraday Efficiency of Agricultural Futures Markets. J. Risk Financ. Manag. 2022, 15, 607. [Google Scholar] [CrossRef]
  25. Mensi, W.; Tiwari, A.K.; Yoon, S.M. Global financial crisis and weak-form efficiency of Islamic sectoral stock markets: An MF-DFA analysis. Phys. A Stat. Mech. 2017, 471, 135–146. [Google Scholar] [CrossRef]
  26. Aslam, F.; Slim, S.; Osman, M.; Tabche, I. The footprints of Russia-Ukraine war on the intraday (in)efficiency of energy markets: A multifractal analysis. J. Risk Financ. 2023, 24, 89–104. [Google Scholar] [CrossRef]
  27. Aslam, F.; Ferreira, P.; Ali, H.; José, A.E. Application of Multifractal Analysis in Estimating the Reaction of Energy Markets to Geopolitical Acts and Threats. Sustainability 2022, 14, 5828. [Google Scholar] [CrossRef]
  28. Fernandez, V. Linear and non-linear causality between price indices and commodity prices. Resour. Policy 2014, 41, 40–51. [Google Scholar] [CrossRef]
  29. Alzahrani, M.; Masih, M.; Al-Titi, O. Linear and non-linear Granger causality between oil spot and futures prices: A wavelet based test. J. Int. Money Financ. 2014, 48, 175–201. [Google Scholar] [CrossRef]
  30. Cao, G.; Han, Y.; Li, Q.; Xu, W. Asymmetric MF-DCCA method based on risk conduction and its application in the Chinese and foreign stock markets. Phys. A Stat. Mech. 2017, 468, 119–130. [Google Scholar] [CrossRef]
  31. Bal, D.P.; Dash, D. Nonlinear Granger causality between oil price and stock returns in India. J. Public Aff. 2020, 21, e2137. [Google Scholar] [CrossRef]
  32. Peng, Y.; Chen, W.; Wei, P.; Yu, G. Spillover effect and Granger causality investigation between China’s stock market and international oil market: A dynamic multiscale approach. J. Comput. Appl. Math. 2020, 367, 112460. [Google Scholar] [CrossRef]
  33. Chu, X.; Wu, C.; Qiu, J. A nonlinear Granger causality test between stock returns and investor sentiment for Chinese stock market: A wavelet-based approach. Appl. Econ. 2015, 48, 1915–1924. [Google Scholar] [CrossRef]
  34. Li, S.; Tu, D.; Yan, Z.; Gong, C.; Yuan, D. Does geopolitical risk matter in crude oil and stock markets? Evidence from disaggregated data. Energy Econ. 2022, 113, 106191. [Google Scholar] [CrossRef]
  35. Bal, D.; Mohanty, S. Sectoral nonlinear causality between stock market volatility and the COVID-19 pandemic: Evidence from India. Asian Econ. Lett. 2021, 2, 1–4. [Google Scholar] [CrossRef]
  36. Ruan, Q.; Zhang, M.; Lv, D.; Yang, H. SAD and stock returns revisited: Nonlinear analysis based on MF-DCCA and Granger test. Phys. A Stat. Mech. 2018, 509, 1009–1022. [Google Scholar] [CrossRef]
  37. Granger, C.W.J. Investigating causal relations by econometric models and cross-spectral methods. Econometrica 1969, 37, 424–438. [Google Scholar] [CrossRef]
  38. Baek, E.; Brock, A. A General Test for Nonlinear Granger Causality: Bivariate Model; Working Paper Iowa State University and University of Wisconsin: Madison, WI, USA, 1992. [Google Scholar]
  39. Hiemstra, J.D.; Jones, J. Testing for linear and nonlinear Granger causality in the stock price-volume relation. J. Financ. 1994, 49, 1639–1664. [Google Scholar] [CrossRef]
Figure 1. Evolution of prices.
Figure 1. Evolution of prices.
Fractalfract 07 00738 g001
Figure 2. Fluctuation functions for US green and sustainable companies.
Figure 2. Fluctuation functions for US green and sustainable companies.
Fractalfract 07 00738 g002
Figure 3. Fluctuation functions for EU green and sustainable companies.
Figure 3. Fluctuation functions for EU green and sustainable companies.
Fractalfract 07 00738 g003
Figure 4. MF-DCCA results.
Figure 4. MF-DCCA results.
Fractalfract 07 00738 g004
Figure 5. US: pre- vs. post-conflict.
Figure 5. US: pre- vs. post-conflict.
Fractalfract 07 00738 g005aFractalfract 07 00738 g005b
Figure 6. EU: pre- vs. post RU-UA conflict.
Figure 6. EU: pre- vs. post RU-UA conflict.
Fractalfract 07 00738 g006
Table 1. Description of data.
Table 1. Description of data.
Tesla (TSLA)A sustainable energy company based in the United States that manufactures electric vehicles and renewable energy products to reduce the world’s dependence on fossil fuels.
First Solar (FSLR)A US-based environmentally friendly technology company that specializes in solar cell manufacturing and is driving the transition to renewable energy sources.
NextEra Energy (NEE)A US-based clean energy company that generates electricity from renewable sources such as wind and solar power and is committed to reducing carbon emissions.
Waste Management (WM)A US-based waste management and environmental services company that works to minimize waste and promote recycling to create a more sustainable future.
Duke Energy (DUK)A US-based energy company that has made significant investments in renewable energy and has a goal of zero carbon emissions by 2050.
Vonovia SE (VNA.DE)A German real estate company focused on sustainable housing, promoting energy efficiency, and reducing the carbon footprint of its buildings.
Vestas Wind Systems A/S (VWS.CO)A Danish wind turbine manufacturer driving the shift to renewable energy sources and reducing carbon emissions.
Schneider Electric (SU.PA)A French multinational company specializing in energy management and automation solutions that help optimize energy efficiency and reduce environmental impact.
GoldMany factors, including global economic trends, political events, and changes in interest rates, influence the price of gold. In recent years, gold prices have fluctuated due to different uncertainties.
Table 2. Descriptive statistics and unit root test results.
Table 2. Descriptive statistics and unit root test results.
Whole PeriodMeanMedMaxMinSDSkeKurJBADFPPObs
TSLA87.5621.68409.979.58107.611.212.95438.99 **−2.03−2.131816
FSLR66.1459.92180.1926.3328.951.746.351764.82 **−1.60−1.631816
NEE47.1139.4490.2518.5322.580.381.64183.64 **−3.26−3.401816
WM98.0491.08175.2944.0537.530.371.94127.08 **−3.45 *3.261816
DUK87.1685.82115.4364.1510.770.452.4188.89 **−4.36 **−4.46 **1816
VNA.DE38.7829.1662.2218.1211.170.141.9687.61 **−0.59−0.421816
VWS.CO129.64110.0231234.0661.320.822.68212.8 **−2.66−2.471816
SU.PA88.4372.34177.8245.9332.680.872.37256.76 **−2.24−2.361816
GOLD1458.491319.42051.51050.8276.230.461.62205.91 **−3.21−2.811816
Before RU-UA conflict
TSLA67.9520.42409.979.5896.311.774.81056.91 **−1.54−1.521604
FSLR59.4155.66121.1426.3318.620.763.17156.54 **−2.69−2.691604
NEE42.8336.2390.2518.5320.410.652.08168.18 **−2.87−2.811604
WM89.9985.81166.8344.0532.140.502.3695.24 **−2.53−2.921604
DUK84.9683.99107.9364.159.170.432.5463.22 **−4.27 **−4.20 **1604
VNA.DE39.9840.4962.2218.1211.020.021.9967.92 **−3.10−2.961604
VWS.CO122.72105.4531234.0661.521.143.29353.37 **−2.10−1.861604
SU.PA82.4470.69177.8245.9329.731.353.77529.71 **−1.85−1.901604
GOLD1412.171293.52051.51050.8258.410.772.11212.21 **−2.73−2.501604
After RU-UA conflict
TSLA235.96236.82381.82108.163.95−0.102.324.41−2.06−2.18212
FSLR117.12122180.1961.440.230.001.4022.51 **−1.93−2.03212
NEE79.5080.1689.7767.025.4−0.222.138.45 *−1.93−2.56212
WM158.93157.95175.29141.917.330.282.574.48−2.88 *−2.91 *212
DUK103.82104.75115.4385.976.56−0.462.728.29 *−2.60−2.84212
VNA.DE26.6927.0747.518.977.520.872.6827.92 **−1.41−1.23212
VWS.CO182.02184.28235.4134.8821.66−0.302.228.51 *−2.13−2.02212
SU.PA133.76133.38156.2811211.590.032.087.53 *−1.79−1.68212
GOLD1808.91811.62040.11626.796.81−0.002.057.95 *−1.45−1.42212
Note: ** indicates significance at 1% level, * indicates significance at 5% level.
Table 3. US: multifractal metrics.
Table 3. US: multifractal metrics.
StockMetricTotalPrePost
TSLA h x y ( 2 ) 1.5161.5091.497
Δ h x y 0.3230.3680.526
Δ α x y 0.4740.5250.714
FSLR h x y ( 2 ) 1.4951.4961.423
Δ h x y 0.3100.2890.215
Δ α x y 0.4460.4150.344
NEE h x y ( 2 ) 1.5371.5221.638
Δ h x y 0.5640.6980.413
Δ α x y 0.7820.9280.559
WM h x y ( 2 ) 1.5301.5201.539
Δ h x y 0.4530.5050.375
Δ α x y 0.6810.7220.525
DUK h x y ( 2 ) 1.5541.5471.652
Δ h x y 0.5910.6810.170
Δ α x y 0.8030.9030.292
Table 4. EU: multifractal metrics.
Table 4. EU: multifractal metrics.
StockMetricTotalPrePost
VNA.DE h x y ( 2 ) 1.5271.5321.486
Δ h x y 0.5790.6160.204
Δ α x y 0.7970.8400.343
VWS.CO h x y ( 2 ) 1.4971.4891.525
Δ h x y 0.2190.2710.203
Δ α x y 0.3460.3960.339
SU.PA h x y ( 2 ) 1.5051.4951.499
Δ h x y 0.4540.6030.294
Δ α x y 0.6430.8310.434
Table 5. Nonlinear Granger causality test results: USA.
Table 5. Nonlinear Granger causality test results: USA.
α TestTSLAGoldFSLRGoldNEEGoldWMGoldDUKGold
=
β GoldTSLAGold FSLRGold NEEGold WMGold DUK
1Statistics1.14565.2201.01165.5741.18068.1391.89574.4303.6113.611
p-value0.325<0.010.431<0.010.299<0.011<0.01<0.01<0.01
2Statistics1.477856.8391.40757.1512.66157.3890.19459.5424.37357.800
p-value0.1<0.010.129<0.01<0.01<0.011<0.01<0.01<0.01
3Statistics0.90736.6390.87736.5041.61436.6130.23337.4712.55736.946
p-value0.585<0.010.627<0.010.036<0.011<0.01<0.01<0.01
4Statistics0.542823.8680.52522.9620.90023.2130.11524.1031.42423.427
p-value0.976<0.010.981<0.010.615<0.011<0.010.070<0.01
Note: α = β denotes the residual series of lag order numbers; TSLA → Gold means the original hypothesis: TSLA is not a Granger causality of gold. The others have the same meaning; e = 1.5 σ , m = 1 .
Table 6. Nonlinear Granger causality test results: EU.
Table 6. Nonlinear Granger causality test results: EU.
α TestVNA.DEGoldVWS.COGoldSU.PAGold
=
β GoldVNA.DEGoldVWS.COGoldSU.PA
1Statistics5.83759.5870.86473.112−0.82874.123
p-value<0.01<0.011<0.011<0.01
2Statistics5.07455.9140.12158.7030.26160.142
p-value<0.01<0.011<0.011<0.01
3Statistics2.98836.3540.17137.0320.275638.078
p-value<0.01<0.011<0.011<0.01
4Statistics1.68823.1100.12323.4640.19324.348
p-value0.014<0.011<0.011<0.01
Note: α = β denotes the residual series of lag order numbers; VNA.DE → Gold means the original hypothesis: VNA.DE is not a Granger causality of gold. The others have the same meaning; e = 1.5 σ , m = 1 .
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kojić, M.; Mitić, P.; Minović, J. Gold and Sustainable Stocks in the US and EU: Nonlinear Analysis Based on Multifractal Detrended Cross-Correlation Analysis and Granger Causality. Fractal Fract. 2023, 7, 738. https://doi.org/10.3390/fractalfract7100738

AMA Style

Kojić M, Mitić P, Minović J. Gold and Sustainable Stocks in the US and EU: Nonlinear Analysis Based on Multifractal Detrended Cross-Correlation Analysis and Granger Causality. Fractal and Fractional. 2023; 7(10):738. https://doi.org/10.3390/fractalfract7100738

Chicago/Turabian Style

Kojić, Milena, Petar Mitić, and Jelena Minović. 2023. "Gold and Sustainable Stocks in the US and EU: Nonlinear Analysis Based on Multifractal Detrended Cross-Correlation Analysis and Granger Causality" Fractal and Fractional 7, no. 10: 738. https://doi.org/10.3390/fractalfract7100738

APA Style

Kojić, M., Mitić, P., & Minović, J. (2023). Gold and Sustainable Stocks in the US and EU: Nonlinear Analysis Based on Multifractal Detrended Cross-Correlation Analysis and Granger Causality. Fractal and Fractional, 7(10), 738. https://doi.org/10.3390/fractalfract7100738

Article Metrics

Back to TopTop