Next Article in Journal
Adaptive Fractional Image Enhancement Algorithm Based on Rough Set and Particle Swarm Optimization
Previous Article in Journal
Generalized Mittag-Leffler Kernel Form Solutions of Free Convection Heat and Mass Transfer Flow of Maxwell Fluid with Newtonian Heating: Prabhakar Fractional Derivative Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Asset Pricing Model Based on Fractional Brownian Motion

1
School of Economics, Peking University, Haidian District, Beijing 100871, China
2
Key Laboratory of Mathematical Economics and Quantitative Finance, Beijing 100871, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Fractal Fract. 2022, 6(2), 99; https://doi.org/10.3390/fractalfract6020099
Submission received: 11 January 2022 / Revised: 2 February 2022 / Accepted: 6 February 2022 / Published: 11 February 2022

Abstract

:
This paper introduces one unique price motion process with fractional Brownian motion. We introduce the imaginary number into the agent’s subjective probability for the reason of convergence; further, the result similar to Ito Lemma is proved. As an application, this result is applied to Merton’s dynamic asset pricing framework. We find that the four order moment of fractional Brownian motion is entered into the agent’s decision-making. The decomposition of variance of economic indexes supports the possibility of the complex number in price movement.

1. Introduction

Many continuous-time stochastic models have been established after researchers realized the role of Brownian motion in describing prices. Merton (1973) [1] and Merton (1975) [2] initiated a new research framework of the asset pricing model under continuous-time and discussed the portfolio theory under this framework. Davis and Norman (1990) introduced the model with transaction cost in [3]. Zhou 1997 introduced the jump process into the option pricing model based on continuous time and used the Poisson process and normal distribution to describe the discontinuities of company value in [4]. Liu 2007 introduced the form of the solution when the rate of return is a quadratic diffusion process in [5].
Many studies consider using fractional Brownian motion, which was first studied by Kolmogorov 1940 [6], as a substitute for Brownian motion.
Definition 1.
For Hurst parameter H ( 0 , 1 ) , fractional Brownian motion (FBM) B H : = { B H ( t ) , t R } , is the Gaussian process. We take B H ( 0 ) = 0 . Its mean for all t is: E [ B H ( t ) ] = 0 and its covariance is:
E [ B H ( t ) B H ( s ) ] = 1 2 { | t | 2 H + | s | 2 H | t s | 2 H } E [ ( B H ( t ) B H ( s ) ) 2 ] = | t s | 2 H .
For H = 1 / 2 , B 1 2 is standard Brownian motion.
The use of fractional Brownian motion to characterise asset prices has attracted the attention of many researchers. Mandelbrot and Van Ness 1968 [7] defined the integral transformation from Brownian motion to fractional Brownian motion. Peters 1994 [8] introduced the fractal efficient market hypothesis to challenge the efficient market hypothesis. Rogers 1997 [9] and Cheridito 2001 [10] introduced the regularizing fractional Brownian motion with long memory but no-arbitrage. Cheridito 2001 [10] showed that if there is an upper limit on trading frequency, the arbitrage example in Peters 1994 [9] may be excluded. Sottinen 2001 discussed an approximation that converges weakly to fractional Brownian motion, then discussed a basic market model under this approximation, and discussed the existence of arbitrage opportunities in [11]. Guasoni 2006 discussed how arbitrage opportunities are eliminated with transaction costs in [12]. Hu and Øksendal 2003 [13] introduced the Ito random integral with different fractional Brownian motion between ( 1 / 2 , 1 ) and showed that there is no arbitrage at this time. Elliott and Van Der Hoek 2003 generalized this result and proposed a framework that combines fractional Brownian motions of different orders between ( 0 , 1 ) in the same process and discussed the option pricing formula under this framework in [14]. Xiao et al. 2010 [15] introduced a currency option pricing model combining the jump process and fractional Brownian motion. Xiao and Yu 2019 [16] developed a class of the Vasicek model with fractional Brownian motion as the random source. Of particular attention is the result made by [17], which is very relevant to us. The core difference is that we use a dynamic optimization approach, the solution process and the final result have more economic implications and we also focus on the possibility of the existence of imaginary numbers. In fact, one of the important reasons why fractional Brownian motion is difficult to shine in economic research is the lack of a matching dynamic optimization framework. A good dynamic optimization framework can help economic researchers understand the agent’s decision-making and the overall movement of the economy. Our main contribution is to give a dynamic optimisation framework that remains feasible under fractional Brownian motion.
In recent years, most of the research on fractional Brownian motion about the economic has focused on option and bond pricing. Miao and Yang 2015 [18] priced convertible bonds involving mixed fractional Brownian motion. Shokrollahi and Kılıçman 2015 [19] priced currency option, when the value of foreign currency option follows the mixed fractional Brownian motion with jumps. Rao 2016 [20] priced the geometric Asian power option. Shokrollahi et al. 2016 priced European options and currency options by time changed mixed fractional Brownian motion in [21]. Zhang et al. 2020 proposed a fuzzy mixed fractional Brownian motion model with jumps in [22]. The article on option pricing is of high applied value. However, the purely pricing equations hardly bring us clearer economic intuition because they are not based on individual decisions. Our core work is to give a solution to the exogenous process as a fractional Brownian motion individual lifetime optimisation solution. A solution based on individual decision making can lead to clearer economic implications.
This paper is arranged as follows. In Section 2, some real data are decomposed to show the necessity of considering H 1 2 . The decomposition of real data also supports the idea of introducing imaginary numbers. In Section 3, we show the classical results and discuss the rationality of our independence hypothesis for different stochastic processes. In Section 4, we will discuss the Merton type single asset and multi-asset continuous-time consumption asset pricing models corresponding to multi fractional Brownian motion process. In Section 5, we summarize our results. The proof in this paper is in the Appendix A, Appendix B and Appendix C.

2. Decomposition of Real Data

We used monthly U. S. consumption data, S & P 500 daily data and tick level data of the Shanghai Composite Index. The monthly U. S. consumption data and S & P 500 daily data are from the CEIC data set, and the Shanghai composite index data is from Resset data. We show descriptive statistics in Table 1.
Among the three data, we find that the variance of the logarithmic growth rate of monthly consumption in the United States has the most obvious nonlinear property with respect to the growth of sampling interval:
It can be seen from Figure 1 that the variance of the logarithmic growth rate of consumption shows an obvious nonlinear growth with respect to the sampling interval, which is contrary to the assumption of random walk.
When an economic stochastic process p satisfies:
d p / p = μ d t + j = 1 σ j d B H j .
When B H j independent of each other, then:
v a r ( l o g ( p ( t + d t ) ) l o g ( p ) ) = j = 1 σ j 2 d t 2 H j .
On the premise that Gaussian processes can be decomposed orthogonally, we find that the variance of multi fractional Brownian motion can be decomposed of different time increments. This formula also implies that the coefficients should be positive when we do a regression of the power exponents for different time intervals using the variance measured at different frequencies. If the coefficients are negative, then it is possible that the presence of imaginary numbers is implied.
Next, we use the empirical method to try to decompose this relationship in Table 2.
From the empirical results, we can see that the coefficients of different power of d t are significant, so we think that the economic data is composed of different fractional Brownian motions. Interestingly, the coefficients of some terms are negative, which means that some σ j 2 in (3) are negative, which is consistent with the introduction of the complex number.
We have used this approach for a wider range of data, seeing the Appendix A, Appendix B and Appendix C.
The real economic time series characteristics we find are derived from autocorrelation. In the past, most studies have been anchored on the Hurst parameter to investigate autocorrelation in economic time series. If an economic time series is a stationary process, its Hurst parameter is 1 / 2 . For the definition of the Hurst parameter, see Definition 1. Eom et al. 2007 [23] demonstrated the role of the Hurst parameter in measuring market efficiency. Many real data support the existence of chaos in real economic laws. Brock and Sayers 1988 [24] tested the U. S. macroeconomic data and confirmed the existence of two-dimensional chaos. Ashley and Patterson 1989 [25] found that nonlinear dynamics is an important feature of the macroeconomic index. Karuppiah and Los 2005 [26] pointed out that the Hurst index of foreign exchange is in the range of ( 0 , 0.5 ) . Wang et al. 2009 found that the Hurst parameters of stocks in the U.S. market obey normal distribution in [27]. Hassani and Thomakos 2010 introduced the successful application of single spectrum analysis in financial time series, which confirmed the non-stationarity of financial time series in [28]. Onali and Goddard 2011 analyzed the European market and found that the Italian stock market has a long memory in [29]. Garnier and Sølna 2017 introduces the option pricing model with volatility correlation in [30]. Caraiani 2012 ([31]) found evidence of the existence of chaos in the central and Eastern European stock markets. Bal and Rath 2015 ([32]) found that there is a nonlinear Granger causality between crude oil and exchange rate in India. Jin 2016 found that during the financial crisis in 2008, the return of the Asian market showed long memory, and the Hurst parameter sometimes was above 0.5 , sometimes below 0.5 in [33]. Bariviera 2017 found that the Hurst parameter results of the bitcoin market measured by different methods are different, but overall, they are greater than 0.5 in [34]. Wei 2018 ([35]) found that among 456 cryptocurrencies, the increase of liquidity will reduce the predictability of prices and found that the mean value of Hurst parameters of cryptocurrencies with stronger predictability is less than 0.5 . Takaishi 2020 thought that the generalized Hurst parameter of logarithmic fluctuation increment of bitcoin price is less than 0.5 in [36].
Of course, some scholars are against it. For example, Shintani and Linton (2003) [37] believed that there may be no chaos in the macroeconomy, or it is difficult to identify because of the small amount of data.
From the past research, we can get the chaotic system from the general equilibrium theory. Brock and Hommes 1997 [38] discussed the local steady-state instability in a heterogeneous belief model with agent learning. Hommes et al. 2005 [39] discussed the chaos of asset prices when heterogeneous traders are allowed to update their beliefs asynchronously. Chiarella et al. 2006 [40] discussed the survival of traders in different chaotic sequences in a heterogeneous belief pricing model that allows traders to update their beliefs. Sommervoll et al. 2010 [41] discussed in a complex dynamic system that if the house is set as the mortgage of the loan, the price will fluctuate sharply.

3. Classical Model and Orthogonalization

In Section 3.1, we recall a classical continuous-time asset pricing model from [2] as a basis for further exploration. In order to give a basis for the assumption that fractional Brownian motions in asset price movements are independent of each other, in Section 3.2 we give an orthogonalization of stochastic processes. The existence of an orthogonalization implies that although stochastic processes in reality may be correlated with each other, we can transform the stochastic processes in reality into stochastic processes that are independent of each other.

3.1. Classical Model and Result

In the classical model, it is assumed that the price of risky asset p satisfies the following conditions: d p / p = r d t + σ d B , and B satisfies the standard Brownian motion: d B N ( 0 , d t ) , E t ( d B d B ) = d t . r is the expected return rate of risk assets, σ is the risk exposure. Risk-free asset price p f satisfies d p f / p f = r f d t , where r f is the return rate of risk-free assets. w is the process of wealth movement. x is the weight of risk assets. c is the consumption. Individual investors need to solve the optimization problem as follows:
m a x c , w E 0 + e ρ t ( c α α ) d t s . t . d w = ( w ( r f ( 1 x ) + r x ) c ) d t + x w σ d B w ( 0 ) = w 0
Finally, it can be transformed into a system of differential Equations (where V is the Bellman Equation):
c α 1 V = 0 V w ( r r f ) + V x ( w σ ) 2 = 0 .
The final solution is obtained: c = 1 1 α ρ α r f α 2 ( 1 α ) r r f σ 2 w , x = r r f ( 1 α ) σ 2 .

3.2. Orthogonalization of Stochastic Processes

In this paper, in order to make the solution easier, some stochastic processes are assumed to be independent. The following describes the orthogonalization process of stochastic processes support our assumption.
Proposition 1.
For a group of stochastic processes: B i , i { 1 , 2 , , I } , which have d B i ( t ) N ( μ B i ( t ) , σ B i ( t ) 2 d t α ) , and c o v ( d B i ( t ) , d B j ( t ) ) = σ i , j , t d t α . Then there exist z i , i { 1 , 2 , , I } which have d z i ( t ) N ( μ z i ( t ) , σ z i ( t ) 2 d t α ) , are independent of each other. d w i , i { 1 , 2 , , I } can be expressed linearly by d z i , i { 1 , 2 , , I } .
Proof. 
We only need to use the Schmidt orthogonalization method and make the coefficient time-varying. Let
d z 1 = d B 1 ,
d z 2 = d B 2 c o v ( d B 2 , d z 1 ) v a r ( d z 1 ) d z 1 = d B 2 c o v ( d B 2 , d B 1 ) v a r ( d B 1 ) d B 1 = d B 2 σ 1 , 2 , t σ B 1 ( t ) 2 d B 1 ,
d z 3 = d B 3 c o v ( d B 3 , d z 1 ) v a r ( d z 1 ) d z 1 c o v ( d B 3 , d z 2 ) v a r ( d z 2 ) d z 2 .
It is easy to verify: d B i , i { 1 , 2 , , I } and d z i , i { 1 , 2 , , I } are equivalent groups of random variables. □

4. Modelling of Fractional Brownian Motion

In this section, we consider a multi fractional Brownian motion process for price p ( μ is the expected return and σ is the risk exposure)
d p = μ ( p ) d t + j = 1 σ j d B H j
where H j [ 1 4 , 1 ) . Our current solution scheme can only cope with the case H j [ 1 4 , 1 ) . If H j [ 1 6 , 1 4 ) , we need to deal with Taylor expansions of order 6. Our current proof for the Ito Lemma type lemma is already very difficult. On the other hand, it is difficult to extend the range of values of H j further to obtain more specific economic implications. Therefore, we choose H j [ 1 4 , 1 ) .
If B H j are independent of each other. We can talk about more concise issues. Since we plan to prove a framework similar to Ito lemma, we only need to focus on H j = 1 2 and H j = 1 4 when a function f ( x ) is expanded by Taylor. The rest will be o ( d t ) , which will affect the expected return of the price but will not affect the higher-order terms in Taylor expansion.
To sum up, we will simplify the process of (4) and further discuss the following one types of price movements:
d p / p = r d t + σ 0 d B 1 2 + ( σ 1 + i σ 1 ) d B 1 4 .
When the fractional Brownian motion with Hurst exponent less than 0.5 is considered directly, its variance is a low order infinitesimal of d t when the time interval d t tends to zero. Therefore, when agents make decisions, they will get some strange asset allocation because of the fear of violent fluctuations in small time intervals (This can be proved by considering the way similar to Ito lemma. After retaining the d t brought by fractional Brownian motion variance, the Bellman equation will change d t into o ( d t ) . ). Therefore, we use the method of introducing imaginary number to let the agent ignore the influence of fractional Brownian motion on the square difference. This scheme is also the best one we think of at present, introducing fractional Brownian motion while maintaining mathematical rigour under the consumer utility framework. The introduction of imaginary numbers is a common scheme in quantum mechanics, such as seeing [42] p. 20. From another point of view, the introduction of an imaginary number shows that the agent ignores the extreme fluctuation caused by fractional Brownian motion in a very short time and has subjectivity.

4.1. Solution of Single Risk Asset Model

In this subsection, we assumed that the price of risky asset p satisfies the following conditions: d p / p = r d t + σ 0 d z 0 + ( σ 1 + i σ 1 ) d z 1 , and z 0 satisfies the standard Brownian motion: d z 0 N ( 0 , d t ) , E t ( d z 0 d z 0 ) = d t . r is the expected return rate of risk assets, σ is the risk exposure. z 0 , z 1 are independent of each other. i is the imaginary root: i 2 = 1 . d z 1 is subject to N ( 0 , d t ) . In fact, it may be problematic to assume that the expectation of fractional Brownian motion increment is zero. However, in reality, traders cannot observe the value of fractional Brownian motion. Therefore, it is reasonable to assume that the increment of fractional Brownian motion in traders’ subjective probability is zero. Risk-free asset price p f satisfies d p f / p f = r f d t , where r f is the return rate of risk-free assets. w is the process of wealth movement. x is the weight of risk assets. c is the consumption. Individual investors need to solve the optimization problem as follows:
m a x E R e ( 0 + e ρ t ( c α α ) d t s . t . d w = ( w ( r f ( 1 x ) + r x ) c ) d t + w x ( σ 0 d z 0 + ( σ 1 + i σ 1 ) d z 1 ) w ( 0 ) = w 0 .
The method of Bellman equation is used to solve the problem:
V ( t , w ) = m a x x , c E t R e t t + d t e ρ t ( c α α ) d t + V ( t + d t , w + d w ) .
It should be noted that Ito lemma (seeing p. 93 in [43]) can not be used directly here, so the following lemma is used:
Lemma 1.
If the random process w satisfies:
d w = f ( t , w ) d t + g 1 ( t , w ) d z 0 + ( g 2 ( t , w ) + i g 2 ( t , w ) ) d z 1 .
Furthermore, z 0 satisfies the standard Brownian motion, there will be d z 0 N ( 0 , d t ) . z 0 , z 1 are independent. i is the imaginary root: i 2 = 1 . z 0 , z 1 obey N ( 0 , d t ) . For u = F ( t , w ) , if F is resolved in a neighbourhood of ( t , w ) , then u satisfies:
d u = F t + f F w + 1 2 g 1 2 F w w 1 2 g 2 4 F w w w w d t + 2 i g 2 2 d t + F w ( g 1 d z 0 + g 2 d z 1 + i g 2 d z 2 ) .
Next, we solve the optimization problem. Bellman equation is,
ρ V = m a x c , x c α α + ( w ( r f ( 1 x ) + r x ) c ) V w + ( w x σ 0 ) 2 2 V w w ( w x σ 1 ) 4 2 V w w w w .
It should be noted that the Bellman equation is a quartic equation with respect to x. According to our later derivation, we can find that the coefficient of x 4 may be positive, which means that there is no suitable solution less than infinity for this equation. For the properties of the fourth derivative of utility function, see [44]. Therefore, when the fourth derivative of the indirect utility function is less than 0, only when there are short selling and borrowing constraints or investors only rely on the first-order conditions for decision-making can there be a finite solution. Next, we will show possible solutions with loan constraints.
The derivation of c yields c α 1 V w = 0 , and the derivative of x can be obtained as follows:
V w w ( r r f ) + V w w x ( w σ 0 ) 2 2 V w w w w x 3 ( w σ 1 ) 4 = 0 .
The guess solution is as follows: V = A w α . (8) can be reduced to:
( r r f ) + ( α 1 ) σ 0 2 w 2 ( α 1 ) ( α 2 ) ( α 3 ) σ 1 4 w 3 = 0 .
Let q 1 = ( α 1 ) σ 0 2 2 ( α 1 ) ( α 2 ) ( α 3 ) σ 1 4 , q 2 = ( r r f ) 2 ( α 1 ) ( α 2 ) ( α 3 ) σ 1 4 , Then (9) can be reduced to x 3 + q 1 x + q 2 = 0 . There must be a real number solution by using the method of Cardano formula (In solving practical problems, we will explore the corresponding properties of the three points and select the best one. This is just an example.):
x * = q 2 2 + q 2 2 2 + q 1 3 3 3 + q 2 2 q 2 2 2 + q 1 3 3 3 .
Put c α 1 = V w and V = A w α in Bellman equation:
A = 1 α 1 1 α ρ α ( r f ( 1 x * ) + r x * ) 1 2 ( x * σ 0 ) 2 α ( α 1 ) + 1 2 ( x * σ 1 ) 4 α ( α 1 ) ( α 2 ) ( α 3 ) α 1 .
Finally, we have c = ( A α ) 1 / ( α 1 ) w . It can be seen that the central result of this subsection is the introduction of σ 1 4 into the decision result, as opposed to the results of Section 3.1. Due to the existence of σ 1 4 , the decision equation for x has the cubic and σ 1 4 appears in the final solution. This means that one of the components of the fourth order moment of the stochastic process of the price of a risky asset, σ 1 4 , affects the proportion of the agent’s investment in the risky asset. The ratio of the agent’s consumption decision to the asset is also only one additional term corresponding to σ 1 4 compared to the result in Section 3.1.
In fact many studies have provided a basis for the importance of higher order moments. Peiro 1999 [45] found that the asymmetry of returns can not be rejected when using other distributions of non-lognormal distribution to characterize returns. Patton 2004 [46] discussed the impact of asset return skewness and asymmetric dependence on the portfolio and used the information of the fourth moment to analyze in the process of analysis. Aït-sahali and Brandt 2001 [47] mentioned the importance of skewness when introducing the selection of investment decision variables corresponding to different utility functions. In relation, Kimball 1993 [48] studied the properties of higher derivative of utility functions.

4.2. Comparison of Results and Numerical Simulations

In this Subsection, we will systematically compare the relationship and differences between our work in Section 4.1 and Merton 1975 [2], Miao and Yang 2015 [18], Shokrollahi and Kılıçman 2015 [19], Rao 2016 [20], Shokrollahi et al. 2016 [21], Zhang et al. [22]. In order to be able to illustrate the value of our work, we conducted numerical simulations using BP500 data and observed whether our decision options improved lifetime utility. Based on the conventions of general economic research, we let 1 α = 3 . Let ρ = 0.08 and r f = 0.03 , see Kyle and Xiong 2001 [49]. Let r, σ 0 and σ 1 be the moment estimates for the past 1000 trading days. After rolling, we computed lifetime utility values under different decisions in Table 3.
An important point about our model is that it has the same microfoundations of individual choice as [2]. The specific form of individual choice allows us to see two important facts. The first aspect is what are the mechanisms that influence individual choice when exogenous shocks occur. The second aspect is the linking of risky asset prices to macro variables such as consumption and savings rates, which can help us to better analyse the impact of macro policies.
Further, our model with its micro-foundations can be extended to a general equilibrium model. Dynamic stochastic general equilibrium models as a central research tool in modern economics allow for the study of economic laws in a broader context. At the meantime, both our model and [2] can be used to study economics under heterogeneous agents. Studying heterogeneous traders can help us to better understand how an economic policy affects different groups of people.
Our central starting point with [2] is the maximisation of investor utility. We therefore consider how our decision making approach and the decision making approach in [2] behave in realistic data.The decision approach in [2] derives a welfare discount of 2.9227 . Our decision approach derives a welfare discount of 2.1076 . The welfare improvement is 0.8151 . Realistic data suggest that our decision scheme does improve relative to the classical work of Merton 1975 ([2]). The reason for the comparison with [2] is that the main aim of our work is to propose a new framework for solving economics problems based on fractional Brownian motion and the example we give is the individual choice problem in [2].

4.3. Multidimensional Situation

Here we assume a variety of risky asset prices of p l satisfies: d p l / p l = r l d t + k = 1 K σ 0 , l , k d z 0 , k + j = 1 J σ 1 , l , j ( 1 + i ) d z 1 , j . d z 0 , k satisfies the standard Brownian motion: d z 0 , k N ( 0 , d t ) , E t ( d z 0 , k d z 0 , k ) = d t . r is the expected return rate of risk assets, σ is the risk exposure. z 0 , k , z 1 , j are independent of each other. i is the imaginary root: i 2 = 1 . d z 1 , j are subject to N ( 0 , d t ) . Risk-free asset price p f satisfies d p f / p f = r f d t , where r f is the return rate of risk-free assets. w is the process of wealth movement. x is the weight of risk assets. c is the consumption. Individual investors need to solve the optimization problem as follows:
m a x E R e 0 + e ρ t ( c α α ) d t s . t . d p l / p l = r l d t + k = 1 K σ 0 , l , k d z 0 , k + j = 1 J σ 1 , l , j ( 1 + i ) d z 1 , j d w = w x d log ( p ) + [ w r f ( 1 x 1 ) c ] d t w ( 0 ) = w 0 .
The method of Bellman equation is used to solve the problem:
V ( t , w ) = m a x x , c E t R e t t + d t e ρ t ( c α α ) d t + V ( t + d t , w + d w ) .
It should be noted that Ito lemma can not be used directly here, so the following lemma is used:
Lemma 2.
If the random process w l satisfies:
d w l = f ( t , w l ) d t + k = 1 K g 1 , l , k ( t , w ) d z 0 , k + j = 1 J g 2 , l , j ( t , w ) ( 1 + i ) d z 1 , j .
Furthermore, d z 0 , k satisfies the standard Brownian motion, there will be
d z 0 , k N ( 0 , d t ) , E t ( d z 0 , k d z 0 , k ) = d t . z 0 , k , z 1 , j are independent. i is the imaginary root: i 2 = 1 . z 1 , j obey N ( 0 , d t ) . For u = F ( t , w ) , if F is resolved in a neighborhood of ( t , w ) , then u satisfies:
R e E t ( d u ) = ( F t + f F w + 1 2 1 g 1 F w w g 1 1 + F w l 1 w l 2 w l 3 w l 4 ( j = 1 J l 1 = 1 L l 2 = 1 L l 3 = 1 L l 4 = 1 L g 2 , l 1 , j g 2 , l 2 , j g 2 , l 3 , j g 2 , l 4 , j 3 2 j 1 = 1 J j 2 = 1 J l 1 = 1 L l 2 = 1 L l 3 = 1 L l 4 = 1 L g 2 , l 1 , j 1 g 2 , l 2 , j 1 g 2 , l 3 , j 2 g 2 , l 4 , j 2 ) ) d t + o ( d t ) .
Next, we solve the optimization problem. Bellman equation is,
ρ V = m a x c , x [ c α α + ( w ( r f ( 1 x 1 ) + x r ) c ) V w + w 2 x σ 0 x 2 V w w w 4 V w w w w ( 3 2 j 1 = 1 J j 2 = 1 J l 1 = 1 L l 2 = 1 L l 3 = 1 L l 4 = 1 L σ 2 , l 1 , j 1 σ 2 , l 2 , j 1 σ 2 , l 3 , j 2 σ 2 , l 4 , j 2 x l 1 x l 2 x l 3 x l 4 j 1 = 1 J l 1 = 1 L l 2 = 1 L l 3 = 1 L l 4 = 1 L σ 2 , l 1 , j 1 σ 2 , l 2 , j 1 σ 2 , l 3 , j 1 σ 2 , l 4 , j 1 x l 1 x l 2 x l 3 x l 4 ) ] .
The derivation of c yields c α 1 V w = 0 , and the derivative of x can be obtained as a cubic equation system of l (the coefficients are all known constants), and the solution can be obtained when there are specific values. After assuming that the solution is obtained, the value of A can be solved. The guess solution is as follows: V = A w α .
Put c α 1 = V w and V = A w α in Bellman equation:
A = 1 α 1 1 α ρ α ( w ( r f ( 1 x * 1 ) + x * r ) c ) 1 2 x * σ 0 x * α ( α 1 ) α ( α 1 ) ( α 2 ) ( α 3 ) 3 2 j 1 = 1 J j 2 = 1 J l 1 = 1 L l 2 = 1 L l 3 = 1 L l 4 = 1 L σ 2 , l 1 , j 1 σ 2 , l 2 , j 1 σ 2 , l 3 , j 2 σ 2 , l 4 , j 2 x l 1 * x l 2 * x l 3 * x l 4 * j 1 = 1 J l 1 = 1 L l 2 = 1 L l 3 = 1 L l 4 = 1 L σ 2 , l 1 , j 1 σ 2 , l 2 , j 1 σ 2 , l 3 , j 1 σ 2 , l 4 , j 1 x l 1 * x l 2 * x l 3 * x l 4 * α 1 .
Finally, we have c = ( A α ) 1 / ( α 1 ) w . In contrast to Section 4.1, this subsection shows the results for multi-asset decisions. We find that the generalisation of the single-asset to multi-asset does not affect the form of the solution. Relative to the results of Section 3.1, the final result for the consumption and investment ratios after the model is multi-asset and is added to the fractional Brownian motion only increases the fourth order moments of the stochastic process for the price of the risky asset.
Aït-Sahalia and Lo 1998 [50] found negative skewness and excess kurtosis for asset returns are important features that affects the pricing of options. Figure 7 in [50] illustrated the effect of skewness and kurtosis on term prices. If we use expectation theory to explain price formation, i.e., investors form expectations which lead to decisions and ultimately determine prices, then we can see that skewness and kurtosis do influence investors’ decisions through. This result from classical observations supports the results of our model. Ahmed and Al Mafrachi 2021 [51] found that kurtosis has a predictive effect on returns in the cryptocurrency market. In our model investors’ decisions are influenced by the kurtosis of an asset. If we still believe in the path: kurtosis affects decisions and hence prices, the kurtosis of an asset will have an impact on future returns. Our model provides an easily understood explanation for [51]. A similar observation is [52]. The role of kurtosis on returns strongly supports our results.

5. Conclusions

We brought the stochastic process with fractional Brownian motion random source into the economic research framework of utility maximization with the complex number. Surprisingly, the decomposition of variance of real economic data in Section 2 supported the scheme of introducing the complex number.
In Section 3, if motions in the random source group are Brownian motion, we proved that the random sources can be set to be independent of each other. In this section, we also reviewed Merton’s classical continuous-time asset pricing model.
Furthermore, we introduced a stochastic process in Section 4 by introducing complex numbers. It had been found that the introduction of fractional Brownian motion affected the decision-making of the agent but did not destroy the original structure. Compared to the results of Section 3.1, the form of the return and variance of the risky asset return did not change in the results of Section 4. The fourth order moments of the stochastic process for the price of the risky asset appeared in the decision result of Section 4.

Author Contributions

Writing—original draft preparation, Y.Y.; writing—review and editing, Y.W. All authors have read and agreed to the published version of the manuscript.

Funding

The APC was funded by he Key Laboratory of Mathematical Economics and Quantitative Finance (Peking University), Ministry of Education of China.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

CEIC Data, an ISI Emerging Markets Group Company. 1976-2020. Global Database, China Premium Database. https://insights.ceicdata.com/login (accessed on 1 January 2022); RESSET DATA.

Acknowledgments

We would like to acknowledge the editor and reviewers for their review and guidance.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Robustness Testing of Data Decomposition Results

Table A1. Robustness Testing of Data Decomposition Results.
Table A1. Robustness Testing of Data Decomposition Results.
Variablesy1y2y3y4y5y6
d t 2 3.95 × 10 5 −1.12 × 10 5 −3.79 × 10 5 1.35 × 10 5 ***1.89 × 10 6 ***1.53 × 10 6 ***
2.74 × 10 5 9.54 × 10 6 2.68 × 10 5 1.18 × 10 6 2.28 × 10 7 2.38 × 10 7
d t 3 2 −0.00203 *0.000470.00146−0.000441 ***−6.98 × 10 5 ***−4.83 × 10 5 ***
0.001030.0003610.001014.44 × 10 5 8.61 × 10 6 8.98 × 10 6
d t 0.0569 **−0.0118−0.03280.00952 ***0.00179 ***0.00117 ***
0.02420.008440.02370.001040.0002010.00021
d t 3 4 −0.184 **0.03730.0972−0.0271 ***−0.00490 ***−0.00307 ***
0.07540.02630.07390.003240.0006280.000655
d t 1 2 0.186 **−0.037−0.08930.0250 ***0.00470 ***0.00299 ***
0.07510.02620.07360.003230.0006250.000652
Constant−0.0540 **0.0166 *0.0355−0.00711 ***−0.00138 ***−0.000865 ***
0.02530.008810.02480.001090.0002110.000219
Observations120120120120120120
R-squared0.9640.0580.089111
y1, y2 and y3 are monthly data from 2005/01 to 2021/09 of Retail sales of consumer goods in China, consumer non-durables in Germany and turnover index in France. y4, y5 and y6 are daily data from 6 June 1996 to 14 December 2021 of Shanghai Stock Exchange Composite Index, German Index HDAX and French Index CAC 40. We can see that the positive and negative alternation is still present, but some of the monthly data is no longer significant due to the quality of the data. The explained variable is the variance of the logarithmic growth rate of several economic indicators. d t is the sampling interval of economic data. v a r i a n c e = β 0 + β 1 d t 2 + β 2 d t 3 2 + β 3 d t + β 4 d t 3 4 + β 5 d t 1 2 . Standard errors in parentheses: *** p < 0.01, ** p < 0.05, * p < 0.1. Source: CEIC dataset and Resset data.

Appendix B. Details of Numerical Simulation

In this section, we describe how the numerical results in Section 4.2 are obtained. We and [2] both derive investment decisions by solving for the maximisation of the investor’s lifetime utility. The core of our numerical simulation is therefore the value of utility under different decision approaches over a very long period of time. We use a lifetime utility in the form of E 0 + e ρ t ( c α α ) d t . The choice of parameters has already been described in the paper.
As we obtain daily returns on asset returns, we discretize the continuous model. At this point the differential equation becomes a difference equation. d t takes 1 / 250 to represent the length of a trading day. Let r, σ 0 and σ 1 be the moment estimates for the past 1000 trading days. We set the initial wealth w 0 to 1.
Investors form their current consumption decision c and investment decision x based on past w, r, σ 0 and σ 1 . Next the wealth value w is iterated according to the wealth value movement equation d w = ( w ( r f ( 1 x ) + r ^ x ) c ) d t ( r ^ is the real rate of return incorporating the stochastic shock). Based on the new w, r, σ 0 and σ 1 , the investor then makes the next decision.
After we obtain two different consumption sequences c based on different decision scenarios, we calculate the lifetime utility corresponding to the different decision options.

Appendix C. Proof

Proof of Lemma 1. 
Let F be expanded at the point of ( t , w ) (We need d z i d z j = 0 , d z 1 2 = d t , d z 1 4 = 3 d t ):
d u = F t d t + F w d w + F w w d w 2 + F w w w d w 3 + F w w w w d w 4 + o ( d t ) = F t d t + F w ( f ( t , w ) d t + g 1 ( t , w ) d z 0 + ( g 2 ( t , w ) + i g 2 ( t , w ) ) d z 1 ) + F w w 2 ( f ( t , w ) d t + g 1 ( t , w ) d z 0 + ( g 2 ( t , w ) + i g 2 ( t , w ) ) d z 1 ) 2 + F w w w 6 ( f ( t , w ) d t + g 1 ( t , w ) d z 0 + ( g 2 ( t , w ) + i g 2 ( t , w ) ) d z 1 ) 3 + F w w w w 24 ( f ( t , w ) d t + g 1 ( t , w ) d z 0 + ( g 2 ( t , w ) + i g 2 ( t , w ) ) d z 1 ) 4 + o ( d t ) = F t d t + F w ( f ( t , w ) d t + g 1 ( t , w ) d z 0 + ( g 2 ( t , w ) + i g 2 ( t , w ) ) d z 1 ) + F w w 2 ( g 1 ( t , w ) 2 d t + g 2 ( t , w ) 2 d t g 2 ( t , w ) 2 d t + 2 i g 2 ( t , w ) 2 d t ) + F w w w w 24 ( 3 ( g 2 4 + 4 i g 2 4 6 g 2 4 + 4 i g 2 4 + g 2 4 ) ) d t + o ( d t ) = F t + f F w + 1 2 g 1 2 F w w 1 2 g 2 4 F w w w w d t + 2 i g 2 2 d t + F w ( g 1 d z 0 + g 2 d z 1 + i g 2 d z 2 )
Proof of Lemma 2. 
Let F be expanded at the point of ( t , w ) :
R e E t ( d u ) = R e E t ( F t d t + F w d w + 1 2 d w F w w d w + 1 6 l 1 = 1 L l 2 = 1 L l 3 = 1 L F w l 1 w l 2 w l 3 d w l 1 d w l 2 d w l 3 + 1 24 l 1 = 1 L l 2 = 1 L l 3 = 1 L l 4 = 1 L F w l 1 w l 2 w l 3 w l 4 d w l 1 d w l 2 d w l 3 d w l 4 + o ( d t ) ) = R e E t ( F t d t + l = 1 L F w l ( f ( t , w l ) d t + k = 1 K g 1 , l , k ( t , w ) d z 0 , k + j = 1 J g 2 , l , j ( t , w ) d z 1 , j + i j = 1 J g 2 , l , j ( t , w ) d z 2 , j ) + l 1 = 1 L l 2 = 1 L F w l 1 w l 2 2 ( f ( t , w l 1 ) d t + k = 1 K g 1 , l 1 , k ( t , w ) d z 0 , k + j = 1 J g 2 , l 1 , j ( t , w ) ( 1 + i ) d z 1 , j ) ( f ( t , w l 2 ) d t + k = 1 K g 1 , l 2 , k ( t , w ) d z 0 , k + j = 1 J g 2 , l 2 , j ( t , w ) ( 1 + i ) d z 1 , j ) + 1 24 l 1 = 1 L l 2 = 1 L l 3 = 1 L l 4 = 1 L F w l 1 w l 2 w l 3 w l 4 ( f ( t , w l 1 ) d t + k = 1 K g 1 , l 1 , k ( t , w ) d z 0 , k + j = 1 J g 2 , l 1 , j ( t , w ) ( 1 + i ) d z 1 , j ) ( f ( t , w l 2 ) d t + k = 1 K g 1 , l 2 , k ( t , w ) d z 0 , k + j = 1 J g 2 , l 2 , j ( t , w ) ( 1 + i ) d z 1 , j ) ( f ( t , w l 3 ) d t + k = 1 K g 1 , l 3 , k ( t , w ) d z 0 , k + j = 1 J g 2 , l 3 , j ( t , w ) ( 1 + i ) d z 1 , j ) ( f ( t , w l 4 ) d t + k = 1 K g 1 , l 4 , k ( t , w ) d z 0 , k + j = 1 J g 2 , l 4 , j ( t , w ) ( 1 + i ) d z 1 , j ) + o ( d t ) = R e E t ( F t d t + l = 1 L F w l f ( t , w l ) d t + 1 2 l 1 = 1 L l 2 = 1 L F w l 1 w l 2 k = 1 K g 1 , l 1 , k g 1 , l 2 , k d t + 1 24 l 1 = 1 L l 2 = 1 L l 3 = 1 L l 4 = 1 L F w l 1 w l 2 w l 3 w l 4 3 ( j 1 = 1 J j 2 = 1 J g 2 , l 1 , j 1 g 2 , l 2 , j 1 g 2 , l 3 , j 2 g 2 , l 4 , j 2 + j 1 = 1 J j 2 = 1 J g 2 , l 1 , j 1 g 2 , l 2 , j 2 g 2 , l 3 , j 1 g 2 , l 4 , j 2 + j 1 = 1 J j 2 = 1 J g 2 , l 1 , j 1 g 2 , l 2 , j 2 g 2 , l 3 , j 2 g 2 , l 4 , j 1 2 j = 1 J g 2 , l 1 , j g 2 , l 2 , j g 2 , l 3 , j g 2 , l 4 , j ) ( 1 + i ) 4 d t ) + o ( d t ) = ( F t + f F w + 1 2 1 g 1 F w w g 1 1 + F w l 1 w l 2 w l 3 w l 4 ( j = 1 J l 1 = 1 L l 2 = 1 L l 3 = 1 L l 4 = 1 L g 2 , l 1 , j g 2 , l 2 , j g 2 , l 3 , j g 2 , l 4 , j 3 2 j 1 = 1 J j 2 = 1 J l 1 = 1 L l 2 = 1 L l 3 = 1 L l 4 = 1 L g 2 , l 1 , j 1 g 2 , l 2 , j 1 g 2 , l 3 , j 2 g 2 , l 4 , j 2 ) ) d t + o ( d t )

References

  1. Merton, R.C. An intertemporal capital asset pricing model. Econom. J. Econom. Soc. 1973, 41, 867–887. [Google Scholar] [CrossRef]
  2. Merton, R.C. Optimum consumption and portfolio rules in a continuous-time model. In Stochastic Optimization Models in Finance; Elsevier: Amsterdam, The Netherlands, 1975; pp. 621–661. [Google Scholar]
  3. Davis, M.H.; Norman, A.R. Portfolio selection with transaction costs. Math. Oper. Res. 1990, 15, 676–713. [Google Scholar] [CrossRef] [Green Version]
  4. Zhou, C. A Jump-Diffusion Approach to Modeling Credit Risk and Valuing Defaultable Securities. Available at SSRN 39800 1997. Available online: https://papers.ssrn.com/sol3/papers.cfm?abstract_id=39800 (accessed on 10 January 2022).
  5. Liu, J. Portfolio selection in stochastic environments. Rev. Financ. Stud. 2007, 20, 1–39. [Google Scholar] [CrossRef]
  6. Kolmogorov, A.N. Wienersche spiralen und einige andere interessante kurven in hilbertscen raum, cr (doklady). Acad. Sci. URSS (NS) 1940, 26, 115–118. [Google Scholar]
  7. Mandelbrot, B.B.; Van Ness, J.W. Fractional Brownian motions, fractional noises and applications. SIAM Rev. 1968, 10, 422–437. [Google Scholar] [CrossRef]
  8. Peters, E.E. Fractal Market Analysis: Applying Chaos Theory to Investment and Economics; John Wiley & Sons: Hoboken, NJ, USA, 1994; Volume 24. [Google Scholar]
  9. Rogers, L.C.G. Arbitrage with fractional Brownian motion. Math. Financ. 1997, 7, 95–105. [Google Scholar] [CrossRef] [Green Version]
  10. Cheridito, P. Regularizing Fractional Brownian Motion with a View Towards Stock Price Modelling. Ph.D. Thesis, ETH Zurich, Zurich, Germany, 2001. [Google Scholar]
  11. Sottinen, T. Fractional Brownian motion, random walks and binary market models. Financ. Stochastics 2001, 5, 343–355. [Google Scholar] [CrossRef]
  12. Guasoni, P. No arbitrage under transaction costs, with fractional Brownian motion and beyond. Math. Financ. 2006, 16, 569–582. [Google Scholar] [CrossRef]
  13. Hu, Y.; Øksendal, B. Fractional white noise calculus and applications to finance. Infin. Dimens. Anal. Quantum Probab. Relat. Top. 2003, 6, 1–32. [Google Scholar] [CrossRef]
  14. Elliott, R.J.; Van Der Hoek, J. A general fractional white noise theory and applications to finance. Math. Financ. 2003, 13, 301–330. [Google Scholar] [CrossRef]
  15. Xiao, W.L.; Zhang, W.G.; Zhang, X.L.; Wang, Y.L. Pricing currency options in a fractional Brownian motion with jumps. Econ. Model. 2010, 27, 935–942. [Google Scholar] [CrossRef]
  16. Xiao, W.; Yu, J. Asymptotic theory for estimating drift parameters in the fractional Vasicek model. Econom. Theory 2019, 35, 198. [Google Scholar] [CrossRef] [Green Version]
  17. Jumarie, G. Merton’s model of optimal portfolio in a Black-Scholes Market driven by a fractional Brownian motion with short-range dependence. Insur. Math. Econ. 2005, 37, 585–598. [Google Scholar] [CrossRef]
  18. Miao, J.; Yang, X. Pricing model for convertible bonds: A mixed fractional Brownian motion with jumps. East Asian J. Appl. Math. 2015, 5, 222–237. [Google Scholar] [CrossRef]
  19. Shokrollahi, F.; Kılıçman, A. Actuarial approach in a mixed fractional Brownian motion with jumps environment for pricing currency option. Adv. Differ. Equs. 2015, 2015, 275. [Google Scholar] [CrossRef]
  20. Rao, B.P. Pricing geometric Asian power options under mixed fractional Brownian motion environment. Phys. A Stat. Mech. Its Appl. 2016, 446, 92–99. [Google Scholar]
  21. Shokrollahi, F.; Kılıçman, A.; Magdziarz, M. Pricing European options and currency options by time changed mixed fractional Brownian motion with transaction costs. Int. J. Financ. Eng. 2016, 3, 1650003. [Google Scholar] [CrossRef]
  22. Zhang, W.G.; Li, Z.; Liu, Y.J.; Zhang, Y. Pricing European Option Under Fuzzy Mixed Fractional Brownian Motion Model with Jumps. Comput. Econ. 2020, 58, 483–515. [Google Scholar] [CrossRef]
  23. Eom, C.; Oh, G.; Jung, W.S. Relationship between degree of efficiency and prediction in stock price changes. arXiv 2007, arXiv:0708.4178. [Google Scholar]
  24. Brock, W.A.; Sayers, C.L. Is the business cycle characterized by deterministic chaos? J. Monet. Econ. 1988, 22, 71–90. [Google Scholar] [CrossRef]
  25. Ashley, R.A.; Patterson, D.M. Linear versus nonlinear macroeconomies: A statistical test. Int. Econ. Rev. 1989, 30, 685–704. [Google Scholar] [CrossRef]
  26. Karuppiah, J.; Los, C.A. Wavelet multiresolution analysis of high-frequency Asian FX rates, Summer 1997. Int. Rev. Financ. Anal. 2005, 14, 211–246. [Google Scholar] [CrossRef]
  27. Wang, F.; Yamasaki, K.; Havlin, S.; Stanley, H.E. Statistical regularities of equity market activity. arXiv 2009, arXiv:0911.4258. [Google Scholar]
  28. Hassani, H.; Thomakos, D. A review on singular spectrum analysis for economic and financial time series. Stat. Its Interface 2010, 3, 377–397. [Google Scholar] [CrossRef]
  29. Onali, E.; Goddard, J. Are European equity markets efficient? New evidence from fractal analysis. Int. Rev. Financ. Anal. 2011, 20, 59–67. [Google Scholar] [CrossRef] [Green Version]
  30. Garnier, J.; Sølna, K. Correction to Black—Scholes formula due to fractional stochastic volatility. SIAM J. Financ. Math. 2017, 8, 560–588. [Google Scholar] [CrossRef] [Green Version]
  31. Caraiani, P. Nonlinear dynamics in CEE stock markets indices. Econ. Lett. 2012, 114, 329–331. [Google Scholar] [CrossRef]
  32. Bal, D.P.; Rath, B.N. Nonlinear causality between crude oil price and exchange rate: A comparative study of China and India. Energy Econ. 2015, 51, 149–156. [Google Scholar]
  33. Jin, X. The impact of 2008 financial crisis on the efficiency and contagion of Asian stock markets: A Hurst exponent approach. Financ. Res. Lett. 2016, 17, 167–175. [Google Scholar] [CrossRef]
  34. Bariviera, A.F. The inefficiency of Bitcoin revisited: A dynamic approach. Econ. Lett. 2017, 161, 1–4. [Google Scholar] [CrossRef] [Green Version]
  35. Wei, W.C. Liquidity and market efficiency in cryptocurrencies. Econ. Lett. 2018, 168, 21–24. [Google Scholar] [CrossRef]
  36. Takaishi, T. Rough volatility of Bitcoin. Financ. Res. Lett. 2020, 32, 101379. [Google Scholar] [CrossRef]
  37. Shintani, M.; Linton, O. Is there chaos in the world economy? A nonparametric test using consistent standard errors. Int. Econ. Rev. 2003, 44, 331–357. [Google Scholar] [CrossRef] [Green Version]
  38. Brock, W.A.; Hommes, C.H. A rational route to randomness. Econom. J. Econom. Soc. 1997, 65, 1059–1095. [Google Scholar] [CrossRef]
  39. Hommes, C.; Huang, H.; Wang, D. A robust rational route to randomness in a simple asset pricing model. J. Econ. Dyn. Control 2005, 29, 1043–1072. [Google Scholar] [CrossRef] [Green Version]
  40. Chiarella, C.; Dieci, R.; Gardini, L. Asset price and wealth dynamics in a financial market with heterogeneous agents. J. Econ. Dyn. Control 2006, 30, 1755–1786. [Google Scholar] [CrossRef]
  41. Sommervoll, D.E.; Borgersen, T.A.; Wennemo, T. Endogenous housing market cycles. J. Bank. Financ. 2010, 34, 557–567. [Google Scholar] [CrossRef] [Green Version]
  42. Landau, L.; Lifshitz, E. Quantum Mechanics: Non-Relativistic Theory; Elsevier: Amsterdam, The Netherlands, 1981; Volume 3. [Google Scholar]
  43. Chung, K.L.; Williams, R.J. Introduction to Stochastic Integration; Springer: Boston, MA, USA, 1990; Volume 2. [Google Scholar]
  44. Gollier, C.; Pratt, J.W. Risk vulnerability and the tempering effect of background risk. Econom. J. Econom. Soc. 1996, 64, 1109–1123. [Google Scholar] [CrossRef]
  45. Peiro, A. Skewness in financial returns. J. Bank. Financ. 1999, 23, 847–862. [Google Scholar] [CrossRef]
  46. Patton, A.J. On the out-of-sample importance of skewness and asymmetric dependence for asset allocation. J. Financ. Econom. 2004, 2, 130–168. [Google Scholar] [CrossRef] [Green Version]
  47. Aït-sahali, Y.; Brandt, M.W. Variable selection for portfolio choice. J. Financ. 2001, 56, 1297–1351. [Google Scholar] [CrossRef] [Green Version]
  48. Kimball, M.S. Standard risk aversion. Econom. J. Econom. Soc. 1993, 61, 589–611. [Google Scholar] [CrossRef]
  49. Kyle, A.S.; Xiong, W. Contagion as a wealth effect. J. Financ. 2001, 56, 1401–1440. [Google Scholar] [CrossRef]
  50. Aït-Sahalia, Y.; Lo, A.W. Nonparametric estimation of state-price densities implicit in financial asset prices. J. Financ. 1998, 53, 499–547. [Google Scholar] [CrossRef] [Green Version]
  51. Ahmed, W.M.; Al Mafrachi, M. Do higher-order realized moments matter for cryptocurrency returns? Int. Rev. Econ. Financ. 2021, 72, 483–499. [Google Scholar] [CrossRef]
  52. Bhandari, R.; Das, S.R. Options on portfolios with higher-order moments. Financ. Res. Lett. 2009, 6, 122–129. [Google Scholar] [CrossRef]
Figure 1. U. S. consumption data increases with sampling intervals.
Figure 1. U. S. consumption data increases with sampling intervals.
Fractalfract 06 00099 g001
Table 1. Descriptive Statistics.
Table 1. Descriptive Statistics.
VariablesConsumptionBP500Tick_Shanghai
Obs12020001200
Mean0.02996770.08732810.0000137
Std. Dev.0.02919150.04921398.78 × 10 6
Min0.00002970.00012192.43 × 10 8
Max0.09880040.18089460.0000294
BeginJan-591 January 19761 December 2018
EndJan-1930 April 202031 December 2018
This table is to show the descriptive statistics of the variance. Source: CEIC dataset and Resset data.
Table 2. Data Decomposition Results.
Table 2. Data Decomposition Results.
VariablesU. S. ConsumptionBP500Tick_Shanghai
d t 2 2.19 × 10 5 ***2.83 × 10 7 ***−5.96 × 10 11 ***
7.11 × 10 8 1.26 × 10 8 0
d t 3 2 −0.000357 ***−3.83 × 10 5 ***5.99 × 10 9 ***
2.69 × 10 6 1.83 × 10 6 2.69 × 10 10
d t 0.00476 ***0.00310 ***−2.91 × 10 7 ***
6.29 × 10 5 0.0001621.85 × 10 8
d t 3 4 −0.0114 ***−0.0166 ***1.37 × 10 6 ***
0.0001960.000979.79 × 10 8
d t 1 2 0.00887 ***0.0287 ***−1.88 × 10 6 ***
0.0001950.001841.64 × 10 7
Constant−0.00188 ***−0.0255 ***1.01 × 10 6 ***
6.57 × 10 5 0.002191.53 × 10 7
Observations12020001200
R-squared10.9971
The explained variable is the variance of the logarithmic growth rate of several economic indicators. d t is the sampling interval of economic data. v a r i a n c e = β 0 + β 1 d t 2 + β 2 d t 3 2 + β 3 d t + β 4 d t 3 4 + β 5 d t 1 2 . Standard errors in parentheses: *** p < 0.01, ** p < 0.05, * p < 0.1. Source: CEIC dataset and Resset data.
Table 3. Comparison of Results and Numerical Simulations.
Table 3. Comparison of Results and Numerical Simulations.
[2][18,19,20,21,22]This Paper
Research Issues
Subject of StudyAny Financial AssetOptionsAny Financial Asset
Consider Personal DecisionsYESNOYES
Number of AssetsMulti-AssetSingle AssetMulti-Asset
Stochastic ProcessBrownian MotionFBMFBM
Solution Method
Solving MethodsStochastic Optimal ControlArbitrage Free PricingStochastic Optimal Control
Core LemmaIto LemmaIto IntegralNew Lemma
Expandability
DSGE ModelYESNOYES
Heterogeneous Agents ModelYESNOYES
Study Results
Difference from RealityNo Consideration of FBMNo MicrofoundationsImaginary Numbers
FBM is IncludedNOYESYES
Economic Variable Correlation is GivenYESNOYES
Agent Welfare
Utility Value−2.9227-−2.1076
The data is daily data from 1 May 1976 to 19 January 2022 of BP500. BM is Brownian Motion. FBM is Fractional Brownian Motion. Source: CEIC dataset and Resset data.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yan, Y.; Wang, Y. Asset Pricing Model Based on Fractional Brownian Motion. Fractal Fract. 2022, 6, 99. https://doi.org/10.3390/fractalfract6020099

AMA Style

Yan Y, Wang Y. Asset Pricing Model Based on Fractional Brownian Motion. Fractal and Fractional. 2022; 6(2):99. https://doi.org/10.3390/fractalfract6020099

Chicago/Turabian Style

Yan, Yu, and Yiming Wang. 2022. "Asset Pricing Model Based on Fractional Brownian Motion" Fractal and Fractional 6, no. 2: 99. https://doi.org/10.3390/fractalfract6020099

APA Style

Yan, Y., & Wang, Y. (2022). Asset Pricing Model Based on Fractional Brownian Motion. Fractal and Fractional, 6(2), 99. https://doi.org/10.3390/fractalfract6020099

Article Metrics

Back to TopTop