key: cord-0067138-la3jcz7f authors: David, S.A.; Inacio Jr., C.M.C.; Nunes, R.; Machado, J.A.T. title: Fractional and fractal processes applied to cryptocurrencies price series date: 2021-01-07 journal: J Adv Res DOI: 10.1016/j.jare.2020.12.012 sha: 6ee9891470acdb581365719af8ce666d1b3b718a doc_id: 67138 cord_uid: la3jcz7f INTRODUCTION: Cryptocurrencies have been attracting the attention from media, investors, regulators and academia during the last years. In spite of some scepticism in the financial area, cryptocurrencies are a relevant subject of academic research. OBJECTIVES: In this paper, several tools are adopted as an instrument that can help market agents and investors to more clearly assess the cryptocurrencies price dynamics and, thus, guide investment decisions more assertively while mitigating risks. METHODS: We consider three methods, namely the Auto-Regressive Integrated Moving Average (ARIMA), Auto-Regressive Fractionally Integrated Moving Average (ARFIMA) and Detrended Fluctuation Analysis, and three indices given by the Hurst and Lyapunov exponents or the Fractal Dimension. This information allows assessing the behaviour of the time series, such as their persistence, randomness, predictability and chaoticity. RESULTS: The results suggest that, except for the Bitcoin, the other cryptocurrencies exhibit the characteristic of mean reverting, showing a lower predictability when compared to the Bitcoin. The results for the Bitcoin also indicate a persistent behavior that is related to the long memory effect. CONCLUSIONS: The ARFIMA reveals better predictive performance than the ARIMA for all cryptocurrencies. Indeed, the obtained residual values for the ARFIMA are smaller for the auto and partial auto correlations functions, as well as for confidence intervals. Automation and cognitive innovation continue to apace, creating opportunities to radically simplify the processes adopted by the human kind. As this transition picks up speed the capacity to add value tend also to be unleashed. In this context, blockchain may play a key role in the future, addressing concepts and tools, such as artificial intelligence, 5G communications, internet of things, and digital finance [1] . Important countries such as China [2, 3] launched an ambitious effort to improve blockchain technology, for issuing digital money and to streamle government services [4] . This strategy evinces that governments understand that to become a high-tech power it is necessary to be positioned at the forefront of the new technologies and that the blockchain is an integral part of this process. Indeed, financial transactions are becoming touchless, as automation and blockchain mature accelerating, therefore, this trend. In this scenario, cryptocurrencies play a major role and become more popular as their use spreads around the world. Indeed, cryptocurrencies can change the procedures adopted for financial transactions [5] [6] [7] and attention must be paid to this novel paradigm. The first and most important digital currency at the present date is the bitcoin (BTC). The BTC first records are dated in 2008 and was introduced with the initial objective of mitigating costs related to electronic transactions (e-commerce) [8] . Nonetheless, the BTC became more popular due to the anonymity that offers in transactions, as well as its independence from traditional financial providers, such as banks or brokers, among others. Baeck and Elbeck [9] point that the BTC is more a speculative commodity rather than a currency. Nonetheless, investors have employed BTC not only as a currency, but also as an investment [10] . Over the year 2013 the BTC skyrocketed to an amazing 8000% of growth of its value, becoming, consequently, the major bull market of that period [8] . Nonetheless, problems appeared due to the increasing amount of users and the popularization of digital transactions, such as the so-called ''double-spending", which is related to the possibility of a same digital currency being involved in two transactions at the same time. Due to this practice, the cryptocurrency became an easy target for hackers [11, 12] . One idea that emerged as a possible solution for the problem was to create a public ledger, where all the information about negotiation would be recorded, including reports about the buyer and the seller of BTC. However, this idea did not completely solved the problem, since it carried to the centralization of the cryptocurrency, and also the lack of anonymity of the transactions [8, 11] . A protection system denoted ''blockchain" was created, consisting of a structure where all the cryptocurrencies are negotiated, serving as a public ledger and being transparent to all transactions. Blockchain is an exceptional anti-fraud system, since every modification in a given transaction (i.e., a block) results in a complete modification of all the precedent blocks posing, therefore, a significant computational load [13, 14] . In recent years, the emergence of the altcoins, that is, cryptocurrencies derivated from the BTC open code, contributed to liquidity and competition in the digital currency market. Among the multitude of present day altcoins, we highlight those explored in this work: Litecoin (LTC), Ripple (XRP), Monero (XMR), Ethereum (ETH), and Ethereum Classic (ETC) released during 2011, 2012, 2014, 2015 and 2016, respectively. How the altcoins will influence the BTC price dynamics is a frequent and still open question. Cryptocurrencies and especially the BTC, have shown to be an innovative and alternative asset that has evolved in the context of absolute returns, volatilities and correlations [15] . At the same time, investors in innovative assets demand for more information to make sure that they are making an appropriate investment while mitigating risks. It is well known that investments in traditional markets require access to data from companies, industries and even the global economy in order to obtain the aforementioned information. Usually, financial advisors and investors rely on analysts to process this data and to provide insight into those assets, that is, to provide an extended and solid investment analysis. Nonetheless, crypto-assets are not related to traditional companies and the processes of evaluating and monitoring their price dynamics is different from those adopted for company stock prices [16] . The main contribution of this work is the use of fractional and fractal mathematical tools, so that market agents and investors have a more refined instrument to assess clearly the cryptocurrencies price dynamics. In this context, recent advances [15] suggest that fractals and fractional calculus are interesting tools to deal with complex behavior with memory effect [17] [18] [19] [20] [21] . Non-integer orders of differentiation and integration describe phenomena such as long-term memory, non-locality, persistence, predictability, chaoticity, thus encouraging their application in oscillatory phenomena [22] , control algorithms [23, 24] , biosystems [25] [26] [27] [28] [29] [30] , economy and finance [31] [32] [33] [34] , among others [35] [36] [37] . The prediction of the cryptocurrencies prices is of key importance since a small advantage obtained through accurate information can lead to significant profits [38] [39] [40] 10] . Lahmiri and Bekiros [38] implemented a novel deep learning model to forecast the Bitcoin, Digital Cash, and Ripple digital currencies prices series. These researchers obseved that a chaotic dynamics, led by long memory effects, can be better understand using fractional operators. Corbet et al. [39] investigated the effect of futures contracts for the Bitcoin market and concluded that the increasing volatility leads to a null effectivity of pricing protection tools, such as, the hedging technique. Thus, speculative approaches are necessary in order to perform in the Bitcoin futures trading. Balcilar et al. [40] examined the causal relation between Bitcoin returns, volatility, and trading volume. It was observed a significant evidence of causality between the returns and trading volume, but no evidence of any causal relation with the volatility was found. This pattern can help models using the return-volume relationship, except during bear and bull market regimes. Despite that, Urquhart [10] found evidence of inefficiency in the Bitcoin market indicating a less random market with memory effects in the prices. In fact, Bitcoin may be moving towards market efficiency. Therefore, for investors to succeed in the digital currency market, it is important to understand the price dynamics of such cryptocurrencies, especially for the fractal dynamics and long-memory effects. Mensi et al. [41] observed that there is a double long-memory process between the BTC and the ETH prices. Aharon and Qadan [42] , and Caporale and Plastum [43] investigated the relationships between the day of the week and the BTC price and, furthermore, the relationship between the prices of the BTC and other altcoins. The so-called ''Monday effect" was found for the returns of the BTC. However, that behavior was not observed for the LTC, XRP and Dash (DASH) cryptocurrencies. Balcilar et al. [40] identified another BTC feature when exploring the relationship between the trading volume, BTC returns, and volatility. The study suggested that the trading volume does not helps in predicting the volatility of the BTC returns. Briere et al. [44] showed that the BTC allows the diversification of the investors benefits, and Cheung et al. [45] unveiled the presence of bubbles in the BTC market. Corbet and Katsiampa [46] noted the asymmetries in the BTC price return series' persistence and found evidences of an higher persistence of positive than negative returns. Corbet et al. [47] explored the effectiveness of technical trading rules in the cryptocurrency markets and obtained results supporting that a variable-length moving average strategy in buy signals outperforms the sell signals in terms of returns. Catania et al. [48, 49] verified that, besides the long-memory, the leverage effect also contributes to the volatility of the price time series (TS), because of the different foreign exchange currencies. Bouri et al. [50] analysed the long-memory process in the structural breaks along with the records of the BTC price. Caporale et al. [51] found the same result for other altcoins, such as the LTC, XRP, and DASH. Since the Efficient-Market hypothesis (EMH) is related to a random-walk type of phenomenon [51, 52] , studies also suggest that the BTC may be in the process of moving towards an efficient market over the last years [10, [53] [54] [55] . This paper explores the dynamics of the price TS for the BTC, LTC, XRP, XMR, ETH and ETC cryptocurrencies. Several mathematical tools are adopted, namely the Auto-Regressive Integrated Moving Average (ARIMA), Auto-Regressive Fractionally Integrated Moving Average (ARFIMA) and Detrended Fluctuation Analysis (DFA) algorithms. Additionally, three metrics corresponding to the Hurst (H) index, Lyapunov (k) exponent and Fractal Dimension (d A ) are also used. The information provided by the distinct approaches captures different characteristics of the price TS, such as, persistence, randomness, predictability and chaoticity. The paper has the following organization. Section 2 introduces the price TS and the proposed methods. Section 3 discusses the results obtained by the mathematical and algorithmic tools. Finally, Section 4 gives the main conclusions. Cryptocurrencies have shown to be highly volatile, allowing significant profit opportunities for experts in this market [56] . Nonetheless, this market, as well as the financial and economic systems, requires accurate mathematical models to understand their complex behavior [57, 58] . The ARIMA and ARFIMA techniques are usually employed in studies involving TS. David et al. [59, 60] applied these tools in the description of energy and agricultural commodities and concluded that the ARFIMA shows a better performance than the ARIMA. Several papers adopt the Hurst exponent H to estimate the market efficiency [61, 51, [62] [63] [64] . The Lyapunov exponent k is also an effective index to indicate the presence of chaos and provides relevant information about the TS predictability [65, 66] . The concept of fractals and multifractals have been also extensively applied to evaluate the market efficiency of cryptocurrencies and other financial TS [67] [68] [69] [70] [71] 53, 72] . Bearing these facts in mind and to achieve an effective forecasting, we adopt the aforementioned mathematical tools to explore the dynamics of cryptocurrencies price TS. The data was obtained in the website https://www.investing.com/crypto/currencies. The TS prices available for the ETC starts in the year of 2016. For this reason, we use daily data (7 values per week and without holidays) of closure prices from 2016 onwards. Furthermore, we consider the same period for all cryptocurrencies investigated in this work, as depicted in Fig. 1 . The time interval explored goes from July/2016 to March/2019 for the prices TS composition, and from April/2019 to Jun/2019 for the future prices prediction. The ARIMA(p,d,q) and ARFIMA(p,d,q) models are applied in this work to predict the prices of the cryptocurrencies. The parameters p and q 2 N and d 2 R stand for the order (number of time lags) of the autoregressive models, the degree of differencing and the order of the moving average, respectively. The ARIMA can be interpreted as a combination of the autoregressive and moving-average models. The ARFIMA [73] generalizes the ARIMA, so that the parameter d can assume non-integer values. The price TS, P t , is integrated and leads to combined process of auto-regressive (AR), integrated (I) and moving average (MA). The ARIMA can be written by means of the discrepancy operator B that, by definition, can be approximated using the Wold decomposition [74] , the autoregressive and moving-average operators, respectively. Therefore, the ARIMA(p,d,q) can be written as: where e t denotes a white noise process, and 1 À B ð Þ d is the differencing operator in the autoregressive model parameters. Similarly, the ARFIMA(p,d,q) is defined by: An asymptotic approximation of n j is given by: where C is the gamma function. The ARFIMA(p,d,q) can grasp the dynamics of a long-range memory process [75, 73, 76] . The Autocorrelation Function (ACF) and the Partial Autocorrelation Function (PACF) correlograms are plotted for obtaining the orders of p and q of the prediction models (see A.7 and A.8). If necessary, the TS are differenced to become stationary. To compare the orders of the models, the Bayesian Information Criterion (BIC) criterion [77] is adopted in the follow-up. Also, the Seasonal and Trend decomposition using the Loess (STL) procedure is applied, where Loess represents an estimating technique for nonlinear relationships [77] . Table 1 lists the BIC values obtained by means of the ARIMA and ARFIMA for the six cryptocurrencies. The model with the smallest BIC value is selected for each technique. The smaller the BIC value, the better the fitting for the cryptocurrency TS [77] . It is known that financial TS may exhibit long-range dependence [78, 79, 77, 59] . The same characteristic can be observed in cryptocurrencies price TS [77, 41, [48] [49] [50] [51] . The long-range dependence is related to the fractional Brownian motion (fBm) [80] that generalizes the classical Brownian motion (Bm) describing a random walk. The properties of the Bm and fBm can be quantified by the Hurst exponent H [62] . Bearing in mind that H 2 0; 1 ð , the main characteristics of the Hurst exponent can be summarized as: a) H ¼ 1=2 for a random walk (Bm) process, i.e., without long-memory behavior, b) H > 1=2 for a persistent process (fBm) and long-memory effect, and c) H < 1=2 for an anti-persistent process related to the short-term memory. The fBm described by a persistent H is a centered Gaussian pro- that has the covariance function [80] given by: The closer the value of H to 1, the higher is the probability for the next step to be positive if the last one was also positive. Several different techniques can be applied to estimate the H exponent of a TS. One of the first methodologies is the rescaled range analysis (R/S), [81, 82] . An approach based on the Fourier analysis, implemented by means of the FFT algorithm, was applied in [83, 84] . The DFA method has also been employed in several works [85, 86] . This technique avoids the false detection of correlation of self-similarities and can be used in the evaluation of nonstationary series, as well as on the detection of long-memory processes. The DFA involves successive steps for calculating the H exponent. Let us consider a stochastic prices TS p j ð Þ; j ¼ 1; . . . ; N, where N 2 N represents the number of observations of the TS. The DFA algorithm involves the three steps: i) Compute the TS mean, , where P n j ð Þ is the ordinary least squares method for removing trend and is subtracted from P j ð Þ, and iii) Repeat the process so that the slope of the straight line relating log F n ð Þ ð Þversus log n ð Þ provides the scaling H exponent. Assuming that the k th order auto-covariance given by: is the k th order autocorrelation, the autocorrelation function q can be determined as: As stated by Peters [87] , there is a relation between H and q described by: Section 3 shows the H exponents determined for the six cryptocurrencies. The Hurst and Lyapunov exponents are important indices for characterizing non-linear and chaotic systems. The H index measures the irregularity of the TS, that is, captures the rate of chaos. On the other hand, the exponent k indicates how the presence of chaos conditions influences the prediction of the future. The value of the k reflects the sensitive dependence on the initial conditions by measuring the exponential divergence of adjacent orbits. Therefore, the evaluation of how distinct trajectories, with nearby initial conditions, diverge, is related to the expansion or contraction of directions in the phase space [88] . Since the values of H and d A of a given TS obey the formula d A ¼ 2 À H, a relation between the Hurst and Lyapunov exponents can be estimated from the global dimension d G . This relationship is used to find the neighboring points in the TS and must be at least equal to or greater than 2d A , and one can write [89, 90] : For the calculation of k, its local dimension D must be determined. The value of D is related to the dimension of Jacobian matrices as stated by Bryant et al. [89] . The values of D must be not much greater than d A and, therefore, an option is to choose the next integer value. Also, if we have D equal to d G then both conditions can be satisfied. The system dynamics maps a D-sphere of states into a Dellipsoid. Consequently, when chaotic motion emerges, some kind of complex dynamics is present. The instabilities and stabilities are associated with the directions where stretching and contraction occur [91] . If we consider e i t ð Þ ¼ e 0 b kt , where e i t ð Þ represents the deformed hyper-volume at time t and b is a given basis, then the Lyapunov exponents can be obtined as [88, 91] , Expression (11) gives a common algorithm for obtaining the Lyapunov spectrum of a system with known equations of motion. For the case of a TS, the Lyapunov exponent can be estimated through the algorithm proposed by Wolf et al. [65] where M and t k À t kÀ1 ¼ D represent the total number of replacement steps and the time step, respectively. The albebraic signs of k i provide information about the system's dynamics and recognize chaotic motion since a positive value indicates that the system is chaotic. Furthermore, the Lyapunov exponent can indicate how far future forecasting can be tried in a TS [92] . Bearing this fact in mind, we applied such technique to calculate the k for the cryptocurrencies price TS. A rolling sample calculation is applied to obtain the values of H and d A , by considering a movable window with fixed length of n ¼ 100 samples, that is, by starting at the first one-hundred observations and rolling until the last group of hundred samples. Indeed, the same methodology is described in SubSection 2.2 to obtain H. With this approach, H is calculated along time and it is possible to analyze the cryptocurrencies behavior from July/2016 to March/2019. Similarly to H, the properties of d A [93, 94] are also related to the memory processes and can be summarized as: (a) 1 < d A 6 2, (b) d A ¼ 3=2, for a random walk (Bm) indicates that the TS does not have a long-memory process and local anti-correlations, (c) d A < 3=2, indicates a persistence process (long-memory or correlated), corresponding to fBm, (d) d A > 3=2, indicates an anti-persistent process (short-term memory, anti-correlated). The index d A is calculated through the Hall-Wood (HW) and Robust Genton (RG) estimators [95] described in the follow-up. The HW estimator [96] is a box-counting algorithm. Let us have a scale l ¼ l=n, where l ¼ 1; 2; 3; . . . ; n. The area of the boxes covers where the operator n=l ½ calculates the integer part of the argument. The HW estimator is given by where L P 2; s l ¼ log l=n ð Þ and s ¼ 1=L ð Þ P L i¼1 s i . If we adopt L ¼ 2, as pointed by Hall-Wood to avoid bias, then we obtain The RG [97] is based on the moments estimator of scale. However, the scheme is not robust and, therefore, the algorithm developed by Genton is adopted. The calculation yieldsV and we obtain the RG estimator as with L P 2; s l ¼ log l=n ð Þand s ¼ 1=L ð Þ P L i¼1 s i . Using L ¼ 2 to mitigate bias, one obtains The ARIMA and ARFIMA prediction results are depicted in Figs. 2 and 3. From Table 1 it is possible to note that an identical order (d ¼ 0:5) is calculated by the ARFIMA for the BTC, ETH, LTC and XMR cryptocurrencies. The range of values higher than d ¼ 0 indicate that these cryptocurrencies have a dynamic fBm. The XRP cryptocurrency also presents a long-memory effect. On the other hand, the ETC price series reveals a short-term memory process. However, other fractional tools, such as the Hurst (H), Fractal Dimension (d A ) and Lyapunov (k) exponents are necessary to determine the dominant process for each price series. Also, Figs. A.7 and A.8 suggest that the ARFIMA has a better fitting for all the cryptocurrencies, since the presence of the residual reveals smaller correlations for the ACF and PACF. Additionally, Figs. 2 and 3 points out smaller confidence intervals values for the ARFIMA that can lead a more precise prediction measurement. The Hurst, Fractal Dimension and Lyapunov exponents, as well as the prediction horizon (1=k) for the cryptocurrencies, are listed in Table 2 . The Hurst exponent and fractal dimension shown in Table 2 and Fig. 4 suggest persistence (long-memory) solely for the BTC, since for the other virtual coins an anti-persistence phenomenon is observed, and are related to a short-memory process. The local and global dimensions, D and d G , with values equal to 3, were identified for all cryptocurrencies, as a way of satisfying the Hurst, Fractal Dimension and Lyapunov conditions. One can note from Table 2 and Fig. 5 , that k is positive for the six cryptocurrencies pointing to a chaotic dynamics. The numerical value 1 k can be viewed as information about the predictability of the future TS price based on its past. The smallest value of k occurs for the BTC and, consequently, the predictability for the BTC is limited to 1 k ¼ 3:020 day. Nonetheless, this value corresponds to the higher horizon prediction among the six cryptocurrencies. The lowest predictability occurs for the LTC with a value of 1 k ¼ 1:769 day. The dynamics in time of the H, d A HW and d A RG for the six criptocurrencies are depicted in Fig. 6 . Fig. 6a shows that, for the most part of the period, the BTC behaves as persistent (long-memory effect) when the indices H, d A HW and d A RG are observed. One can also note from Fig. 6 that the other five cryptocurrencies alternate between persistent and anti-persitent behavior. Nonetheless, a predominant anti-persistent (short-memory effect) phenomenon can be observed (Fig. 6b, d and e) for the ETH, XRP and LTC, respectively. The ETC and XMR (Figs. 6c and f) also fluctuate between persistent and anti-persistent behavior, but it seems clear that they remain preferably with anti-persistence. The results achieved for the BTC are consistent with the one mentioned in the literature [41, 46, [48] [49] [50] . We highlighted that despite the classical fixed window technique showing evidence for an specific effect for a TS (for example, a long-memory effect), the dynamic behavior of an TS can significantly bias its analysis. For this reason, it is possible that the volatility of an certain asset can eventually impair the efficiency of the prediction, when only the aforesaid technique is used. In our study, besides the adoption of a fixed window, we also consider the rolling window approach, that is, we take into account its important dynamic behavior during the analysis, mitigating possible bias phenomena. Fig. 6 shows the time evolution of the three indices H, d A HW and d A RG . Differently of our findings, previous studies [48, 41, 50, 49] point to a persistent behavior for the cryptocurrencies. Nonetheless, such papers do not employ rolling windows and important information may have been overlooked. Despite the price volatility revealed by the oscillation between persistence and anti-persistence along the whole period, it is possible to observe that the measures for H and d A point to similar behaviors as listed in Table 2 , evincing the consistency of the results. This study employed a variety of mathematical tools for the analysis of six significant cryptocurrencies. Its main contribution is related to the use of fractional and fractal mathematical tools as an instrument that can help market agents and investors to more clearly assess the cryptocurrencies price dynamics and, thus, guide investment decisions more assertively while mitigating risks. Classical and fractional integration were explored by means of the ARIMA and ARFIMA processes. ARFIMA performed better than the ARIMA for all cryptocurrencies, since the residual values revealed smaller correlations for the ACF and PACF. Moreover, the smaller confidence interval values for the ARFIMA indicated a more precise prediction measurement. The DFA method and the Wolf algorithm were used for obtaining the Hurst index and the Lyapunov exponent, respectively. The fractal dimension was computed by means of HW and RG estimators. The Hurst exponent and fractal dimension versus time were calculated using sliding windows of constant width, that is, the so-called ''rolling sample approach". The BTC was the only cryptocurrency that presented more consistent long-memory behavior and the smallest value of the Lyapunov exponent. The LTC exhibited the lowest predictable horizon compared to the other cryptocurrencies, pointing to a chaotic behavior and presenting the highest Lyapunov exponent. The ETH and the XMR presented values of H near to the random walk phenomenon. Nevertheless, they behaved mainly as an anti-persistent process showing a short memory effect. With exception of the BTC, the other five cryptocurrencies TS are mean reverting, showing a lower predictability than the BTC, revealing a behavior that was verified to be persistent. Following the results, a future research topic is the multivariate analysis about the influence or the price transmission between the altcoins and the BTC. This article does not contain any studies with human or animal subjects. whereP t represents the forecast value and e t ¼ P t ÀP t , which is the error value of t. The results for the six cryptocurrencies are summarized in Table A Blockchain technology in the energy sector: A systematic review of challenges and opportunities A brief overview on China and cryptocurrency An empirical investigation on the adoption of cryptocurrencies among the people of mainland China Taming the blockchain beast? Regulatory implications for the cryptocurrency market The effects of US-China trade war and Trumponomics, MPRA Paper 93682. Germany: University Library of Munich Mastering Bitcoin: unlocking digital cryptocurrencies Understanding bitcoin: cryptography, engineering and economics The Oxford handbook of innovation Bitcoins as an investment or speculative vehicle? A first look The inefficiency of bitcoin Understanding bitcoins: facts and questions Bitcoin and beyond: A technical survey on decentralized digital currencies The future of cryptocurrencies: Bitcoin and beyond Sok: Research perspectives and challenges for bitcoin and cryptocurrencies Multiscale characteristics of the emerging global cryptocurrency market Exploring the dynamic relationships between cryptocurrencies and other financial assets Fractional calculus view of complexity: A tutorial Recent history of fractional calculus Fractional order calculus: historical apologia, basic concepts and some applications The fractional calculus: theory and applications of differentiation and integration to arbitrary order Fractional calculus: an introduction for physicists The fractional features of a harmonic oscillator with position-dependent mass A new adaptive synchronization and hyperchaos control of a biological snap oscillator A new fractional HRSV model and its optimal control: A non-singular operator approach Can fractional calculus help improve tumor growth models? The ACF and PACF residuals for the ARFIMA for the six cryptocurrencies On multistep tumor growth models of fractional variable-order The role of fractional calculus in modeling biological phenomena: A review Numerical solutions for time-fractional cancer invasion system with nonlocal diffusion A hybrid fractional optimal control for a novel coronavirus (2019-ncov) mathematical model A delay fractional order model for the co-infection of malaria and HIV/AIDS On history of mathematical economics: Application of fractional calculus Fractional econophysics: Market price dynamics with memory effects The recovery of global stock markets indices after impacts due to pandemics Fractional dynamics and pseudophase space of country economic processes Analysis of a fractional model of the Ambartsumian equation Dynamic stability analysis of fractional order leaky integrator echo state neural networks Non-commensurate fractional linear systems: New results Cryptocurrency forecasting with deep learning chaotic neural networks Bitcoin futures-what use are they? Can volume predict bitcoin returns and volatility? A quantiles-based approach Structural breaks and double long memory of cryptocurrency prices: A comparative analysis from bitcoin and ethereum Bitcoin and the day-of-the-week effect The day of the week effect in the cryptocurrency market Virtual currency, tangible return: Portfolio diversification with bitcoin Crypto-currency bubbles: an application of the Phillips-Shi-Yu (2013) methodology on Mt. Gox bitcoin prices Asymmetric mean reversion of bitcoin price returns The effectiveness of technical trading rules in cryptocurrency markets Modelling crypto-currencies financial time-series Mathematical and statistical methods for actuarial sciences and finance Modelling long memory volatility in the bitcoin market: Evidence of persistence and structural breaks Persistence in the cryptocurrency market Efficient capital markets: a review of theory and empirical work On bitcoin markets (in)efficiency and its evolution The inefficiency of bitcoin revisited: A dynamic approach On the inefficiency of bitcoin Long-and short-term cryptocurrency volatility components: A GARCH-MIDAS analysis Partial chaos suppression in a fractional order macroeconomic model Fractional state space analysis of economic systems Fractional dynamic behavior in ethanol prices series Dynamics of commodities prices: integer and fractional models Time-varying long-term memory in bitcoin market The long-term dependence in stock returns Use and misuse of some Hurst parameter estimators applied to stationary and non-stationary financial time series Fractional dynamic behavior in Forcados oil price series: An application of detrended fluctuation analysis Determining Lyapunov exponents from a time series Chaos and time-series analysis Cryptocurrencies market efficiency ranking: Not so straightforward Chaos, randomness and multi-fractality in bitcoin market Does bitcoin exhibit the same asymmetric multifractal cross-correlations with crude oil, gold and DJIA as the Euro, Great British Pound and Yen? Measuring capital market efficiency: long-term memory, fractal dimension and approximate entropy currencies and market efficiency Multifractal behavior of price and volume changes in the cryptocurrency market Time series analysis: forecasting & control A study in the analysis of stationary time series Long memory processes and fractional integration in econometrics Bootstrap approaches and confidence intervals for stationary and non-stationary long-range dependence processes An introduction to long-memory time series models and fractional differencing Non-gaussian state-space modeling of nonstationary time series Consumption and fractional differencing: Old and new anomalies Fractional calculus and fractional processes with applications to financial economics Statistical methodology for nonperiodic cycles: from the covariance to R/S analysis Robustness of the rescaled range R/S in the measurement of noncyclic long run statistical dependence The use of fast Fourier transform for the estimation of power spectra: A method based on time averaging over short, modified periodograms Reconstructing cloudfree NDVI composites using Fourier analysis of time series Mosaic organization of DNA nucleotides Long memory and sampling frequencies: Evidence in stock index futures markets Chaos and order in the capital markets: A new view of cycles, prices, and market volatility Quantifying chaos with Lyapunov exponents Lyapunov exponents from observed time series Chaotic analysis of gold price index Nonlinear dynamics, chaos, control, and their applications to engineering sciences Nonlinear prediction of chaotic time series Measuring capital market efficiency: Global and local correlations structure How are rescaled range analyses affected by different memory and distributional properties? A Monte Carlo study Estimators of fractal dimension: Assessing the roughness of time series and spatial data On the performance of box-counting estimators of fractal dimension Variogram fitting by generalized least squares using an explicit formula for the covariance structure UNESP-Brazil) in 2012. Ph.D and Master's degree in Mechanical Engineering at University of Campinas (UNICAMP-Brazil) in 2003 and 1996, respectively. Executive MBA in Finance at Fundação Getúlio Vargas (FGV-RJ-Brazil)/Ohio University (OHIO-USA) in 2002. Bachelor degree in Physics at São Paulo State University (UNESP-Brazil) in 1993. He works on Mathematical Methods Applied to Science and Engineering. Main research interests: mathematical modeling & simulation in mechanical, financial, economic and biomedical sciences, fractional calculus, systems dynamics, nonlinear phenomena and energy economics The authors wish to acknowledge the FAPESP (São Paulo Research Foundation), grants 2017/13815-3 and 2017/15517-0, for funding support. We also acknowledge Paloma O. Santos for partial contribution in numerical simulations. To measure the accurary of the prediction models, several criteria are adopted for comparison. We consider the Entropy-Theil's