Computational technique for simulating variable-order fractional Heston model with application in US stock market

In this paper, a numerical technique is developed to discretize variable-order fractional Heston differential equation. The proposed strategy is followed by an optimization technology, genetic algorithm, for tuning the unknown parameters in the proposed model. The performance of the model is analyzed to profit and loss 500 close index from the US stock markets. Simulations illustrate the application of the proposed technique.


Introduction
The Black-Scholes model has been introduced based on the geometric Brownian motion to providing a closed-form pricing formula for the European options and describing the behavior of underlying asset prices [1][2][3]. The assumptions of the Black-Scholes model are unrealistic due to its inability to generate volatility satisfying the market observations, being nonnegative and mean-reverting [4][5][6]. To overcome this limitation, in 1993, Heston suggested Heston's stochastic volatility (HSV) model [6]. The HSV model has been extended in finance for modeling the dynamics of implied volatilities and providing their user with simple breakeven accounting conditions for the profit and loss (S&P) of a hedged position [2,[7][8][9][10][11][12][13][14][15][16][17].
The HSV model, in the risk-neutral measure, is presented as follows is the integer-order derivative operator, S(t) denotes the stock price, and V(t) denotes it's return variance at time t that referred to in the literature as volatility [18]. The parameters r and show the risk-free rate and the volatility long-run mean, respectively. Moreover, the parameter reveals the mean reversion rate toward that controls the speed volatility going back to its average value. The volatility of volatility, , is a significant risk factor that depicts the kurtosis in the stock return rate distribution. The Brownian motions, 1 (t) and 2 (t) with correlation coefficient, ∈ [−1, 1] , are in the risk-neutral measure.
Despite its popularity, research on efficient discretizations of the continuous time dynamics of the Heston model has generated less attention by researchers. In 2006, an exact simulation technique was suggested by Broadie and Kaya for the HSV model [19]. Subsequently, in [20], a simulation scheme including the quadratic-exponential technique has been developed, and in [21], a second-order discretization scheme has been discussed. Along the same line of thought, the drift-implicit Milstein algorithm for the volatility [22], a Euler discretization for the log-Heston price [16] and an analytic approach for degenerate parabolic problem associated with the HSV model [17] were presented.
In the last decades, the theory of fractional calculus was motivated as a useful mathematical tool to handle application of associated concepts in the areas of physics, chemistry, economics, finance and engineering sciences [23][24][25][26][27][28][29][30][31][32][33][34]. To the author's best knowledge, while numerous stochastic volatility models driven by fractional Brownian motion have been considered [35][36][37][38], stochastic volatility models with fixed-order and variable-order fractional derivative operators have not been carried out and this topic is far from being fully explored. More recently, just a numerical discretization technique has been proposed for approximation of a class of fractional stochastic differential equations driven by Brownian motion [39].
In this work, we assume that Brownian motions, 1 (t) and 2 (t) , t ≥ 0 , be real continuous functions that defined on a filtered probability space ( ,  , ℙ) with a normal filtration ( t ) t≥0 . We consider the variable-order fractional HSV (VOF-HSV) model as (t) ∈ ℝ + and (t) ∈ ℝ + , respectively. In general, variableorder fractional derivative operator was defined by Lorenzo and Hartley as [40] where u(t) is (m − 1) , m ∈ ℕ , times continuously differentiable and u (m) (t) is once integrable, is an auxiliary variable that belongs to the interval [0, t], and (⋅) denotes the Gamma function. The concept and properties of variableorder fractional derivative operators and their application have been discussed in [41][42][43][44][45].
The outline of the paper is organized as follows. In "Optimal VOF-HSV model" section, the tuning technique is described for designing optimal VOF-HSV model for US stock market. For this propose, a discretization algorithm for the VOF-HSV model is proposed and followed by an optimization technique, genetic algorithm (GA), for tuning the unknown parameters in the VOF-HSV model. The advantages of using the VOF-HSV model are shown and verified by considering different performance criteria in "Numerical results" section. Finally, "Conclusion" section draws the conclusions.

Optimal VOF-HSV model
Throughout this paper, we let = [0, T] with a uniform grid t j = jh , j = 0, 1, … , n , n ∈ ℤ + , such that n + 1 = T h , and n = (t n ) , n = (t n ) , S n = S(t n ) and V n = V(t n ) . In this section, an optimization strategy is formulated based on discretizing the equations in (2) and minimizing performance index under the following propositions.

e n o t e s t h e u n i f o r m s t e p s i z e ,
where (⋅) denotes convergence order of approximation.
and ‖u �� (t)‖ ≤ , j = 0, 1, … , n . Then, there exists an n+1 dependent constant, i.e., C n+1 > 0 , such that the truncated error of the variable-order fractional derivative operator obtained by the M-algorithm satisfies Let ℰ(t) be the error function in the interval (0, t n ] ; therefore, □ Using Propositions 2.1, the discretized equations of VOF-HSV model (2) are derived as follows There are different performance criteria for measuring goodness of fit typically summarizing discrepancy between observed values and estimated values by the prediction model. We analyze the performance of proposed models, in the viewpoint of root-mean-square error (RMSE) or the residual standard deviation defined as where RSS denotes the square root of residual sum of squares, N and p are the number of observations and parameters in the equation, respectively [47].
For estimating the Bayesian information criterion (BIC), we apply the following formula [47] Moreover, the relative quality of statistical models for a given set of data can be analyzed in the perspective of the Akaike information criterion (AIC) defined as [48] A lower RMSE , BIC or AIC value indicates a better fit.

Numerical results
In this section, fixed-order and VOF-HSV models are used to estimate the behavior of the price of S&P 500 index of American options. The S&P 500 is a US stock market index that covers all the large market public companies recorded on the New York Stock Exchange (NYSE) or National Association of Securities Dealers Automated Quotations (NAS-DAQ) [49]. For practical applicability, in this work, the parameters of the financial market model are chosen based on [50]. In [50], authors obtained estimate values of parameters for the typical Heston model by considering S&P 500 index returns. Throughout the numerical analysis, the main parameters satisfy that the mean reversion rate converges to = 2.75 , the long-run mean of the volatility converges to = 0.035 , the volatility of volatility converges to = 0.425 , and the correlation coefficient converges to = −0.4644.
For constructing VOF-HSV model, (t) = c 1 + c 2 t and (t) = c 3 + c 4 t functions where 0 ≤ c i ≤ 1 , i = 1, … , 4 , with four unknown parameters being considered. Next, an optimization method, GA, is used to find the optimal values of the unknowns c i , i = 1, … , 4 , by minimizing the mean absolute error (MAE) at the discretized points goodness, that is, by minimizing MAE = 1 where S(t n ) and Ŝ n are the discretized value of the model and experimental value, respectively. The optimization algorithm is based on discretizing (6) by using the M-algorithm formulated in Propositions 2.1 with n = 1227 equal mesh points, corresponding to observed data, so that approximation error can be determined by using Proposition 2.2. All the computations are performed under Maple v18 on an Intel (R) Core (TM) i7-7500U CPU @ 2.70 GHz machine. Table 1 shows the obtained values for the unknown parameters of the integer-order (i.e., c 1 = c 3 = 1 and c 2 = c 4 = 0 ), the fixed-order fractional (i.e., c 2 = c 4 = 0 ), and the VOF-HSV models when t ≤ 1 . In Table 2, the results in the perspective of the indices RMSE , AIC and BIC reveal that there is evidence of positive serial correlation in all considered cases of HSV models and the variable-order fractional model gives the best fitting out of the all proposed models. Moreover, there is a significant improvement in the response estimation by using the fractional Table 1 The MAEs and optimal parameters of the integer-order, the fixed-order fractional and VOF-HSV models for (t) = c 1 + Since proposed models cannot be solved analytically using the standard mathematical techniques, the upper 95% confidence interval (CI) was utilized to predict the behavior of sample trajectories of solutions by means of the approximated HSV models. Moreover, the statistical analysis of the process is provided in Table 3. In Table 3, the approximation values of the mean, median, first and third quartiles, kurtosis, skewness, standard deviation (STD) and 95% CI of the 50 simulated trajectories in time for several variation of orders for stock prices are listed.
The comparison of experimental stock price and simulation results of 50 trajectories of the stock prices of proposed HSV models is shown in Fig. 3. In Fig. 3, the red circles indicate the trajectories corresponding to the observed stock price. The blue lines represent the simulation results of 50 trajectories of the approximated stock prices (Fig. 3, left column). The black and green lines indicate the 95% CI area of the 50 trajectories and the empirical mean of the process (Fig. 3. right column),  . 1 The experimental stock price and the variation of the orders (t) and (t) in the approximated HSV models for the experimental stock price tests Fig. 2 The experimental stock price volatility and the variation of the orders (t) and (t) in the approximated HSV models for the experimental volatility tests respectively, for proposed HSV models. The purple lines are first and third quartiles of these models.

Conclusion
It has been shown that using fixed-order and variable-order fractional derivative operators instead of an integer-order derivative operator in models can apparently lead to better results since one has an extra freedom degree. In this paper, the objective of introducing modified Heston model was to capture the complexities in simulating the scenarios of equity movement in finance. The proposed Heston's stochastic volatility models, due to random volatility and memory or hereditary properties, have more flexibility with market data than the constant volatility models. In this study, an optimization technique, genetic algorithm, was utilized to find an optimal set of unknown parameters of models. For this propose, an explicit numerical technique has been developed for solving Fig. 3 The red circles indicate the trajectories corresponding to the observed stock price. The blue lines represent the simulation results of 50 trajectories of the approximated stock prices for proposed HSV models (left column). The black and the purple lines indicate the 95% CI and first and third quartiles area of the 50 trajectories, respectively, and green line is the empirical mean of the process (right column), for proposed HSV models Heston's stochastic volatility models. Furthermore, experimental stock price tests of a US stock market were used to estimate fixed-order and variable-order fractional Heston's stochastic volatility models. It was also shown that proposed fractional Heston's stochastic volatility models are significantly better in estimating the stock price in comparison to integer-order Heston's stochastic volatility models. Moreover, variable-order fractional Heston's stochastic volatility model is superior to fixed-order fractional Heston's stochastic volatility model.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creat iveco mmons .org/licen ses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.