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

Evolutionary correlation, regime switching, spectral dynamics and optimal trading strategies for cryptocurrencies and equities

Nick James nick.james@unimelb.edu.au School of Mathematics and Statistics, University of Melbourne, Victoria, Australia
Abstract

This paper uses new and recently established methodologies to study the evolutionary dynamics of the cryptocurrency market, and compares the findings with that of the equity market. We begin by applying random matrix theory and principal components analysis (PCA) to correlation matrices of both collections, highlighting clear differences in the eigenspectra exhibited. We then explore the heterogeneity of both asset classes, studying the time-varying dynamics of underlying sector behaviours, and determine the collective similarity within each collection. We then turn to a study of structural break dynamics and evolutionary power spectra, where we quantify the collective affinity in structural breaks and evolutionary behaviours of underlying sector time series. Finally, we implement two algorithms simulating ‘portfolio choice’ dynamics to compare the effectiveness of stock selection and sector allocation in cryptocurrency portfolios. There, we highlight the importance of both endeavours and comment on noteworthy implications for cryptocurrency portfolio management.

keywords:
Time series analysis , Random Matrix Theory , Regime switching , Trading strategies , Financial markets , Cryptocurrency
journal: Elsevier

1 Introduction

Since its inception, the cryptocurrency market has seen rapid growth in its total market capitalisation and gradual changes in its underlying dynamics. The increase in the cryptocurrency market’s maturity has been accompanied by rapid growth in the number of cryptocurrencies available to trade, and a general rise in the level of sophistication across the asset class. In this paper, we seek to explore the evolution of the cryptocurrency market, and compare its behaviours to the larger, and more sophisticated equity market.

A variety of communities have been interested in the study of financial market dynamics, and the evolutionary changes in structural behaviours, including researchers in applied mathematics, econometrics and econophysics [1, 2, 3]. A variety of mathematical techniques have been used to analyze the evolutionary dynamics exhibited by such financial assets including principal components analysis (PCA) [2, 4, 5, 6], random matrix theory [7, 8, 9, 10, 11, 1], various types of clustering [12, 13], structural break detection [14, 15] and a variety of mathematical and statistical modelling techniques [16, 17]. These techniques have been used on a wide variety of asset classes such as equities [6], foreign exchange [18] and fixed income [19]. In more recent years, researchers have explored the application of such techniques as a means of understanding the burgeoning cryptocurrency market. Various aspects of the market have been studied including time-varying behaviours [14, 15, 20, 21, 22, 23, 24, 25], Bitcoin’s price behaviours [26, 27, 28, 29, 30], fractal-type dynamics and pattern formation, [31, 32, 33, 17], and trading strategies [34].

The study of nonstationarity behaviours and the time-varying nature of (autoregressive) model parameters for phenomena such as volatility, has been of great interest to financial markets researchers for many years. Such research dates back to Autoregressive conditionally heteroskedastic (ARCH) models [35], Generalized ARCH (GARCH) [36] and stochastic volatility models [37, 38, 39]. In more recent years there have been various extensions and adaptations to the core models, allowing for the specific consideration of various time series features. A non-exhaustive collection of model extensions may mention EGARCH [40], GJR-GARCH [41], T-GARCH [42] and T-SV [43], Markov Switching GARCH [44, 45, 46] and MS-SV [47]. Bayesian estimation has also been prominent in this area [48, 49, 50]. Such modelling procedures and various adaptations have been widely applied to study the properties of cryptocurrencies and other asset classes [51, 52, 53, 54, 55, 56]. In this work specifically, we build on the concept of local stationarity in stochastic processes [57, 58] and time-varying spectral density estimation [59, 60, 61, 62], to study the propensity for regime switching, and the associated changes in periodic behaviour, in cryptocurrency and equity sectors.

Portfolio optimization, and more generally the concept of portfolio construction, has been a topic of considerable interest to the mathematics and quantitative finance communities [63, 64, 65, 66, 67, 68, 69, 70]. Although there is literature exploring the effectiveness of cryptocurrencies as safe-haven portfolio assets, we are unaware of any work explicitly testing for comparative effectiveness in alpha generation when contrasting security selection, and sector allocation. In this work, we introduce two algorithmic procedures to test this hypothesis in the cryptocurrency market.

This paper is structured as follows. In Section 3 we apply random matrix theory in stationary and time-varying contexts, to reveal differences in the impact of the ‘market factor’ and other factors on the correlations of both collections’ security returns. In Section 4 we extend our study of asset class heterogeneity, and compare the evolution in collective strength among underlying cryptocurrency and equity sectors. In Section 5 we contrast the similarity in underlying sectors’ structural break dynamics, and further distinguish between nonstationarity and long-range dependence by quantifying the distance between representative time series’ time-varying power spectra. Finally in Section 6 we introduce two algorithms to determine the relative importance of security selection and sector allocation when assembling cryptocurrency portfolios. In Section 7 we conclude.

2 Data

In this paper we study cryptocurrency and equity prices between 01-01-2019 to 01-12-2021. We study a collection of 45 unique cryptocurrency securities and 77 unique equity securities. Cryptocurrencies are partitioned into 8 sectors (Centralized exchange, Collectibles/NFTs, Decentralized finance, Platform, Smart contracts, Store of value, Wallet) and equities are partitioned into 11 sectors (Communication services, Consumer discretionary, Consumer staples, Energy, Financials, Healthcare, Industrials, Information technology, Materials, Real estate, Utilities). Equity securities are classified based on their Global Industry Classification Standard (GICS) sector. Cryptocurrency data is sourced from Coinmarketcap.com (https://coinmarketcap.com/) and equity data is sourced from Yahoo Finance (https://au.finance.yahoo.com/). Some cryptocurrencies may be classified in more than one sector (as is consistent with the Coinmarketcap.com classification process).

3 Time-varying random matrix theory and spectral dynamics

Throughout this paper we deal with two central objects of study, a collection of cryptocurrency time series and a collection of equity time series, where the underlying securities are uniformly sampled from a variety of respective sectors (see B for the complete list of securities). Given that cryptocurrencies trade 7 days per week (unlike equities), although both collections span 01-01-2019 to 01-12-2021, the length of time series within each collection will be different. Let our cryptocurrency collection be indexed t1=1,,T1t_{1}=1,...,T_{1} where T1=1065T_{1}=1065 and our equity collection be indexed t2=1,,T2t_{2}=1,...,T_{2} where T2=733T_{2}=733. Let pic(t1)p^{c}_{i}(t_{1}) and pje(t2)p^{e}_{j}(t_{2}) refer to the multivariate time series of cryptocurrency daily closing prices and equity adjusted daily closing prices respectively i=1,,Ni=1,...,N and j=1,,Kj=1,...,K. A collection of representative price time series for our cryptocurrency and equity collections are shown in Figures 1 and 2 respectively. We generate two multivariate time series of log returns Ric(t1)R^{c}_{i}(t_{1}) and Rje(t2)R^{e}_{j}(t_{2}) as follows:

Ric(t1)\displaystyle R^{c}_{i}{(t_{1})} =log(pic(t1)pic(t11)),\displaystyle=\log\left(\frac{p^{c}_{i}{(t_{1})}}{p^{c}_{i}{(t_{1}-1)}}\right), (1)
Rje(t2)\displaystyle R^{e}_{j}{(t_{2})} =log(pje(t2)pje(t21)).\displaystyle=\log\left(\frac{p^{e}_{j}{(t_{2})}}{p^{e}_{j}{(t_{2}-1)}}\right). (2)

We standardise cryptocurrency and equity log returns as follows. For our cryptocurrency collection, let R~ic(t1)=[Ric(t1)Ric]/σ(Ric)\tilde{R}^{c}_{i}(t_{1})=[R^{c}_{i}(t_{1})-\langle R^{c}_{i}\rangle]/\sigma(R^{c}_{i}), where .\langle.\rangle is an average over a candidate period of time and σ(.)\sigma(.) is a standard deviation operator, computed over the same interval. Equity log returns are standardised analogously. Over any candidate smoothing window SS, the cryptocurrency and equity correlation matrices, Ωc\Omega^{c} and Ωe\Omega^{e} are defined as follows:

Ωc(t1)=1SR~c(t1)R~cT(t1),t1S,,T1,\displaystyle\Omega^{c}(t_{1})=\frac{1}{S}\tilde{R}^{c}(t_{1})\tilde{R}^{cT}(t_{1}),\forall t_{1}\in S,...,T_{1}, (3)
Ωe(t2)=1SR~e(t2)R~eT(t2),t2S,,T2.\displaystyle\Omega^{e}(t_{2})=\frac{1}{S}\tilde{R}^{e}(t_{2})\tilde{R}^{eT}(t_{2}),\forall t_{2}\in S,...,T_{2}. (4)

In our proceeding experiments, we set S=150S=150, and study the evolution in structure of the cryptocurrency and equity correlation matrix time series. Previous literature has outlined potential sensitivity in findings to the parameter SS [1]. If SS is too small, although correlation structure can be computed dynamically, results may be noisy. Alternatively if SS is too large, our time-evolving model may be insensitive to abrupt changes in behaviours (such as the COVID-19 market crisis). We also wish to ensure that Q=S/N1Q=S/N\geq 1 in subsequent experiments.

Refer to caption
(a)
Refer to caption
(b)
Refer to caption
(c)
Refer to caption
(d)
Figure 1: A sample of cryptocurrency price time series for (a) Bitcoin (b) Ethereum, (c) Binance Coin, (d) NEO.
Refer to caption
(a)
Refer to caption
(b)
Refer to caption
(c)
Refer to caption
(d)
Figure 2: A sample of equity price time series for (a) Amgen (b) Amazon, (c) JP Morgan, (d) Nike.

We apply principal components analysis (PCA) to our two time-varying correlation matrices. In particular, we wish to study the evolution in the magnitude and explanatory variance of the correlation matrix’s first eigenvalue. At every point in time each correlation matrix generates a sequence of ordered eigenvalues λ1c(t1),,λNc(t1)\lambda^{c}_{1}(t_{1}),...,\lambda^{c}_{N}(t_{1}) and λ1e(t2),,λKe(t2)\lambda^{e}_{1}(t_{2}),...,\lambda^{e}_{K}(t_{2}). At each point in time we normalize the eigenvalues in each collection by defining λc~i(t1)=λic(t1)z=1Nλzc(t1)\tilde{\lambda^{c}}_{i}(t_{1})=\frac{\lambda^{c}_{i}(t_{1})}{\sum^{N}_{z=1}\lambda^{c}_{z}(t_{1})} and λe~j(t2)=λje(t2)z=1Kλze(t2)\tilde{\lambda^{e}}_{j}(t_{2})=\frac{\lambda^{e}_{j}(t_{2})}{\sum^{K}_{z=1}\lambda^{e}_{z}(t_{2})}.

We now wish to compare the eigenspectrum and properties of our equity and cryptocurrency correlation matrices with that of random time series. Before proceeding, we provide some background on this problem from the literature. As pointed out in [71], decoupling correlation matrix behaviours into actual dependencies among the degrees of freedom and random correlation structure is a problem of significant difficulty. Many authors have employed random matrix theory to address this problem [72, 73]. The Marchenko-Pastur law is widely applied in the econophysics literature. We provide a basic overview of the law, following [74], however more detailed exposition can be found in the original paper [75]. For any random matrix XX of dimension p×mp\times m where matrix elements are randomly samples from a standard normal distribution, the limiting distribution of the eigenspectrum, which we denote λj(W)\lambda_{j}(W), of the sample covariance matrix W=1mXXTW=\frac{1}{m}XX^{T} is determined by the Marchenko-Pastur law. This result is shown to be valid where the dimensions of our matrix XX increase to infinity, but the aspect ratio (the ratio of an objects size in different dimensions) remains constant. That is, as pp\rightarrow\infty and p/mλ(0,)p/m\rightarrow\lambda\in(0,\infty). From here, one can state with absolute certainty that the empirical spectral distribution of the p×pp\times p matrix WW will converge to a deterministic distribution, referred to as the Marchenko-Pastur law, with a single parameter λ\lambda. In the case λ\lambda lies in (0,1)(0,1) the following will hold for all xx\in\mathbb{R}:

FW(x):=1p#{1jp:λj(W)x}xfλ(t)𝑑t\displaystyle F^{W}(x):=\frac{1}{p}\#\{1\leq j\leq p:\lambda_{j}(W)\leq x\}\rightarrow\int^{x}_{\infty}f_{\lambda}(t)\,dt\ (5)

where fλf_{\lambda} is the Marchenko-Pastur density

fλ(x)=12πλx[(λ+x)(xλ)]+,\displaystyle f_{\lambda}(x)=\frac{1}{2\pi\lambda x}\sqrt{[(\lambda_{+}-x)(x-\lambda_{-})]_{+}}, (6)

where λ±=(1±λ)2\lambda_{\pm}=(1\pm\sqrt{\lambda})^{2}. We now turn to our application. The resulting correlation matrix among a collection of NN uncorrelated time series each of length SS where elements are drawn randomly from a Gaussian distribution, belongs to a class referred to as Wishart matrices [76]. Assuming the respective constraints Qc=S/N1Q^{c}=S/N\geq 1 and Qe=S/K1Q^{e}=S/K\geq 1 are satisfied for both collections, in the limit where N,KN,K\to\infty and SS\to\infty, the probability density functions of our two eigenvalue collections are given by:

p(λc(t1))=Qc2πσ2(R~c(t1))(λ+c(t1)λc(t1))(λc(t1)λc(t1))λc(t1)\displaystyle p(\lambda^{c}(t_{1}))=\frac{Q^{c}}{2\pi\sigma^{2}(\tilde{R}^{c}(t_{1}))}\frac{\sqrt{(\lambda^{c}_{+}(t_{1})-\lambda^{c}(t_{1}))(\lambda^{c}_{-}(t_{1})-\lambda^{c}(t_{1}))}}{\lambda^{c}(t_{1})} (7)
p(λe(t2))=Qe2πσ2(R~e(t2))(λ+e(t2)λe(t2))(λe(t2)λe(t2))λe(t2).\displaystyle p(\lambda^{e}(t_{2}))=\frac{Q^{e}}{2\pi\sigma^{2}(\tilde{R}^{e}(t_{2}))}\frac{\sqrt{(\lambda^{e}_{+}(t_{2})-\lambda^{e}(t_{2}))(\lambda^{e}_{-}(t_{2})-\lambda^{e}(t_{2}))}}{\lambda^{e}(t_{2})}. (8)

where σ2(R~c)\sigma^{2}(\tilde{R}^{c}) and σ2(R~e)\sigma^{2}(\tilde{R}^{e}) refer to the variance among the elements of R~c\tilde{R}^{c} and R~e\tilde{R}^{e} at each point in time t1t_{1} and t2t_{2}. λ+\lambda_{+} and λ\lambda_{-} refer to the maximum and minimum eigenvalues of the respective correlation matrices, and are computed as follows:

λ±c(t1)=σ2(R~c(t1))(1+1Qc±21Qc),\displaystyle\lambda^{c}_{\pm}(t_{1})=\sigma^{2}(\tilde{R}^{c}(t_{1}))\bigg{(}1+\frac{1}{Q^{c}}\pm 2\sqrt{\frac{1}{Q^{c}}}\bigg{)}, (9)
λ±e(t2)=σ2(R~e(t2))(1+1Qe±21Qe).\displaystyle\lambda^{e}_{\pm}(t_{2})=\sigma^{2}(\tilde{R}^{e}(t_{2}))\bigg{(}1+\frac{1}{Q^{e}}\pm 2\sqrt{\frac{1}{Q^{e}}}\bigg{)}. (10)

As we apply our random matrix analysis in a rolling capacity, we must ensure that QcQ^{c} and QeQ^{e} are greater than 1 with S=150S=150. Given that our cryptocurrency collection consists of N=45N=45 unique cryptocurrencies and our equity collection consists of K=77K=77 unique equities, Qc=3.3Q^{c}=3.3 and Qe=1.95Q^{e}=1.95. It is well understood that when Q=1Q=1, the minimum eigenvalue is λ=0\lambda_{-}=0 and the maximum eigenvalue is λ+=4σ2(R~)\lambda_{+}=4\sigma^{2}(\tilde{R}). For the two collections of securities we study, λ+c=1.45\lambda^{c}_{+}=1.45 and λ+e=1.75\lambda^{e}_{+}=1.75.

Refer to caption
(a)
Refer to caption
(b)
Figure 3: Plot of theoretical density and empirical density for (a) cryptocurrencies and (b) equities between 01-01-2019 to 01-12-2021. The equity collection clearly exhibits a greater number of non-random eigenvalues.

Figure 3 shows the theoretical and empirical probability density functions of each collections’ eigenvalues. Several noteable differences between cryptocurrency and equity behaviours are observed. First, the cryptocurrency eigenvalues, which are shown in Figure 3(a), exhibit one dominant eigenvalue, λ1c\lambda^{c}_{1}, which is significantly greater than the theoretical upper bound, λ+c\lambda_{+}^{c}. This confirms the results demonstrated by [24] in Figure 3 of their work. By contrast, the equity collection (shown in Figure 3(b)) consists of four eigenvalues (λ1e,λ2e,λ3e\lambda_{1}^{e},\lambda_{2}^{e},\lambda_{3}^{e} and λ4e\lambda_{4}^{e}) all of which are larger in magnitude than the theoretical largest eigenvalue, λ+e\lambda^{e}_{+}. Our finding for the equity market is consistent with the results shown in [71]. The authors highlight that the presence of several significant eigenvalues (beyond the standard Wishart regime) corresponds to the effect of grouping stocks within various market sectors. Further discussion on this phenomenon, and similar phenomena, can be found in their extensive review. There is also significant literature detailing more complex models surrounding Wishart matrices [77, 78, 79, 80]. The difference in the number of non-random eigenvalues between our collection is striking, and may be suggestive of a more profound difference in the impact of exogenous variables on both collections. Given that the first (dominant) eigenvalue represents the collective behaviour of the market [1], the existence of meaningful subsequent eigenvalues in equities may highlight the sophistication of the market (particularly in comparison to cryptocurrencies), where non-market related collective behaviours influence the correlation structure between underlying securities. By contrast, the cryptocurrency market’s exhibiting one non-random eigenvalue draws attention to the overwhelming importance of the ‘market’ factor among the collection.

Refer to caption
(a)
Refer to caption
(b)
Figure 4: Time-varying computation of leading eigenvalue magnitude and number of non-random eigenvalues for (a) cryptocurrencies and (b) equities. The equity collection consistently produces a larger number of non-random eigenvalues.

Having studied the distribution of correlation matrix eigenvalues in the stationary case, we now investigate how this structure evolves over time in both collections. Figure 4 displays the evolution in the magnitude of the dominant eigenvalue (λ1c(t1)\lambda_{1}^{c}(t_{1}) and λ1e(t2)\lambda_{1}^{e}(t_{2})) and the number of eigenvalues greater than the theoretical maximum eigenvalue λ+c(t1)\lambda^{c}_{+}(t_{1}) and λ+e(t2)\lambda^{e}_{+}(t_{2}), respectively. There are several significant findings. The number of non-random eigenvalues in the cryptocurrency collection, shown in Figure 4(a), ranges between 1 and 2. By contrast, the equity collection exhibits between 2 and 4 non-random eigenvalues over the entirety of our analysis window. This is shown in Figure 4(b). This confirms the persistence of the findings first highlighted in our static analysis, with the equity collection consistently exhibiting a larger number of non-random eigenvalues than the cryptocurrency collection. Both figures reveal a consistently inverse relationship between the number of non-random eigenvalues and the magnitude of the first eigenvalue. The intuition behind this relationship is interesting and open to a variety of interpretations. One possible argument is that when there is an increase in the magnitude of the dominant eigenvalue, the collective strength of the market increases. Given that this phenomenon is shown to coincide with a smaller number of non-random eigenvalues, one could assume that when the collective strength of the market is high - there is less opportunity for alternative factors to influence correlation behaviours of underlying securities. Given that spikes in the first eigenvalue tend to coincide with bear markets and the indiscriminate selling of assets, one can see that during market crises, exogenous factors other than the collective strength of the market (represented by eigenvalues λ2\lambda_{2},…,λm\lambda_{m} such that λm>λ+\lambda_{m}>\lambda_{+}) have less influence on correlation dynamics. This pattern appears to exist in both the cryptocurrency and equity markets.

4 Sector evolutionary dynamics

In this section we wish to explore the potential heterogeneity of each time series collection, and compare the evolution in the strength of market dynamics, represented by the correlation matrix’s dominant eigenvalue. We begin by partitioning our collection of cryptocurrencies into 8 sectors, and computing a daily rolling correlation matrix for each sector, again setting S=150S=150. Let our cryptocurrency correlation matrices at time t1t_{1} be Ωc,ψl(t1)\Omega^{c,\psi_{l}}(t_{1}), where ψl:l=1,,8\psi_{l}:l=1,...,8 refers to any one of L=8L=8 potential cryptocurrency sectors. Similarly, let Ωe,ψp(t2)\Omega^{e,\psi_{p}}(t_{2}) refer to one of ψp:p=1,,11\psi_{p}:p=1,...,11 be any one of P=11P=11 equity sector correlation matrices at time t2t_{2}. We are particularly interested in studying the evolution of collective behaviour within candidate sectors, and how much this evolution varies between different sectors. Let λ~1c,ψl(t1)\tilde{\lambda}^{c,\psi_{l}}_{1}(t_{1}) and λ~1e,ψp(t2)\tilde{\lambda}^{e,\psi_{p}}_{1}(t_{2}) be the explanatory variance exhibited by the first eigenvalue of each sector’s correlation matrix at any point in time. This normalization is computed as described in Section 3. We determine the affinity between sectors’ evolutionary collective similarity by computing an L1L^{1} norm between any two paths of the first eigenvalue’s explanatory variance. Our two distance matrices are computed as follows:

Refer to caption
(a)
Refer to caption
(b)
Figure 5: Time-varying explanatory variance exhibited by the first eigenvalue of each sector’s correlation matrix among (a) cryptocurrencies and (b) equities. The cryptocurrency market is partitioned into 8 sub-groups (sectors) including: centralized exchange, collectibles/NFTs, Decentralized finance, Platforms, Privacy, Smart contracts, Store of value and Wallet. Equities are classified based on their GICS classification: Communication services, Consumer discretionary, Consumer staples, Energy, Financials, Healthcare, Industrials, Information technology, Materials, Real estate and Utilities.
Dijc,λ1~=1T1St1=ST1|λ~1c,ψi(t1)λ~1c,ψj(t1)|,\displaystyle D_{ij}^{c,\tilde{\lambda_{1}}}=\frac{1}{T_{1}-S}\sum^{T_{1}}_{t_{1}=S}\bigg{|}\tilde{\lambda}^{c,\psi_{i}}_{1}(t_{1})-\tilde{\lambda}^{c,\psi_{j}}_{1}(t_{1})\bigg{|}, (11)
Dije,λ1~=1T2St2=ST2|λ~1e,ψi(t2)λ~1e,ψj(t2)|.\displaystyle D_{ij}^{e,\tilde{\lambda_{1}}}=\frac{1}{T_{2}-S}\sum^{T_{2}}_{t_{2}=S}\bigg{|}\tilde{\lambda}^{e,\psi_{i}}_{1}(t_{2})-\tilde{\lambda}^{e,\psi_{j}}_{1}(t_{2})\bigg{|}. (12)
Refer to caption
(a)
Refer to caption
(b)
Figure 6: Hierarchical clustering applied to our two distance matrices (a) Dc,λ1~D^{c,\tilde{\lambda_{1}}} and (b) De,λ1~D^{e,\tilde{\lambda_{1}}}. The dendrograms reveal some notable differences in the heterogeneity of evolutionary dynamics of cryptocurrency and equity sectors. The cryptocurrency collection exhibits a propensity for coupling behaviour between sectors, while the equity collection produces a predominant cluster of high affinity (and a small collection of outlier sectors).

Figure 5 shows the time-varying explanatory variance exhibited by respective cryptocurrency and equity sectors. There are some broad similarities in the evolution of respective sectors and the time-varying affinity among both collections. Both the cryptocurrency and equity collections, shown in Figures 5(a) and 5(b), display surprising consistency in the ranking of each sector’s first eigenvalue explanatory variance over time. In the cryptocurrency sector, smart contract-themed cryptocurrencies consistently display the highest collective market strength, while wallet-themed cryptocurrencies consistently display the lowest. Within our equity collection, utilities display the strongest collective market dynamics, while our sample of communications stocks exhibits the weakest collective market strength. Both collections produce a similar range in λ~1\tilde{\lambda}_{1} scores over time, suggesting that each asset class would present similarly difficult problems for portfolio managers seeking diversification among security returns. One notable difference between the collections however, is the range in λ~1\tilde{\lambda}_{1} values during the COVID-19 market crisis. The range among cryptocurrency sectors is approximately 0.3, while the range among respective equity sectors is approximately 0.2. This may suggest that for portfolio managers focused solely on individual asset classes, cryptocurrencies could provide more opportunity for diversifying across collective behaviours during market crises.

We further investigate the similarity in evolutionary market dynamics by applying hierarchical clustering to our two distance matrices, Dc,λ1~D^{c,\tilde{\lambda_{1}}} and De,λ1~D^{e,\tilde{\lambda_{1}}}. Figures 6(a) and 6(b) show the cryptocurrency and equity evolutionary collective behaviour dendrograms respectively. Both figures displays broad similarity in the collective affinity between respective sectors, with the scale of both dendrograms ranging from 0 to 0.32 and 0 to 0.35 respectively. Figure 6(a) reveals the propensity for coupling among cryptocurrency sectors. Three notable pairs are revealed between Store of value and privacy, platform and decentralized finance and collectibles/nfts and centralized exchange sectors. The smart contracts and wallet-themed cryptocurrencies exhibit more anomalous evolutionary collective dynamics with respect to the rest of the collection. By contrast the collection of equities, which are clustered together in Figure 6(a), displays a significantly different cluster structure. There is a predominant cluster of similarity among 8 sectors (real estate, materials, information technology, consumer staples, consumer discretionary, industrials, healthcare and communications), with a small collection of oultier sectors which includes utilities, energy and financials. The predominant similarity between constituent sectors of the primary cluster are consistent with the λ1e~\tilde{\lambda^{e}_{1}} paths shown in Figure 5(b), where collective behaviours are shown to be of closer proximity than that of the cryptocurrency market.

5 Structural breaks and spectral regime switching

5.1 Structural break dynamics

In this section we investigate the structural break dynamics and regime switching behaviour of cryptocurrencies and equities. We extend a recently introduced framework [81], which marries a Bayesian change point detection algorithm [62, 60] with a new distance metric which quantifies distance between sets of probability distributions. The Bayesian changepoint detection framework utilizes a reversible jump Markov chain Monte Carlo (RJMCMC) algorithm to produce a distribution over possible models, where each model consists of a varying number of change points, mm. In all RJMCMC experiments we run 10,000 iterations and discard the first 5,000 of these iterations (burn-in). We select the maximum a posteriori model for each time series, and accordingly determine an appropriate number of change points for each cryptocurrency and equity security respectively. In subsequent experiments we elect one representative cryptocurrency/equity from each sector and apply our respective methodology. This is primarily for computational savings, given the high cost in running the Reversible Jump Markov Chain Monte Carlo algorithm.

Refer to caption
(a)
Refer to caption
(b)
Refer to caption
(c)
Refer to caption
(d)
Figure 7: A sample of cryptocurrency sector log returns with annotated change points (a) Centralized exchange (b) Decentralized finance, (c) Wallet, (d) Platform. The transparency of the black vertical lines indicate the probability of a changepoint at any candidate location (darker colours are consistent with greater probability mass). Although the four sectors above all display structural breaks corresponding to the beginning and end of the COVID-19 market crash, there is significantly greater variability in the number and location of structural breaks than that exhibited by the equity collection.
Refer to caption
(a)
Refer to caption
(b)
Refer to caption
(c)
Refer to caption
(d)
Figure 8: A sample of equity sector log returns with annotated change points (a) Energy (b) Financials, (c) Industrials, (d) Materials. The transparency of the black vertical lines indicate the probability of a changepoint at any candidate location (darker colours are consistent with greater probability mass). All equity sectors exhibit two clear structural breaks indicating the beginning and end of the COVID-19 market crash.

The location of each changepoint is associated with a probability density function. Let S~,T~\tilde{S},\tilde{T} refer to any two sets with uncertainty within our respective collections. For any candidate time series pair, let S~={f1,,fa}\tilde{S}=\{f_{1},...,f_{a}\} and T~={g1,,gb}\tilde{T}=\{g_{1},...,g_{b}\}, be our two sets with uncertainty which contain aa and bb probability distributions respectively. Within each set, the probability distributions 1,,a1,...,a and 1,,b1,...,b associate a collection of time series indices with a changepoint probability. In the proceeding experiments, we set the polynomial order o=1o=1, and compute two MJ-Wasserstein distance matrices as follows:

Dc,MJW(S~ic,T~jc)\displaystyle D^{c,MJW}(\tilde{S}^{c}_{i},\tilde{T}^{c}_{j}) =(gT~jcdW(g,S~ic)o2|T~jc|+fS~icdW(f,T~jc)o2|S~ic|)1o\displaystyle=\Bigg{(}\frac{\sum_{g\in\tilde{T}^{c}_{j}}d_{W}(g,\tilde{S}^{c}_{i})^{o}}{2|\tilde{T}^{c}_{j}|}+\frac{\sum_{{f}\in\tilde{S}^{c}_{i}}d_{W}(f,\tilde{T}^{c}_{j})^{o}}{2|\tilde{S}^{c}_{i}|}\Bigg{)}^{\frac{1}{o}} (13)
De,MJW(S~ie,T~je)\displaystyle D^{e,MJW}(\tilde{S}^{e}_{i},\tilde{T}^{e}_{j}) =(gT~jedW(g,S~ie)o2|T~je|+fS~iedW(f,T~je)o2|S~ie|)1o.\displaystyle=\Bigg{(}\frac{\sum_{g\in\tilde{T}^{e}_{j}}d_{W}(g,\tilde{S}^{e}_{i})^{o}}{2|\tilde{T}^{e}_{j}|}+\frac{\sum_{{f}\in\tilde{S}^{e}_{i}}d_{W}(f,\tilde{T}^{e}_{j})^{o}}{2|\tilde{S}^{e}_{i}|}\Bigg{)}^{\frac{1}{o}}. (14)

In the two equations above (outlining the distance among our cryptocurrency and equity collections), for each distribution fif_{i} in S~\tilde{S}, the nearest distribution gjg_{j} in T~\tilde{T} is identified according to the Wasserstein distance metric, and the minimal distance between any two sets is computed. This allows us to determine the similarity between any two time series’ collection of changepoints. Importantly, this metric allows us to capture the uncertainty associated with the location of each change point within any time series’ changepoint profile.

Estimated number of regimes by asset class sector
Asset class Sector # Regimes
Cryptocurrency Centralized exchange 6
Cryptocurrency Collectibles/NFTs 6
Cryptocurrency Decentralized finance 2
Cryptocurrency Platform 3
Cryptocurrency Privacy 5
Cryptocurrency Smart contracts 5
Cryptocurrency Store of value 2
Cryptocurrency Wallet 5
Equity Communication services 4
Equity Consumer discretionary 3
Equity Consumer staples 4
Equity Energy 3
Equity Financials 4
Equity Healthcare 4
Equity Industrials 4
Equity Information Technology 4
Equity Materials 3
Equity Real estate 4
Equity Utilities 4
Table 1: Estimated number of sector regimes based on Reversible Jump MCMC algorithm’s determination.

The methodology we use highlights marked differences between cryptocurrency and equity sectors’ structural breaks. Figure 7 shows the log returns and estimated structural breaks of four cryptocurrency sectors; centralized exchange (Figure 7(a)), decentralized finance (Figure 7(b)), wallet (Figure 7(c)) and platform (Figure 7(d)). There is clearly great variability in the number, location and associated uncertainty of each sector’s structural breaks. The Centralized exchange sector exhibits just one change point, while decentralized finance and wallet both exhibit four and platform exhibits five. Turning to equities, Figure 8 shows four equity sectors’ empirical log returns with annotated change points determined by the change point algorithm. It is clear from the Figure that energy (Figure 8(a)), financials (Figure 8(b)), industrials (Figure 8(c)) and materials (Figure 8(d)) all display highly similar change point dynamics around the COVID-19 market crisis. The energy sector and materials sector both exhibit two structural breaks, while financials and industrials exhibit three change points. Although both the financials and materials sectors are determined to exhibit an additional (similarly located) structural break corresponding to a shift in these sectors’ log returns behaviours, in general, equity sectors appear to exhibit greater affinity in their structural break profile than that of the cryptocurrency sectors. However, the financials and industrials sectors both exhibit one additional changepoint located in early 2021 - corresponding to a shift in the behaviours of these sectors’ log returns.

We further investigate differences in the propensity for various sectors’ structural breaks in Table 1, which indicates the determined number of locally stationary segments. Among cryptocurrency sectors, the average number of regimes is 4.25 with a modal value of 5. The average number of equity sector regimes by contrast is 3.7, with a modal value of 4. The stark difference in structural break similarity is evidenced in Figure 9, which highlights significantly greater similarity in the structural break profile of equities. Such analysis considers both the number and location of such change points. An indication of this, is the comparative scale of both dendrograms. The cryptocurrency structural break dendrogram, shown in Figure 9(a), has a range of 0.75, while the equity structural break dendrogram, shown in Figure 9(b), has a range of just 0.13. One may conclude from this analysis, that cryptocurrency sectors are more prone to regime switching behaviours, and display less similarity in their structural break profile than equity sectors.

Refer to caption
(a)
Refer to caption
(b)
Figure 9: Hierarchical clustering applied to cryptocurrency and equity structural break distance matrices, Dc,MJWD^{c,MJW} and De,MJWD^{e,MJW}. The two dendrograms highlight significantly greater similarity between equity sector structural breaks than that exhibited by cryptocurrency sectors.

5.2 Spectral regime switching

Next, we seek to further investigate the collective similarity in evolutionary power spectra among cryptocurrency and equity sectors. In the previous section, we focused solely on the measurement between sets of structural breaks, ignoring all the data within any candidate segment. In this section, we study the evolutionary power spectral density of each sector and determine respective distances. This measure may provide a better indication of the similarity in the behaviours between any two sectors, where regime switching (or abrupt changes in behaviour) may exist.

Unlike many regime switching methods which rely on strong parametric assumptions, our nonparametric Bayesian framework determines the number of regimes associated with each time series. Structural breaks are determined based on changes in the time series’ evolutionary power spectral density. Accordingly, we can study the similarity in the number of regimes between various sectors, and the similarity in their evolutionary power spectral density. Let the time-varying power spectrum for each cryptocurrency sector be fψlc(ν,t1)f^{c}_{\psi_{l}}(\nu,t_{1}) where the pair (ν,t1)(\nu,t_{1}) refers to a vector of respective frequency components at points in time t1t_{1}, and ψl:l=1,,8\psi_{l}:l=1,...,8 indexes one of L=8L=8 possible cryptocurrency sectors. Similarly, let fψpe(ν,t2)f^{e}_{\psi_{p}}(\nu,t_{2}) refer to each equity sector’s evolutionary power spectrum, where again (ν,t2)(\nu,t_{2}) refers to a vector of frequency components at various point in time t2t_{2}, and ψp:p=1,,11\psi_{p}:p=1,...,11 indexes one of P=11P=11 possible equity sectors.

Dijc,(ν,t1)=1T1×Kt1=1T1k=1K|fψic(νk,t1)fψjc(νk,t1)|,\displaystyle D_{ij}^{c,({\nu,t_{1}})}=\frac{1}{T_{1}\times K}\sum^{T_{1}}_{t_{1}=1}\sum^{K}_{k=1}\bigg{|}f_{\psi_{i}}^{c}(\nu_{k},t_{1})-f_{\psi_{j}}^{c}(\nu_{k},t_{1})\bigg{|}, (15)
Dije,(ν,t2)=1T2×Kt2=1T2k=1K|fψie(νk,t2)fψje(νk,t2)|.\displaystyle D_{ij}^{e,({\nu,t_{2}})}=\frac{1}{T_{2}\times K}\sum^{T_{2}}_{t_{2}=1}\sum^{K}_{k=1}\bigg{|}f_{\psi_{i}}^{e}(\nu_{k},t_{2})-f_{\psi_{j}}^{e}(\nu_{k},t_{2})\bigg{|}. (16)
Refer to caption
(a)
Refer to caption
(b)
Figure 10: Hierarchical clustering applied to cryptocurrency and equity time-varying spectrum distance matrices, Dc,(ν,t1)D^{c,(\nu,t_{1})} and De,(ν,t2)D^{e,(\nu,t_{2})}.
Refer to caption
(a)
Refer to caption
(b)
Refer to caption
(c)
Refer to caption
(d)
Figure 11: A sample of cryptocurrency evolutionary power spectra (a) Privacy (b) Wallet, (c) Smart Contracts, (d) Platform.
Refer to caption
(a)
Refer to caption
(b)
Refer to caption
(c)
Refer to caption
(d)
Figure 12: A sample of equity evolutionary power spectra (a) Energy (b) Financials, (c) Industrials, (d) Materials.

The time-varying power spectral densities for the cryptocurrency and equity collections are shown in Figures 11 and 12. Their surfaces reveal some interesting differences in the evolutionary behaviour of cryptocurrency and equity sectors. Among the four cryptocurrency sectors shown, all exhibit highly varying power spectra over time. In general, there is an overwhelming tendency for higher frequency components to dominate, indicating the abrupt changes in behaviours exhibited by cryptocurrencies. This is particularly pronounced in early 2020, with all sectors showing a significant translation upward in the level of the power spectrum, and a clear shift toward higher frequency components exhibiting most power. This is indicative of the volatility experienced in the market during the COVID-19 market crash. These trends are relatively consistent in Figures 11(a), 11(b), 11(c) and 11(d) representing privacy, wallet, smart contract and platform focused cryptocurrencies.

Next we turn to our collection of equities. Equity sectors are shown in Figure 12, with energy, financials, industrials and materials’ time-varying power spectra shown in Figures 12(a), 12(b), 12(c) and 12(d) respectively. The power spectra exhibited by equities is significantly smoother, and there are clearly less significant re-allocations of power among frequency components when transitioning to adjacent segments. Like the cryptocurrency surfaces, the equity spectra exhibit a pronounced translation during the COVID-19 market crash, with significant power moving to higher frequency components among all equity sectors.

The clustered distance matrix between the time-varying spectra of respective sectors for the cryptocurrency and equity collections are shown in Figure 10. Turning first to cryptocurrencies, shown in Figure 10(a), two characteristic clusters of sector behaviours emerge. Cluster one consists of smart contracts, privacy centralized exchange, platform and store of value themed cryptocurrency sectors. The second cluster consists of wallet, decentralized finance and collectibles/NFT-themed cryptocurrency sectors. The overall cluster structure in evolutionary spectra among cryptocurrency sectors is shown to be relatively diffuse, without the emergence of a predominant collection of sectors with characteristic spectral properties. By contrast, the similarity in time-varying power spectra of equity sectors shown in Figure 10(b), highlights the emergence of a dominant cluster, consisting of most equity sectors. The primary cluster consists of: Utilities, Healthcare, Consumer staples, Information technology, Financials, Consumer discretionary, Communications, Materials and Real estate. All sectors exhibit power spectra with similar levels of amplitude, and clear peaks in the power spectra at high frequency components. In the outlier cluster, which consists of the Industrials and Energy sectors, spectral peaks are exhibited at slightly lower frequency components - differentiating the surfaces from the rest of the collection. Viewing both dendrograms in conjunction with the time-varying power spectra figures, more consistent similarity is seen among equity sector behaviours. The time-varying power spectra of cryptocurrencies and the increased frequency of structural breaks reflects the higher levels of volatility exhibited by the underlying cryptocurrency sectors.

6 Cryptocurrency time-varying optimization

In this section we introduce two new simulation frameworks to test the relative importance of security selection and sector allocation in the context of cryptocurrency portfolio management. There is a litany of research in the quantitative finance and optimization fields focused on untangling this question for more traditional asset classes. However, given the nascence of the cryptocurrency market, it is worth investigating the most purposeful decisions for candidate portfolio managers.

Specifically, we wish to determine whether there is more portfolio alpha available at the stock selection, or sector allocation level. To test our hypothesis, we implement two iterative algorithms which aim to maximize portfolio Sharpe ratio subject to various portfolio constraints.

Algorithm 1 Optimal security selection
1:Choose a smoothing rate SS over time.
2:Choose a number of ‘best’ securities to sample from each sector, BB.
3:Let the cryptocurrency sectors be indexed ψl:l=1,,L\psi_{l}:l=1,...,L.
4:Let the number of securities per cryptocurrency sector be J=7J=7.
5:for t1=S to T1t_{1}=S\,\text{ to }\,T_{1} do
6:     for ψl:l=1 to L\psi_{l}:l=1\,\text{ to }\,L do
7:         for j=1 to Jj=1\,\text{ to }\,J do
8:              Compute total log returns Zjψl(t1)=d=t1St1Rjψl(d)Z^{\psi_{l}}_{j}(t_{1})=\sum^{t_{1}}_{d=t_{1}-S}R^{\psi_{l}}_{j}(d).
9:              Compute volatility σj2,ψl(t1)\sigma_{j}^{2,\psi_{l}}(t_{1}).
10:              Compute risk-adjusted return RARjψl(t1)=Zjψl(t1)/σj2,ψl(t1)RAR^{\psi_{l}}_{j}(t_{1})={Z^{\psi_{l}}_{j}(t_{1})}/{\sigma_{j}^{2,\psi_{l}}(t_{1})}          
11:         Select BψlB_{\psi_{l}} best securities within sector ψl\psi_{l} ranked by RARjψl(t1)RAR_{j}^{\psi_{l}}(t_{1}).
12:         Compute sector return Zψl(t1)=1Bψlj=1BψlZjψl(t1)Z^{\psi_{l}}(t_{1})=\frac{1}{B_{\psi_{l}}}\sum^{B_{\psi_{l}}}_{j=1}Z^{\psi_{l}}_{j}(t_{1})
13:         Append sector return Zψl(t1)Z^{\psi_{l}}(t_{1}) .      
14:     Compute portfolio return Z(t1)=1Ll=1LZψl(t1)Z(t_{1})=\frac{1}{L}\sum^{L}_{l=1}Z^{\psi_{l}}(t_{1}).
15:     Append total portfolio return Z(t1)Z(t_{1}).
16:Compute optimal cryptocurrency portfolio return over all time Rpo,c=t1=ST1Z(t1)R^{o,c}_{p}=\sum^{T_{1}}_{t_{1}=S}Z(t_{1}).

In Algorithm 1, we implement a time-varying simulation where we select a collection of cryptocurrency securities from each sector that exhibit a maximal Sharpe ratio over some trailing window of size SS. In doing so, we allocate uniform weight across our 8 candidate cryptocurrency sectors, and choose a sample of the BB best securities within each sector, measured by their trailing volatility-adjusted log returns. To simulate the dynamics of a portfolio management operation, there is no optimization routine performed at the sector-specific level, and all sectors are assumed to have the same amount of allocated capital at all points in time t1t_{1}. This is consistent with the concept of daily portfolio rebalancing.

Algorithm 2 Optimal sector allocation
1:Choose a smoothing rate SS over time.
2:Choose a ‘best’ number of sectors, BB, to allocate capital to.
3:Let the cryptocurrency sectors be indexed ψl:l=1,,L\psi_{l}:l=1,...,L.
4:Let the number of securities per cryptocurrency sector be J=7J=7.
5:for t1=S to T1t_{1}=S\,\text{ to }\,T_{1} do
6:     for ψl:l=1 to L\psi_{l}:l=1\,\text{ to }\,L do
7:         for j=1 to Jj=1\,\text{ to }\,J do
8:              Compute total log returns for each security Zjψl(t1)=d=t1St1Rjψl(d)Z_{j}^{\psi_{l}}(t_{1})=\sum^{t_{1}}_{d=t_{1}-S}R^{\psi_{l}}_{j}(d)          
9:         Assume uniform weight for each security wjψl=1/Jj{1,,J}w_{j}^{\psi_{l}}=1/J\quad\forall j\in\{1,...,J\}
10:         Compute sector returns Zψl(t1)=j=1JZjψl(t1)wjψlZ^{\psi_{l}}(t_{1})=\sum^{J}_{j=1}Z^{\psi_{l}}_{j}(t_{1})w_{j}^{\psi_{l}}
11:         Compute sector volatility σ2,ψl(t1)=𝐰ψlTΣ(t1)𝐰ψl\sigma^{2,\psi_{l}}(t_{1})=\mathbf{w}^{\psi_{l}T}\Sigma(t_{1})\mathbf{w}^{\psi_{l}}
12:         Compute sector risk-adjusted return RARψl(t1)=Zψl(t1)/σ2,ψl(t1)RAR^{\psi_{l}}(t_{1})=Z^{\psi_{l}}(t_{1})/\sigma^{2,\psi_{l}}(t_{1})      
13:     Select BB best sectors ranked by RARψl(t1)RAR^{\psi_{l}}(t_{1})
14:     Compute total portfolio return Z(t1)=1Bq=1BZqψl(t1)Z(t_{1})=\frac{1}{B}\sum^{B}_{q=1}Z^{\psi_{l}}_{q}(t_{1})
15:     Append total portfolio return Z(t1)Z(t_{1})
16:Compute portfolio return over time Rpo,s=t1=ST1Z(t1)R^{o,s}_{p}=\sum^{T_{1}}_{t_{1}=S}Z(t_{1}).

In Algorithm 2, we seek to determine the BB best performing cryptocurrency sectors over some trailing window SS, and uniformly allocate capital across these sectors. To judge the rolling performance of each sector, we compute the sector risk-adjusted returns at each point in time t1t_{1} (assuming uniform weight among constituent securities), and again select the BB best sectors to allocate capital toward. In this simulation, we have a constant assumption of uniform weight allocated across cryptocurrency securities within any candidate sector, and uniform weight allocated across the BB best sectors. This is done to explicitly test for alpha generation at the level of sector allocation.

Cryptocurrency portfolio return simulations
SS BB Rpo,cR^{o,c}_{p} Rpo,sR^{o,s}_{p}
120 2 215% 177%
120 3 210% 201%
120 4 210% 219%
120 5 214% 235%
150 2 155% 130%
150 3 167% 162%
150 4 177% 207%
150 5 177% 208%
180 2 166% 168%
180 3 149% 179%
180 4 179% 175%
180 5 187% 180%
Table 2: Simulated cryptocurrency portfolio return simulations. In various iterations of the simulation, BB and SS are varied. Our results are inconclusive as to whether there is more value in cryptocurrency sector allocation or security selection.

We run Algorithm 1 and Algorithm 2 for a range of values SS and BB, and compute the portfolio returns for optimal cryptocurrency selection Rpo,cR_{p}^{o,c} and optimal security selection Rpo,sR_{p}^{o,s} for all possible parameter combinations. The results are shown in Table 2, and are relatively stable and consistent among all possible parameter values. The table highlights several interesting findings. The most notable observation however, is that there is no clear winner when comparing the realised performance of the security selection (Algorithm 1) and sector allocation (Algorithm 2) algorithms. The second notable observation is the significantly improved performance when S=120S=120. In both Algorithm 1 and Algorithm 2, the highest average portfolio returns are exhibited with a smaller rolling window. This reflects the dynamic nature of cryptocurrency behaviours, and demonstrates the need for continued monitoring of market dynamics such as correlation. The final interesting observation, is that there is typically an improvement in portfolio performance (subject to a constant level SS) as BB increases. This phenomenon occurs in both portfolios, but is particularly pronounced in the sector allocation simulation. There are two takeaways from this observation. First, there is clearly a meaningful benefit when cryptocurrency portfolios are well diversified. The second point is more subtle. Given that the improvement (subject to a constant SS) when BB increases is more pronounced in Algorithm 2, it highlights a potential risk in cryptocurrency characteristic behaviours. There appears to be a tendency for certain cryptocurrency sectors to trend, and holding more exposure to a variety of cryptocurrency sectors could be of great benefit to any investor in the cryptocurrency markets. This finding may question the notion that holding a major (dominant) cryptocurrency such as Bitcoin or Ethereum, provides satisfactory exposure to the entire cryptocurrency market.

7 Conclusion

In Section 3 we use random matrix theory in a time-varying capacity to explore the evolutionary properties of cryptocurrency and equity correlation matrices. Our methodology produces several interesting findings. First, the equity collection consistently exhibits a larger number of non-random eigenvalues over our analysis window. This is consistent with the findings shown in [24]. By contrast with the equity market, we show that the cryptocurrency market consistently displays a smaller number of non-random eigenvalues. When interpreting this distinction between the two asset classes, it is worth clarifying the possible interpretation of such findings. Large outliers in the eigenspectrum between the bulk of the distribution of eigenvalues and the largest eigenvalue (representing the strength of collective behaviour, or the market) measures the degree of collectivity in groups of sectors. The origins of such collectivity in the sectors could stem from exogenous or endogenous factors. In the case of the former, this may represent new laws or regulations which only impact one industry for example, and in the case of the latter, investor market expectations may shift with respect to only one sector for some sentiment-driven investor reasoning. Our analysis certainly implies that there is a greater degree of complexity in the eigenspectrum exhibited by equities, and this is consistent over time. The second key finding is the notable inverse relationship observed between the magnitude of the first correlation matrix eigenvalue and the number of non-random eigenvalues based on the random matrix theory. This relationship implies that when the collective strength of market behaviours increases (associated with the magnitude of the first eigenvalue), the effect of the market itself becomes more pronounced on underlying financial market dynamics.

In Section 4, we demonstrate marked differences in the similarity of cryptocurrency and equity sectors’ evolutionary dynamics. In particular, cryptocurrencies exhibit quite an obvious ‘coupling’ behaviour, where sectors with similar underlying dynamics are revealed. By contrast, hierarchical clustering on our equity collection reveals one predominant cluster of similar behaviours, and a small group of outlier sectors (Utilities, Energy and Financials). As the market continues to develop, such techniques could be used to identify sub-groups within the cryptocurrency space that behave most similarly. Given the nascence and rapid changes in the market, this is worth monitoring closely.

In Section 5 we study the collective similarity in structural break propagation and time-varying power spectra in both collections. Our frameworks reveal significantly greater similarity in the structural break dynamics of equity sectors than that of cryptocurrency sectors. This is consistent with prior research which demonstrates the existence of erratic behaviours in cryptocurrency returns and volatility [14]. Our experiments also reveal that cryptocurrency sectors tend to consist of a greater number locally-stationary segments than equities, further highlighting the tendency for changes in time series behaviour (structural breaks) in cryptocurrency sectors. To better understand the subtle difference in nonstationarity and long range dependence, we study the similarity between sectors’ time-varying power spectra. Our experiments reveal that equities produce one predominant cluster of evolutionary spectra with just two outliers (Industrials and Energy), while cryptocurrencies are partitioned into two, (approximately) equally-sized clusters.

Finally in Section 6 we focus exclusively on the burgeoning cryptocurrency market, and investigate the relative importance of security selection and sector allocation in cryptocurrency portfolio management. We implement two algorithms, with a range of values for our smoothing window parameter SS, and best number of samples BB (where we sample underlying securities in Algorithm 1 and cryptocurrency sectors in Algorithm 2). Our experiments reveal that both decisions are crucial in cryptocurrency portfolio management, and neither piece is of greater importance than the other. Optimal portfolio returns are realised with the smallest possible smoothing window (S=120S=120), which demonstrates the dynamic nature of cryptocurrency behaviours and the need for adaptability in portfolio management decision-making. Finally, we demonstrate (in both algorithms), the benefit of diversifying across securities and sectors. Our experiments reveal that subject to a candidate level of SS, as BB increases portfolio returns also tend to increase.

We hope that this paper may serve as a current accompaniment to academics and practitioners seeking to understand the evolution of individual and collective cryptocurrency behaviours, and the implications for managing a portfolio of cryptocurrency investments. However, there are several avenues to further develop upon the frameworks introduced in this paper. First, one could extend all the analysis in this paper to other asset classes such as bonds, currencies and commodities. The study of time-varying sector explanatory variance could be accompanied by a study of sector equity returns. One could test the ‘resilience’ or ‘elasticity’ of sector returns to spikes in collective movements (typically associated with market crises). The study of sector structural breaks and evolutionary power spectra could be extended. In this work, we choose one representative security as a proxy for each sector. A focused paper on this particular issue could study intra-sector similarity in addition to inter-sector similarity. This may reveal specific sectors that exhibit more self-similarity in their structural break profile. One could also test for the consistency in collective similarity between evolutionary spectra and structural breaks, among various asset classes. Such analysis could identify anomalous sectors with regards to their evolutionary power spectra, or structural break profile. Finally, one could further develop the trading strategy simulations by incorporating other portfolio management decisions, or treating the problem as a more explicit, higher dimensional matter - where one could learn model parameters such as SS and BB subject to some candidate loss function.

8 Acknowledgments

The author would like to thank Kevin Chin from Arowana & Co. who helped provide some of the key motivating questions for the author.

Appendix A Mathematical objects

Table 3: Mathematical objects
Mathematical object Description
t1=1,,T1t_{1}=1,...,T_{1} Cryptocurrency time index, T1=1065T_{1}=1065
t2=1,,T2t_{2}=1,...,T_{2} Equity time index, T2=733T_{2}=733
picp^{c}_{i} Multivariate time series of cryptocurrency prices.
pjep^{e}_{j} Multivariate time series of equity prices.
RicR^{c}_{i} Multivariate time series of cryptocurrency log returns.
RjeR^{e}_{j} Multivariate time series of equity log returns.
NN Number of cryptocurrency time series.
KK Number of equity time series.
Ωc\Omega^{c} Cryptocurrency correlation matrix.
Ωe\Omega^{e} Equity correlation matrix.
λ1c,,λNc\lambda_{1}^{c},...,\lambda^{c}_{N} Cryptocurrency eigenvalues ordered by magnitude.
λ1e,,λKe\lambda_{1}^{e},...,\lambda^{e}_{K} Equity eigenvalues ordered by magnitude.
λ~1c,,λ~Nc\tilde{\lambda}_{1}^{c},...,\tilde{\lambda}^{c}_{N} Normalized cryptocurrency eigenvalues ordered by magnitude.
λ~1e,,λ~Ke\tilde{\lambda}_{1}^{e},...,\tilde{\lambda}^{e}_{K} Normalized equity eigenvalues ordered by magnitude.
p(λc)p(\lambda^{c}) Probability density function of cryptocurrency eigenvalues.
p(λe)p(\lambda^{e}) Probability density function of equity eigenvalues.
ll Cryptocurrency sector time series index, L=8L=8.
pp Equity sector time series index, P=11P=11.
Ωc,ψl\Omega^{c,\psi_{l}} Cryptocurrency correlation matrix for sector ll.
Ωe,ψp\Omega^{e,\psi_{p}} Equity correlation matrix for sector pp.
λ1~c,ψl\tilde{\lambda_{1}}^{c,\psi_{l}} Cryptocurrency dominant eigenvalue explanatory variance for sector ll.
λ1~e,ψp\tilde{\lambda_{1}}^{e,\psi_{p}} Equity dominant eigenvalue explanatory variance for sector pp.
Dc,λ1~D^{c,\tilde{\lambda_{1}}} Cryptocurrency dominant eigenvalue distance matrix (collective behaviour).
De,λ1~D^{e,\tilde{\lambda_{1}}} Equity dominant eigenvalue distance matrix (collective behaviour).
Dc,MJWD^{c,MJW} Cryptocurrency MJ-Wasserstein distance matrix (structural breaks).
De,MJWD^{e,MJW} Equity MJ-Wasserstein distance matrix (structural breaks).
fψlc(ν,t1)f_{\psi_{l}}^{c}(\nu,t_{1}) Cryptocurrency time-varying power spectrum for sector ll.
fψpe(ν,t2)f_{\psi_{p}}^{e}(\nu,t_{2}) Equity time-varying power spectrum for sector pp.
Dc,(ν,t1)D^{c,(\nu,t_{1})} Cryptocurrency time-varying power spectrum distance matrix.
De,(ν,t2)D^{e,(\nu,t_{2})} Equity time-varying power spectrum distance matrix.

Appendix B Securities

Table 4: Cryptocurrency securities.
Security Ticker Sector
Binance Coin BNB Centralized exchange
Bankera BNK Centralized exchange
Crypto.com CRO Centralized exchange
Huobi Token HT Centralized exchange
KuCoin Token KCS Centralized exchange
LATOKEN LA Centralized exchange
ZB Token ZB Centralized exchange
Enjin Coin ENJ Collectibles/NFTs
Decentraland MANA Collectibles/NFTs
Phantasma SOUL Collectibles/NFTs
Syscoin SYS Collectibles/NFTs
THETA THETA Collectibles/NFTs
WAXP WAXP Collectibles/NFTs
Tezos XTZ Collectibles/NFTs
Basic Attention Token BAT Decentralized finance
Bancor BNT Decentralized finance
Fantom FTM Decentralized finance
Chainlink LINK Decentralized finance
Loopring LRC Decentralized finance
Maker MKR Decentralized finance
Tezos XTZ Decentralized finance
Cardano ADA Platform
Ethereum Classic ETC Platform
Fantom FTM Platform
Chinlink LINK Platform
Decentraland MANA Platform
Neo NEO Platform
Tezos XTZ Platform
Decred DCR Privacy
Flux FLUX Privacy
iExec RLC RLC Privacy
Monero XMR Privacy
Verge XVG Privacy
Zcash ZEC Privacy
Horizen ZEN Privacy
Cardano ADA Smart Contracts
Binance Coin BNB Smart Contracts
Ethereum Classic ETC Smart Contracts
Ethereum ETH Smart Contracts
Chainlink LINK Smart Contracts
VeChain VET Smart Contracts
Stellar XLM Smart Contracts
Bitcoin Cash BCH Store of Value
Bitcoin SV BSV Store of Value
Bitcoin BTC Store of Value
Decred DCR Store of Value
Maker MKR Store of Value
Bitcoin Plus XBC Store of Value
Nano XNO Store of Value
Bread BRD Wallet
Divi DIVI Wallet
Electroneum ETN Wallet
STASIS EURO EURS Wallet
Metronome MET Wallet
Pillar PLR Wallet
Voyager Token VGX Wallet
Table 5: Equity securities.
Security Ticker Sector
The Walt Disney Co. DIS Communication Services
Alphabet Inc. GOOG Communication Services
The Interpublic Group IPG Communication Services
Lumen Technologies LUMN Communication Services
Netflix Inc. NFLX Communication Services
News Corp. NWSA Communication Services
Verizon Communications VZ Communication Services
Amazon.com AMZN Consumer Discretionary
AutoZone Inc. AZO Consumer Discretionary
Best Buy Co. BBY Consumer Discretionary
General Motors GM Consumer Discretionary
Home Depot HD Consumer Discretionary
NIKE Inc. NKE Consumer Discretionary
Starbucks Corp. SBUX Consumer Discretionary
Kellogg Co. K Consumer Staples
The Kraft Heinz Co. KHC Consumer Staples
The Coca-Cola Co. KO Consumer Staples
PepsiCo. Inc PEP Consumer Staples
Procter & Gamble Co. PG Consumer Staples
Tyson Foods Inc. TSN Consumer Staples
Walmart Inc. WMT Consumer Staples
Chevron cO. CVX Energy
Haliburton Co. HAL Energy
Hess Corp. HES Energy
Kinder Morgan Inc. KMI Energy
Schlumberger NV SLB Energy
Valero Energy Corp. VLO Energy
Exxon Mobil Corp. XOM Energy
American Express Co. AXP Financials
Bank of America Corp. BAC Financials
Citigroup C Financials
JPMorgan Chase & Co. JPM Financials
MSCI Inc. MSCI Financials
Wells Fargo & Co. WFC Financials
Willis Towers Watson PLC WLTW Financials
Amgen Inc. AMGN Healthcare
Baxter International Inc. BAX Healthcare
Danaher Corp. DHR Healthcare
Johnson & Johnson JNJ Healthcare
McKesson Corp. MCK Healthcare
Merck & Co. Inc. MRK Healthcare
Pfizer Inc. PFE Healthcare
The Boeing Co. BA Industrials
Caterpillar Inc. CAT Industrials
Deere & Co. DE Industrials
Equifax Inc. EFX Industrials
FedEx Corp. FDX Industrials
General Electric Co. GE Industrials
3M Co. MMM Industrials
Autodesk Inc. ADSK Information Technology
Broadcom Inc. AVGO Information Technology
Cisco Systems Inc. CSCO Information Technology
Citrix Systems CTXS Information Technology
Hewlett Packard Enterprise Co. HPE Information Technology
Microsoft Corp. MSFT Information Technology
NortonLifeLock Inc. NLOK Information Technology
Ball Corp BLL Materials
DuPont de Nemours Inc. DD Materials
FMC Corp. FMC Materials
International Paper Co. IP Materials
LyondellBasell Industries NV LYB Materials
Vulcan Materials Co. VMC Materials
Westrock Co. WRK Materials
American Tower Corp. AMT Real estate
Boston Properties Inc. BXP Real estate
CBRE Group Inc. CBRE Real estate
Essex Property Trust Inc. ESS Real estate
Extra Space Storage Inc. EXR Real estate
Public Storage PSA Real estate
Welltower Inc. WELL Real estate
American Electric Power Co. Inc. AEP Utilities
American Water Works Co. AWK Utilities
Duke Energy Corp. DUK Utilities
Exelon Corp. EXC Utilities
Alliant Energy Corp. LNT Utilities
Pinnacle West Capital Corp. PNW Utilities
WEC Energy Group Inc. WEC Utilities

Appendix C Bayesian change point detection algorithm

In this section, we describe the Bayesian change point detection algorithm that was used in Section 5 of this paper. We use a reversible jump Markov chain Monte Carlo (RJMCMC) algorithm [61] that continues to partition a (presumably nonstationary) time series into mm segments, where mm also changes with time. The sampling scheme produces a distribution over possible models mm and also produces a set of mm change points.

The method identifies structural breaks/change points based on changes in the evolutionary power spectral density [60]. Many change point detection methodologies that exist in the time domain fail to identify suitable change points when time series exhibit dependence. This is commonly the case with autoregressive processes, which often display similar patterns to that exhibited by financial time series. Previous research [61] has highlighted that the methodology used in this paper can detect changes in a sequence of autoregressive processes, highlighting the algorithm’s ability to partition data that exhibits dependence.

We follow Rosen et al. and James and Menzies [60, 62] in our implementation of the RJMCMC sampling scheme. In particular, James and Menzies introduce a framework where an optimal number of basis functions is used when estimating a candidate power spectral density. We denote a varying partition of the time series by 𝝃m=(ξ0,m,,ξm,m)\bm{\xi}_{m}=(\xi_{0,m},...,\xi_{m,m}); these are our mm change points (excluding ξm,m\xi_{m,m}, which by convention is always going to be the final time point). The algorithm must consider a vector of amplitude parameters 𝝉m2=(τ1,m2,,τm,m2)\bm{\tau}_{m}^{2}=(\tau_{1,m}^{2},...,\tau_{m,m}^{2})^{\prime} and regression coefficients 𝜷m=(𝜷1,m,,𝜷m,m)\bm{\beta}_{m}=(\bm{\beta^{\prime}}_{1,m},...,\bm{\beta^{\prime}}_{m,m}) that we wish to estimate, for the jjth component within a partition of mm segments, j=1,,m.j=1,...,m. For notational convenience, 𝜷j,m,j=1,,m,\bm{\beta}_{j,m},j=1,...,m, is assumed to include the first entry, α0j,m.\alpha_{0j,m}. In the sections that follows, the superscripts cc and pp refer to the current and proposed values in the sampling scheme, respectively. This is consistent with most notation used in Markov chain Monte Carlo algorithms.

We begin by describing the between-model moves, where the number of change points mm may be changed. Let 𝜽m=(𝝃m,𝝉m2,𝜷m)\bm{\theta}_{m}=(\bm{\xi}^{\prime}_{m},\bm{\tau}^{2^{\prime}}_{m},\bm{\beta^{\prime}}_{m}) be the model parameters at some point in the RJMCMC sampling scheme and assume that the chain starts at (mc,𝜽mcc)(m^{c},\bm{\theta}_{m^{c}}^{c}). The algorithm proposes a move to (mp,𝜽mpp)(m^{p},\bm{\theta}_{m^{p}}^{p}), by drawing (mp,𝜽mpp)(m^{p},\bm{\theta}_{m^{p}}^{p}) from a proposal distribution q(mp,𝜽mpp|mc,𝜽mcc)q(m^{p},\bm{\theta}_{m^{p}}^{p}|m^{c},\bm{\theta}_{m^{c}}^{c}). This draw is accepted with probability

α= min {1,p(mp,𝜽mpp|𝒙)q(mc,𝜽mcc|mp,𝜽mpp)p(mc,𝜽mcc|𝒙)q(mp,𝜽mpp|mc,𝜽mcc)},\alpha=\text{ min }\Bigg{\{}1,\frac{p(m^{p},\bm{\theta}_{m^{p}}^{p}|\bm{x})q(m^{c},\bm{\theta}_{m^{c}}^{c}|m^{p},\bm{\theta}_{m^{p}}^{p})}{p(m^{c},\bm{\theta}_{m^{c}}^{c}|\bm{x})q(m^{p},\bm{\theta}_{m^{p}}^{p}|m^{c},\bm{\theta}_{m^{c}}^{c})}\Bigg{\}},

where p()p(\cdot) is a target distribution, namely the product of the likelihood and the prior. The target and proposal distributions vary based on the type of move taken in the sampling scheme. The distribution q(mp,𝜽mpp|mc,𝜽mcc)q(m^{p},\bm{\theta}_{m^{p}}^{p}|m^{c},\bm{\theta}_{m^{c}}^{c}) is described as follows:

q(mp,𝜽mpp|mc,𝜽mcc)=q(mp|mc)q(𝜽mpp|mp,mc,𝜽mcc)=q(mp|mc)q(𝝃𝒎𝒑𝒑,𝝉mp2p,𝜷mpp|mp,mc,𝜽mcc)\displaystyle q(m^{p},\bm{\theta}_{m^{p}}^{p}|m^{c},\bm{\theta}_{m^{c}}^{c})=q(m^{p}|m^{c})q(\bm{\theta}_{m^{p}}^{p}|m^{p},m^{c},\bm{\theta}_{m^{c}}^{c})=q(m^{p}|m^{c})q(\bm{\xi^{p}_{m^{p}}},\bm{\tau}_{m^{p}}^{2p},\bm{\beta}_{m^{p}}^{p}|m^{p},m^{c},\bm{\theta}_{m^{c}}^{c})
=q(mp|mc)q(𝝃mpp|mp,mc,𝜽mcc)q(𝝉mp2p|𝝃mpp,mp,mc,𝜽mcc)q(𝜷mpp|𝝉mp2p,𝝃mpp,mp,mc,𝜽mcc).\displaystyle=q(m^{p}|m^{c})q(\bm{\xi}_{m^{p}}^{p}|m^{p},m^{c},\bm{\theta}_{m^{c}}^{c})q(\bm{\tau}_{m^{p}}^{2p}|\bm{\xi}_{m^{p}}^{p},m^{p},m^{c},\bm{\theta}_{m^{c}}^{c})q(\bm{\beta}_{m^{p}}^{p}|\bm{\tau}_{m^{p}}^{2p},\bm{\xi}_{m^{p}}^{p},m^{p},m^{c},\bm{\theta}_{m^{c}}^{c}).

To draw (mp,𝜽mpp)(m^{p},\bm{\theta}_{m^{p}}^{p}), one must first draw mpm^{p}, and then 𝝃mpp\bm{\xi}_{m^{p}}^{p}, τmp2p, and 𝜷mpp\tau_{m^{p}}^{2p},\text{ and }\bm{\beta}_{m^{p}}^{p}. The number of segments mpm^{p} is drawn from the proposal distribution q(mp|mc)q(m^{p}|m^{c}). Let MM be the maximum number of segments and m2,mincm^{c}_{2,\text{min}} be the number of current segments containing at least 2tmin2t_{min} data points. The proposal is as follows:

q(mp=k|mc)={1/2 if k=mc1,mc+1 and mc1,M,m2,minc01 if k=mc1 and mc=M or m2,minc=01 if k=mc+1 and mc=1\displaystyle q(m^{p}=k|m^{c})=\left\{\begin{array}[]{ll}1/2\text{ if }k=m^{c}-1,m^{c}+1\text{ and }m^{c}\neq 1,M,m_{2,\text{min}}^{c}\neq 0\\ 1\text{ if }k=m^{c}-1\text{ and }m^{c}=M\text{ or }m_{2,\text{min}}^{c}=0\\ 1\text{ if }k=m^{c}+1\text{ and }m^{c}=1\end{array}\right.

Conditional on the proposed number of change points mpm^{p}, a new partition 𝝃mpp\bm{\xi}_{m^{p}}^{p}, a new vector of covariance amplitude parameters 𝝉mp2p\bm{\tau}_{m^{p}}^{2p}, and a new vector of regression coefficients, 𝜷mpp\bm{\beta}_{m^{p}}^{p} are proposed. τ2\tau^{2} is described as a smoothing parameter [60] or amplitude parameter.

Now, we describe the process of the birth of new segments, where the number of proposed change points increases. Suppose that mp=mc+1m^{p}=m^{c}+1. As defined before, a time series partition,

𝝃mpp=(ξ0,mcc,,ξk1,mcc,ξk,mpp,ξk,mcc,,ξmc,mcc)\displaystyle\bm{\xi}^{p}_{m^{p}}=(\xi^{c}_{0,m^{c}},...,\xi^{c}_{k^{*}-1,m^{c}},\xi_{k^{*},m^{p}}^{p},\xi_{k^{*},m^{c}}^{c},...,\xi_{m^{c},m^{c}}^{c})

is drawn from the proposal distribution q(𝝃mpp|mp,mc,𝜽mcc)q(\bm{\xi}_{m^{p}}^{p}|m^{p},m^{c},\bm{\theta}_{m^{c}}^{c}). The algorithm first proposes a partition by selecting a random segment j=kj=k^{*} to split. Then, a random point tt^{*} within the segment j=kj=k^{*} is selected to be the proposed partition point. This is subject to a constraint, ξk1,mcc+tmintξk,mcctmin\xi_{k^{*}-1,m^{c}}^{c}+t_{\text{min}}\leq t*\leq\xi_{k^{*},m^{c}}^{c}-t_{\text{min}}. The proposal distribution is then computed as follows:

q(ξj,mpp=t|mp,mc,𝝃mcc)=p(j=k|mp,mc,𝝃mcc);\displaystyle q(\xi_{j,m^{p}}^{p}=t^{*}|m^{p},m^{c},\bm{\xi}_{m^{c}}^{c})=p(j=k^{*}|m^{p},m^{c},\bm{\xi}_{m^{c}}^{c});
p(ξk,mpp=t|j=k,mp,mc,𝝃mcc)=1m2minc(nk,mc2tmin+1).\displaystyle p(\xi_{k^{*},m^{p}}^{p}=t^{*}|j=k^{*},m^{p},m^{c},\bm{\xi}_{m^{c}}^{c})=\frac{1}{m_{2\text{min}}^{c}(n_{k^{*},m^{c}}-2t_{\text{min}}+1)}.

Then, the vector of amplitude parameters

τmp2p=(τ1,mc2c,,τk1,mc2c,τk,mp2p,τk+1,mp2p,τk+1,mc2c,,τmc,mc2c)\displaystyle\tau_{m^{p}}^{2p}=(\tau_{1,m^{c}}^{2c},...,\tau_{k^{*}-1,m^{c}}^{2c},\tau_{k^{*},m^{p}}^{2p},\tau_{k^{*}+1,m^{p}}^{2p},\tau_{k^{*}+1,m^{c}}^{2c},...,\tau_{m^{c},m^{c}}^{2c})

is drawn from the proposal distribution q(𝝉mp2p|mp,𝝃mpp,mc,𝜽mcc)=q(𝝉mp2p|mp,𝝉mc2c).q(\bm{\tau}_{m^{p}}^{2p}|m^{p},\bm{\xi}_{m^{p}}^{p},m^{c},\bm{\theta}_{m^{c}}^{c})=q(\bm{\tau}_{m^{p}}^{2p}|m^{p},\bm{\tau}_{m^{c}}^{2c}). The algorithm is based on the RJMCMC algorithm of [82]. This draws from a uniform distribution uU[0,1]u\sim U[0,1] and defines τk,mp2p\tau_{k^{*},m^{p}}^{2p} and τk+1,mp2p\tau_{k^{*}+1,m^{p}}^{2p} in terms of uu and τk,mc2c\tau_{k^{*},m^{c}}^{2c} as follows:

τk,mp2p=u1uτk,mc2c;\displaystyle\tau_{k^{*,m^{p}}}^{2p}=\frac{u}{1-u}\tau_{k^{*},m^{c}}^{2c}; (17)
τk+1,mp2p=1uuτk,mc2c.\displaystyle\tau_{k^{*}+1,m^{p}}^{2p}=\frac{1-u}{u}\tau_{k^{*},m^{c}}^{2c}. (18)

Then, the vector of coefficients

𝜷mpp=(𝜷1,mcc,,𝜷k1,mcc,𝜷k,mpp,𝜷k+1,mpp,𝜷k+1,mcc,,𝜷mc,mcc)\displaystyle\bm{\beta}_{m^{p}}^{p}=(\bm{\beta}_{1,m^{c}}^{c},...,\bm{\beta}_{k^{*}-1,m^{c}}^{c},\bm{\beta}_{k^{*},m^{p}}^{p},\bm{\beta}_{k^{*}+1,m^{p}}^{p},\bm{\beta}_{k^{*}+1,m^{c}}^{c},...,\bm{\beta}_{m^{c},m^{c}}^{c})

is drawn from the proposal distribution q(𝜷mpp|𝝉mp2p,𝝃mp2p,mp,mc,𝜽mcc)=q(𝜷mpp|𝝉mp2p,𝝃mpp,mp)q(\bm{\beta}_{m^{p}}^{p}|\bm{\tau}_{m^{p}}^{2p},\bm{\xi}_{m^{p}}^{2p},m^{p},m^{c},\bm{\theta}_{m^{c}}^{c})=q(\bm{\beta}_{m^{p}}^{p}|\bm{\tau}_{m^{p}}^{2p},\bm{\xi}_{m^{p}}^{p},m^{p}). The vectors 𝜷k,mpp\bm{\beta}_{k^{*},m^{p}}^{p} and 𝜷k+1,mpp\bm{\beta}_{k^{*}+1,m^{p}}^{p} are drawn from Gaussian approximations to the respective posterior conditional distributions p(𝜷k,mpp|𝒙kp,τk,mp2p,mp)p(\bm{\beta}_{k^{*},m^{p}}^{p}|\bm{x}_{k^{*}}^{p},\tau_{k^{*},m^{p}}^{2p},m^{p}) and p(𝜷k+1,mpp|𝒙k+1p,τk+1,mp2p,mp)p(\bm{\beta}_{k^{*}+1,m^{p}}^{p}|\bm{x}_{k^{*}+1}^{p},\tau_{k^{*}+1,m^{p}}^{2p},m^{p}), respectively. Here, 𝒙kp\bm{x}_{k^{*}}^{p} and 𝒙k+1p\bm{x}_{k^{*}+1}^{p} refer to the subsets of the time series across segments kk^{*} and k+1k^{*}+1, respectively. Then, 𝝃mpp\bm{\xi}_{m^{p}}^{p} determines 𝒙p=(𝒙kp,𝒙k+1p)\bm{x_{*}}^{p}=(\bm{x}_{k^{*}}^{p^{\prime}},\bm{x}_{k^{*}+1}^{p^{\prime}})^{\prime}. We provide an example for illustration: the coefficient 𝜷k,mpp\bm{\beta}_{k^{*},m^{p}}^{p} is drawn from the Gaussian distribution N(𝜷kmax,Σkmax)N(\bm{\beta}_{k^{*}}^{\text{max}},\Sigma_{k^{*}}^{\text{max}}), where

𝜷kmax=argmax𝜷k,mppp(𝜷k,mpp|𝒙kp,τk,mp2p,mp)\displaystyle\bm{\beta}_{k^{*}}^{\text{max}}=\operatorname*{argmax}_{\bm{\beta}_{k^{*},m^{p}}^{p}}p(\bm{\beta}^{p}_{k^{*},m^{p}}|\bm{x}_{k^{*}}^{p},\tau_{k^{*},m^{p}}^{2p},m^{p})

and

Σkmax={2logp(𝜷k,mpp|𝒙kp,τk,mp2p,mp)𝜷k,mpp𝜷k,mpp|𝜷k,mpp=𝜷kmax}.1\displaystyle\Sigma_{k^{*}}^{\text{max}}=-\Bigg{\{}\partialderivative{\log p(\bm{\beta}_{k^{*},m^{p}}^{p}|\bm{x}_{k^{*}}^{p},\tau_{k^{*},m^{p}}^{2p},m^{p})}{\bm{\beta}_{k^{*},m^{p}}^{p}}{\bm{\beta}_{k^{*},m^{p}}^{p^{\prime}}}\Bigg{|}_{\bm{\beta}_{k^{*},m^{p}}^{p}=\bm{\beta}_{k^{*}}^{\text{max}}}\Bigg{\}}_{.}^{-1}

For the birth move (that is the increase in mm), the probability of acceptance is α=min{1,A}\alpha=\min\{1,A\}, where AA is equal to

|(τk,mp2p,τk+1,mp2p)(τk,mc,u2c)|p(𝜽mpp|𝒙,mp)p(𝜽mpp|mp)p(mp)p(𝜽mpp|𝒙,mp)p(𝜽mcc|mc)p(mc)p(mc|mp)p(𝜷k,mcc)p(mp|mc)p(ξk,mpmp|mp,mc)p(u)p(𝜷k,mpp)p(𝜷k+1,mpp).\displaystyle\Bigg{|}\partialderivative{(\tau_{k^{*},m^{p}}^{2p},\tau_{k^{*}+1,m^{p}}^{2p})}{(\tau_{k^{*},m^{c},u}^{2c})}\Bigg{|}\frac{p(\bm{\theta}_{m^{p}}^{p}|\bm{x},m^{p})p(\bm{\theta}_{m^{p}}^{p}|m^{p})p(m^{p})}{p(\bm{\theta}_{m^{p}}^{p}|\bm{x},m^{p})p(\bm{\theta}_{m^{c}}^{c}|m^{c})p(m^{c})}\frac{p(m^{c}|m^{p})p(\bm{\beta}_{k^{*},m^{c}}^{c})}{p(m^{p}|m^{c})p(\xi_{k^{*},m^{p}}^{m^{p}}|m^{p},m^{c})p(u)p(\bm{\beta}^{p}_{k^{*},m^{p}})p(\bm{\beta}_{k^{*}+1,m^{p}}^{p})}.

Above, p(u)=1,0u1,p(u)=1,0\leq u\leq 1, while p(𝜷k,mpp)p(\bm{\beta}_{k^{*},m^{p}}^{p}) and p(𝜷k+1,mpp)p(\bm{\beta}_{k^{*}+1,m^{p}}^{p}) are the density functions of Gaussian proposal distributions N(𝜷kmax,Σkmax)N(\bm{\beta}_{k^{*}}^{\text{max}},\Sigma_{k^{*}}^{\text{max}}) and N(𝜷k+1max,Σk+1max)N(\bm{\beta}_{k^{*}+1}^{\text{max}},\Sigma_{k^{*}+1}^{\text{max}}), respectively. The above Jacobian is computed as

|(τk,mp2p,τk+1,mp2p)(τk,mc2c,u)|=2τkmc2cu(1u)=2(τk,mpp+τk+1,mpp)2.\displaystyle\bigg{|}\partialderivative{(\tau_{k^{*},m^{p}}^{2p},\tau_{k^{*}+1,m^{p}}^{2p})}{(\tau^{2c}_{k^{*},m^{c}},u)}\bigg{|}=\frac{2\tau_{k^{*}m^{c}}^{2c}}{u(1-u)}=2(\tau_{k^{*},m^{p}}^{p}+\tau_{k^{*}+1,m^{p}}^{p})^{2}.

Now, we move to describe the process of the death of new segments. This is where the number of proposed change points will decrease, or mp=mc1m^{p}=m^{c}-1. A time series partition is denoted

𝝃mpp=(ξ0,mcc,,ξk1,mcc,ξk+1,mcc,,ξmc,mcc),\displaystyle\bm{\xi}_{m^{p}}^{p}=(\xi_{0,m^{c}}^{c},...,\xi_{k^{*}-1,m^{c}}^{c},\xi_{k^{*}+1,m^{c}}^{c},...,\xi_{m^{c},m^{c}}^{c}),

is proposed by randomly selecting a single change point from mc1m^{c}-1 candidates, and then removing it. The change point selected for the removal is denoted j=kj=k^{*}. There are mc1m^{c}-1 possible change points available for removal among the mcm^{c} segments currently in existence in the model. The proposal can choose each change point with equal probability, that is,

q(ξj,mpp|mp,mc,𝝃mcc)=1mc1.\displaystyle q(\xi_{j,m^{p}}^{p}|m^{p},m^{c},\bm{\xi}_{m^{c}}^{c})=\frac{1}{m^{c}-1}.

The updated vector of amplitude parameters

𝝉mp2p=(τ1,mc2c,,τk1,mc2c,τk,mp2c,τk+2,mc2c,,τmc,mc2c)\displaystyle\bm{\tau}_{m^{p}}^{2p}=(\tau_{1,m^{c}}^{2c},...,\tau_{k^{*}-1,m^{c}}^{2c},\tau_{k^{*},m^{p}}^{2c},\tau_{k^{*}+2,m^{c}}^{2c},...,\tau_{m^{c},m^{c}}^{2c})

is then drawn from the proposal distribution q(𝝉mp2p|mp,𝝃mpp,mc,𝜽mcc)=q(𝝉mp2p|mp,𝝉mc2c)q(\bm{\tau}_{m^{p}}^{2p}|m^{p},\bm{\xi}_{m^{p}}^{p},m^{c},\bm{\theta}_{m^{c}}^{c})=q(\bm{\tau}_{m^{p}}^{2p}|m^{p},\bm{\tau}_{m^{c}}^{2c}). Then, one amplitude parameter τk,mp2p\tau_{k^{*},m^{p}}^{2p} is formed from two candidate amplitude parameters, τk,mc2c\tau_{k^{*},m^{c}}^{2c} and τk+1,mc2c\tau_{k^{*}+1,m^{c}}^{2c},. by reversing the equations (17) and (18). Specifically,

τk,mp2p=τk,mc2cτk+1,mc2c.\displaystyle\tau_{k^{*},m^{p}}^{2p}=\sqrt{\tau_{k^{*},m^{c}}^{2c}\tau_{k^{*}+1,m^{c}}^{2c}}.

Finally, the updated vector of regression coefficients,

𝜷mpp=(β1,mcc,,βk1,mcc,βk,mpp,βk+2,mcc,,βmc,mcc)\displaystyle\bm{\beta}_{m^{p}}^{p}=(\beta_{1,m^{c}}^{c},...,\beta_{k^{*}-1,m^{c}}^{c},\beta_{k^{*},m^{p}}^{p},\beta_{k^{*}+2,m^{c}}^{c},...,\beta_{m^{c},m^{c}}^{c})

is drawn from the proposal distribution q(𝜷mpp|𝝉mp2p,𝝃mpp,mp,mc,θmcc)=q(𝜷mpp|𝝉mp2p,𝝃mpp,mp)q(\bm{\beta}_{m^{p}}^{p}|\bm{\tau}_{m^{p}}^{2p},\bm{\xi}_{m^{p}}^{p},m^{p},m^{c},\theta_{m^{c}}^{c})=q(\bm{\beta}_{m^{p}}^{p}|\bm{\tau}_{m^{p}}^{2p},\bm{\xi}_{m^{p}}^{p},m^{p}). The vector of regression coefficients is drawn from a Gaussian approximation to the posterior distribution p(βk,mp|𝒙,τk,mp2p,𝝃mpp,mp)p(\beta_{k^{*},m^{p}}|\bm{x},\tau_{k^{*},m^{p}}^{2p},\bm{\xi}^{p}_{m^{p}},m^{p}) with the same underlying procedure for the vector of coefficients within the birth step. The probability of an acceptance is the inverse of the corresponding birth step. If the move is accepted then the following updates of the current values occur: mc=mpm^{c}=m^{p} and 𝜽mcc=𝜽mpp\bm{\theta}_{m^{c}}^{c}=\bm{\theta}_{m^{p}}^{p}.

Finally, we describe the within-model moves: from this point on, the number of change points mm will be fixed and the notation describing the dependence on the number of change points (explicitly) will be removed. There are two sub-components of the within-model move. In the first piece, the change points could be relocated, and conditional on this relocation, the coefficients of the candidate basis functions will be updated. These steps are jointly accepted or rejected within a Metropolis-Hastings step, and the amplitude parameters are also updated in a separate Gibbs sampling step.

We assume the chain is located at 𝜽c=(𝝃c,𝜷c)\bm{\theta}^{c}=(\bm{\xi}^{c},\bm{\beta}^{c}). The proposal 𝜽p=(𝝃p,𝜷p)\bm{\theta}^{p}=(\bm{\xi}^{p},\bm{\beta}^{p}) is as follows: first, a change point ξk\xi_{k^{*}} is selected for potential relocation from m1m-1 candidate change points. Next, a position within the interval [ξk1,ξk+1][\xi_{k^{*}-1},\xi_{k^{*}+1}] is chosen, subject to the constraint that the new location will be at least tmint_{\text{min}} data points away from ξk1\xi_{k^{*}-1} and ξk+1\xi_{k^{*}+1}, so that

Pr(ξkp=t)=Pr(j=k)Pr(ξkp=t|j=k),\displaystyle\Pr(\xi^{p}_{k^{*}}=t)=\Pr(j=k^{*})\Pr(\xi_{k^{*}}^{p}=t|j=k^{*}),

where Pr(j=k)=(m1)1\Pr(j=k^{*})=(m-1)^{-1}. A mixture distribution for Pr(ξkp=t|j=k)\Pr(\xi_{k^{*}}^{p}=t|j=k^{*}) is constructed to explore the space efficiently, so

Pr(ξkp=t|j=k)=πq1(ξkp=t|ξkc)+(1π)q2(ξkp=t|ξkc),\displaystyle\Pr(\xi_{k^{*}}^{p}=t|j=k^{*})=\pi q_{1}(\xi_{k^{*}}^{p}=t|\xi_{k^{*}}^{c})+(1-\pi)q_{2}(\xi_{k^{*}}^{p}=t|\xi_{k^{*}}^{c}),

where q1(ξkp=t|ξkc)=(nk+nk+12tmin+1)1,ξk1+tmintξk+1tminq_{1}(\xi_{k^{*}}^{p}=t|\xi_{k^{*}}^{c})=(n_{k^{*}}+n_{k^{*}+1}-2t_{\text{min}}+1)^{-1},\xi_{k^{*}-1}+t_{\text{min}}\leq t\leq\xi_{k^{*}+1}-t_{\text{min}} and

q2(ξkp=t|ξkc)={0 if |tξkc|>1;1/3 if |tξkc|1,nktmin and nk+1tmin;1/2 if tξkc1,nk=tmin and nk+1tmin;1/2 if ξkct1,nktmin and nk+1=tmin;1 if t=ξkc,nk=tmin and nk+1=tmin.\displaystyle q_{2}(\xi_{k^{*}}^{p}=t|\xi_{k^{*}}^{c})=\left\{\begin{array}[]{ll}0\text{ if }|t-\xi^{c}_{k^{*}}|>1;\\ 1/3\text{ if }|t-\xi^{c}_{k^{*}}|\leq 1,n_{k^{*}}\neq t_{\text{min}}\text{ and }n_{k^{*}+1}\neq t_{\text{min}};\\ 1/2\text{ if }t-\xi_{k^{*}}^{c}\leq 1,n_{k^{*}}=t_{\text{min}}\text{ and }n_{k^{*}+1}\neq t_{\text{min}};\\ 1/2\text{ if }\xi_{k^{*}}^{c}-t\leq 1,n_{k^{*}}\neq t_{\text{min}}\text{ and }n_{k^{*}+1}=t_{\text{min}};\\ 1\text{ if }t=\xi_{k^{*}}^{c},n_{k^{*}}=t_{\text{min}}\text{ and }n_{k^{*}+1}=t_{\text{min}}.\end{array}\right.

The support of q1q_{1} has nk+nk+12tmin+1n_{k^{*}}+n_{k^{*}+1}-2t_{\text{min}}+1 points while q2q_{2} has at most three. The term q2q_{2} alone may result in a meaningfully higher acceptance rate for the M-H step, but it will slowly explore the parameter space. This issue is exacerbated in higher dimensional settings. The q1q_{1} component will allow for larger steps, and can produce a compromise between higher acceptance rates and a deep exploration of parameter space.

Then, 𝜷𝒋𝒑,j=k,k+1\bm{\beta^{p}_{j}},j=k^{*},k^{*}+1 is drawn from an approximation to j=kk+1p(𝜷j|𝒙jp,τj2)\prod^{k^{*}+1}_{j=k^{*}}p(\bm{\beta}_{j}|\bm{x}_{j}^{p},\tau_{j}^{2}), following the corresponding step in the between-model move. The proposal distribution, evaluated at 𝜷jp,j=k,k=1\bm{\beta}^{p}_{j},j=k^{*},k^{*}=1, is

q(𝜷p|𝒙p,𝝉2)=j=kk+1q(𝜷jp|𝒙jp,τj2),\displaystyle q(\bm{\beta}_{*}^{p}|\bm{x}_{*}^{p},\bm{\tau}_{*}^{2})=\prod^{k^{*}+1}_{j=k^{*}}q(\bm{\beta}_{j}^{p}|\bm{x}_{j}^{p},\tau_{j}^{2}),

where 𝜷p=(𝜷kp,𝜷k+1p)\bm{\beta}_{*}^{p}=(\bm{\beta}^{p^{\prime}}_{k^{*}},\bm{\beta}^{p^{\prime}}_{k^{*}+1})^{\prime} and 𝝉2=(τk2,τk+12)\bm{\tau}_{*}^{2}=(\tau^{2}_{k^{*}},\tau^{2}_{k^{*}+1})^{\prime}. This proposal distribution is also evaluated at current values of 𝜷c=(βkc,βk+1c)\bm{\beta}_{*}^{c}=(\beta^{c^{\prime}}_{k^{*}},\beta^{c^{\prime}}_{k^{*}+1})^{\prime}. βp\beta_{*}^{p} is then accepted with probability

α=min{1,p(𝒙p|𝜷p)p(𝜷p|𝝉2)q(𝜷c|𝒙c,𝝉2)p(𝒙c|𝜷c)p(𝜷c|𝝉2)q(𝜷p|𝒙p,𝝉2)},\alpha=\min\Bigg{\{}1,\frac{p(\bm{x}_{*}^{p}|\bm{\beta}_{*}^{p})p(\bm{\beta}_{*}^{p}|\bm{\tau}_{*}^{2})q(\bm{\beta}_{*}^{c}|\bm{x}^{c}_{*},\bm{\tau}_{*}^{2})}{p(\bm{x}_{*}^{c}|\bm{\beta}_{*}^{c})p(\bm{\beta}_{*}^{c}|\bm{\tau}_{*}^{2})q(\bm{\beta}_{*}^{p}|\bm{x}^{p}_{*},\bm{\tau}_{*}^{2})}\Bigg{\}},

where 𝒙c=(𝒙kc,𝒙k+1c)\bm{x}_{*}^{c}=(\bm{x}^{c^{\prime}}_{k^{*}},\bm{x}^{c^{\prime}}_{k^{*}+1}). When the draw is accepted, we update the partition and regression coefficients (ξkc,βc)=(ξkp,βp)(\xi^{c}_{k^{*}},\beta_{*}^{c})=(\xi^{p}_{k^{*}},\beta_{*}^{p}). Finally, we draw τ2p\tau^{2p} from

p(τ2|𝜷)=j=kk+1p(τj2|βj).\displaystyle p(\tau_{*}^{2}|\bm{\beta}_{*})=\prod^{k^{*}+1}_{j=k^{*}}p(\tau_{j}^{2}|\beta_{j}).

This is a Gibbs sampling step, and as such the draw is accepted with probability 1.

References

  • 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.
  • Miceli and Susinno [2004] M. Miceli, G. Susinno, Ultrametricity in fund of funds diversification, Physica A: Statistical Mechanics and its Applications 344 (2004) 95–99. doi:10.1016/j.physa.2004.06.094.
  • Conlon et al. [2007] T. Conlon, H. Ruskin, M. Crane, Random matrix theory and fund of funds portfolio optimisation, Physica A: Statistical Mechanics and its Applications 382 (2007) 565–576. doi:10.1016/j.physa.2007.04.039.
  • Bouchaud and Potters [2003] J.-P. Bouchaud, M. Potters, Theory of Financial Risk and Derivative Pricing, Cambridge University Press, 2003. doi:10.1017/cbo9780511753893.
  • Burda et al. [2004] Z. Burda, A. Görlich, A. Jarosz, J. Jurkiewicz, Signal and noise in correlation matrix, Physica A: Statistical Mechanics and its Applications 343 (2004) 295–310. doi:10.1016/j.physa.2004.05.048.
  • Sharifi et al. [2004] S. Sharifi, M. Crane, A. Shamaie, H. Ruskin, Random matrix theory for portfolio optimization: a stability approach, Physica A: Statistical Mechanics and its Applications 335 (2004) 629–643. doi:10.1016/j.physa.2003.12.016.
  • 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.
  • Tilfani et al. [2022] O. Tilfani, L. Kristoufek, P. Ferreira, M. Y. E. Boukfaoui, Heterogeneity in economic relationships: Scale dependence through the multivariate fractal regression, Physica A: Statistical Mechanics and its Applications 588 (2022) 126530. doi:10.1016/j.physa.2021.126530.
  • 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.
  • Kwapień et al. [2021] J. Kwapień, M. Wątorek, S. Drożdż, Cryptocurrency market consolidation in 2020–2021, Entropy 23 (2021) 1674. doi:10.3390/e23121674.
  • 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 901 (2020) 1–82. doi:10.1016/j.physrep.2020.10.005.
  • 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. [2020a] 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 (2020a) 023122. doi:10.1063/1.5139634.
  • Drożdż et al. [2020b] S. Drożdż, J. Kwapień, P. Oświęcimka, T. Stanisz, M. Wątorek, Complexity in economic and social systems: Cryptocurrency market at around COVID-19, Entropy 22 (2020b) 1043. doi:10.3390/e22091043.
  • 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.
  • 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.
  • Engle [1982] R. Engle, Autoregressive Conditional Heteroskedasticity with Estimates of the Variance of United Kingdom Inflation, Econometrica 50 (1982) 987–1007.
  • Bollerslev [1986] T. Bollerslev, Generalized Autoregressive Conditional Heteroskedasticity, Journal of Econometrics 31 (1986) 307–327.
  • Taylor [1982] S. Taylor, Financial Returns Modelled by the product of two stochastic processes, a study of daily sugar prices 1961-1979, Time Series Analysis: Theory and Practice 1 (1982) 203–226.
  • Taylor [1986] S. Taylor, Modelling Financial Time Series, New York, Wiley, 1986.
  • Taylor [1994] S. J. Taylor, Modelling Stochastic Volatility: A review and comparative study, Mathematical Finance 4 (1994) 183–204. doi:10.1111/j.1467-9965.1994.tb00057.x.
  • Nelson [1991] D. Nelson, Conditional heteroskedasticity in asset returns: A new approach, Econometrica 59 (1991) 347–370.
  • Glosten et al. [1993] L. R. Glosten, R. Jagannathan, D. E. Runkle, On the relation between the expected value and the volatility of the nominal excess return on stocks, The Journal of Finance 48 (1993) 1779–1801. doi:10.1111/j.1540-6261.1993.tb05128.x.
  • Zakoian [1994] J. Zakoian, Threshold heteroskedastic models, Journal of Economic Dynamics and Control 18 (1994) 931–955.
  • So et al. [2002] M. So, W. Li, K. Lam, On a Threshold Stochastic Volatility Model, Journal of Forecasting 22 (2002) 473–500.
  • Cai [1994] J. Cai, A markov model of switching-regime arch, Journal of Business & Economics Statistics 12 (1994) 309–316.
  • Hamilton and Susmel [1994] J. Hamilton, R. Susmel, Autoregressive conditional heteroskedasticity and changes in regime, Journal of Econometrics 64 (1994) 307–333.
  • Gray [1996] S. Gray, Modelling the conditional distribution of interest rates as a regime-switching process, Journal of Financial Econometrics 2 (1996) 211–250.
  • So et al. [1998] M. So, K. Lam, W. Li, A stochastic volatility model with markov switching, Journal of Business & Economic Statistics 16 (1998) 244–253.
  • Andersen et al. [2003] T. Andersen, T. Bollerslev, F. Diebold, P. Labys, Modeling and forecasting realized volatility, Econometrica 71 (2003) 579–625.
  • Hwang and Satchell [2007] S. Hwang, S. Satchell, GARCH Model with Cross-sectional Volatility: GARCHX Models, Applied Financial Economics 15 (2007) 203–216.
  • Hansen et al. [2011] P. R. Hansen, Z. Huang, H. H. Shek, Realized GARCH: a joint model for returns and realized measures of volatility, Journal of Applied Econometrics 27 (2011) 877–906. doi:10.1002/jae.1234.
  • Cerqueti et al. [2020] R. Cerqueti, M. Giacalone, R. Mattera, Skewed non-gaussian GARCH models for cryptocurrencies volatility modelling, Information Sciences 527 (2020) 1–26. doi:10.1016/j.ins.2020.03.075.
  • Wan and Si [2017] Y. Wan, Y.-W. Si, A formal approach to chart patterns classification in financial time series, Information Sciences 411 (2017) 151–175. doi:10.1016/j.ins.2017.05.028.
  • Stehlík et al. [2017] M. Stehlík, C. Helperstorfer, P. Hermann, J. Šupina, L. Grilo, J. Maidana, F. Fuders, S. Stehlíková, Financial and risk modelling with semicontinuous covariances, Information Sciences 394-395 (2017) 246–272. doi:10.1016/j.ins.2017.02.002.
  • Chu et al. [1996] C.-S. J. Chu, G. J. Santoni, T. Liu, Stock market volatility and regime shifts in returns, Information Sciences 94 (1996) 179–190. doi:10.1016/0020-0255(96)00117-x.
  • yong Chen et al. [2018] G. yong Chen, M. Gan, G. long Chen, Generalized exponential autoregressive models for nonlinear time series: Stationarity, estimation and applications, Information Sciences 438 (2018) 46–57. doi:10.1016/j.ins.2018.01.029.
  • Cerqueti et al. [2019] R. Cerqueti, M. Giacalone, D. Panarello, A generalized error distribution copula-based method for portfolios risk assessment, Physica A: Statistical Mechanics and its Applications 524 (2019) 687–695. doi:10.1016/j.physa.2019.04.077.
  • Dahlhaus [1997] R. Dahlhaus, Fitting time series models to nonstationary processes, The Annals of Statistics 25 (1997) 1–37. doi:10.1214/aos/1034276620.
  • Adak [1998] S. Adak, Time-dependent spectral analysis of nonstationary time series, Journal of the American Statistical Association 93 (1998) 1488–1501. doi:10.1080/01621459.1998.10473808.
  • Rosen et al. [2009] O. Rosen, D. S. Stoffer, S. Wood, Local spectral analysis via a Bayesian mixture of smoothing splines, Journal of the American Statistical Association 104 (2009) 249–262. doi:10.1198/jasa.2009.0118.
  • Rosen et al. [2012] O. Rosen, S. Wood, D. S. Stoffer, AdaptSPEC: Adaptive spectral estimation for nonstationary time series, Journal of the American Statistical Association 107 (2012) 1575–1589. doi:10.1080/01621459.2012.716340.
  • Rosen et al. [2017] O. Rosen, S. Wood, D. Stoffer, BayesSpec: Bayesian Spectral Analysis Techniques, 2017. URL: https://CRAN.R-project.org/package=BayesSpec, R package version 0.5.3.
  • James and Menzies [2021] N. James, M. Menzies, Optimally adaptive Bayesian spectral density estimation for stationary and nonstationary processes, 2021. arXiv:2003.02367.
  • 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.
  • Kwapień and Drożdż [2012] J. Kwapień, S. Drożdż, Physical approach to complex systems, Physics Reports 515 (2012) 115–226. doi:10.1016/j.physrep.2012.01.007.
  • Brody et al. [1981] T. A. Brody, J. Flores, J. B. French, P. A. Mello, A. Pandey, S. S. M. Wong, Random-matrix physics: spectrum and strength fluctuations, Reviews of Modern Physics 53 (1981) 385–479. doi:10.1103/revmodphys.53.385.
  • Guhr et al. [1998] T. Guhr, A. Müller–Groeling, H. A. Weidenmüller, Random-matrix theories in quantum physics: common concepts, Physics Reports 299 (1998) 189–425. doi:10.1016/s0370-1573(97)00088-4.
  • Bryson et al. [2021] J. Bryson, R. Vershynin, H. Zhao, Marchenko–pastur law with relaxed independence conditions, Random Matrices: Theory and Applications (2021) 2150040. doi:10.1142/s2010326321500404.
  • Marčenko and Pastur [1967] V. A. Marčenko, L. A. Pastur, Distribution of eigenvalues for some sets of random matrices, Mathematics of the USSR-Sbornik 1 (1967) 457–483. doi:10.1070/sm1967v001n04abeh001994.
  • Wishart [1928] J. Wishart, The Generalised product moment distribution in samples from a normal multivariate population, Biometrika 20A (1928) 32–52. doi:10.1093/biomet/20a.1-2.32.
  • Recher et al. [2010] C. Recher, M. Kieburg, T. Guhr, Eigenvalue densities of real and complex wishart correlation matrices, Physical Review Letters 105 (2010). doi:10.1103/physrevlett.105.244101.
  • Hachem et al. [2016] W. Hachem, A. Hardy, J. Najim, Large complex correlated wishart matrices: the pearcey kernel and expansion at the hard edge, Electronic Journal of Probability 21 (2016). doi:10.1214/15-ejp4441.
  • Wirtz and Guhr [2013] T. Wirtz, T. Guhr, Distribution of the smallest eigenvalue in the correlated wishart model, Physical Review Letters 111 (2013). doi:10.1103/physrevlett.111.094101.
  • Waltner et al. [2015] D. Waltner, T. Wirtz, T. Guhr, Eigenvalue density of the doubly correlated wishart model: exact results, Journal of Physics A: Mathematical and Theoretical 48 (2015) 175204. doi:10.1088/1751-8113/48/17/175204.
  • James and Menzies [2021] N. James, M. Menzies, A new measure between sets of probability distributions with applications to erratic financial behavior, Journal of Statistical Mechanics: Theory and Experiment 2021 (2021) 123404. doi:10.1088/1742-5468/ac3d91.
  • Green [1995] P. J. Green, Reversible jump Markov chain Monte Carlo computation and Bayesian model determination, Biometrika 82 (1995) 711–732. doi:10.1093/biomet/82.4.711.