Next Article in Journal
Estimating Bargaining Power in Real Estate Pricing Models: Conceptual and Empirical Issues
Next Article in Special Issue
Robust Inference in the Capital Asset Pricing Model Using the Multivariate t-distribution
Previous Article in Journal
Corporate Governance and Agency Cost: Empirical Evidence from Vietnam
Previous Article in Special Issue
Forecasting the Term Structure of Interest Rates with Dynamic Constrained Smoothing B-Splines
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Hypothesis Test Method for Detecting Multifractal Scaling, Applied to Bitcoin Prices

School of Finance & Applied Statistics, Australian National University, Canberra ACT 0200, Australia
*
Author to whom correspondence should be addressed.
J. Risk Financial Manag. 2020, 13(5), 104; https://doi.org/10.3390/jrfm13050104
Submission received: 2 March 2020 / Revised: 13 May 2020 / Accepted: 14 May 2020 / Published: 20 May 2020
(This article belongs to the Special Issue Financial Statistics and Data Analytics)

Abstract

:
Multifractal processes reproduce some of the stylised features observed in financial time series, namely heavy tails found in asset returns distributions, and long-memory found in volatility. Multifractal scaling cannot be assumed, it should be established; however, this is not a straightforward task, particularly in the presence of heavy tails. We develop an empirical hypothesis test to identify whether a time series is likely to exhibit multifractal scaling in the presence of heavy tails. The test is constructed by comparing estimated scaling functions of financial time series to simulated scaling functions of both an iid Student t-distributed process and a Brownian Motion in Multifractal Time (BMMT), a multifractal processes constructed in Mandelbrot et al. (1997). Concavity measures of the respective scaling functions are estimated, and it is observed that the concavity measures form different distributions which allow us to construct a hypothesis test. We apply this method to test for multifractal scaling across several financial time series including Bitcoin. We observe that multifractal scaling cannot be ruled out for Bitcoin or the Nasdaq Composite Index, both technology driven assets.

1. Introduction

The cryptocurrency market is an emerging market comprised of thousands of digital assets including Bitcoin. A cryptocurrency is a digital asset that uses cryptography and decentralised governance to secure a ledger of transactions. Cryptocurrencies can also be used as a medium of exchange; for background, refer to the review article by Milutinović (2018). Bitcoin, introduced in Nakamoto (2009), is an example of the world’s first cryptocurrency; its core innovation established digital scarcity, regulated and audited via a novel decentralised governance mechanism. Financial analysts, chartists and the news media often report on Bitcoin’s fractal features (Chambers 2019); however, academic research has shown mixed evidence for multifractal or monofractal scaling with Mensi et al. (2019); Lahmiri and Bekiros (2018); Stavroyiannis et al. (2019) presenting a case for multifractal scaling, while Nadarajah and Chu (2017); Bariviera (2017); Bariviera et al. (2017); Zhang et al. (2018) presenting a case for monofractal scaling. The inconclusive evidence of fractal scaling could be a result of the different detection methods along with varied assumptions. Salat et al. (2017) emphasise that implementation of existing multifractal detection methods can be a perilous undertaking. Moreover, Sly (2006) and Grahovac and Leonenko (2014) demonstrate that attempts to use the scaling function to detect multifractality may mistakenly give rise to false positive results, due to heavy-tailed effects. This calls into question the reliability of the above-mentioned results. Therefore, a detailed examination of scaling functions in the presence of heavy tails is required to bridge the gap between the theory and empirical realities. As an application of the new method we present here, we seek to understand Bitcoin’s scaling properties to assist practitioners with model selection for the purpose of risk management and volatility forecasting. Given the novel nature of digital assets, it is also of interest to observe similarities in scaling behaviour between Bitcoin and other asset classes. In this paper, we address the central research question: Do Bitcoin prices exhibit multifractal scaling? To answer this question, we investigate the effect of heavy-tails on the scaling function and study scenarios in which false positive1 detection occurs. If multifractal scaling is plausible, it could indicate the need to incorporate more complex scaling behaviour into pricing models. In the case of Bitcoin, capturing complex scaling behaviour could result in better volatility forecasting. This paper develops an empirical method that tests for multifractal scaling in the presence of heavy tails. Its application is not limited to financial time series, however. This paper is structured as follows. In the rest of this section, we introduce multifractal stochastic processes and their simulation methods. We then explain how the multifractal scaling function is distorted by the presence of heavy tails. In Section 2, we propose a new empirical hypothesis test and construct the look-up table to form rejection and acceptance regions. Then, in Section 3, we apply our hypothesis test to detect multifractal scaling in Bitcoin prices and compare the scaling behaviour of Bitcoin to other financial assets, including the S&P500 index, the Nasdaq Composite Index, the USD/JPY exchange rate, and Gold Futures. Our method detects more complex scaling behaviour in emerging technology asset classes such as Bitcoin and the Nasdaq Composite Index.

1.1. Monofractal vs. Multifractal Processes

The multifractal system was introduced as a generalisation of the fractal system to describe more complicated dynamics in time series. Multifractals are common in nature and have enjoyed great application in finance and science, e.g., modelling the turbulence in fluid dynamics (Sreenivasan 1991). In the finance field, multifractal processes are able to capture many stylised characteristics of high-volatility financial assets.
Multifractal processes are defined based on the scaling property of their moments, when they are finite, in Mandelbrot et al. (1997).
Definition 1.
(Multifractal Stochastic Process) A stochastic process { X ( t ) , t 0 } is multifractal if it has stationary increments and there exist functions c ( q ) > 0 and τ ( q ) > 0 and positive constants Q and T such that
E | X ( t ) | q = c ( q ) t τ ( q ) , for all q [ 0 , Q ] and t [ 0 , T ] .
The function τ ( q ) is called the scaling function.
The multifractal processes are said to be multiscaling when the scaling function is nonlinear. For a multifractal process, the scaling function τ ( q ) is always concave when q > 0 on a bounded interval as is shown in Sly (2006).
In contrast, the multifractal process is called uniscaling or a monofractal process when the scaling function is linear, τ ( q ) = H q . Monofractal processes enjoy a self-similarity.2
Theorem 1.
(Lamperti (1962)). If { X ( t ) , t 0 } is self-similar and stochastically continuous at t = 0, then there exists a unique H > 0 such that for all a > 0 ,
{ X ( a t ) } = d { a H X ( t ) }
where equality is of finite dimensional distributions.
Here, H 0 is a constant known as the Hurst parameter. H measures in some way the persistence of a process. Given all moments of the process are finite, the process is long-range dependent if H ( 1 2 , 1 ) , while anti-persistence occurs if H ( 0 , 1 2 ) (Embrechts and Maejima 2000). Typical examples of monofractal processes are the fractional Brownian motions and the stable Lévy processes. The fractional Brownian motion B H ( t ) introduced in Mandelbrot et al. (1997) is a unique class of self-similar Gaussian process with stationary increments possessing a dependence structure described by the following covariance equation. When H = 1 2 , the fractional Brownian motion becomes Brownian motion with independent increments.

1.2. Examination of Multifractality

Copious methods have been proposed to examine the multifractality in a certain process { X ( t ) , 0 t T } , including the multifractal detrended fluctuation analysis method (MF-DFA) by Kantelhardt et al. (2002) and the wavelet transform modulus maxima method (WTMM). In this paper, we employ the standard partition method introduced by Mandelbrot et al. (1997). This method attempts to detect multifractality based on concavity in the scaling function. From Definition 1, one can easily derive
log E | X ( t ) | q = log c ( q ) + τ ( q ) log t , for all q [ 0 , Q ] and t [ 0 , T ] ,
where the q-th moment of the process X ( t ) , E | X ( t ) | q , can be estimated by the sample statistic
S q ( T , Δ t ) = 1 N i = 1 N | X ( ( i 1 ) Δ t ) X ( i Δ t ) | q ,
with N = T / Δ t .
The scaling function can be derived by fitting a linear regression between log S q ( T , Δ t ) and log t with various values of q [ 0 , Q ] . For a fixed value of q, we obtain an estimate τ ^ ( q ) , the slope of our linear regression. A plot of τ ^ ( q ) against q then provides visual display of the scaling function. The function of ln S q ( T , Δ t ) for various values of ln Δ t is called the partition function. If the scaling function shows concavity, we have evidence that multifractality is present in the process of interest.

1.3. Simulation of Multifractal Processes

The Brownian motion in multifractal time (BMMT) B H ( θ ( t ) ) introduced in Mandelbrot et al. (1997) is employed to simulate a multifractal process. The BMMT model shows appealing features that coincide with some stylised features of financial time series. It displays heavy-tails while not necessarily implying infinite variance. It also implies long-term dependence in absolute values of returns while the price increments themselves can remain uncorrelated. By construction, the BMMT is defined as the subordinate model of a fractional Brownian motion B H ( t ) with a multifractal process θ ( t ) .
Definition 2.
(Brownian Motion in Multifractal Time). Brownian motion in multifractal time is defined as
X ( t ) = B H ( θ ( t ) ) , t [ 0 , T ] ,
where θ ( t ) is a positive multifractal stochastic process and B H is an independent fractional Brownian motion with Hurst parameter H.
θ ( t ) can also be seen as a trading time where the index t denotes clock time. We assume the activity time θ ( t ) to be the cumulative density function of a random multifractal measure μ defined on [ 0 , T ] . That is, θ ( t ) is a multifractal process with continuous, non-decreasing paths, and stationary increments.
The multiplicative multifractal measure is implemented for simulation. Without loss of generality, we assume that the time series of interest is defined on a compact interval [ 0 , 1 ] . A multiplicative measure is constructed as follows:
  • In Stage 1, divide the time interval [ 0 , 1 ] into b > 1 non-overlapping subintervals with equal length 1 / b . Assign multipliers M 1 , β ( 0 β b 1 ) to each subinterval, where the M 1 , β are random variables with distributions that are not necessarily discrete. For computational convenience, we assume the M 1 , β to be identically distributed with a common distribution M.
  • In Stage 2, each of the b intervals is further divided into b subintervals of length 1 / b 2 . Again, we assign multipliers M 2 , β ( 0 β b 1 ) to each subinterval. The M 2 , β are assumed to be identically distributed with distribution M. Thus, after the second stage, the mass on an interval, for example [ 0 , 1 / b 2 ] , will be μ 2 [ 0 , 1 / b 2 ] = m 0 m 1 if M 1 , 0 = m 0 and M 2 , 0 = m 1 with probability
    P ( μ 2 [ 0 , 1 / b 2 ] = m 0 m 1 ) = P ( M = m 0 ) P ( M = m 1 ) ,
    since the multipliers at different stages are independent. The measure μ 2 represents the multiplicative measure at Stage 2.
  • Repetition of this scheme generates a sequence of measures ( μ k ) k N which converges to our desired multiplicative measure μ as k .
Remark 1.
To preserve the mass at each stage, some restrictions are required on the values of the M l , β , 1 l k . If we strictly assume that β = 0 b 1 M l , β = 1 at each stage, the resulting measure is called microcanonical or micro-conservative. If we loosen the assumption so that the mass at each stage is only conserved “on average", that is, E ( β = 0 b 1 M l , β ) = 1 , 1 l k , the resulting measure is called canonical.
In this paper, we only consider canonical measures, as they impose less restriction on the distribution of M. Let
t = 0 . η 1 η k = i = 1 k η i b i
be a b-adic number and set Δ t = b k . The mass on the b-adic cell [ t , t + Δ t ] for a canonical measure at Stage k is
μ ( Δ t ) = μ [ t , t + Δ t ] = Ω ( η 1 , , η k ) M ( η 1 ) M ( η 1 , η 2 ) M ( η 1 , , η k ) ,
where the random variable Ω represents the total mass. The high-frequency component Ω ( η 1 , , η k ) is assumed to have the same distribution as Ω . It captures changes in total mass of the interval caused by stages beyond k. As μ here is a canonical measure, we assume the multipliers M l , β ( 0 β b 1 ) satisfy E ( M l , β ) = 1 or equivalently E M = 1 / b .
Since multipliers at different stages of subdivision are independent, E [ μ ( Δ t ) q ] = E ( Ω q ) [ E ( M q ) ] k , for all q . Defining τ ( q ) = log b E ( M q ) and recalling Δ t = b k , we have
E [ μ ( Δ t ) q ] = E ( Ω q ) ( Δ t ) τ ( q ) .
Consequently, the constructed multiplicative measure μ is multifractal according to Definition 1. Because Equation (5) only holds for Δ t = b k and t = 0 . η 1 η k = i = 1 k η i b i being a b-adic number, these multiplicative measures are said to be grid-bound multifractal measures.
The relation between the distribution of M and the scaling function can be inferred from the following theorem.
Theorem 2.
(Calvet et al. 1997). Define p ( α ) to be the continuous density of V = log b M , and p k ( α ) to be the density of the k-th convolution product of p. The scaling function of a multiplicative measure satisfies
τ ( q ) = inf α [ α q lim k 1 k log b [ k p k ( k α ) ] .
In this paper, we have estimated sample scaling functions using the BMMT with a log-normal multiplicative measure. We have chosen this measure because its multifractal spectrum is a better fit for financial time series3. When V follows a normal distribution, N ( λ , σ 2 ) where λ > 1 , then M = b V follows a log-normal distribution that is, log b M = d N ( λ , σ 2 ) . We then have
E M = E ( e V log b ) = exp λ log b + ( log b ) 2 σ 2 2 ,
As M ( 0 , ) , we need the constraint E ( M ) = 1 / b to make M a canonical measure. This constraint gives a relationship between λ and σ 2 , namely log b = 2 ( λ 1 ) σ 2 .
Based on Theorem 2, the scaling function for the log-normal multiplicative cascade is derived as follows:
τ ( q ) = inf α α q lim k 1 k log b [ k p α ( k α ) ] = λ q σ 2 log b 2 q 2 .
This quadratic form reaches its maximum at q = λ σ 2 log b .
The cumulative distribution of the simulated multiplicative measures μ gives a multifractal multiplicative cascade, which is taken to be the trading time θ ( t ) = 0 t μ ( [ 0 , s ] ) d s in a BMMT model B H ( θ ( t ) ) .

1.4. Multifractality and Heavy Tails

While a concave scaling function is necessary for multifractality, it is not sufficient. Several papers have explored other potential drivers of concavity in the scaling function. One is heavy-tailedness (Sly 2006; Heyde 2009 and Grahovac and Leonenko 2014). Sly (2006) has shown that, after removing the extreme values of the S&P500 price process, the scaling function approaches linearity. He proceeded further to study the asymptotic behaviour of the estimated scaling functions and proved the following theorem.
Theorem 3.
(Sly 2006). Let X be a random variable with E ( X ) = 0 such that the distribution function of | X | has a regularly varying tail of order α where α > 2 ; that is,
P ( | X | > x ) = x α L ( x ) ,
where L ( x ) is slowly varying. Then, for an i.i.d. sequence with distribution X and q > α , for each s ( 0 , 1 ) , with S q as defined in Equation (4) and the time increment Δ t = n s ,
ln S q ( n , n s ) ln n p max s + q α 1 , s q 2
as n , where p stands for convergence in probability.
Grahovac and Leonenko (2014) generalised this result to a type E set of stochastic processes and extended Sly’s (2006) work by investigating the relationship between the asymptotic behaviour of scaling functions and tail indices.
Definition 3
(Type E stochastic processes). A stochastic process { X ( t ) , t 0 } is said to be of type E if Y ( t ) = X ( t ) X ( t 1 ) , t N , is a strictly stationary sequence having heavy-tailed marginal distribution with index α, satisfying the strong mixing property with an exponentially decaying rate and such that E Y ( t ) = 0 when α > 1 .
Theorem 4
(Grahovac and Leonenko 2014). Suppose { X ( t ) , 0 t T } is of type E and suppose Δ t i is of the form T i N for i = 1 , , N . Then, for every q > 0 and s ( 0 , 1 )
lim N plim T τ ^ N , T ( q ) = τ ( q ) ,
where plim stands for limit in probability and
τ ( q ) = q α i f 0 < q α a n d α 2 1 i f q > α a n d α 2 q 2 i f 0 < q α a n d α > 2 q 2 + 2 ( α q ) 2 ( 2 α + 4 q 3 α q ) α 3 ( q 2 ) 2 i f q > α a n d α > 2 .
Besides the heavy-tailed effects, other factors are discussed in a number of papers. Matia et al. (2003) concludes that price fluctuations can contribute to concavity of scaling function. Bouchaud et al. (2000) states that another factor is the long range nature of the volatility correlations. Nevertheless, neither of these scenarios excludes the effect of heavy-tailedness. Therefore, we regard the heavy-tailedness as a main contributor to the possible misdetection of multifractality in this paper.
Figure 1 displays the situation considered in Theorem 4. To illustrate the heavy-tailed effect, the estimated scaling functions for both S&P500 Index and 20 Student t 4 -distributed processes before and after truncating heavy tails are shown in Figure 2.4

2. Methodology

Our hypothesis test aims to identify multifractal scaling in the presence of heavy tails. It is constructed based on the different distributions of concavity measures for various tail indices:
H 0 : The process of interest is not multifractal . vs. H A : The process of interest is multifractal .
In this section, we first introduce the measure of concavity and tail index implemented in this paper. Then, we examine the distributions of concavity measures with regard to various tail indices for different processes. The following four instances are chosen: monofractal process (Brownian motion and fractional Brownian motion); Student t-distributed processes with degrees of freedom between 1 and 20; and a multifractal process simulated using BMMT with a log-Normal multiplicative measure.

2.1. Measure of Concavity

The global and localised simplex statistics proposed in Abrevaya and Jiang (2005) are employed. These loosen the assumption of normality and homoscedasticity of the error disturbances, which is required for most other concavity measures.
Definition 4
(Global Simplex Statistics). Assume that
1. 
there is an iid sample { ( x i , ϵ i ) } i = 1 n drawn from the joint distribution of the random variables ( x , ϵ ) , where ϵ is symmetrically distributed about 0 (conditional on x), so that E ( ϵ i | x i ) = 0 , i = 1 , , n ;
2. 
the observed sample is { ( y i , x i ) } i = 1 n , where y i is generated by y i = f ( x i ) + ϵ i , i = 1 , , n , and the functional form of f is left unspecified.
Then, for a sample of size n, the global simplex statistic is defined as
U n = n 3 1 1 t 1 t 2 t 3 n sign ( a 1 y [ t 1 ] + a 2 y [ t 2 ] y [ t 3 ] ) ,
where t 1 , t 2 , t 3 { 1 , , n } ,
sign ( q ) = 1 i f q > 0 , 1 i f q < 0 , 0 i f q = 0 .
and a 1 , a 2 are non-negative numbers with a 1 + a 2 = 1 .
Abrevaya and Jiang (2005) proves the asymptotic normality for the global simplex statistic. An interval is considered to be a concave interval if its standardised global simplex statistics U n ˜ is less than Z α / 2 , where Z α / 2 is the ( α / 2 ) -th percentile of the standard normal distribution, vice versa for a convex interval.
For the localised simplex statistic, calculation starts with dividing the time interval [ 0 , T ] into G = T / 2 h subintervals, where 2 h is the window width. Let x 1 , x 2 , , x G denote evaluation points. In this case, they are assumed to be the mid-points of each sub-interval. Then, the sub-population for which the x values fall in a local window is defined to be
V h ( x ) = { ( y , x ) : x h < x < x + h } .
Definition 5.
Let p h ( x ) be the number of observations in the set V h ( x ) . The localised simplex statistic at a given x is then defined as
U n , h ( x ) = p h ( x ) 3 1 t 1 < t 2 < t 3 sign ( a 1 y [ t 1 ] + a 2 y [ t 2 ] y [ t 3 ] ) × k = 1 3 K h ( x t k x ) ,
where a 1 + a 2 = 1 and K h ( v ) = K ( v / h ) where the kernel function K is to be specified.
In this paper, we assume a uniform kernel function since no addition information is given. Asymptotic normality also holds for the localised simplex statistics. A sub-interval is defined to be concave if U n , h < Z α / 2 , vice versa for the convex case. We define a concavity measure (localised) as the difference in the proportions of the concave intervals and the convex intervals:
L n , h = proportion of convex intervals proportion of concave intervals .
Here, L n , h [ 1 , 1 ] measures the level of concavity (convexity). Convexity is concluded for the overall interval when L n , h [ 0 , 1 ] , while we report strict concavity when L n , h [ 1 , 0 ) .
As a measure of tail index, the classic Hill’s estimator is implemented. If X 1 , X 2 , , X n is a sequence of independent and identically distributed random variables with distribution function F and n is the sample size, let X ( i , n ) be the i-th order statistic of X 1 , X 2 , , X n . The Hill’s estimator is
h ( k ( n ) , n ) Hill = 1 k ( n ) i = n k ( n ) + 1 n log ( X ( i , n ) ) log ( X ( n k ( n ) + 1 , n ) ) 1 ,
where k ( n ) { 1 , 2 , , n 1 } is a sequence to be specified. Under these assumptions, h ( k ( n ) , n ) Hill converges in probability to the true tail index, and is asymptotically normal when k ( n ) .

2.2. Simulation Results

We explore the possible range of behaviours for scaling functions for simulated multifractal processes. The findings are compared with some simulated monofractal processes (e.g., Brownian motions and fractional Brownian motions) and Student t-distributed processes. Both Brownian motion and fractional Brownian motion yield convex scaling functions only, which puts them out of consideration in the case of multifractality detection. The distributions of the concavity measures of the simulated multifractal processes and Student t-distributed processes are found to differ significantly in terms of their ranges and skewness. This provides us with a way to design a hypothesis test. As a result of these simulations, we construct a look-up table for a multifractality hypothesis test.
The simulated distributions of both the global and localised5 simplex statistics are quite similar. Here, we only present the simulation results based on the global simplex statistics. The simulation results are shown in Figure 3 and Figure 4. The distributions for the Student t-distributed process and a simulated multifractal process both span [ 1 , 1 ] but display different characteristics. A detailed comparison for the specific tail index = 3.06 is shown in Figure 4. For the simulated BMMTs, more than 20% of the scaling functions have global simplex statistics of exactly 1 indicating strict concavity. For the Student t-distributed process, the concavity statistic spans [ 1 , 1 ] with its median at around 0.36 in comparison to the one for the simulated BMMTs, 0.45 .
If the concavity of the observed scaling function for a given process is more severe than for the Student t-distributed processes, we are led to conclude that the concavity is driven by multifractality instead of heavy tails. Therefore, hypothesis tests can be constructed using the concavity measures as the test statistic.

2.3. Look-Up Table

To perform this hypothesis test on a given set of data, we firstly estimate the heaviness of tail using Hill’s estimator h Hill . Then, the null distribution of concavity measures is constructed by simulating from Student t-distributed processes corresponding to tail indices ranging over the interval [ h Hill 0.5 , h Hill + 0.5 ] . A larger range of tail index is used to allow for variation in our tail index estimation, while allowing more observations for constructing the concavity measure distribution.6 The 5th percentile of this distribution is the critical value for the hypothesis test.
The null hypothesis is rejected if the observed test statistic is less than the calculated critical value. To carry out the hypothesis test, a look-up table is generated to summarise the tail indices and the corresponding critical values based on both the localised and global concavity measures. A simplified version of the look-up table is presented in Appendix A. The critical values are presented for tail indices 0.5, 1.0, 1.5, …, with the largest tail index being 10.

3. Results

3.1. Application to Bitcoin

In this section, we apply our proposed hypothesis test to Bitcoin prices in order to address our research question: Do Bitcoin prices exhibit multifractal scaling? Both the global simplex statistic and the localised concavity measures are employed and reported. Two Bitcoin data sets are used in this paper:
  • the daily Bitcoin open price, daily data (in USD) from 28 April 2013 to 3 September 2019 with 2,320 observations, retrieved from CoinMarketCap (2019); and
  • the high frequency Bitcoin open price data minute-by-minute (in USDT) covering the period 22 May 2018 14:00 to 1 March 2019 11:00 with 406,089 observations, retrieved from Binance (2019).
Let { P ( t ) , t = 0 , 1 , , T } denote the open prices of Bitcoin. We define X ( t ) to be the mean-centered log-prices
X ( t ) = ln P ( t ) ln P ( 0 ) μ t ,
where μ = 1 T i = 1 T ln P ( i ) P ( i 1 ) . The mean-centered log-returns are then defined as
r ( t ) = X ( t ) X ( t 1 ) .

3.1.1. Daily Bitcoin Price Data

Mean-centered log-prices X ( t ) and the log-returns for daily Bitcoin price are displayed in Figure 5. Two bull runs followed by two bear runs are observed in late 2013 and late 2017 represented by two prominent spikes. These two price surges are believed to be results of the increases in Bitcoin’s popularity and media coverage. A price drop is spotted in early 2018 followed by the price falls of most cryptocurrencies. This period is known as ‘the 2018 cryptocurrency crash’ (Popken 2018). The Bitcoin price collapsed by 80% from January to September 2018, which is reported to be worse than the dot-com bubble’s 78% collapse (Patterson 2018). Subsequently, a price hike on a smaller scale is observed between early 2019 and 3 September 2019.
A stylised feature check gives uncorrelated and stationary return increments but correlations in the absolute value of returns.
The estimation of τ ( q ) exhibits a relatively linear scaling function, as is shown in Figure 6. However, both concavity measures are less than 0, indicating the presence of concavity. A global simplex statistic of value 0.2222 and a localised concavity measure of value 0.6256 are obtained. The hypothesis tests based on these two measures indicate that the null hypothesis should be rejected in favour of the alternative; see Table 1.
Thus, despite appearances in Figure 6, we have preliminary evidence that Bitcoin may follow a multifractal process.
To observe whether different results appear when the sample changes, the log-price time series is broken down into smaller time periods. The same hypothesis tests are performed to examine multifractality. Some literature indicates that 2017 is a turning point for Bitcoin’s price behaviour. Zhang et al. (2018) concluded an increase in market inefficiency from late 2016, which is believed to be a result of an increase in speculation. Consequently, we break the daily price process into two parts,
  • X 1 ( t ) from 28 April 2013 to 16 July 2017 with 1541 observations; and,
  • X 2 ( t ) from 17 July 2017 to 3 September 2019 with 779 observations.
The corresponding scaling functions are shown in Figure 7 with hypothesis test results in Table 2 and Table 3. The hypotheses tests for X 1 ( t ) give various results. This may be a result of a smaller data set or it could indicate a change to the way Bitcoin prices scale.

3.1.2. High-Frequency Bitcoin Price Data

Multifractality tests are performed on the high-frequency minute-by-minute data of Bitcoin open prices from 22 May 2018 14:00 to 1 March 2019 11:00 with 406,089 observations. Its mean-centered log-prices and log-returns are displayed in Figure 8. Similar to other financial assets, a lack of correlation in the return increments but high correlation in the absolute returns are found. However, compared with the daily open Bitcoin price data, the high-frequency Bitcoin data display stronger auto-correlation in the absolute return increments.
This data set covers approximately a one year period from May 2018 to March 2019. The two corresponding annual breakdown tests for daily Bitcoin open prices covering this period in Appendix B.1 in Appendix B do not indicate multifractal scaling, despite some evidence in Table 1 and Table 3 for multifractality in daily data, generally.
The analysis of the high-frequency minute-by-minute data are broadly consistent with this result. With the estimated scaling function τ ^ ( q ) shown in Figure 9 and the test results in Table 4, there appears to be little evidence of multifractal scaling for the high frequency data. Nevertheless, as shown in Appendix B.2, multifractality can be detected in some cases if we partition the data into smaller subperiods.
Given these contradictory results and that the high-frequency data only spans one year, we are reluctant to infer at this stage that the Bitcoin time series scales in the same way over longer or shorter time intervals.
Overall, it is possible that multifractal scaling exists for Bitcoin prices; however, a longer time series is necessary to substantiate the use of a multifractal model and further tests are needed to discern between monofractal models as a parsimonious alternative. In the following section, we test for the presence of multifractal scaling in other financial assets. These assets have more established markets with a long history of price data. It is thus useful to compare their scaling properties to that of the Bitcoin and to identify any similarities.

3.2. Bitcoin Compared to Other Financial Assets

Cryptocurrencies are a new financial asset whose value depends upon the evolution of its underlying technology as well as the design of its economic model. Bitcoin in particular appears to have characteristics that span a commodity, medium of exchange and technology. As a result, it is unclear which asset class Bitcoin most resembles. In the previous section, we identified that Bitcoin’s price dynamics could be characterised by multifractal scaling laws after accounting for heavy tails. In this section, we compare Bitcoin’s scaling properties to other financial assets. Included in our set of comparable assets are stock indices, a foreign exchange rate series and gold futures. The S&P500 index7 has been chosen to represent a large class of global equities; the Nasdaq Composite Index8 illustrates the price behaviour of technology-stocks; the USD/JPY exchange rate9 represents foreign currency; and gold futures prices10 are also included in the comparison, as Bitcoin is often referred to as the digital version of gold.
The results of the hypothesis tests are summarised in Table 5, Table 6 and Table 7 and the associated scaling functions are depicted in Figure 10, Figure 11 and Figure 12. Table 5 compares Bitcoin daily prices to other assets, over a short time period spanning Bitcoin’s limited but available price history, from 2013–2019. We observe that the local and global results can produce conflicting evidence. The local test for the S&P500 marginally accepts the null hypothesis, where the global test result marginally accepts the alternative in favour of multifractality for BTC Daily. The global test statistic reveals that the greatest evidence for multifractal scaling is for the USD/JPY exchange series. The scaling functions for each asset class look visibly concave for all assets except the S&P 500; however, multifractality in the presence of heavy tails is only indicated for Bitcoin and the the USD/JPY, indicating on prima facie evidence that Bitcoin could share similar scaling properties to foreign exchange. However, Table 6 and Table 7 reveal more compelling evidence given the larger data set spanning longer time periods. The lengthy time period is able to capture scaling behaviour for a wide range of time scales, as is necessary to establish a multifractal scaling relationship for an arbitrary set of rescaling factors. The results in Table 6 and Table 7 indicate that both the NASDAQ and the USD/JPY foreign exchange series could share multifractal scaling properties along with Bitcoin. The scaling function of the Gold Futures series looks visibly concave; however, the hypothesis test accepts the null in the presence of heavy tails indicating that the concavity is a result of heavy tails as opposed to multifractal scaling.

4. Conclusions

Multifractal scaling cannot be assumed to exist a priori, but instead should be established after accounting for the heavy tail effect on scaling functions. In this paper, we demonstrate how one can develop a multifractal scaling hypothesis test that accounts for the heavy tail effect on scaling functions. The test outlined in this paper doesn’t incorporate the possibility of monofractal scaling, but it can be incorporated into the test in the future and we will undertake such an approach in further research. The hypothesis test presented in this paper distinguishes between the heavy tail effect that distorts a linear scaling function to look concave, and true multifractal scaling. To implement the test, a look-up table is employed to simplify the hypothesis test procedure. This makes it easy to implement the test on various time series with marginal distributions of varying tail indexes. Our test results are of course contingent upon the validity of the underlying distributional assumptions. While a sound statistical theory awaits further development, a thorough examination of the test properties by means of Monte Carlo simulation can be carried out straightforwardly, in future research.
We apply this hypothesis test to Bitcoin prices and reveal that Bitcoin exhibits scaling behaviour more similar to a multifractal model than to a heavy tail process. We then extend the test to a set of other financial assets to ascertain whether Bitcoin prices are likely to share multifractal scaling relationships akin to other financial time series. Our results show that Bitcoin, USD/JPY exchange rates, and the technology heavy NASDAQ could all share multifractal scaling properties, after accounting for heavy tails. The findings of this paper are that, while Bitcoin prices span a relatively short period, the hypothesis test indicates that multifractal scaling is plausible and such scaling could be a feature of foreign exchange markets and technology stocks as well. This indicates some helpful methodology for model selection in risk analysis. Furthermore, the research suggests that financial time series may be classified by their statistical scaling properties in addition to the asset class they belong to. This additional type of classification could allow practitioners to construct statistically diverse portfolios based on assets grouped by their scaling dynamics.

Author Contributions

Research and analysis, C.J.; Supervision, P.D. and R.A.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by ARC grant DP160104737.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A. Look-Up Table

Table A1. Look-up tables for multifractality hypothesis test.
Table A1. Look-up tables for multifractality hypothesis test.
Localised Concavity MeasureGlobal Concavity Measure
Tail IndexCritical ValueTail IndexCritical Value
0.50−0.88890.50−0.9979
0.75−0.88890.75−0.9987
1.00−0.83331.00−0.9987
1.25−0.77781.25−0.9986
1.50−0.66671.50−0.9963
1.75−0.55561.75−0.8543
2.00−0.38892.00−0.7503
2.25−0.33332.25−0.5482
2.50−0.22222.50−0.5264
2.75−0.11112.75−0.4992
3.00−0.05563.00−0.6231
3.25−0.11113.25−0.7011
3.50−0.16673.50−0.7562
3.75−0.22223.75−0.8330
4.00−0.27784.00−0.8791
4.25−0.33334.25−0.9274
4.50−0.38894.50−0.9291
4.75−0.50004.75−0.9228
5.00−0.55565.00−0.9259
5.25−0.55565.25−0.9149
5.50−0.55565.50−0.9133
5.75−0.55565.75−0.9086
6.00−0.50836.00−0.9006
6.25−0.50006.25−0.8835
6.50−0.44446.50−0.8548
6.75−0.44446.75−0.8233
7.00−0.44447.00−0.8047
7.25−0.33337.25−0.7321
7.50−0.33337.50−0.7018
7.75−0.25287.75−0.6804
8.00−0.22228.00−0.6545
8.25−0.22228.25−0.6203
8.50−0.25288.50−0.5517
8.75−0.22228.75−0.4959
9.00−0.22229.00−0.4449
9.25−0.22229.25−0.4065
9.50−0.22229.50−0.3928
9.75−0.11679.75−0.3969
10.00−0.111110.00−0.3293

Appendix B. Multifractality Test Results

Appendix B.1. Hypothesis Test Results on Daily Bitcoin Open Price Data—Annual Breakdown

We explore for the existence of multifractality for each calendar year. The results are shown in Table A2 and Table A3.
Table A2. Local hypothesis test results on daily Bitcoin open prices.
Table A2. Local hypothesis test results on daily Bitcoin open prices.
Time PeriodTail IndexLocalised MeasureSample SizeTest Result
28/04/2013–31/12/20132.0886−0.55562026Multifractal
01/01/2014–31/12/20141.4785−0.33331904Non-Multifractal
01/01/2015–31/12/20150.94640.2222980Non-Multifractal
01/01/2016–31/12/20161.7584−0.77782036Multifractal
01/01/2017–31/12/20172.2187−0.61111997Multifractal
01/01/2018–31/12/20181.61450.16671977Non-Multifractal
01/01/2019–03/09/20191.44110.61111867Non-Multifractal
Table A3. Global Hypothesis test results on daily Bitcoin open prices.
Table A3. Global Hypothesis test results on daily Bitcoin open prices.
Time PeriodTail IndexGlobal MeasureSample SizeTest Result
28/04/2013–31/12/20132.0886−0.926072Multifractal
01/01/2014–31/12/20141.4785−0.730077Non-Multifractal
01/01/2015–31/12/20150.9464−0.051150Non-Multifractal
01/01/2016–31/12/20161.7584−0.947380Multifractal
01/01/2017–31/12/20172.2187−0.955170Multifractal
01/01/2018–31/12/20181.61450.270674Non-Multifractal
01/01/2019–03/09/20191.44110.843980Non-Multifractal

Appendix B.2. Hypothesis Test Results on High Frequency Bitcoin Price Data

The high frequency Bitcoin price data are divided into 6 equal-length time periods with 67,681 observations in each time period. The results are presented as follows:
Table A4. Local hypothesis test results on high frequency Bitcoin open prices.
Table A4. Local hypothesis test results on high frequency Bitcoin open prices.
Time PeriodTail IndexLocalisedSample SizeTest Result
22/05/2018 14:01–08/07/2018 23:232.12260.00002015Non-Multifractal
08/07/2018 23:24–24/08/2018 23:242.9599−0.33331572Multifractal
24/08/2018 23:25–11/10/2018 00:252.5649−0.44441889Multifractal
11/10/2018 00:26–27/11/2018 10:562.5984−0.66671847Multifractal
27/11/2018 10:57–13/01/2019 10:572.8381−0.27781654Multifractal
13/01/2019 10:58–01/03/2019 10:582.6245−0.50001845Multifractal
Table A5. Global Hypothesis test results on high frequency Bitcoin open prices.
Table A5. Global Hypothesis test results on high frequency Bitcoin open prices.
Time PeriodTail IndexGlobalSample SizeTest Result
22/05/2018 14:01–08/07/2018 23:232.1226−0.481172Non-Multifractal
08/07/2018 23:24–24/08/2018 23:242.9599−0.387296Non-Multifractal
24/08/2018 23:25–11/10/2018 00:252.5649−0.558774Multifractal
11/10/2018 00:26–27/11/2018 10:562.5984−0.828774Multifractal
27/11/2018 10:57–13/01/2019 10:572.8381−0.669088Multifractal
13/01/2019 10:58–01/03/2019 10:582.6245−0.412575Non-Multifractal

References

  1. Abrevaya, Jason, and Wei Jiang. 2005. A nonparametric approach to measuring and testing curvature. Journal of Business & Economic Statistics 23: 1–19. [Google Scholar]
  2. Bariviera, Aurelio F. 2017. The inefficiency of bitcoin revisited: A dynamic approach. Economics Letters 161: 1–4. [Google Scholar] [CrossRef] [Green Version]
  3. Bariviera, Aurelio F., María José Basgall, Waldo Hasperué, and Marcelo Naiouf. 2017. Some stylized facts of the bitcoin market. Physica A: Statistical Mechanics and Its Applications 484: 82–90. [Google Scholar] [CrossRef] [Green Version]
  4. Binance. 2019. BTC/USDT. Available online: https://www.binance.com/en/trade/BTC_USDT (accessed on 20 March 2019).
  5. Bouchaud, Jean-Philippe, Marc Potters, and Martin Meyer. 2000. Apparent multifractality in financial time series. The European Physical Journal B-Condensed Matter and Complex Systems 13: 595–99. [Google Scholar] [CrossRef] [Green Version]
  6. Calvet, Laurent E., Adlai J. Fisher, and Benoit B. Mandelbrot. 1997. Large Deviations and the Distribution of Price Changes. Cowles Foundation Discussion Paper No. 1165. New Haven: Cowles Foundation for Research in Economics Yale University. [Google Scholar]
  7. Chambers, Clem. 2019. Bitcoin Is Fractal. Available online: https://www.forbes.com/sites/investor/2019/01/23/bitcoin-is-fractal/#27bbcc49208d (accessed on 30 April 2019).
  8. CoinMarketCap. 2019. Bitcoin (BTC) Price, Charts, Market Cap, and Other Metrics. Available online: https://coinmarketcap.com/currencies/bitcoin/ (accessed on 5 September 2019).
  9. Embrechts, Paul, and Makoto Maejima. 2000. An introduction to the theory of self-similar stochastic processes. International Journal of Modern Physics B 14: 1399–420. [Google Scholar] [CrossRef]
  10. Grahovac, Danijel, and Nikolai N. Leonenko. 2014. Detecting multifractal stochastic processes under heavy-tailed effects. Chaos, Solitons & Fractals 65: 78–89. [Google Scholar]
  11. Heyde, Chris C. 2009. Scaling issues for risky asset modelling. Mathematical Methods of Operations Research 69: 593–603. [Google Scholar] [CrossRef]
  12. investing.com Australia. 2020a. Gold Futures—Apr 20 (GCJ0). Available online: https://au.investing.com/commodities/gold-historical-data (accessed on 27 February 2020).
  13. investing.com Australia. 2020b. USD JPY Historical Data—Investing.com AU. Available online: https://au.investing.com/currencies/usd-jpy-historical-data (accessed on 27 February 2020).
  14. Kantelhardt, Jan W., Stephan A. Zschiegner, Eva Koscielny-Bunde, Shlomo Havlin, Armin Bunde, and H. Eugene Stanley. 2002. Multifractal detrended fluctuation analysis of nonstationary time series. Physica A: Statistical Mechanics and Its Applications 316: 87–114. [Google Scholar] [CrossRef] [Green Version]
  15. Lahmiri, Salim, and Stelios Bekiros. 2018. Chaos, randomness and multi-fractality in bitcoin market. Chaos, Solitons & Fractals 106: 28–34. [Google Scholar]
  16. Mandelbrot, Benoit B., Adlai J. Fisher, and Laurent E. Calvet. 1997. A Multifractal Model of Asset Returns. Cowles Foundation Discussion Paper No. 1164. New Haven: Cowles Foundation for Research in Economics Yale University. [Google Scholar]
  17. Matia, Kaushik, Yosef Ashkenazy, and H. Eugene Stanley. 2003. Multifractal properties of price fluctuations of stocks and commodities. EPL (Europhysics Letters) 61: 422. [Google Scholar] [CrossRef] [Green Version]
  18. Mensi, Walid, Yun-Jung Lee, Khamis Hamed Al-Yahyaee, Ahmet Sensoy, and Seong-Min Yoon. 2019. Intraday downward/upward multifractality and long memory in bitcoin and ethereum markets: An asymmetric multifractal detrended fluctuation analysis. Finance Research Letters 31: 19–25. [Google Scholar] [CrossRef]
  19. Milutinović, Monia. 2018. Cryptocurrency. Екoнoмика-Часoпис за екoнoмску теoрију и праксу и друштвена питањаа 1: 105–122. [Google Scholar] [CrossRef] [Green Version]
  20. Nadarajah, Saralees, and Jeffrey Chu. 2017. On the inefficiency of bitcoin. Economics Letters 150: 6–9. [Google Scholar] [CrossRef] [Green Version]
  21. Nakamoto, Satoshi. 2009. Bitcoin: A Peer-to-Peer Electronic Cash System. Available online: https://bitcoin.org/bitcoin.pdf (accessed on 20 November 2018).
  22. Patterson, Michael. 2018. Crypto’s 80% Plunge Is Now Worse Than the Dot-Com Crash. Available online: https://www.bloomberg.com/news/articles/2018-09-12/crypto-s-crash-just-surpassed-dot-com-levels-as-losses-reach-80 (accessed on 7 August 2019).
  23. Platen, Eckhard, and Renata Rendek. 2008. Empirical evidence on student-t log-returns of diversified world stock indices. Journal of Statistical Theory and Practice 2: 233–251. [Google Scholar] [CrossRef] [Green Version]
  24. Popken, Ben. 2018. Bitcoin Loses More Than Half Its Value Amid Crypto Crash. Available online: https://www.nbcnews.com/tech/internet/bitcoin-loses-more-half-its-value-amid-crypto-crash-n844056 (accessed on 7 August 2019).
  25. Salat, Hadrien, Roberto Murcio, and Elsa Arcaute. 2017. Multifractal methodology. Physica A: Statistical Mechanics and Its Applications 473: 467–487. [Google Scholar] [CrossRef]
  26. Sly, Allan. 2006. Self-similarity, multifractionality and multifractality. Master’s thesis, Australian National University, Canberra, Australia. [Google Scholar]
  27. Sreenivasan, Katepalli. 1991. Fractals and multifractals in fluid turbulence. Annual Review of Fluid Mechanics 23: 539–604. [Google Scholar] [CrossRef]
  28. Stavroyiannis, Stavros, Vassilios Babalos, Stelios Bekiros, Salim Lahmiri, and Gazi Salah Uddin. 2019. The high frequency multifractal properties of bitcoin. Physica A: Statistical Mechanics and Its Applications 520: 62–71. [Google Scholar] [CrossRef]
  29. Yahoo Finance. 2019. S&P 500 (^GSPC) Historical Data. Available online: https://au.finance.yahoo.com/quote/^GSPC/history?p=^GSPC (accessed on 27 February 2020).
  30. Yahoo Finance. 2020. NASDAQ Composite (^IXIC) Historical Data. Available online: https://finance.yahoo.com/quote/%5Eixic/history?ltr=1 (accessed on 27 February 2020).
  31. Zhang, Wei, Pengfei Wang, Xiao Li, and Dehua Shen. 2018. The inefficiency of cryptocurrency and its cross-correlation with dow jones industrial average. Physica A: Statistical Mechanics and Its Applications 510: 658–670. [Google Scholar] [CrossRef]
1.
“False positive” corresponds to the scenario that we mistakenly detect multifractality when the underlying process does not possess the multifractal property.
2.
For convenience, in the rest of this paper, we call uniscaling multifractal processes monofractal processes while referring to multiscaling multifractal processes as multifractal processes.
3.
We compared the scaling functions of the Bitcoin and other financial assets with the ones of BMMTs simulated using multiplicative cascades with Poisson distribution, Gamma distribution and Normal distribution. The scaling function of the BMMT simulated through log-Normal multiplicative cascade displays the most similar behaviour.
4.
The choice of Student’s t-distribution with 4 degrees of freedom is suggested by the findings of Platen and Rendek (2008).
5.
Knots were chosen so that the function is evaluated on 18 equal sub-intervals. This gave the best approximation considering computational efficiency.
6.
For example, if Hill’s estimator takes value h Hill = 2 , we generate simulated student t-distributed processes with 2 degrees of freedom, select those with tail indices in the interval [ 1.5 , 2.5 ] , and construct the null distribution using the empirical distribution of their concavity measures.
7.
Retrieved from Yahoo Finance (2019) for the period 30 December 1927 to 26 February 2020 with 23,147 observations.
8.
Retrieved from Yahoo Finance (2020) for the period 5 February 1971 to 25 February 2020 with 12,372 observations.
9.
Retrieved from investing.com Australia (2020b) for the period 4 March 1988 to 28 February 2020 with 8337 observations.
10.
Retrieved from investing.com Australia (2020a) for the period 27 December 1979 to 28 February 2020 with 10,190 observations.
Figure 1. Asymptotic scaling functions for various values of α . The grey line is the q / 2 reference line.
Figure 1. Asymptotic scaling functions for various values of α . The grey line is the q / 2 reference line.
Jrfm 13 00104 g001
Figure 2. Left—Scaling Functions of the original S&P 500 open price data (black solid line), the S&P 500 open price data after removing 19 October 1987 (dotted line), the S&P 500 open price data after truncating any return larger than 4 standard deviations from the mean (grey solid line). Right—Scaling Functions of the Student t-distributed processes before (black lines) and after truncation (grey lines). For comparison, the reference lines q / 2 (black dashed line) are included.
Figure 2. Left—Scaling Functions of the original S&P 500 open price data (black solid line), the S&P 500 open price data after removing 19 October 1987 (dotted line), the S&P 500 open price data after truncating any return larger than 4 standard deviations from the mean (grey solid line). Right—Scaling Functions of the Student t-distributed processes before (black lines) and after truncation (grey lines). For comparison, the reference lines q / 2 (black dashed line) are included.
Jrfm 13 00104 g002
Figure 3. Distribution of global simplex statistics for Brownian motion (left) and fractional Brownian motion (right).
Figure 3. Distribution of global simplex statistics for Brownian motion (left) and fractional Brownian motion (right).
Jrfm 13 00104 g003
Figure 4. Distributions of global simplex statistics for Student t-distributed process (black) and the simulated BMMT with log-normal multiplicative cascade (grey) around tail index 3.06. The two vertical lines represent the 95th percentile for the simulated BMMTs (left) and 5th percentiles for the Student t-distributed processes, respectively (right).
Figure 4. Distributions of global simplex statistics for Student t-distributed process (black) and the simulated BMMT with log-normal multiplicative cascade (grey) around tail index 3.06. The two vertical lines represent the 95th percentile for the simulated BMMTs (left) and 5th percentiles for the Student t-distributed processes, respectively (right).
Jrfm 13 00104 g004
Figure 5. Daily mean-centered log Bitcoin open price process (top) and mean-centered log returns process (bottom).
Figure 5. Daily mean-centered log Bitcoin open price process (top) and mean-centered log returns process (bottom).
Jrfm 13 00104 g005
Figure 6. Left—Autocorrelations of mean-centered log-returns (black) and absolute returns (grey) for daily Bitcoin open price (28 April 2013–3 April 2019). Right—Scaling function for daily Bitcoin open price (28 April 2013–3 September 2019). The dashed line is the q / 2 reference line (same for all the scaling function figures).
Figure 6. Left—Autocorrelations of mean-centered log-returns (black) and absolute returns (grey) for daily Bitcoin open price (28 April 2013–3 April 2019). Right—Scaling function for daily Bitcoin open price (28 April 2013–3 September 2019). The dashed line is the q / 2 reference line (same for all the scaling function figures).
Jrfm 13 00104 g006
Figure 7. Scaling functions for Daily Bitcoin open price from 28 April 2013 to 16 July 2017 (left) vs. from 17 July 2017 to 3 September 2019 (right).
Figure 7. Scaling functions for Daily Bitcoin open price from 28 April 2013 to 16 July 2017 (left) vs. from 17 July 2017 to 3 September 2019 (right).
Jrfm 13 00104 g007
Figure 8. High frequency mean-centered log Bitcoin open price process (top) and mean-centered log returns process (bottom).
Figure 8. High frequency mean-centered log Bitcoin open price process (top) and mean-centered log returns process (bottom).
Jrfm 13 00104 g008
Figure 9. Left—Autocorrelations of mean-centered log returns (black) and absolute returns (grey) for high frequency Bitcoin open prices (22 May 2018–1 March 2019). Right—Scaling function for high frequency Bitcoin open prices (22 May 2018–01 March 2019).
Figure 9. Left—Autocorrelations of mean-centered log returns (black) and absolute returns (grey) for high frequency Bitcoin open prices (22 May 2018–1 March 2019). Right—Scaling function for high frequency Bitcoin open prices (22 May 2018–01 March 2019).
Jrfm 13 00104 g009
Figure 10. Scaling functions for S&P500, NASDAQ Index, USD/JPY Exchange and Gold Futures Prices the period 28 April 2013 to 3 September 2019.
Figure 10. Scaling functions for S&P500, NASDAQ Index, USD/JPY Exchange and Gold Futures Prices the period 28 April 2013 to 3 September 2019.
Jrfm 13 00104 g010
Figure 11. Scaling functions for S&P500, NASDAQ Index, USD/JPY Exchange and Gold Futures Prices during the dot-com bubble (3 January 1994–8 October 2004). The time period from 3 January 1994 to 8 October 2004 corresponds to the dot-com bubble and the dot-com crash.
Figure 11. Scaling functions for S&P500, NASDAQ Index, USD/JPY Exchange and Gold Futures Prices during the dot-com bubble (3 January 1994–8 October 2004). The time period from 3 January 1994 to 8 October 2004 corresponds to the dot-com bubble and the dot-com crash.
Jrfm 13 00104 g011
Figure 12. Scaling functions for S&P500, NASDAQ Index, USD/JPY Exchange and Gold Futures Prices the period 4 March 1988 to 3 September 2019.
Figure 12. Scaling functions for S&P500, NASDAQ Index, USD/JPY Exchange and Gold Futures Prices the period 4 March 1988 to 3 September 2019.
Jrfm 13 00104 g012
Table 1. Hypothesis test results on daily Bitcoin open prices (28 April 2013–3 September 2019).
Table 1. Hypothesis test results on daily Bitcoin open prices (28 April 2013–3 September 2019).
Tail Index3.0630
Sample Size 1Test StatisticRejection RegionTest Result
Localised Test1479 0.2222 [ 1 , 0.0556 ) Multifractal
Global Test99 0.6256 [ 1 , 0.6148 ) Multifractal
1Sample size refers to the number of student t-distributed processes used in estimating the distribution of concavity measures under H 0 . In the case of Bitcoin open prices from 28/04/2013 to 03/09/2019, there are 1479 student t-distributed processes with tail indices ranging in the interval [2.5630, 3.5630] when constructing the distribution of localised concavity measures. 99 student t-distributed processes are used in the construction of the distributions of global simplex statistics.
Table 2. Hypothesis test results on daily Bitcoin open prices (28 April 2013–16 July 2017).
Table 2. Hypothesis test results on daily Bitcoin open prices (28 April 2013–16 July 2017).
Tail Index3.2791
Sample SizeTest StatisticRejection RegionTest Result
Localised Test1403 0.4444 [ 1 , 0.1111 ) Multifractal
Global Test98 0.5378 [ 1 , 0.7014 ) Non-Multifractal
Table 3. Hypothesis test results on daily Bitcoin open prices (17 July 2017–3 September 2019).
Table 3. Hypothesis test results on daily Bitcoin open prices (17 July 2017–3 September 2019).
Tail Index1.6528
Sample SizeTest StatisticRejection RegionTest Result
Localised Test1992 0.8333 [ 1 , 0.6111 ) Multifractal
Global Test75 0.9950 [ 1 , 0.8703 ) Multifractal
Table 4. Hypothesis test results on high frequency Bitcoin open prices (22 May 2018–1 March 2019).
Table 4. Hypothesis test results on high frequency Bitcoin open prices (22 May 2018–1 March 2019).
Tail Index2.7759
Sample SizeTest StatisticRejection RegionTest Result
Localised Test1682 0.2778 [ 1 , 0.0556 ) Non-multifractal
Global Test83 0.2714 [ 1 , 0.5049 ) Non-multifractal
Table 5. Test results among different financial assets (28 April 2013–3 September 2019).
Table 5. Test results among different financial assets (28 April 2013–3 September 2019).
Financial AssetTail IndexTestTest StatisticsRejection RegionTest Result (Local)
BTC Daily3.06Local−0.22[−1, −0.06)Multifractal
Global−0.63[−1, −0.61)Multifractal
S&P5003.13Local−0.11[−1, −0.11)Non-Multifractal
Global−0.14[−1, −0.66)Non-Multifractal
NASDAQ3.19Local0.11[−1, −0.11)Non-Multifractal
Global−0.44[−1, −0.66)Non-Multifractal
USD/JPY2.69Local−0.28[−1, −0.17)Multifractal
Global−0.86[−1, −0.49)Multifractal
Gold Futures1.35Local−0.67[−1, −0.72)Non-Multifractal
Global−0.89[−1, −0.9979)Non-Multifractal
Table 6. Test results among different financial assets during the dot-com bubble (3 January 1994–8 October 2004).
Table 6. Test results among different financial assets during the dot-com bubble (3 January 1994–8 October 2004).
Financial AssetTail IndexTestTest StatisticsRejection RegionTest Result (Local)
S&P5003.41Local−0.22[−1, −0.17)Multifractal
Global−0.40[−1, −0.71)Non-Multifractal
NASDAQ2.81Local−0.44[−1, −0.06)Multifractal
Global−0.79[−1, −0.52)Multifractal
USD/JPY3.50Local−0.33[−1, −0.17)Multifractal
Global−0.80[−1, −0.76)Multifractal
Gold Futures3.48Local0.22[−1, −0.17)Non-Multifractal
Global−0.48[−1, −0.71)Non-Multifractal
Table 7. Test results among different financial assets (4 March 1988–3 September 2019).
Table 7. Test results among different financial assets (4 March 1988–3 September 2019).
Financial AssetTail IndexTestTest StatisticsRejection RegionTest Result (Local)
S&P5003.20Local0.78[−1, −0.11)Non-Multifractal
Global0.93[−1, −0.66)Non-Multifractal
NASDAQ3.84Local−0.61[−1, −0.22)Multifractal
Global−0.93[−1, −0.84)Multifractal
USD/JPY3.33Local−0.28[−1, −0.17)Multifractal
Global−0.72[−1, −0.70)Multifractal
Gold Futures1.36Local−0.11[−1, −0.72)Non-Multifractal
Global0.11[−1, −0.9979)Non-Multifractal

Share and Cite

MDPI and ACS Style

Jiang, C.; Dev, P.; Maller, R.A. A Hypothesis Test Method for Detecting Multifractal Scaling, Applied to Bitcoin Prices. J. Risk Financial Manag. 2020, 13, 104. https://doi.org/10.3390/jrfm13050104

AMA Style

Jiang C, Dev P, Maller RA. A Hypothesis Test Method for Detecting Multifractal Scaling, Applied to Bitcoin Prices. Journal of Risk and Financial Management. 2020; 13(5):104. https://doi.org/10.3390/jrfm13050104

Chicago/Turabian Style

Jiang, Chuxuan, Priya Dev, and Ross A. Maller. 2020. "A Hypothesis Test Method for Detecting Multifractal Scaling, Applied to Bitcoin Prices" Journal of Risk and Financial Management 13, no. 5: 104. https://doi.org/10.3390/jrfm13050104

APA Style

Jiang, C., Dev, P., & Maller, R. A. (2020). A Hypothesis Test Method for Detecting Multifractal Scaling, Applied to Bitcoin Prices. Journal of Risk and Financial Management, 13(5), 104. https://doi.org/10.3390/jrfm13050104

Article Metrics

Back to TopTop