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

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.


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 positive 1 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.

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 (1) 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) = Hq. 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, where equality is of finite dimensional distributions.
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.
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.

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 where the q-th moment of the process X(t), E|X(t)| q , can be estimated by the sample statistic 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.

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 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: 1. 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. 2. 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 since the multipliers at different stages are independent. The measure µ 2 represents the multiplicative measure at Stage 2. 3. 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 ∑ b−1 β=0 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(∑ b−1 β=0 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 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 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, 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 = ∑ k i=1 η 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 . 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 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 series 3 . When V follows a normal distribution, 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: 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

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 andGrahovac andLeonenko 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, 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 , as n → ∞, where p → stands for convergence in probability. 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. 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.
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 EY(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) where plim stands for limit in probability and 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

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.

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 )} n i=1 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; .., n, and the functional form of f is left unspecified.
Then, for a sample of size n, the global simplex statistic is defined as where t 1 , t 2 , t 3 ∈ {1, . . . , n}, 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Ũ 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/2h subintervals, where 2h 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 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 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.
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 where k(n) ∈ {1, 2, ..., n − 1} is a sequence to be specified. Under these assumptions, h Hill (k(n),n) converges in probability to the true tail index, and is asymptotically normal when k(n) → ∞.

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 localised 5 simplex statistics are quite similar. Here, we only present the simulation results based on the global simplex statistics. The simulation results are shown in Figures 3 and 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.

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 5 Knots were chosen so that the function is evaluated on 18 equal sub-intervals. This gave the best approximation considering computational efficiency. 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.

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: We define X(t) to be the mean-centered log-prices . The mean-centered log-returns are then defined as r(t) = X(t) − X(t − 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. 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.   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 Tables 2 and 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.

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.  Tables 1  and 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.

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 index 7 has been chosen to represent a large class of global equities; the Nasdaq Composite Index 8 illustrates the price behaviour of technology-stocks; the USD/JPY exchange rate 9 represents foreign currency; and gold futures prices 10 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 Tables 5-7 and the associated scaling functions are depicted in Figures 10-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, Tables 6 and 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 Tables 6 and 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.

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.

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. Table   Table A1. Look-up tables for multifractality hypothesis test. We explore for the existence of multifractality for each calendar year. The results are shown in Tables A2 and A3.