This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

On the systemic nature of global inflation, its association with equity markets and financial portfolio implications

Nick James [email protected] Kevin Chin School of Mathematics and Statistics, University of Melbourne, Victoria, Australia Arowana & Co.
Abstract

This paper uses new and recently introduced mathematical techniques to undertake a data-driven study on the systemic nature of global inflation. We start by investigating country CPI inflation over the past 70 years. There, we highlight the systemic nature of global inflation with a judicious application of eigenvalue analysis and determine which countries exhibit most "centrality" with an inner-product based optimization method. We then turn to inflationary impacts on financial market securities, where we explore country equity indices’ equity robustness and the varied performance of equity sectors during periods of significant inflationary pressure. Finally, we implement a time-varying portfolio optimization to determine which asset classes were most beneficial in increasing portfolio Sharpe ratio when an investor must hold a core (and constant) allocation to equities.

keywords:
Nonlinear time series analysis , Nonlinear financial market dynamics , Inflation , Clustering , Portfolio optimization , Evolutionary correlation
journal: Physica A

1 Introduction

Inflation has been dormant for the best part of 20 years in most developed countries globally, providing an environment for increasing business prosperity and consistent valuation multiple expansion across a number of asset classes including equities. In the wake of the COVID-19 pandemic and the associated economic and financial market crises, we are left with a new question of great pertinence - is the current elevation in global inflation levels transitory, or symptomatic of something more pathological? Rather than seeking to predict what is likely to happen next, we analyze previous periods of inflationary pressure and the associated impacts on economies, equity indices, equity sectors and investor portfolios. Our analysis takes a data-driven and mathematical approach, seeking to learn from the results of bespoke experiments, rather than approaching this problem with (often rigid) econometric theory.

Inflation is a well researched problem in the quantitative finance and economics literature [1, 2, 3, 4]. A significant amount of prior work on the topic has applied well-established econometric theory [5] to gain insights into the time-varying nature of inflationary behaviours. In the academic literature, however, there is limited work exploring the explicit impact of inflation on equity markets. Anecdotally, low-to-moderate inflation is typically associated with buoyant periods for equity markets. Inflationary extremes, where inflation may be particularly high or low, is generally considered to be harmful to equity market returns. In this paper, we test this hypothesis in a variety of contexts. Importantly, we study periods of time consistent with significant changes in inflation, rather than simply studying periods of high or low inflation. This way, we may ensure that our findings are relevant and topical in a period that is characterised by such pronounced economic uncertainty.

More generally, however, the study of financial market dynamics has been a topic of great interest to the applied mathematics, econometrics and econophysics communities. In more recent years the study of time-varying dynamics and correlation structures [6, 7, 8] has attracted the interest of researchers. A wide variety of techniques have been used to study these evolutionary dynamics including principal components analysis (PCA) [7, 9, 10, 11], clustering [12, 13], change point detection [14, 15] and various statistical modelling frameworks [Tilfani2022, 16]. Such analysis has been applied to a wide variety of underlying security types including equities [11], foreign exchange [17], and fixed income instruments [18].

There has been a recent surge in interest among researchers studying various aspects of cryptocurrencies, the blockchain and other aspects of decentralised finance. Recent work that has attracted research interest include analyses of Bitcoin’s behaviour [19, 20, 21, 22, 23], fractal-type behaviour and fractal patterns, [24, 25, 26, 16], cross-correlations and various scaling effects [27, 28, 29, 30, 31], and trading strategies [32]. Quite separately, the topic of portfolio optimization has been of significant interest to the quantitative finance community for many years [33, 34], and the core problem has been extended upon considerably since its inception [35, 36, 37, 38, 39, 40]. Although there is literature examining the effectiveness of various asset classes’ role as a portfolio ‘safe-haven’, we are unaware of any work explicitly studying the benefit of cryptocurrencies in the context of extreme inflationary periods. This work studies the effect of cryptocurrencies in multi-asset portfolios, exploring their time-varying optimal portfolio allocation.

Many of these studies are concerned with the time-varying nature of such dynamics, or the behaviour of cryptocurrencies during various market regimes. Quite naturally, the impact of COVID-19 on cryptocurrency behaviours has been widely studied [41, 42, 43, 44]. For recent work on time-varying similarity and its implications for the cryptocurrency market, and the significant impact of societal physics on the economy, readers should see [45, 46, 47].

In Section 3, we examine the similarity in country CPI trajectories and the associated collective similarity in global inflation. We then introduce a new optimization framework to identify countries that are most "central" in global inflation, and which countries’ equity indices are most robust during periods of significant inflationary change. Then in Section 4 we turn to inflationary impacts on equity market performance, where we study the evolutionary correlation structure of various equity sectors. Finally in Section 5, we implement a time-varying portfolio optimization over the past 5.5 years to determine which asset class would serve as the optimal complement to a core equity exposure, with the purpose of maximizing an investor’s Sharpe ratio. Our methods demonstrate that, despite the associated volatility, cryptocurrency was the most essential additional asset class.

2 Data

In Section 3 we study country CPI data between 01/01/1955-01/09/2021 and equity index data 01/01/1990 - 01/09/2021. This data comes from a variety of sources including The International Monetary Fund (https://www.imf.org/en/Home), Bank of Canada (https://www.bankofcanada.ca/), Organisation for Economic Cooperation and Development (https://www.oecd.org/australia/), Bank of Japan (https://www.boj.or.jp/en/index.htm/) and Bloomberg. In Section 4 we study equity sector behaviours based on underlying equity data between 01/01/2005 - 31/12/2009, this data is sourced from Bloomberg. In Section 5 we study various asset class performance (in the context of portfolio optimization) between 01/01/2016-30/06/2021. This data is sourced from Bloomberg and Coinmarketcap.com (https://coinmarketcap.com/).

3 Evolutionary country inflation dynamics

3.1 CPI Inflation trajectories

In this section of the paper, we study a collection of n=8n=8 countries listed alphabetically and indexed i=1,,ni=1,...,n. The countries in this section include: Australia, Canada, France, Germany, Italy, Japan, the United Kingdom (UK) and the United States of America (USA). We study monthly CPI inflation data between 01/01/1955 to 01/09/2021, a total of T=801T=801 months. Let xi(t)x_{i}(t)\in\mathbb{R} be the multivariate time series of CPI inflation for country ii in month tt. We wish to study the evolutionary change in inflation over time, and accordingly we compute the log returns of each country’s CPI inflation and contrast these behaviours over time. Each country’s CPI log returns are computed as follows

Δix(t)\displaystyle\Delta^{x}_{i}{(t)} =log(xi(t)xi(t1)),\displaystyle=\log\left(\frac{x_{i}{(t)}}{x_{i}{(t-1)}}\right), (1)

while Δix(t)\Delta^{x}_{i}{(t)}, references country ii’s CPI log returns at time tt. To compare inflation trajectories of countries where raw inflation varies significantly, the CPI returns time series are normalized. Each country generates an inflation returns time series, 𝚫ixT\mathbf{\Delta}^{x}_{i}\in\mathbb{R}^{T}. We let 𝚫ix1=t=1T|Δix(t)|\|\mathbf{\Delta}^{x}_{i}\|_{1}=\sum^{T}_{t=1}|\Delta^{x}_{i}(t)| be the L1L^{1} norm of the inflation returns function, and define a normalized inflation returns trajectory by 𝐓ix=𝚫ix𝚫ix1\mathbf{T}^{x}_{i}=\frac{\mathbf{\Delta}^{x}_{i}}{\|\mathbf{\Delta}^{x}_{i}\|_{1}}. Any distance between two countries’ trajectory vectors would highlight the relative similarity in CPI changes over time. We compute an inflationary trajectory matrix Dijx=𝐓ix𝐓jx1D^{x}_{ij}=\|\mathbf{T}^{x}_{i}-\mathbf{T}^{x}_{j}\|_{1}, and apply hierarchical clustering to the resulting distance matrix.

Refer to caption
(a)
Refer to caption
(b)
Refer to caption
(c)
Refer to caption
(d)
Refer to caption
(e)
Refer to caption
(f)
Refer to caption
(g)
Refer to caption
(h)
Figure 1: CPI log return time series (a) Australia (b) Canada, (c) France, (d) Germany, (e) Italy, (f) Japan, (g) UK, (h) USA. There are some notable similarities in the evolutionary inflation differences including spikes around 1980. Australia and Japan appear to exhibit a spike in inflationary differences slightly earlier than other countries.

Figure 1 reveals several noticeable similarities between various countries’ inflation between 1955-2021. The most prominent similarity is the inflationary spike exhibited by all countries during the 1970s and early 1980s. Most countries (including Canada, France, Germany, Italy, UK and the USA shown in Figures 1(b), 1(c), 1(d), 1(e), 1(g) and 1(h) respectively) display similarly-timed peaks in inflationary changes. The two notable outliers in this regard are Australia and Japan (Figures 1(a) and 1(f)) who exhibit earlier inflationary spikes than other countries. The second theme that is broadly consistent among countries are sharp drops in inflation corresponding to the GFC and COVID-19 market crises. Although there is some variability among all countries, the most interesting trajectory is that of the USA which exhibits the most pronounced shock during the GFC, followed by relatively buoyant inflationary levels. Other countries such as France, Germany, Italy and the UK exhibit negative inflationary differences more frequently than other countries.

Refer to caption
Figure 2: Hierarchical clustering applied to the distance matrix DxD^{x}. One predominant country cluster is formed with two outliers detected - Australia and Japan. This is most likely due to the earlier onset of their inflationary pressures.

To further investigate overall country inflationary structure, we apply hierarchical clustering to our inflationary trajectory matrix DxD^{x}. dendrogram displayed in Figure 2 reveals a cluster structure consistent with the CPI log return time series shown in Figure 1. The dendrogram reveals one predominant cluster consisting of the USA, Canada, France, Italy, Germany and the UK. The two outlier countries identified are Japan and Australia. Their separate classification with respect to the remaining countries is most likely due to the spike in inflation levels which precedes the bulk of countries.

3.2 Collective country inflation behaviours

Having identified two outlier countries with our trajectory analysis, we further explore the structural similarity between countries’ CPI inflationary patterns with a judicious application of eigenvalue analysis. Given the results from our methods in preceding sections, we expect the majority of inflation time series to display significant similarity over time.

We start with the following contrived scenario. Consider a large collection of financial time series, 5050, for example. If 4040 out of the 5050 time series display a high degree of affinity, then the distance matrix DxD^{x} should have the following structure:

(c1c2c3c40r10r20r30r400 r1  r2  r40     r1Tr2Tr3Tr40T    000)\begin{pmatrix}\begin{matrix}&c_{1}&c_{2}&c_{3}&\ldots&c_{40}\\ r_{1}&0&*&*&\ldots&*\\ r_{2}&*&0&*&\ldots&*\\ r_{3}&*&*&0&\ldots&*\\ &\vdots&\vdots&\vdots&\ddots\\ r_{40}&*&*&*&*&0\\ \end{matrix}&\vline&\begin{array}[]{ccc}\rule[2.15277pt]{10.76385pt}{0.5pt}&r_{1}&\rule[2.15277pt]{10.76385pt}{0.5pt}\\ \rule[2.15277pt]{10.76385pt}{0.5pt}&r_{2}&\rule[2.15277pt]{10.76385pt}{0.5pt}\\ &\vdots&\\ \rule[2.15277pt]{10.76385pt}{0.5pt}&r_{40}&\rule[2.15277pt]{10.76385pt}{0.5pt}\end{array}&\\ \hline\cr\begin{array}[]{ccccc}\rule[-4.30554pt]{0.5pt}{10.76385pt}&\rule[-4.30554pt]{0.5pt}{10.76385pt}&\rule[-4.30554pt]{0.5pt}{10.76385pt}&&\rule[-4.30554pt]{0.5pt}{10.76385pt}\\ r^{T}_{1}&r^{T}_{2}&r^{T}_{3}&\ldots&r^{T}_{40}\\ \rule[-4.30554pt]{0.5pt}{10.76385pt}&\rule[-4.30554pt]{0.5pt}{10.76385pt}&\rule[-4.30554pt]{0.5pt}{10.76385pt}&&\rule[-4.30554pt]{0.5pt}{10.76385pt}\end{array}&\vline&\begin{matrix}0&\\ &0\\ &&\ddots\\ &&&0\end{matrix}\end{pmatrix}

where rows r1,,r40r_{1},\ldots,r_{40} are time series that display a high degree of affinity and elements * are close to zero. This interpretation of this is that small deformations in elements of the matrix make the first 4040 rows identical. Therefore, the matrix is a small perturbation from a rank 1111 matrix, with 3939 eigenvalues that are equal to zero. So in this hypothetical case, if 4040 out of our 5050 time series exhibit similar behaviours, then 3939 of our eigenvalues would be close to zero.

We now return to our current data. Given a candidate threshold δ\delta, we are able to order the absolute values of the distance matrix eigenvalues |λ1||λn||\lambda_{1}|\leq...\leq|\lambda_{n}|. If |λ1|,,|λk|<δ|\lambda_{1}|,...,|\lambda_{k}|<\delta then we can deduce k+1k+1 of our country CPI inflation time series are similar with respect to their evolutionary behaviours. Eigenvalue analysis such as this provide a fast and intuitive way of assessing the scale of any distance matrix. Since our distance matrix is symmetric, it can be conjugated by an orthogonal matrix which yields a diagonal matrix of eigenvalues (this is also known as the spectral theorem). The result is that the operator norm of the matrix coincides with |λn||\lambda_{n}|, or,

maxynDxyy=Dxoper=|λn|\max_{y\in\mathbb{R}^{n}}\frac{||D^{x}y||}{||y||}=||D^{x}||_{oper}=|\lambda_{n}|

.

Refer to caption
Figure 3: Eigenspectrum of the distance matrix DxD^{x}. This is perhaps the most concise and expressive way of detailing the collective similarity in evolutionary inflation patterns. Eigenvalue analysis suggests that there are 1 or 2 countries that are dissimilar to the rest of the collection.

Figure 3 shows the magnitude of eigenvalues associated with the distance matrix, DxD^{x}. It is clear from the eigenspectrum that there is a large collection of countries with highly similar CPI inflation trajectories. For our purposes, we set δ=2.5\delta=2.5, and identify two outlier countries (consistent with both findings in the preceding section 3.1. Again, Australia and Japan are identified as the most dissimilar countries based on our eigenvalue analysis.

3.3 Inner product-based "Centrality-identification"

Given the predominant structural similarity in country CPI inflation that we have now identified, we wish to further study the systemic nature of such inflation and determine whether some countries are more ‘central’ to global inflation than others. To examine this phenomenon, we implement a time-varying linear regression on inflation log returns, and study the evolution of each model’s coefficients. Individually, this methodology allows us to see the evolution of inflationary trends over time. When viewed in conjunction with other countries’ results, this framework provides us with a general overview of the time-varying global similarity in inflationary trends. To explore the dynamic nature of each country’s CPI inflation, we study the evolution of model coefficients in two settings. Model 1 uses a 60-month rolling window, and model 2 uses a 120-month rolling window. In each case, we restrict our model to specific intervals of either 60 months or 120 months, and study the evolution of of model coefficients over time. We now compute:

Δix[t60:t]=β0,i[t60:t]+β1,i[t60:t]+ϵi[t60:t], and\displaystyle\Delta^{x}_{i}[t-60:t]=\beta_{0,i}[t-60:t]+\beta_{1,i}[t-60:t]+\epsilon_{i}[t-60:t],\text{ and} (2)
Δix[t120:t]=β0,i[t120:t]+β1,i[t120:t]+ϵi[t120:t],\displaystyle\Delta^{x}_{i}[t-120:t]=\beta_{0,i}[t-120:t]+\beta_{1,i}[t-120:t]+\epsilon_{i}[t-120:t], (3)

for Model 1 and Model 2 respectively. Implementing our model for all ii countries in the data yields a vector of evolutionary trend coefficients which we label β1,iM1(t)t{61,,T}\mathbb{\beta}^{M_{1}}_{1,i}(t)\forall t\in\{61,...,T\} and β1,iM2(t)t{121,,T}\mathbb{\beta}^{M_{2}}_{1,i}(t)\forall t\in\{121,...,T\}, for models 1 and 2 respectively. The choice of our smoothing windows is quite deliberate. For the 60-month rolling window, we wish to study inflationary trends in a more dynamic manner, using 5 year windows. For the study of longer term country inflationary patterns, we use a 120-month rolling window, as is reflected in Model 2. Although there is literature to suggest that macroeconomic cycles typically last somewhere in the order of 5-7 years, we do not wish for excessive transience in our regression parameters.

Although we study an evolutionary trend coefficient with a time-varying linear regression implementation, other modelling approaches could be used. Some possibilities of alternative models include moving average (MA) and autoregressive moving average (ARMA) models. Such analysis could be the foundation of future research, where one could study the capacity for change in the two polynomial terms among country CPI time series. Furthermore, one could explore the existence of change points and structural breaks based on changes in the behaviour of the autoregressive model parameters.

Refer to caption
(a)
Refer to caption
(b)
Figure 4: Evolutionary slope coefficients for all 8 countries studied for (a) Model 1 (60-month rolling window) and (b) Model 2 (120-month rolling window). Most countries exhibit a broadly similar trend of "up-down-up" early in their trajectories (corresponding to the 1970s, 1980s and 1990s, respectively) This is followed by a prolonged period of slope coefficients of smaller magnitude.

The general findings for both models are relatively consistent. Both Model 1 and Model 2, shown in Figures 4(a) and 4(b) respectively display a general "up-down-up" pattern of inflation trends. That is they display a positive inflation trend during the 1970s, followed by a negative trend in the 1980s and a reversion to normal levels by the early 2000s. From the early 2000s onwards, all country slope coefficients are of significantly lower magnitude, and display greater affinity between each other. One observation that can be made from both figures, is that although the general shape of inflationary trends are highly similar between countries, some countries appear to precede others in reaching local maxima and minima in slope coefficient magnitudes. This raises several interesting questions. First, which countries display the greatest differences to the rest of the collection with respect to their inflation behaviours? Second, which countries are most "central" in systemic global inflation.

To explore these questions, we introduce an inner product-based optimisation framework to determine a suitable translation (between all candidate country pairings) that would align country inflation trends most appropriately. For each country, we seek to learn an offset ϕ\phi, such that we maximize a normalized inner product. In the proceeding section where we describe our optimization methodology, we drop the superscripts M1M_{1} and M2M_{2} denoting Models 1 and 2, for the sake of more accessible exposition.

For any two countries that we compare, we allow ϕ\phi to be either positive or negative, so the general approach is flexible to the possibility that any country inflation time series may be leading or lagging another. We now define our normalized inner product between any two evolutionary trend coefficient trajectories. For any two countries’ slope trajectories β1,i\beta_{1,i} and β1,j\beta_{1,j} we define two possible inner product expressions, based on the direction of our offset ϕ\phi. If ϕ>0\phi>0, we may express our objective function

<β1,i(0:Tϕ),β1,j(ϕ:T)>n=β1,i(0)β1,j(ϕ)++β1,i(Tϕ)β1,j(T)(β1,i(0)2++β1,i(Tϕ)2)12(β1,j(ϕ)2++β1,j(T)2)12,\displaystyle<\beta_{1,i}(0:T-\phi),\beta_{1,j}(\phi:T)>_{n}=\frac{\beta_{1,i}(0)\beta_{1,j}(\phi)+...+\beta_{1,i}(T-\phi)\beta_{1,j}(T)}{(\beta_{1,i}(0)^{2}+...+\beta_{1,i}(T-\phi)^{2})^{\frac{1}{2}}(\beta_{1,j}(\phi)^{2}+...+\beta_{1,j}(T)^{2})^{\frac{1}{2}}}, (4)

If ϕ<0\phi<0 we let σ=ϕ\sigma=-\phi, and this can now be expressed

<β1,i(σ:T),β1,j(0:Tσ)>n=β1,i(σ)β1,j(0)++β1,i(T)β1,j(Tσ)(β1,i(σ)2++β1,i(T)2)12(β1,j(0)2++β1,j(Tσ)2)12.\displaystyle<\beta_{1,i}(\sigma:T),\beta_{1,j}(0:T-\sigma)>_{n}=\frac{\beta_{1,i}(\sigma)\beta_{1,j}(0)+...+\beta_{1,i}(T)\beta_{1,j}(T-\sigma)}{(\beta_{1,i}(\sigma)^{2}+...+\beta_{1,i}(T)^{2})^{\frac{1}{2}}(\beta_{1,j}(0)^{2}+...+\beta_{1,j}(T-\sigma)^{2})^{\frac{1}{2}}}. (5)

We record the optimal offset between all country slope trajectories in an inflation offset matrix which we denote DϕM1D^{M_{1}}_{\phi} and DϕM2D^{M_{2}}_{\phi}, for Models 1 and 2 respectively. To determine countries that are most and least central in the evolution of global inflation, we compute an inflation centrality score for each matrix as follows centralityjM1=iDϕ,ijM1centrality^{M_{1}}_{j}=\sum_{i}D^{M_{1}}_{\phi,ij}, and centralityjM2=iDϕ,ijM2centrality^{M_{2}}_{j}=\sum_{i}D^{M_{2}}_{\phi,ij}. The interpretation of this measure is that the lower the centrality score, the closer any country is to the rest of the collection.

Country inflation centrality scores
Country centralityM1centrality^{M_{1}} centralityM2centrality^{M_{2}}
Australia 144 94
Canada 38 49
France 74 146
Germany 78 86
Italy 14 64
Japan 130 199
UK 43 54
USA 81 56
Table 1: Inflation centrality scores for all countries under Model 1 and Model 2. In both cases, there is a majority collection of central countries. Under both models, Japan is most consistently the least central (most anomalous) country.
Refer to caption
(a)
Refer to caption
(b)
Figure 5: Hierarchical clustering applied to the two distance matrices (a) DϕM1D^{M_{1}}_{\phi} and (b) DϕM2D^{M_{2}}_{\phi}. The dendrogram shows that under both models, there is a majority cluster of countries who display similar timing in global inflationary trends. Japan is most consistently the "least central" country with respect to its inflation behaviours.

We further investigate the evolution of each country’s inflation centrality by partitioning our analysis into two discrete windows. First, we study inflation centrality between January 1960 to December 1989, and then study the period between January 1990 and December 2019. It is clear from Figure 4 that the period until 1990 exhibits significantly higher inflation - and greater volatility in country inflation trends. The period between 1990 and 2020 exhibits significantly lower inflation, and less variability in country inflation trends. Given that we have halved the length of our analysis window, we introduce a new model (Model 3) and use a 30-month rolling window.

Country inflation centrality scores: partitioned analysis
Country centralityM3centrality^{M_{3}}: 1960-1990 centralityM3centrality^{M_{3}}: 1990:2020
Australia 31 47
Canada 14 22
France 23 13
Germany 35 66
Italy 10 22
Japan 14 73
UK 24 44
USA 23 3
Table 2: Inflation centrality scores for all countries under Model 3 (30-month rolling window) between 1960-1990 and 1990-2020.

Table 2 captures some interesting insights. In our first partition, studying the period between January 1960 to December 1989, there is a material difference in inflation centrality structure. Most significantly, Australia and Germany are identified to be most anomalous (with respective scores of 31 and 35), while Japan is much more central in its inflationary behaviours. During our second analysis window which studies the period between January 1990 and December 2019, Japan is identified as the most anomalous (least central) country with respect to its inflationary behaviours. It is worth noting that Japan’s identification as an inflationary anomaly is most pronounced during a period of relative stability in global inflation.

3.4 Country equity index robustness

For our 8 countries examined in previous sections, we now wish to study the robustness of each country’s equity market during extreme inflationary periods. To do so, we introduce a new multivariate time series yi(t)y_{i}(t), which refers to each country’s respective equity index. We study data between 01/01/1990 - 01/09/2021, a total of T1=381T_{1}=381 months. We index this window of time t1=1,,T1t_{1}=1,...,T_{1}.

We compute the log returns of each equity index as follows

Δiy(t1)\displaystyle\Delta^{y}_{i}{(t_{1})} =log(yi(t1)yi(t11)).\displaystyle=\log\left(\frac{y_{i}{(t_{1})}}{y_{i}{(t_{1}-1)}}\right). (6)

From each country’s log returns time series, we construct two truncated probability density functions (PDFs). The first consists of the country equity index returns during the top decile and bottom decile of inflationary returns. That is, we extract equity index returns consistent with a country’s most extreme inflationary and deflationary pressures. The second distribution consists of equity index returns during the remaining months. Let the first pdf containing extreme equity returns be giex(x)g^{ex}_{i}(x), and the stable period pdf be gist(x)g^{st}_{i}(x). For each country ii, we compute the distance between the two pdfs. A small distance would indicate less extreme equity behaviours during inflationary crises, while a large distance may show countries that are more prone to pronounced equity market fluctuations during strongly inflationary and deflationary periods. For each country, we measure the distance between these two distributions with a Wasserstein distance[48], and term the resulting measure an equity robustness distance

ERi=dW(gist(x),giex(x)),\displaystyle ER_{i}=d^{W}(g^{st}_{i}(x),g_{i}^{ex}(x)), (7)
Refer to caption
(a)
Refer to caption
(b)
Refer to caption
(c)
Refer to caption
(d)
Refer to caption
(e)
Refer to caption
(f)
Refer to caption
(g)
Figure 6: Country equity distributions during extreme and stable inflationary periods for (a) Canada, (b) France, (c) Germany, (d) Italy, (e) Japan, (f) UK, (g) USA. For the purposes of visualisation, we have applied kernel density estimation to all countries’ two conditional distributions. Australian data is only available on a sparser periodicity, and Australia has been excluded for this reason.
Country equity robustness
Country ERER
Canada .012
France .009
Germany .010
Italy .012
Japan .008
UK .006
USA .006
Table 3: Country equity robustness scores computed using Wasserstein distance between two conditional distributions of log returns. The USA and the UK have equity indices that are most robust during extreme inflationary periods.

where dWd^{W} is a Wasserstein distance between our two distributions. Results are shown in Table 3. The table reveals that the countries exhibiting the smallest distance between the two distributions are the USA and UK. This implies that these countries’ equity indices exhibit the greatest robustness to inflationary pressures. By contrast, Italy and Canada display the largest Wasserstein distance between the extreme and stable periods of inflationary market dynamics. This suggests that these countries’ equity indices may display the least robustness during periods of extreme inflationary pressure.

As a subsequent experiment, we explicitly compare the difference in equity returns in high and low inflationary periods for all countries. The results highlight that the country with the smallest discrepancy in equity returns is the United States, with a difference in average monthly log returns of -0.003 between the respective inflationary periods (high and low). The country with the largest discrepancy in average log returns is Canada, which has a difference in average monthly log returns of 0.029 in periods of extremely high and low inflation. This suggests that the United States’ equity index exhibits the most similar behaviours during high and low inflationary periods. By contrast, Canada exhibits the greatest difference in equity returns during high and low inflationary periods. Examining the rest of the collection, there is significant variability as to whether average log returns perform better in anomalously high or low inflationary periods. In Australia, Japan, Germany, the UK and Canada, equity indices outperform during extremely high inflationary periods (relative to extremely low inflationary periods), while in France, Italy and the USA the opposite occurs. Further research could test this phenomenon more explicitly, on a larger collection of countries and their respective equity indices.

4 Underlying equity sector dynamics

Having explored equity index robustness, we now focus on the performance of underlying equities, grouped by sector, during a more recent period of time. Here, we study correlation behaviours of S&P 500 equities to determine differences in sector behaviours during periods of significant inflationary pressure. We partition these equities into S=11S=11 sectors, and study the evolutionary correlation dynamics of each sector. Each underlying equity is grouped into one sector from Communication Services, Consumer Discretionary, Consumer Staples, Energy, Financials, Healthcare, Industrials, Information Technology, Materials, Real Estate and Utilities. We study the evolution of intra-sector correlation dynamics over a relatively recent period that exhibited extreme inflationary behaviours, 1/1/2005-31/12/2009. Let t2=1,,T2t_{2}=1,...,T_{2} be an index in time reflecting this date range. Let each sector’s collection of underlying equity log returns be the vector 𝐑S(t2)\mathbf{R}^{S}(t_{2}), where SS indexes one of 11 sectors, and t2t_{2} denotes a specific point in time for the vector of log return values. We use a rolling window of 120 days. For each sector, we implement a time-varying correlation matrix as follows:

ΩijS(t2)=k=t2120t2(RiS(k)R¯iS)(RjS(k)R¯jS))(k=t2120t2(RiS(k)R¯iS)2k=t2120t2(RjS(k)R¯jS)2)1/2,t2{121,,T2}\displaystyle\Omega_{ij}^{S}(t_{2})=\frac{\sum_{k=t_{2}-120}^{t_{2}}(R^{S}_{i}(k)-\bar{R}^{S}_{i})(R^{S}_{j}(k)-\bar{R}^{S}_{j}))}{\left(\sum_{k=t_{2}-120}^{t_{2}}(R^{S}_{i}(k)-\bar{R}^{S}_{i})^{2}\sum_{k=t_{2}-120}^{t_{2}}(R^{S}_{j}(k)-\bar{R}^{S}_{j})^{2}\right)^{1/2}},\forall{t_{2}}\in\{121,...,T_{2}\} (8)

and study the evolution of each sector’s correlation behaviours. For each evolutionary correlation matrix ΩijS(t)\Omega_{ij}^{S}(t), we wish to compute an average sector correlation μ~S\tilde{\mu}^{S} during our window of interest. To determine the average correlation of each sector during the period 1/1/2005-31/12/2009, we compute

μ~S=2Sn(Sn1)×1303t2=130726101i<jNΩijS(t2)\displaystyle\tilde{\mu}^{S}=\frac{2}{S_{n}(S_{n}-1)\times 1303}\sum^{2610}_{t_{2}=1307}\sum_{1\leq i<j\leq N}\Omega^{S}_{ij}(t_{2}) (9)

where SnS_{n} is the number of equities in each sector, t2=1307t_{2}=1307 and t2=2610t_{2}=2610 correspond to 1/1/2005 and 31/12/2009 respectively.

Average sector correlation: 1/1/2005: 31/12/2009
Country μ~S\tilde{\mu}^{S}
Communication Services 0.40
Consumer Discretionary 0.42
Consumer Staples 0.34
Energy 0.70
Financials 0.50
Healthcare 0.32
Industrials 0.45
Information Technology 0.42
Materials 0.50
Real estate 0.66
Utilities 0.62
Table 4: Average sector correlation during the period 1/1/2005 - 31/12/2009. There is significant variability in the correlation of underlying equities based on the sector. The sectors with highest and lowest average correlation among constituent equities are energy and healthcare.

There are several reasons why we elect to study this period of time. First, rather than simply studying periods of high inflation, we wish to study equity performance during periods of significant changes in inflation (we are indifferent to the direction of such changes). In the current economic climate, the presence of inflation may be transient - and we could see a sharp spike in CPI followed by a reversion to more stable levels. Alternatively, we could see a structural increase in global inflation, and these levels may remain elevated for a prolonged period of time (such as that seen across the world in the 1970s and 1980s). Given the vast range of possibilities, we deliberately choose a country (USA) which experienced profound inflationary changes (both up and down) for a substantial period of time. Second, given the wealth of market commentary currently speculating on an imminent market crisis, we wish to study the behaviour of various equity sectors during a more prolonged crash such as the Global Financial Crisis (GFC). We avoid examining equity behaviours during the COVID-19 market crisis due to its brevity and the immediate nature of this market shock.

Figure 7 displays the rolling average correlation exhibited by each equity sector. Some general patterns are clear. First, the ranking of various sectors’ average correlation over time is remarkably consistent. That is, most sectors’ average correlation trajectory appear to deviate from any other by some translation τ\tau, which is approximately constant over time. The second most notable insight is the obvious spike in average correlation levels during the GFC market crisis. This is consistent with a significant body of prior research, demonstrating that equity behaviours become more homogenous during market crises, with average correlations increasing. That being said, the degree of this increase varies significantly based on the sector - and this warrants further investigation.

The findings of our analysis are shown in Table 4. There is a striking difference in the average correlation between difference sectors. During our period of analysis, the sector with the lowest average correlation is healthcare and the sector with the highest average correlation is energy. This is consistent with prior portfolio research, which has highlighted the typically low portfolio beta displayed by healthcare, and the high portfolio beta displayed by more volatile sectors such as energy. One interpretation of this analysis for those seeking to make stock selection decisions during periods of meaningful inflationary pressure, is the degree of homogeneity within candidate equity sectors. Sectors such as: Utilities, Real estate and Energy, all of which exhibit average correlation levels above 0.60.6, display more homogenous behaviours and will provide less opportunities for relative alpha generation. By contrast, sectors such as Consumer Staples and Healthcare, both of which have average correlation below 0.350.35, exhibit greater heterogeneity in equity behaviours. These sectors are likely to provide more opportunity for alpha generation during an environment of extreme inflationary pressure.

Refer to caption
Figure 7: Average time-varying correlation for 11 sectors analyzed. Sectors considered include (a) Communication services, (b) Consumer Discretionary, (c) Consumer Staples, (d) Energy, (e) Financials, (f) Healthcare, (g) Industrials, (h) Information Technology, (i) Materials, (j) Real Estate, (k) Utilities.

5 Dynamic multi-asset optimization

We now explore a similar theme within the context of managing a portfolio of different asset classes. One notable development over the past several years has been the proliferation of digital currencies such as Bitcoin. Although investments in many coins such as Bitcoin and Ethereum provided very strong returns over the past several years, this was accompanied by significant volatility. But just how beneficial was cryptocurrency exposure to an investor focused on maximizing their risk-adjusted returns? We now try to address this pertinent and topical question. In particular, we study (in the context of managing a multi-asset financial portfolio), which asset classes would have provided the greatest utility to a core equity portfolio over the past \sim 5 years? We take a collection of a=1,,Aa=1,...,A daily asset class returns where we let A=7A=7, and study their behaviour over a date range from 01/01/2016 - 30/06/2021. We index these date t3=1,,T3t_{3}=1,...,T_{3}. The 7 securities studied include S&P 500 equities, the Bloomberg LNG Index, London Metals Exchange, S&P 500 energy index, CRB Commodity Index, Bloomberg REIT index and Bitcoin (as a representative of the cryptocurrency asset class). Let each asset classes’ daily log returns be Ra(t3)R_{a}(t_{3}), where t3=1,,T3t_{3}=1,...,T_{3} references a point in time t3t_{3} during our window of analysis. We implement a time-varying portfolio optimization, where we use a rolling window of 250250 days, and study the resulting functions of time-varying optimal portfolio weights. We carefully select a 250250 day rolling window to simulate a portfolio manager who examines data over the past year to inform trading decisions. In the optimization problem below, we seek to learn portfolio weights w1,,wAw_{1},...,w_{A} such that we maximize the following

𝔼(Rp(t3))Rfσp2(t3),t3{251,,T3}, where\displaystyle\frac{\operatorname{\mathbb{E}}(R_{p}(t_{3}))-R_{f}}{\sigma^{2}_{p}(t_{3})},\forall t_{3}\in\{251,...,T_{3}\},\text{ where} (10)
𝔼(Rp)=a=1AwaRa(t3), and\displaystyle\operatorname{\mathbb{E}}(R_{p})=\sum^{A}_{a=1}w_{a}R_{a}(t_{3}),\text{ and} (11)
σp2(t3)=𝒘(t3)TΣ(t3)𝒘(t3).\displaystyle\sigma^{2}_{p}(t_{3})=\bm{w}(t_{3})^{T}\Sigma(t_{3})\bm{w}(t_{3}). (12)

In the above expression, 𝒘(t3)\bm{w}(t_{3}) is a vector of portfolio weights at each point in time, Σij(t3)=cov(Ra,i(t3),Ra,j(t3))\Sigma_{ij}(t_{3})=\text{cov}(R_{a,i}(t_{3}),R_{a,j}(t_{3})) is the time-varying covariance matrix between rolling log returns and RfR_{f} is the current risk free rate which we set to 0.25%. The constraints in this optimisation are as follows. We hold the S&P 500 weight constant a1=0.4a_{1}=0.4, and restrict weights for the remaining portfolio assets 0.025wa0.3,a=2,,A0.025\leq w_{a}\leq 0.3,a=2,...,A and assume that the portfolio is always fully invested in a long-only capacity a=1Awa=1\sum^{A}_{a=1}w_{a}=1. At each point in time, the optimization generates a vector of optimal portfolio weights which we denote wa(t3)w_{a}^{*}(t_{3}). For each asset class, we compute the mean and variance of the portfolio weights which we denote μaw\mu_{a}^{w^{*}} and σa2,w\sigma_{a}^{2,w^{*}}. These results are shown in table 8 which displays the mean and variance of the time-varying optimal portfolio weights among our collection of A=7A=7 representative securities/indices. In conjunction with Figure 8 which displays the optimal weight trajectories, there are several interesting findings.

Dynamic optimal portfolio weights
Securities μaw\mu_{a}^{w^{*}} σa2,w\sigma_{a}^{2,w^{*}}
S&P 500 equities 0.4 n/a
LNG Index 0.045 0.004
London Metals Exchange 0.13 0.009
Energy index 0.049 0.0045
CRB Commodity index 0.116 0.0056
REIT index 0.108 0.011
Bitcoin 0.149 0.0088
Table 5: Mean and variance of time-varying optimal portfolio weight trajectories. Among our collection of more traditional hedging instruments, Bitcoin (cryptocurrency) was the most effective complement to a core equity portfolio in maximizing our portfolio’s Sharpe Ratio.

First, Bitcoin has the highest average optimal weight of 14.9%. This suggests that for an investor with a core allocation to equities over the past 5.5 years, cryptocurrency was the most important asset in increasing the overall portfolio Sharpe Ratio. This was followed closely by the London Metals Exchange index with a total weight of 13%. Although having an average optimal weight of 10.8%, the REIT index exhibited the highest variance in its optimal weight over time. This suggests that REITs exposure may have to be monitored more closely than other asset classes, and may require more frequent rebalancing. Finally, although the LNG Index’s average optimal allocation was only 4.5%, it exhibited the lowest variance in its weight trajectory. This indicates that for a security exposure such as this, portfolio managers may require less trading and rebalancing.

Refer to caption
Figure 8: Time-varying optimal weight for each of our 7 candidate assets: (a) Bloomberg LNG Index, (b) London Metal Exchange, (c) S&P 500 Energy Index, (d) CRB Commodity Index, (e) Bloomberg REIT Index, (f) Bitcoin, (g) S&P 500.

6 Conclusion

This paper uses various data, mathematical techniques and frameworks to provide a holistic view of inflation, and its impact on equity markets and (predominantly) equity market investors. We believe this is the first work to explore this topic from a macroeconomic, financial market dynamics and portfolio optimization perspective.

In Section 3 we study the similarity of evolutionary CPI inflation patterns between a collection of 8 critically important economies. We demonstrate that Australia and Japan exhibited the most anomalous CPI trajectories, and this was primarily due to their hyper-inflationary onset preceding that of other major economies in the 1970s and 1980s. Our eigenvalue analysis demonstrated the strength of collective behaviours among country inflation patterns, reinforcing the systemic and largely ubiquitous nature of global inflation. We then introduce an inner-product based optimization framework to identify countries that are most central in global inflation. Japan is shown to be the least central (or most anomalous) country with respect to the timing of its inflationary patterns, with the majority of other countries forming one predominant cluster. Finally, we apply the Wasserstein distance between two conditional distributions of country equity index log returns to demonstrate that the USA and UK equity indices possess the most robustness during times of immense inflationary pressure.

In Section 4 we partition S&P 500 equities into 11 disjoint sectors, and study the evolutionary correlation behaviours during a period of extreme inflationary behaviours. Our analysis reveals significant differences between various equity sectors’ correlation during times of inflationary crisis. Consumer staples and Healthcare exhibited the lowest average correlation of 0.340.34 and 0.320.32 respectively, while energy and real estate exhibited the highest average correlation among constituents with 0.70.7 and 0.660.66 respectively. Sectors with lower average correlation levels may produce more opportunities for alpha generation, relative to the remaining collection of equities, given the heterogeneity in dynamics. By contrast, sectors with higher average levels of correlation may provide less opportunities for portfolio alpha, with more homogeneity among the collection.

Finally in Section 5 we turn our attention to core equity investors over the past several years who are seeking favourable portfolio diversification. Among our collection of candidate asset classes one could use to diversify portfolio returns, cryptocurrency (Bitcoin) was identified as the most effective complement to a core exposure of equities. The second most effective asset class was the London Metals Exchange. In studying the resulting time-varying optimal weight trajectories, we computed the variance of each function as an indication as to how frequently asset classes may have to be rebalanced. Our analysis indicates that REITS exhibit the highest variance, and may have to be monitored most closely by portfolio managers.

There are several avenues for further work. First, one could explore similar questions using different types of mathematical and statistical analysis. One possible methodology for the study of country inflation trajectories could be the application of multiplex network analysis, for example. Second, one could apply the techniques introduced in this manuscript to answer more general macroeconomic questions, not focused on inflation specifically. Specific modelling techniques, such as our time-varying linear regression, could be replaced with moving average or autoregressive moving average models. Finally, our equity sector and portfolio optimization analysis could be applied on a more varied collection of underlying securities. Equity securities could include equities outside the USA, and the optimization section could include more asset classes, such as fixed income, for example.

7 Appendix

7.1 Optimization assets

Below we include a list of the assets used in the optimization simulations.

Asset definition
Asset Definition
S&P 500 equity index (SPX) Collection of 500 US equities.
LNG Index (BCOMNG) Bloomberg Natural Gas subindex.
London Metals Exchange (LMEX) Narrowly focused basket of industrial metals (aluminium, copper, nickel, zinc, tin, lead).
Energy Index (S5ENRS) S&P 500 Energy index.
CRB Commodity Index (CRB CMDT) Index consisting of 19 commodity futures (eg. Oil, hogs, Cattle, Silver, Soybeans, etc.).
REIT Index (BBREIT) Collection of US REITs.
Bitcoin Bitcoin share price quoted in USD.
Table 6: Assets used in Section 5 and their definition.

Below we test the robustness of our optimization finding under a variety of constant weights for the S&P 500 equity index allocation. We set this weight equal to 0.2, 0.3, 0.4 and 0.5 in 4 different simulations. In all scenarios, cryptocurrency was the largest position in maximizing portfolio Sharpe Ratio.

Dynamic optimal portfolio weights
Securities μaw\mu_{a}^{w^{*}} μaw\mu_{a}^{w^{*}} μaw\mu_{a}^{w^{*}} μaw\mu_{a}^{w^{*}}
S&P 500 equities 0.2 0.3 0.4 0.5
LNG Index 0.05 0.048 0.045 0.042
London Metals Exchange 0.18 0.16 0.13 0.11
Energy index 0.06 0.05 0.048 0.04
CRB Commodity index 0.175 0.148 0.12 0.075
REIT index 0.14 0.12 0.10 0.09
Bitcoin 0.184 0.164 0.148 0.14
Table 7: Mean weight of portfolio assets under various assumptions for fixed S&P 500 equity weight. Under all scenarios cryptocurrency (Bitcoin) was the most effective portfolio complement in maximizing the portfolio Sharpe Ratio.
Refer to caption
(a)
Refer to caption
(b)
Refer to caption
(c)
Refer to caption
(d)
Figure 9: Time-varying optimal weight for all 7 portfolio assets with varying levels of the constant S&P 500 equity allocation. The S&P 500 equity weight ranges in (a) 0.2, (b) 0.3, (c) 0.4 and (d) 0.5.

References

  • Siu [2011] T. K. Siu, Long-term strategic asset allocation with inflation risk and regime switching, Quantitative Finance 11 (2011) 1565–1580. doi:10.1080/14697680903055588.
  • Mishkin [1990] F. S. Mishkin, The information in the longer maturity term structure about future inflation, The Quarterly Journal of Economics 105 (1990) 815. doi:10.2307/2937901.
  • Fama [1990] E. F. Fama, Term-structure forecasts of interest rates, inflation and real returns, Journal of Monetary Economics 25 (1990) 59–76. doi:10.1016/0304-3932(90)90045-6.
  • Mundell [1963] R. Mundell, Inflation and real interest, Journal of Political Economy 71 (1963) 280–283. doi:10.1086/258771.
  • Tobin [1965] J. Tobin, Money and economic growth, Econometrica 33 (1965) 671. doi:10.2307/1910352.
  • Fenn et al. [2011] D. J. Fenn, M. A. Porter, S. Williams, M. McDonald, N. F. Johnson, N. S. Jones, Temporal evolution of financial-market correlations, Physical Review E 84 (2011). doi:10.1103/physreve.84.026109.
  • Laloux et al. [1999] L. Laloux, P. Cizeau, J.-P. Bouchaud, M. Potters, Noise dressing of financial correlation matrices, Physical Review Letters 83 (1999) 1467–1470. doi:10.1103/physrevlett.83.1467.
  • Münnix et al. [2012] M. C. Münnix, T. Shimada, R. Schäfer, F. Leyvraz, T. H. Seligman, T. Guhr, H. E. Stanley, Identifying states of a financial market, Scientific Reports 2 (2012). doi:10.1038/srep00644.
  • Kim and Jeong [2005] D.-H. Kim, H. Jeong, Systematic analysis of group identification in stock markets, Physical Review E 72 (2005). doi:10.1103/physreve.72.046133.
  • Pan and Sinha [2007] R. K. Pan, S. Sinha, Collective behavior of stock price movements in an emerging market, Physical Review E 76 (2007). doi:10.1103/physreve.76.046116.
  • Wilcox and Gebbie [2007] D. Wilcox, T. Gebbie, An analysis of cross-correlations in an emerging market, Physica A: Statistical Mechanics and its Applications 375 (2007) 584–598. doi:10.1016/j.physa.2006.10.030.
  • Heckens et al. [2020] A. J. Heckens, S. M. Krause, T. Guhr, Uncovering the dynamics of correlation structures relative to the collective market motion, Journal of Statistical Mechanics: Theory and Experiment 2020 (2020) 103402. doi:10.1088/1742-5468/abb6e2.
  • James and Menzies [2021] N. James, M. Menzies, Association between COVID-19 cases and international equity indices, Physica D: Nonlinear Phenomena 417 (2021) 132809. doi:10.1016/j.physd.2020.132809.
  • James et al. [2021] N. James, M. Menzies, J. Chan, Changes to the extreme and erratic behaviour of cryptocurrencies during COVID-19, Physica A: Statistical Mechanics and its Applications 565 (2021) 125581. URL: https://doi.org/10.1016/j.physa.2020.125581. doi:10.1016/j.physa.2020.125581.
  • James [2021] N. James, Dynamics, behaviours, and anomaly persistence in cryptocurrencies and equities surrounding COVID-19, Physica A: Statistical Mechanics and its Applications 570 (2021) 125831. URL: https://doi.org/10.1016/j.physa.2021.125831. doi:10.1016/j.physa.2021.125831.
  • Ferreira et al. [2020] P. Ferreira, L. Kristoufek, E. J. de Area Leão Pereira, DCCA and DMCA correlations of cryptocurrency markets, Physica A: Statistical Mechanics and its Applications 545 (2020) 123803. doi:10.1016/j.physa.2019.123803.
  • Ausloos [2000] M. Ausloos, Statistical physics in foreign exchange currency and stock markets, Physica A: Statistical Mechanics and its Applications 285 (2000) 48–65. doi:10.1016/s0378-4371(00)00271-5.
  • Driessen et al. [2003] J. Driessen, B. Melenberg, T. Nijman, Common factors in international bond returns, Journal of International Money and Finance 22 (2003) 629–656. doi:10.1016/s0261-5606(03)00046-9.
  • Chu et al. [2015] J. Chu, S. Nadarajah, S. Chan, Statistical analysis of the exchange rate of Bitcoin, PLOS ONE 10 (2015) e0133678. doi:10.1371/journal.pone.0133678.
  • Lahmiri and Bekiros [2018] S. Lahmiri, S. Bekiros, Chaos, randomness and multi-fractality in Bitcoin market, Chaos, Solitons & Fractals 106 (2018) 28–34. doi:10.1016/j.chaos.2017.11.005.
  • Kondor et al. [2014] D. Kondor, M. Pósfai, I. Csabai, G. Vattay, Do the rich get richer? an empirical analysis of the Bitcoin transaction network, PLoS ONE 9 (2014) e86197. doi:10.1371/journal.pone.0086197.
  • Bariviera et al. [2017] A. F. Bariviera, M. J. Basgall, W. Hasperué, M. Naiouf, Some stylized facts of the Bitcoin market, Physica A: Statistical Mechanics and its Applications 484 (2017) 82–90. doi:10.1016/j.physa.2017.04.159.
  • Alvarez-Ramirez et al. [2018] J. Alvarez-Ramirez, E. Rodriguez, C. Ibarra-Valdez, Long-range correlations and asymmetry in the Bitcoin market, Physica A: Statistical Mechanics and its Applications 492 (2018) 948–955. doi:10.1016/j.physa.2017.11.025.
  • Stosic et al. [2019a] D. Stosic, D. Stosic, T. B. Ludermir, T. Stosic, Multifractal behavior of price and volume changes in the cryptocurrency market, Physica A: Statistical Mechanics and its Applications 520 (2019a) 54–61. doi:10.1016/j.physa.2018.12.038.
  • Stosic et al. [2019b] D. Stosic, D. Stosic, T. B. Ludermir, T. Stosic, Exploring disorder and complexity in the cryptocurrency space, Physica A: Statistical Mechanics and its Applications 525 (2019b) 548–556. doi:10.1016/j.physa.2019.03.091.
  • Manavi et al. [2020] S. A. Manavi, G. Jafari, S. Rouhani, M. Ausloos, Demythifying the belief in cryptocurrencies decentralized aspects. A study of cryptocurrencies time cross-correlations with common currencies, commodities and financial indices, Physica A: Statistical Mechanics and its Applications 556 (2020) 124759. doi:10.1016/j.physa.2020.124759.
  • Drożdż et al. [2018] S. Drożdż, R. Gębarowski, L. Minati, P. Oświęcimka, M. Wątorek, Bitcoin market route to maturity? Evidence from return fluctuations, temporal correlations and multiscaling effects, Chaos: An Interdisciplinary Journal of Nonlinear Science 28 (2018) 071101. doi:10.1063/1.5036517.
  • Drożdż et al. [2019] S. Drożdż, L. Minati, P. Oświęcimka, M. Stanuszek, M. Wątorek, Signatures of the crypto-currency market decoupling from the forex, Future Internet 11 (2019) 154. doi:10.3390/fi11070154.
  • Drożdż et al. [2020] S. Drożdż, L. Minati, P. Oświęcimka, M. Stanuszek, M. Wątorek, Competition of noise and collectivity in global cryptocurrency trading: Route to a self-contained market, Chaos: An Interdisciplinary Journal of Nonlinear Science 30 (2020) 023122. doi:10.1063/1.5139634.
  • Gębarowski et al. [2019] R. Gębarowski, P. Oświęcimka, M. Wątorek, S. Drożdż, Detecting correlations and triangular arbitrage opportunities in the forex by means of multifractal detrended cross-correlations analysis, Nonlinear Dynamics 98 (2019) 2349–2364. doi:10.1007/s11071-019-05335-5.
  • Wątorek et al. [2020] M. Wątorek, S. Drożdż, J. Kwapień, L. Minati, P. Oświęcimka, M. Stanuszek, Multiscale characteristics of the emerging global cryptocurrency market, Physics Reports (2020). doi:10.1016/j.physrep.2020.10.005.
  • Fil and Kristoufek [2020] M. Fil, L. Kristoufek, Pairs trading in cryptocurrency markets, IEEE Access 8 (2020) 172644–172651. doi:10.1109/access.2020.3024619.
  • Markowitz [1952] H. Markowitz, Portfolio selection, The Journal of Finance 7 (1952) 77. URL: https://doi.org/10.2307/2975974. doi:10.2307/2975974.
  • Sharpe [1966] W. F. Sharpe, Mutual fund performance, The Journal of Business 39 (1966) 119–138. URL: https://doi.org/10.1086/294846. doi:10.1086/294846.
  • Almahdi and Yang [2017] S. Almahdi, S. Y. Yang, An adaptive portfolio trading system: A risk-return portfolio optimization using recurrent reinforcement learning with expected maximum drawdown, Expert Systems with Applications 87 (2017) 267–279. doi:10.1016/j.eswa.2017.06.023.
  • Calvo et al. [2014] C. Calvo, C. Ivorra, V. Liern, Fuzzy portfolio selection with non-financial goals: exploring the efficient frontier, Annals of Operations Research 245 (2014) 31–46. doi:10.1007/s10479-014-1561-2.
  • Soleimani et al. [2009] H. Soleimani, H. R. Golmakani, M. H. Salimi, Markowitz-based portfolio selection with minimum transaction lots, cardinality constraints and regarding sector capitalization using genetic algorithm, Expert Systems with Applications 36 (2009) 5058–5063. doi:10.1016/j.eswa.2008.06.007.
  • Vercher et al. [2007] E. Vercher, J. D. Bermúdez, J. V. Segura, Fuzzy portfolio optimization under downside risk measures, Fuzzy Sets and Systems 158 (2007) 769–782. doi:10.1016/j.fss.2006.10.026.
  • Bhansali [2007] V. Bhansali, Putting economics (back) into quantitative models, The Journal of Portfolio Management 33 (2007) 63–76. doi:10.3905/jpm.2007.684755.
  • Moody and Saffell [2001] J. Moody, M. Saffell, Learning to trade via direct reinforcement, IEEE Transactions on Neural Networks 12 (2001) 875–889. doi:10.1109/72.935097.
  • Corbet et al. [2020] S. Corbet, C. Larkin, B. Lucey, The contagion effects of the COVID-19 pandemic: Evidence from gold and cryptocurrencies, Finance Research Letters 35 (2020) 101554. doi:10.1016/j.frl.2020.101554.
  • Conlon and McGee [2020] T. Conlon, R. McGee, Safe haven or risky hazard? bitcoin during the Covid-19 bear market, Finance Research Letters 35 (2020) 101607. doi:10.1016/j.frl.2020.101607.
  • Conlon et al. [2020] T. Conlon, S. Corbet, R. J. McGee, Are cryptocurrencies a safe haven for equity markets? an international perspective from the COVID-19 pandemic, Research in International Business and Finance 54 (2020) 101248. doi:10.1016/j.ribaf.2020.101248.
  • Ji et al. [2020] Q. Ji, D. Zhang, Y. Zhao, Searching for safe-haven assets during the COVID-19 pandemic, International Review of Financial Analysis 71 (2020) 101526. doi:10.1016/j.irfa.2020.101526.
  • Sigaki et al. [2019] H. Y. D. Sigaki, M. Perc, H. V. Ribeiro, Clustering patterns in efficiency and the coming-of-age of the cryptocurrency market, Scientific Reports 9 (2019). doi:10.1038/s41598-018-37773-3.
  • Jusup et al. [2022] M. Jusup, P. Holme, K. Kanazawa, M. Takayasu, I. Romic, Z. Wang, S. Gecek, T. Lipic, B. Podobnik, L. Wang, W. Luo, T. Klanjscek, J. Fan, S. Boccaletti, M. Perc, Social physics, 2022. arXiv:2110.01866.
  • Perc [2019] M. Perc, The social physics collective, Scientific Reports 9 (2019). doi:10.1038/s41598-019-53300-4.
  • Kantorovich [1960] L. V. Kantorovich, Mathematical methods of organizing and planning production, Management Science 6 (1960) 366–422. doi:10.1287/mnsc.6.4.366.