A VECM analysis of Bitcoin price using time-varying cointegration approach

Yong Lee (Soongsil University, Seoul, Republic of Korea)
Joon Hee Rhee (Soongsil University, Seoul, Republic of Korea)

Journal of Derivatives and Quantitative Studies: 선물연구

ISSN: 1229-988X

Article publication date: 1 July 2022

Issue publication date: 25 July 2022

2316

Abstract

This study proposed an optimal model to examine the relationship between the Bitcoin price and six macroeconomic variables – the Bitcoin price, Standard and Poor's 500 volatility index, US treasury 10-year yield, US consumer price index, gold price and dollar index. It also examined the effectiveness of the vector error correction model (VECM) in analyzing the interrelationship among these variables. The authors employed the following approach: first, the authors sampled the period August 2010–February 2022. This is because Bitcoin achieved a market capitalization of more than US$1 tn over this period, gaining market attention and acceptance from retail, corporate and institutional investors. Second, the authors employed a VECM with the six macroeconomic variables. Finally, the authors expanded the long-run equilibrium relationship (time-invariant cointegration)-based VECM to develop a time-varying cointegration (TVC) VECM. The authors estimated the TVC VECM using the Chebyshev polynomial specification based on various information criteria. The results showed that the Bitcoin price can be modeled with the VECM (p = 1, r = 1). The TVC approach generated more explanatory power for Bitcoin pricing, indicating the effectiveness of the approach for modeling the long-run relationship between Bitcoin price and macroeconomic variables.

Keywords

Citation

Lee, Y. and Rhee, J.H. (2022), "A VECM analysis of Bitcoin price using time-varying cointegration approach", Journal of Derivatives and Quantitative Studies: 선물연구, Vol. 30 No. 3, pp. 197-218. https://doi.org/10.1108/JDQS-01-2022-0001

Publisher

:

Emerald Publishing Limited

Copyright © 2022, Yong Lee and Joon Hee Rhee

License

Published in Journal of Derivatives and Quantitative Studies: 선물연구. Published by Emerald Publishing Limited. This article is published under the Creative Commons Attribution (CC BY 4.0) licence. Anyone may reproduce, distribute, translate and create derivative works of this article (for both commercial and non-commercial purposes), subject to full attribution to the original publication and authors. The full terms of this licence maybe seen at http://creativecommons.org/licences/by/4.0/legalcode.


1. Introduction

Bitcoin exceeded a market capitalization of US$1 tn in February 2021, thereby attracting the interests of institutional investors and listed corporations in the cryptocurrency market (Coinmarketcap https://coinmarketcap.com/). Its market capitalization can be compared to the market capitalization rate of gold, silver, Apple and Samsung Electronics at 6.7, 58.6, 29.7, and 208%, respectively, as of March 2022 (8marketcap https://www.8marketcap.com/). In addition to this appreciation in the Bitcoin's market capitalization since 2010, there has been growing debate on Bitcoin's price, which is fundamental to the analysis of Bitcoin's price model.

In this context, while some studies use the security market analysis theory, such as the capital asset pricing model (CAPM) and the Fama–French multi-factor model, to predict the Bitcoin Price, some others use historical time series-based models, such as the univariate autoregressive integrated moving average (ARIMA) model and the multivariate vector autoregressive (VAR) models to explain Bitcoin's price dynamics.

From the perspective of the Bitcoin price metric, there has been a rise in the Bitcoin price since the mining of the genesis block in early 2009. Specifically, the price increased from $1 to $10 to $100 to below $1000 by April 2011, August 2012, April 2013, and the end of 2016, respectively, and it has been hovering at $10,000 since July 2020. From the market capitalization perspective, Bitcoin has staged only 10% of the time horizon from our study period. This study captures the period when Bitcoin exceeded a market capitalization of US$100 bn computed based on the total number of mined coins, which holds significance for the financial market.

Despite Bitcoin's price and market capitalization history and its higher return/volatility characteristics as an emerging asset class (Figure 1), previous studies have failed to analyze Bitcoin price in relation to the macro, financial and economic variables. Most of the Bitcoin price analysis, in mid-2010, have reported the significant relationship between the Bitcoin price and the factors (e.g. search volume and media exposure) extracted from the social network service-based statistics. However, they have paid little attention to the causality between the macroeconomic variables. Despite functioning as an alternative investment asset class facilitating portfolio diversification and yield enhancement, studies have produced inconsistent outcomes on the relationship between the Bitcoin price and the macroeconomic variables explaining the rapid price appreciation and market adoption (see Figure 2).

Based on the previous Bitcoin research, we adopt the following approach to analyze the Bitcoin price. First, we extend the data observation period to early 2022 when Bitcoin was more exposed to financial market interactions. Second, we interpret data using the VAR and vector error correction model (VECM) and the impulse response function (IRF) to analyze the explanatory power among the macroeconomic variables. Given that most of the explanatory variables are integrated of order one (I(1)), we conduct a cointegration test to determine the applicability of VECM. The results lead to the selection of the VECM over the VAR. Third, under the general VECM structure, the cointegration equation shows a long-run equilibrium with a constant coefficient (termed “time-invariant cointegration” [TIC]). However, this study presents a VECM model based on a time-varying coefficient (termed “time-varying cointegration” [TVC]). To check the significance of the TVC VECM model against the general TIC VECM, we test the null hypothesis using the likelihood ratio (LR) test across multiple TVC parameters. Based on the results, we propose an optimal TVC parameter based on information criteria.

The remaining paper is organized as follows. Chapter 2 reviews the literature review, and Chapter 3 presents the data and methodology. Chapter 4 estimates the TIC VECM and TVC VECM and explains the hypothesis used to report the optimal model parameter. Chapter 5 concludes the study.

2. Literature review

Based on the correlation between the Bitcoin price and the macroeconomic variables, Son and Kim (2019) presented Bitcoin as a safe-haven asset during high-volatility periods. The study considered the period from January 2019 to June 2019. It sourced data on the Bitcoin price (Korean Won) from Bithumb and on the macroeconomic variables from the Federal Reserve Economic Data of the Federal Reserve Bank of St. Louis. The results showed a low correlation between gold and the Bitcoin price. They did not show Granger causality between the price and the six macroeconomic variables – the Bitcoin price, Standard and Poor's 500 volatility index, US treasury 10-year yield, US consumer price index, gold price and dollar index.

Lee et al. (2019a, b) empirically analyzed the determinants of the Bitcoin price in the Korean market for the period July 2015–May 2018. They classified the pricing factors into three groups composed of 12 variables – supply/demand (7), real economic (4), and psychological factors (1) – and examined the effect of these factors on the Bitcoin price. They employed five regression analysis models. The results showed a significant relationship of the price only with one psychological variable extracted from the Naver Trend Index. They also showed a significant negative relationship between the net stock purchase from the retail investors and the price. This finding implied the role of Bitcoin as an alternative investment from retail investors' standpoint.

Kim et al. (2019) analyzed the mutual relationship between Bitcoin and real, speculative and currency assets using VECM and daily data from July 2010 to April 2018. They found that Bitcoin is more closely associated with the speculative than with real or currency assets.

Bouoiyour and Selmi (2015) analyzed the Bitcoin pricing factors using the autoregressive distributed lag (ARDL) bounds test with data spanning from May 2010 to July 2014. They showed Bitcoin's strong speculative characteristics and its weak role as a safe-haven asset.

Zhu et al. (2017) performed Bitcoin analysis using VECM. From 2011 to 2016, they used monthly data on macroeconomic variables such as the price index, dollar index, Dow Jones industrial average, index, Fed fund rate, gold price and bitcoin price. They showed that the economic variables have explanatory power to Bitcoin price on a long-term basis. They also that the dollar index most significantly influenced the Bitcoin price movement, while gold had the least influence. They concluded that the Bitcoin as a speculative asset may not be a stable currency.

Thaker and Mand (2021) examined the relationship between the Bitcoin price and the Asian stock market index using the VECM, generalized autoregressive conditional heteroskedasticity in the mean-GARCH and wavelet analysis. They found a long-term significant relationship between the Bitcoin price and the stock market index.

Syafiqah and Mohamad (2021) included various cryptocurrencies (e.g. Bitcoin, Monero and Stella) to model the long-term relationship between cryptocurrencies and economic variables, using VECM and monthly data from January 2016 to December 2020. They found that only Monero had a significant long-run relationship with macroeconomic variables.

Yang et al. (2020) analyzed the price discovery relationship between the Bitcoin spot and futures price employing the time-varying Granger causality, cointegration and information sharing approaches. They found that the Bitcoin Futures price granger causes the Bitcoin spot price and that the former contributes more to the price discovery function than the latter.

This study is different from the previous research in the following aspects. First, it takes a longer data observation period of 10 years covering the US$1 trillion market capitalization period. Second, it expands the general VECM used in the initial stage of the research to include the TVC VECM in order to capture the dynamic nature of the long-term equilibrium relationship embedded in the VECM.

In summary, until 2015, the research used the GARCH and ARDL to determine the relationship between the Bitcoin price and the macroeconomic variables. Since 2017, studies have been using the VECM and TVC for exploring the long-term relationship between the variables of interest. However, this study expands the Bitcoin price model based on the TVC VECM to explain the dynamic long-run equilibrium between the Bitcoin price and the macroeconomic variables.

3. Methodology

3.1 Analysis flow

Despite the short history and idiosyncratic nature of its functioning, Bitcoin has been significantly influencing the financial markets. This background motivates our Bitcoin price analysis based on the VECM approach. In this context, we consider the dynamic nature of the macroeconomic variables and examine the long-run relationship by expanding the VECM to include the TVC approach. To test the significance of the TVC approach, we set a null hypothesis. This is because of the same cointegration coefficient relationship between the TIC VECM and the TVC VECM. We estimate the coefficient of the TVC VECM using the Chebyshev polynomials. We also propose an optimal TVC VECM under any Chebyshev polynomial specification, based on various information criteria outcomes (Bierens and Martins, 2010).

We employ the following steps to apply the VECM. First, we check the stationarity of the VAR variables using the augmented Dickey–Fuller (ADF) tests. Based on the test, we configure the variables as an I(0) or I(1) variable. For an appropriate lag selection for the ADF test, we use the Akaike information criterion (AIC) test. Second, we apply the VAR estimation with numerous information criteria outcomes to select a VAR lag model including the intercept adoption. Third, we estimate the VAR (p = lag) model. Fourth, given that most of the macroeconomic variables are I(1), we conduct Johansen tests. This helps us to set up the model correctly in the presence of cointegration among I(1) variables. Fifth, after confirming the presence of cointegration, we set the VECM model as the analysis model. We validate the model from various standpoints by employing additional tests such as the test for the significance of the model coefficient and the test for the characteristics of the residuals generated from the model. The test for residuals included the tests for heteroscedasticity, autocorrelation of residuals and normality of residuals. Sixth, based on the estimated VECM model IRF and the forecast error variance decomposition (FEVD), we explain the relationship among the selected variables. Finally, we adjust for the time-varying long-run relationship of the cointegration equation by using the Chebyshev polynomials. This included the following steps. First, we estimate the TVC VECM model (Bierens and Martins, 2010). Second, we run a hypothesis of the TIC VECM over the TVC VECM using an LR test. Third, we choose optimal the TVC VECM model with some Chebyshev polynomials of (m) to illustrate the trend of each TVC VECM long-run cointegration parameter alongside the polynomial (m).

3.2 Data

The dataset comprised six macroeconomic variables (https://data.nasdaq.com) sampled monthly from August 2010 to February 2022. Table 1 presents the summary statistics of these variables.

The six economic variables are the Bitcoin price (in US dollars; BTCUSD), Standard and Poor's (S&P) 500 volatility index (VIX), US treasury 10-year yield (US 10Y), US consumer price index (US CPI), gold price and dollar index. VIX derived from option price of the S&P500 index is a volatility index representing broad market uncertainties moving in the opposite direction of the stock market or risky asset. US10Y standing for the US treasury's 10-year yield is regarded as the constant maturity treasury (CMT). Rate operates as a safe asset during high volatility periods and absorbs the inflationary pressure with rising yields. BTCUSD is the log of the Bitcoin price, USCPI stands for the US CPI, GOLD stands for the gold price in US$, and DXY represents the dollar index.

3.3 VECM and cointegration

As the starting point for the multivariate analysis, we use VAR, where this model assumes the stationarity of the times series dataset. If there are k variables for the analysis with time lag = p (VAR (lags = p)), the VAR model can be defined as below (Sims, 1980).

(1)Yt=c+Γ1Yt1+Γ2Yt2++ΓpYtp+εt where t=1,,T

From Equation (1), Yt denotes the value of Y at t as the vector of kx1; Yt1 depicts the value of Yt with the time lag = 1; c denotes the constant value with vector kx1; Γi represents a parameter vector of kxk with tune lag of i = 1, …,p and εt is a white noise with vector of kx1.

Given that VAR(p) model assumes the stationarity of the time series variables, if these variables fail to achieve stationarity at the raw level during the analysis, we differentiate the time series until the time series achieves stationarity and facilitates the application of VAR(p) model. However, the stationarity of the underlying variables can be obtained by differentiating the nonstationary time series at the expense of the cointegration information among variables. Hence, we use the Johansen test to check for the presence of cointegration among variables. This helps in validating the eligibility of the VECM.

Using the test results in Table 2, we conduct the ADF to test the time series stationarity (ADF lags were chosen based on the AIC).

According to ADF test result, all variables, except for USCPI, were reported as nonstationary, under no intercept condition (test regression: none). When we included the Trend term, the ADF test reported USCPI and DXY as nonstationary time series. The first differencing of all variables led to stationarity in the time series, except for USCPI with the Trend condition.

From Equation (1), if the solution of det(IkΓ1 z   Γpzp) ≠0 for |z|1 with z = 1 case have an unit root, all or part of variables of VAR(p) model can be I(1) variables. Under this case, the VECM could be appropriate (Pfaff, 2008).

The long-run relationship between the Bitcoin and the macroeconomic variables under VECM can be explained by the linear cointegration vector. In this case (transitory), the VECM can be represented as follows (Engle and Granger, 1987).

(2)ΔYt=ΠYt1+j=1p1ΓjΔYtj+εt where ΔYt=YtYt1
(3)Π=αβ=(IΓ1... Γp)

From Equation (2), Γj stands for the transitory relationship, under the reduced rank of αβ matrix in the presence of cointegration. α means loading matrix (or weights), and β of k x r dimension matrix implies the long-run relationship among variables. r represents the number of cointegration equations.

3.4 Time varying cointegration

Utilizing the VECM from Equation (2), the TVC VECM can be modeled as follows.

(4)ΔYt=ΠtYt1+j=1p1ΓjΔYtj+εt where Πt=αβt
Πt is subjected to βt , contrary to Π from Equation (2), where it depends on the time(t) and is defined as βt=f(βVIX,t,βUS10Y,t,βBTCUSD,t,βUSCPI,t,βGOLD,t,βDXY,t). From this expression, we know that β presents TIC, while βt includes the TVC information.

Our study uses the Chebyshev polynomials with different m (dimension) to estimate smoothly changing βt, in line with Bierens and Martins (2010). Specifically, we approximate TVC βt by βt(m), and this can be presented as follows:

(5)βt(m)=i=0mξiPi,T(t)
where ξi stands for the Fourier coefficient,  m denotes the order of the Chebyshev polynomials and βt with i(=0,,) are estimated as smoothed values [1]. Substituting the above information into Equation (4), we obtain Equation (6).
(6)ΔYt=α(i=0mξiPi,T(t))Yt1+j=1p1ΓjΔYtj+εt

We use the LR test to test for the effectiveness of the TVC VECM over the TIC VECM. This test is expressed as below in Equation (7).

(7)LR TVC=LRtvc=2[l^T(r,0)l^T(r,m)]

We use the following values for the hypothesis test: l^T(r,0) from the TIC VECM(p), the log-likelihood value [2] from Equation (2), l^T(r,m) from the TVC VECM(p) and the log-likelihood value from Equation (4) [3].

4. Empirical analysis

4.1 VECM and cointegration

To study the Bitcoin pricing model using macro-economic variables, we test the VAR lag selection procedure without and with the constant + Trend term. We also use the following information criteria: the AIC, final prediction error (FPE), Schwarz information criterion (SIC) and Hannan–Quinn (HQ) information criterion. The AIC and FPE produce Lag = 2, while the HQ and SC produce Lag = 1 (see Table 3).

Our analysis adopted the AIC of Lag = 2 and VAR (p = 2) with Trend model. According to the ADF test, most of variables were I(1). Additionally, we use the Johansen test to check for the presence of cointegration. This helps in the selection of a suitable model between VAR and VECM. The Johansen test utilized the trace and Eigen values. The test confirmed the existence of a single cointegration equation; this presence is regardless of Trend term being under the 5% significance level. Table 4 summaries the outcomes of the Johansen test.

Based on the results of the Johansen test, we select the VECM (p = 1, r = 1) [4]. In the following sections, we present the weights (or loading matrix) and coefficient of cointegration from the VECM (see Table 5).

Concerning the significance of the loading matrix (or weights), VIX, BTCUSD and GOLD were significant at the 5% level with the Trend term. In the absence of the constant term, only VIX and GOLD were significant at the 5% level (see Table 6).

With the loading matrix and the cointegrating coefficient, we present ΠYt1=αβYt1 part of the Bitcoin (ΔBTCUSDt) equation using Equation (2) (with Trend model).

ΠYt1 of ΔBTCUSDt=0.0110|0.2281,1.0000,4.5776,2.7270,0.3617,0.0443,0.6762||Trendt1VIXt1US10Yt1BTCUSDt1USCPIt1GOLDt1DXYt1|

F-statistics of ΔBTCUSDt equation was estimated to be 4.322 (P-value: 0.0001). It implies the significance of the model equation for the Bitcoin price is significant. However, we noted a variance in the individual coefficients. For example, the constant, 1st differenced BTCUSD (lag = 1) and 1st differenced DXY (lag = 1) were significant at the 5, 1 and 10% levels, respectively (see Table 7).

Table 8 summarizes the residual analysis for the VECM (p = 1, r = 1). The heteroscedasticity test failed to reject the null hypothesis, implying the absence of an ARCH effect. The Portmanteau test (PT) for serial correlation also failed to reject the null hypothesis, implying the absence of a serial correlation. The Jarque–Bera (JB) test also rejected the null hypothesis, rejecting the normality of residuals.

4.2 Granger causality, impulse response function and forecast error variance decomposition analyses

To understand the economic relationship between variables, we employed the Granger causality test. According to the analysis, the USCPI and DXY granger causes the Bitcoin price (10 and 5% significance levels, respectively). When analyzing the opposite direction, we find that only DXY granger causes Bitcoin at a 5% significance level. This outcome contradicts the previous results in Kim et al. (2019). Specifically, they reported a reciprocal Granger causality between VIX and the Bitcoin, while we found a mutual Granger causality between DXY and Bitcoin [5] (Table 9).

The IRF results with lag (n = 12) are summarized by n = 1, 3, 5 and 12 (left axis = impulse variable). In case of the impulse from the VIX, the Bitcoin price first moved negatively and gradually mitigated negative shocks over time. Concerning the impulse from the USCPI and US10Y, they provided a (+) positive impulse to the Bitcoin price. Concerning GOLD and DXY, these safe assets exerted a persistent (−) negative impact until n = 5 (roughly 5 months). In summary, Bitcoin has moved positively toward the inflation linked variables, like USCPI and USCPI, while it moved against the GOLD and DXY impulse shock (Table 10).

Concerning the outcomes of FEVD with lag (n = 12), they are summarized as follows. In case of VIX, GOLD and Bitcoin accounted for 19.3 and 2.2%, respectively, each at n = 12 (at the 12th month). For Bitcoin at n = 12, Bitcoin, DXY and GOLD accounted for 89.6, 3.7 and 2.1%, respectively. For DXY, when n increases from 1 to 12, the explanatory weights of GOLD decreased from 21 to 12.8%. However, the explanatory power of DXY, USCPI and GOLD increased from 52 to 55.9, 2.0–2.9, and 0.2–1.7%, respectively. The explanatory power of GOLD declined from 85.6 to 57.5%, when n moved from 1 to 12. However, US10Y to GOLD increased from 8.3 to 16.6%, USCPI to GOLD increased from 4.6 to 11.0% and Bitcoin to GOLD also increased from 1.0 to 2.7% (Table 11).

4.3 Time-varying cointegration test

Given the dynamic relationship between the Bitcoin price and the macroeconomic time series, it would be appropriate to consider the varying degree of their long-run relationship. This means that a time-varying approach should replace a constant coefficient equation from the general VECM. Hence, we develop the TVC VECM and test the feasibility of this expansion using the LR test. The cointegration rank for both the TIC and TVC VECM was set at r = 1 (see Table 12).

  1. Null hypothesis: The TIC VECM and TVC VECM share the same cointegration coefficient relationship across the Chebyshev polynomials m = 1, 2, 3 and 4

The tests yield the following results. First, when Chebyshev polynomials (m) = 1, the null hypothesis was not rejected at a 10% significance level (P value = 0.18). Second, when Chebyshev polynomial of m = 2 or m > 3, the null hypothesis was rejected (P value = 0.00152 for m = 2, P value = 0.00233 for m = 3, P value = 0.00019 for m = 4), with TVC revealed visually. When we applied the HQ information criteria, the optimal Chebyshev polynomials m were 2 (Bierens and Martins, 2010).

When we used the Chebyshev polynomial of m = 1, the null hypothesis was not rejected. As shown in Figure 3, the TVC parameters have shown like constant numbers across the time horizon. At this point, we believe that they are similar to the TIC VECM cointegration parameters (Beta chart with m = 1). Hence, the null hypothesis holds.

When we use the Chebyshev polynomial of m = 2, the null hypothesis was rejected, implying the significance of the TVC parameters. Figure 4 shows the smooth and gradual evolution of the long-run relationship coefficient of three economic variables, including US10Y.

While we have shown the coefficient of the long-run relationship using the Chebyshev polynomials of m (see Appendix 1 for m = 3 and 4), the following Figures 5–7 exhibit the ΠYt1=αβYt1 (cointegration relationship) under the VECM. Figure 5 (cointegration relation with nondeterministic term) shows a cointegration relationship under the TIC VECM, and Figures 6 and 7 show a cointegration relationship (m = 1,2) under the TVC VECM [6].

5. Conclusion

This study analyzes the Bitcoin price in relation to six macroeconomic variables studied over a 10-year time horizon, using the VECM. Most studies have reported the speculative characteristics of Bitcoin owing to its short history and association with the financial market. Hence, we use a period when Bitcoin was more exposed to the financial market interactions, that is, when it achieved and crossed a market capitalization of US$1 tn. We also construct a VAR/VECM model suitable for capturing the interrelationship among the sampled macroeconomic variables. We validate these relationships using the time series stationary test and an optimal lag selection. We also test for the long-run cointegration both under the TIC VECM (p = 1, r = 1) and TVC VECM. Given the dynamic nature of the macroeconomic variables, based on the results of the LR test, we use a TVC VECM for Bitcoin modeling.

This study has implications in demonstrating that the emerging digital assets, including Bitcoin, would prevail and consolidate. Given this, it is important to analyze the dynamic long-run relationship between digital assets and macroeconomic variables from a modeling perspective. This study has used the time-varying parameter approach for long-run cointegration from the VECM. However, the future studies may apply the time-varying approach to loading factor.

Figures

Bitcoin market capitalization (left) and price (right). Chart (Y-axis: log scale)

Figure 1

Bitcoin market capitalization (left) and price (right). Chart (Y-axis: log scale)

Underlying variables – historical chart (data period: August 2010 to February 2022; monthly data: 139 observations). Bitcoin price (BTCUSD) is denoted by log (BTCUSD)

Figure 2

Underlying variables – historical chart (data period: August 2010 to February 2022; monthly data: 139 observations). Bitcoin price (BTCUSD) is denoted by log (BTCUSD)

Chebyshev polynomials with m = 1

Figure 3

Chebyshev polynomials with m = 1

Chebyshev polynomials with m = 2

Figure 4

Chebyshev polynomials with m = 2

Cointegration relationship of Chebyshev polynomials with m = 0

Figure 5

Cointegration relationship of Chebyshev polynomials with m = 0

Cointegration relationship of Chebyshev polynomials with m = 1

Figure 6

Cointegration relationship of Chebyshev polynomials with m = 1

Cointegration relationship of Chebyshev polynomials with m = 2

Figure 7

Cointegration relationship of Chebyshev polynomials with m = 2

Underlying variables – correlation statistics (Data period: Aug-2010 to Feb-2022, monthly data, 139 observations)

Figure A1

Underlying variables – correlation statistics (Data period: Aug-2010 to Feb-2022, monthly data, 139 observations)

VECM (p = 1, r = 1) impulse response function (Response = BTCUSD)

Figure A2

VECM (p = 1, r = 1) impulse response function (Response = BTCUSD)

VECM (p = 1, r = 1) impulse response function (Impulse = BTCUSD)

Figure A3

VECM (p = 1, r = 1) impulse response function (Impulse = BTCUSD)

VECM (p = 1, r = 1) FEVD chart representation from the main paper

Figure A4

VECM (p = 1, r = 1) FEVD chart representation from the main paper

Cointegration relation (with Chebyshev polynomials of m = 3)

Figure A5

Cointegration relation (with Chebyshev polynomials of m = 3)

Cointegration relation (with Chebyshev polynomials of m = 4)

Figure A6

Cointegration relation (with Chebyshev polynomials of m = 4)

Cointegration beta (with Chebyshev polynomials of m = 3)

Figure A7

Cointegration beta (with Chebyshev polynomials of m = 3)

Cointegration beta – normalized on BTCUSD (with Chebyshev polynomials of m = 3)

Figure A8

Cointegration beta – normalized on BTCUSD (with Chebyshev polynomials of m = 3)

Cointegration beta (with Chebyshev polynomials of m = 4)

Figure A9

Cointegration beta (with Chebyshev polynomials of m = 4)

Cointegration beta – normalized on BTCUSD (with Chebyshev polynomials of m = 4)

Figure A10

Cointegration beta – normalized on BTCUSD (with Chebyshev polynomials of m = 4)

Summary statistics

VIXUS10YBTCUSDUSCPIGOLDDXY
# of Observation139139139139139139
Max9.51000.5500−2.7806218.31201060.000085.5999
Min53.54003.470011.0305283.71601964.9000122.8165
Mean18.30712.10236.3133243.54921438.6662105.3285
Median16.48002.16006.4841240.23601326.5000110.0848
Stdev6.90790.64483.287314.7723238.307210.9168
Skewness1.9509−0.2706−0.70360.56480.5271−0.3695
Kurtosis5.2782−0.2170−0.2192−0.2483−1.0416−1.4860

Results of augmented Dickey–Fuller and test regression with none condition

VIXUS10YBTCUSDUSCPIGOLDDXY
A. Level
Test regression, None
ADF(t-value)−0.09−0.490.742.530.391.00
P-value0.9287400.6250000.4631800.012767**0.6980000.322400
Lag (from AIC)72104211
Test regression, Trend
ADF(t-value)−4.462.08−4.440.53−2.82−1.64
P-value0.0000212***0.0399**0.000024***0.5976980.005723***0.1046
Lag (from AIC)229428
VIXUS10YBTCUSDUSCPIGOLDDXY
B. 1st Difference
Test regression, None
ADF(t-value)−6.69−6.53−2.830.60−7.64−2.83
P-value0.0000000013***0.0000000024***0.005598***0.5533000.00000000001***0.00571***
Lag (from AIC)62811210
Test regression, Trend
ADF(t-value)−6.65−6.47−3.10−4.44−8.03−3.34
P-value0.0000000016***0.0000000033***0.00253***0.000022***0.00000000000***0.00121***
Lag (from AIC)6293210

Note(s): *** significant at 1% level, ** significant at 5% level and * significant at 10% level

Results of the VAR order selection criteria

None Both (=Constant + Trend)
LagAlC(n)HQ(n)SC(n)FPE(n)LagAIC(n)HQ(n)SC(n)FPE(n)
15.89996.22106.6901365.159015.77236.20046.8258321.5622
25.76916.41127.3494321.242425.63846.38757.4820282.3844
35.96456.92778.3349393.476335.77896.84918.4126328.0954
46.27957.56389.4400547.085146.06307.45439.4869443.6370
56.61748.222710.5680786.014256.40248.114710.6164640.8920
66.73998.666311.4807922.337466.61858.652011.6227829.9974
76.78509.032412.31591,017.978176.54028.894612.3345814.7621
86.83099.399513.15201,147.420186.60519.280713.1896943.1980

Results of the Johansen test

None (lag = 2) Critical valueTrend (lag = 2) Critical value
Null hypothesisTest statistics1%5% Null hypothesisTest statistics1%5%
Tracer ≤ 231.7755.4348.28Tracer ≤ 237.6270.0562.99
r ≤ 152.0778.8770.6 r ≤ 171.1896.5887.31
r = 096.37104.290.39 r = 0115.77124.75114.9
Eigenr ≤ 215.8732.1432.14Eigenr ≤ 216.6736.6531.46
r ≤ 120.3138.7833.32 r ≤ 133.5642.3637.52
r = 044.344.5939.43 r = 044.5949.5143.97

VECM loading matrix (or weights) and statistical significance test

Trend*None
Coefficientt-valuePr(>|t|)Coefficientt-valuePr(>|t|)
VIX0.316944.132200.00010−0.36320−4.423700.00003
US10Y0.001010.326000.752000.000760.229200.82400
ECT1**BTCU SD0.011082.290800.027960.008941.699300.10160
USCPI−0.00804−0.841600.41600−0.00961−0.932800.36700
GOLD2.886302.802900.007433.224302.907700.00554
DXY−0.00893−0.347500.73600−0.01685−0.608100.55600

Note(s): * VECM with Trend adopted from the analysis ** ECT1: error correction term (α)

VECM cointegration relationship: the coefficient (β) of cointegrating equation estimated under Trend and none assumption

Trend*None
Coefficientt-valuePr(>|t|)Coefficientt-valuePr(>|t|)
VIX1.0000 1.0000
US10Y4.57762.03680.0420−3.6978−1.80670.0710
BTCUSD2.72702.42170.0150−1.4801−1.99760.0460
ßetaUSCPI0.36170.95940.33700.51222.34080.2190
GOLD−0.0443−5.02380.0000−0.0376−5.11660.0070
DXY0.67621.75100.0800−0.3418−1.97130.1730
Trend0.22810.91210.3620

Note(s): * VECM with Trend adopted from the analysis

VECM (p = 1, r = 1) representation

D.BTCUSDTrend
Coefficientt-valuePr(>|t|)
ECT10.01112.22300.0280
Const0.50232.59600.0l05
D.VIX(l1)−0.0081−1.31300.1915
D.US10Y(l1)0.05710.37500.7085
D.BTCUSD(l1)0.22792.68900.0081Adj R square0.1625
D.USCPI(l1)0.04691.25300.2126F-statistics4.3220
D.GOLD(l1)−0.0007−1.46200.1463Pr(>|t|)0.0001
D.DXY(l1)0.04221.97400.0505

Note(s): ECT1: error correction term (α); coefficient: short-term relationship;

adjusted R square: 0.1625; F-statistics: 4.322; p-value: 0.001 (significance of ΔBTCUSD equation)

Model diagnostic statistics using ARCH, PT and JB

VECM(r = 1)Trend (Y/N)Multivariate ARCH (lags = 5)P valuePortmanteau test (asymptotic, lags = 16)P valueJB normality testP value
p = 1Y2251.60.23990527.330.28860370.22.20E-16

Granger causality test based on F-statistics with lag = 2 (Wald test comparing the unrestricted with the restricted model)

CausedBITCUSD CausingBITCUSD
F-valuePr(>|t|) F-valuePr(>|t|)
VIX1.383100.25440 VIX0.344700.70910
US10Y1.213600.30040 US10Y1.110400.33250
CausingUSCPI2.389900.09559CausedUSCPI0.258400.77260
GOLD0.798800.45200 GOLD0.072600.93000
DXY3.393600.03654 DXY2.525900.08384

Impulse response function

Response VariableResponse Variable
n = 1VIXUS10YBTCUSDUSCPIGOLDDXYn = 3VIXUS10YBTCUSDUSCPIGOLDDXY
VIX4.945−0.0520.059−0.0644.8010.818VIX2.586−0.0710.087−0.40816.1310.883
US10Y0.1920.0200.101−19.110−0.097US10Y−0.4860.2200.0290.310−25.838−0.077
BTCUSD0.3060.0666.6210.081BTCUSD0.2020.0210.3850.1748.9160.237
USCPI0.60014.244−0.234USCPI0.113−0.0130.0421.14820.474−0.275
GOLD61.415−0.759GOLD1.039−0.0230.0480.11848.074−0.612
DXY1.195DXY0.1330.0040.072−0.111−11.6561.276
IRF outcome (vertical: impulse variables), n = number of time lapse
Response VariableResponse Variable
n = 5VIXUS10YBTCUSDUSCPIGOLDDXYn = 12VIXUS10YBTCUSDUSCPIGOLDDXY
VIX1.809−0.0580.054−0.50519.5290.891VIX1.433−0.0500.030−0.53721.3590.882
US10Y−0.3830.2200.0370.384−26.297−0.088US10Y−0.2930.2180.0360.416−26.499−0.091
BTCUSD0.3970.0210.3780.20110.930.239BTCUSD0.4930.0190.3700.21411.740.234
USCPI0.226−0.0170.0551.31321.969−0.298USCPI0.337−0.0190.0571.38222.184−0.306
GOLD1.240−0.0270.0600.14047.007−0.601GOLD1.304−0.0290.0650.14646.743−0.599
DXY0.1730.0040.080−0.143−12.6471.291DXY0.1610.0040.082−0.157−12.7481.293

Forecast error variance decomposition by dependent variables

Dependent variableVIX Dependent variableUS10Y Dependent variableBTCUSD
n = VIXUS10YBTCUSDUSCPIGOLDDXYVIXUS10YBTCUSDUSCPIGOLDDXYVIXUS10YBTCUSDUSCPIGOLDDXY
1100.0%0.0%0.0%0.0%0.0%0.0%6.8%93.2%0.0%0.0%0.0%0.0%3.5%0.4%96.0%0.0%0.0%0.0%
296.6%0.5%0.0%0.0%2.8%0.1%9.7%89.6%0.3%0.1%0.3%0.0%5.7%0.2%92.0%0.3%0.5%1.4%
394.1%1.0%0.1%0.0%4.7%0.1%9.5%89.3%0.5%0.2%0.6%0.0%5.3%0.3%90.9%0.6%0.8%2.1%
491.4%1.2%0.3%0.1%6.9%0.2%9.0%89.5%0.6%0.2%0.7%0.0%4.6%0.4%90.4%0.9%1.1%2.5%
588.8%1.3%0.5%0.2%8.9%0.2%8.4%89.7%0.6%0.3%0.9%0.0%4.0%0.5%90.1%1.1%1.4%2.8%
686.4%1.4%0.8%0.3%10.9%0.2%8.0%90.0%0.7%0.3%1.0%0.0%3.5%0.6%90.0%1.3%1.6%3.1%
784.1%1.5%1.1%0.4%12.7%0.3%7.6%90.3%0.7%0.4%1.1%0.0%3.1%0.6%89.9%1.4%1.7%3.3%
882.0%1.5%1.4%0.5%14.3%0.3%7.3%90.5%0.7%0.4%1.1%0.0%2.8%0.7%89.8%1.5%1.8%3.4%
980.2%1.6%1.6%0.6%15.7%0.3%7.0%90.6%0.7%0.5%1.2%0.0%2.6%0.7%89.7%1.5%1.9%3.5%
1078.5%1.6%1.8%0.7%17.1%0.3%6.8%90.8%0.7%0.5%1.2%0.0%2.4%0.7%89.7%1.6%2.0%3.6%
1177.0%1.6%2.0%0.8%18.2%0.3%6.6%90.9%0.7%0.5%1.3%0.0%2.2%0.7%89.7%1.6%2.1%3.7%
1275.6%1.6%2.2%0.9%19.3%0.3%6.5%91.0%0.7%0.5%1.3%0.0%2.1%0.7%89.6%1.7%2.1%3.7%
Dependent variableUSCPI Dependent variableGOLD Dependent variableDXY
n = VIXUS10YBTCUSDUSCPIGOLDDXYVIXUS10YBTCUSDUSCPIGOLDDXYVIXUS10YBTCUSDUSCPIGOLDDXY
11.1%2.7%1.2%95.1%0.0%0.0%0.5%8.3%1.0%4.6%85.6%0.0%24.4%0.3%0.2%2.0%21.0%52.0%
25.8%4.4%1.7%87.1%0.7%0.3%1.0%14.3%1.0%6.7%74.9%2.2%26.1%0.2%1.0%2.0%16.2%54.4%
38.1%5.2%1.8%83.6%0.8%0.5%2.9%15.3%1.4%8.0%69.8%2.6%26.3%0.2%1.3%2.2%15.1%54.9%
49.3%5.7%1.8%81.8%0.8%0.7%4.1%15.8%1.7%9.0%66.4%3.0%26.5%0.2%1.5%2.4%14.3%55.1%
510.0%6.0%1.8%80.7%0.8%0.7%5.2%16.1%1.9%9.6%64.0%3.2%26.5%0.2%1.6%2.5%13.9%55.3%
610.4%6.2%1.8%80.0%0.8%0.8%6.1%16.3%2.1%10.0%62.2%3.3%26.5%0.2%1.6%2.6%13.6%55.4%
710.7%6.3%1.8%79.5%0.9%0.8%6.7%16.4%2.3%10.3%60.9%3.4%26.5%0.2%1.7%2.7%13.4%55.6%
810.8%6.4%1.8%79.2%0.9%0.9%7.3%16.4%2.4%10.5%59.9%3.5%26.5%0.2%1.7%2.7%13.2%55.7%
911.0%6.5%1.8%78.9%0.9%0.9%7.7%16.5%2.5%10.7%59.1%3.5%26.5%0.2%1.7%2.8%13.1%55.7%
1011.1%6.6%1.8%78.7%0.9%0.9%8.0%16.5%2.6%10.8%58.5%3.6%26.4%0.2%1.7%2.8%13.0%55.8%
1111.1%6.6%1.8%78.6%0.9%0.9%8.3%16.6%2.7%10.9%57.9%3.6%26.4%0.2%1.7%2.8%12.9%55.9%
1211.2%6.7%1.8%78.5%0.9%0.9%8.5%16.6%2.7%11.0%57.5%3.6%26.4%0.3%1.7%2.9%12.8%55.9%

Summary results on time-varying cointegration

ModelVECM(p = 1, r = 1)
m = 1m = 2m = 3m = 4
LR TVC statistics8.8739531.7388739.6544956.63646
P-value0.180790.00152***0.00233***0.00019***
Log likelihood−1511.23162−1499.79916−1495.84136−1487.35037
AIC22.7041122.5956122.8298022.76424
BIC23.6419223.6186724.2791324.29883
HQC23.0852123.0113523.4187723.38786

Note(s): *** significant at 1% level, ** significant at 5% level and * significant at 10% level

Notes

2.

This is the case when Chebyshev polynomials m = 0

3.

For the maximum likelihood estimates and LR TVC, refer to Bierens and Martins (2010)

4.

As the VAR model was estimated with lag = 2 (VAR (p = 2)), the VECM model is presented as lag = 1 (VECM (p = 1)).

5.

Kim et al. (2019) used the daily data from July 2010 to April 2018, while used monthly data from August 2010 to Feb 2022.

6.

Cointegration Relation with m = 3,4 were reported at Appendix 3 and normalization in charts was based on the coefficient of BTCUSD.

Appendix 1

Figure A1-A10

Appendix 2

Orthonormal Chebyshev Polynomials Pi,T(t) can be presented as below,

(A)Pi,T(t)=2cos(iπ(t0.5)T), P0,T(t)=1,  t=1,2,,T,  i=1,2,3,..

Using and substituting βt  βt(m)=i=0mξiPi,T(t) into TVC VECM equation (6) from main article results in

(B)ΔYt=α(i=0mξiPi,T(t))Yt1+j=1p1ΓjΔYtj+εt
With respect to arbitrary kxr matrix ξi (B) equation can be arranged as follows,
(C)ΔYt=αξYt1(m)+ΓXt+εt
where ξ=(ξ0,ξ1,,ξm) is an r(m+1)k matrix of rank r 

Let  Yt1(m)= (Yt1,P1,T(t)Yt1,P2,T(t)Yt1,,Pm,T(t)Yt1) and Xt= (ΔYt1,,ΔYtp+1) then,

Coefficient β of TIC (time invariant cointegration) for null hypothesis is derived from  ξYt1(m)= βYt1(0) with  Yt1(0)= Yt1  which means ξ=(β, Or,k,m)

Appendix 3

References

Bierens, H.J. and Martins, L.F. (2010), “Time-varying cointegration. econometric theory 143–1490”, Econometric Theory, Vol. 26 No. 5, pp. 1453-1490, doi: 10.1017/S0266466609990648.

Bouoiyour, J. and Selmi, R. (2015), “What does bitcoin look like?”, Analysis of Economics and Finance, Vol. 16 No. 2, pp. 449-492.

Engle, R.F. and Granger, C.W. (1987), “Co-integration and error correction: representation, estimation, and testing”, Econometrica: Journal of the Econometric Society, Vol. 55 No. 2m, pp. 21-276, doi: 10.2307/1913236.

Kim, J.M. Kang, N.Y. and Park, Y.J. (2019), “A study on cross – effects of prices of bitcoin, traditional assets, and traditional currencies”, Korean Business Education Review, Vol. 34 No. 5, pp. 151-169, doi: 10.23839/kabe.2019.34.5.151.

Lee, G.S., Joe, Y.M. and Jeong, J.H. (2019a), “An investigation of dynamic price movements of the cryptocurrency coin in Korea”, Asian Review of Financial Research, Vol. 32 No. 3, pp. 383-400.

Lee, K.K., Cho, S.J., Min, G.S. and Yang, W.C. (2019b), “The determinant of bitcoin prices in Korea”, Korean Journal of Financial Studies, Vol. 48 No. 4, pp. 393-415, doi: 10.26845/KJFS.2019.08.48.4.393.

Pfaff, B. (2008), “VAR, SVAR and SVEC models: implementation within r package vars”, Journal of Statistical Software, Vol. 27 No. 4, doi: 10.18637/jss.v027.i04.

Sims, C.A. (1980), “Macroeconomics and reality”, Econometrica, Vol. 48 No. 1, pp. 1-48, doi: 10.2307/1912017.

Son, J.H. and Kim, J.Y. (2019), “A study in bitcoin volatility through economic factors”, The Journal of Society for E-Business Studies, Vol. 24 No. 4, pp. 109-118, doi: 10.7838/jsebs.2019.24.4.109.

Syafiqah, I. and Mohamad, Y. (2021), “An analysis on cryptocurrencies and macroeconomic variables using VECM, ASEAN”, Journal of Management and Business Studies, Vol. 3 No. 1, pp. 08-15, doi: 10.26666/rmp.ajmbs.2021.1.2.

Thaker, H.M.T. and Mand, A.A. (2021), “Bitcoin and stock markets: a revisit of relationship”, Journal of Derivatives and Quantitative Studies, Vol. 29 No. 3.

Yang, H., Yang, G.H. and Oxley, L. (2020), “What role do futures markets play in bitcoin pricing? Causality, cointegration and price discovery from a time-varying perspective?”, International Review of Financial Analysis, Vol. 72, doi: 10.1016/j.irfa.2020.101569.

Zhu, Y., Dickinson, D. and Li, J. (2017), “Analysis on the influence factors of bitcoin's price based on VEC model”, Financial Innovation, Vol. 3 No. 3, pp. 1-13, doi: 10.1186/s40854-017-0054-0.

Further reading

Badev, A. and Chen, M. (2014), “Bitcoin: technical background and data analysis”, Working Papers, No. 2014-104, Finance and Economics Discussion Series (FEDS), Divisions of Research & Statistics and Monetary Affairs Federal Reserve Board, Washington, DC, pp. 1-39.

Barai, M. and Kundu, S. (2019), “The role of the federal reserve in the U.S housing crisis: a VAR analysis with endogenous structural breaks”, Journal of Risk Financial Management, Vol. 12 No. 3, doi: 10.3390/jrfm12030125.

Berentsen, A. and Schar, F. (2018), “A short introduction to the world of cryptocurrencies. Federal Reserve Bank of St. Louis Review”, First Quarter, Vol. 100 No. 1, doi: 10.20955/r.2018.1-16.

Casas, I. and Fernandez-Casal, R. (2019), “tvReg: time-varying coefficient linear regression for single and multi-equations in R”, available at: https://ssrn.com/abstract=3363526.

Dyhrberg, A.H. (2015), “Bitcoin, gold and the dollar – a GARCH volatility analysis”, Finance Research Letters, Vol. 16, pp. 85-92, doi: 10.1016/j.frl.2015.10.008.

Choi, S.Y. and Shin, J.S. (2019), “Analysis of cryptocurrency volatility”, The Korean Journal of Financial Management, Vol. 36 No. 2, pp. 65-82.

D'Agostino, A., Gambeti, L. and Giannone, D. (2011), “Macroeconomic forecasting and structural change”, Journal of Applied Econometrics, Vol. 28 No. 1, pp. 82-101, doi: 10.1002/jae.1257.

Glaser, F., Zimmermann, K., Haferkorn, M., Weber, M. and Siering, M. (2014), “Bitcoin – asset or currency? revealing user's hidden intentions”, Twenty Second European Conference on Information Systems, Tel Aviv 2014, Electronic copy, available at: http://ssrn.com/abstract=2425247.

Gottschalk, J. (2001), “An introduction into the SVAR methodology: identification, interpretation and limitations of SVAR models”, Kiel Working Paper No. 1072, Kiel Institute for the World Economy (IfW Kiel), available at: http://hdl.handle.net/10419/17887.

Haslbeck, J., Bringmann, L. and Waldorp, L. (2020), “A tutorial on estimating time-varying vector autoregressive models”, Multivariate Behavioral Research, Vol. 56 No. 1, pp. 120-149, doi: 10.1080/00273171.2020.1743630.

Jang, S.I. and Kim, J.Y. (2017), “A study on the asset characterization of bitcoin”, The Journal of Society for E-Business Studies, Vol. 22 No. 4, pp. 117-128, doi: 10.7838/jsebs.2017.22.4.117.

Johansen, S. (1988), “Statistical analysis of cointegration vectors”, Journal of Economic Dynamics and Control, Vol. 12 Nos 2-3, pp. 231-254, doi: 10.1016/0165-1889(88)90041-3.

Kjaerland, F., Khazal, A., Krogstad, E.A., Nordstrom, F.B.G. and Oust, A. (2018), “An analysis of bitcoin's price dynamics”, Journal of Risk and Financial Management, Vol. 11, doi: 10.3390/jrfm11040063.

Kuschnig, N. and Vashold, L. (2020), “BVAR: Bayesian vector autoregressions with hierarchical prior selection in R”, Working paper, 296, WU Vienna University of Economics and Business, available at: https://epub.wu.ac.at/id/eprint/7216.

Lange, A., Dalheimer, B., Herwartz, H. and Maxand, S. (2021), “An r package for data-driven identification in multivariate time series analysis”, Journal of Statistical Software, Vol. 97 No. 5, doi: 10.18637/jss.v097.i05.

Lee, J.S., Kim, K.W. and Park, D.H. (2018), “Empirical analysis on bitcoin price change by consumer, industry and macro-economy variables”, Journal of Intelligence and Information Systems, Vol. 24 No. 2, pp.195-220, doi: 10.13088/jiis.2018.24.2.195.

Martins, L.F. (2018), “Bootstrap tests for time varying cointegration”, Econometric Reviews, Vol. 37 No. 5, pp. 466-483, doi: 10.1080/07474938.2015.1092830.

Meynkhard, A. (2019), “Fair market value of bitcoin: halving effect”, Investment Management and Financial Innovations, Vol. 16 No. 4, doi: 10.21511/imfi.16(4).2019.07.

Nam, Y.S. and Kang, S.H. (2022), “Effect of quantitative growth of the cryptocurrency market on herding phenomenon”, The Korean Journal of Financial Management, Vol. 39 No. 1, pp. 63-96, doi: 10.22510/kjofm.2022.39.1.003.

Thaker, H.M.T. and Mand, A.A. (2021), “Bitcoin and stock markets: a revisit of relationship”, Journal of Derivatives and Quantitative Studies, Vol. 29 No. 3, pp. 234-256, doi: 10.1108/JDQS-07-2020-0016.

Yoo, J.H., Kang, J.Y. and Park, S.U. (2018), “Measuring return and volatility spillovers across major virtual currency market”, The Journal of Information Systems, Vol. 27 No. 3, pp. 43-62, doi: 10.5859/KAIS.2018.27.3.43.

Corresponding author

Joon Hee Rhee can be contacted at: joonrh@ssu.ac.kr

Related articles