Pricing arithmetic Asian option under a two-factor stochastic volatility model with jumps

This paper presents an efficient Monte Carlo simulation scheme based on the variance reduction methods to evaluate arithmetic average Asian options in the context of the double Heston's stochastic volatility model with jumps. This paper consists of two essential parts. The first part presents a new flexible stochastic volatility model, namely, the double Heston model with jumps. In the second part, by combining two variance reduction procedures via Monte Carlo simulation, we propose an efficient Monte Carlo simulation scheme for pricing arithmetic average Asian options under the double Heston model with jumps. Numerical results illustrate the efficiency of our method.     

Analytical and numerical studies on the second-order asymptotic expansion method for European option pricing under two-factor stochastic volatilities

The celebrated Black–Scholes model made the assumption of constant volatility but empirical studies on implied volatility and asset dynamics motivated the use of stochastic volatilities. Christoffersen in 2009 showed that multi-factor stochastic volatilities models capture the asset dynamics more realistically. Fouque in 2012 used it to price European options. In 2013, Chiarella and Ziveyi considered Christoffersen’s ideas and introduced an asset dynamics where the two volatilities of the Heston type act separately and independently on the asset price, and using Fourier transform for the asset price process and double Laplace transform for the two volatilities processes, solved a pricing problem for American options. This paper considers the Chiarella and Ziveyi model and parameterizes it so that the volatilities revert to the long-run-mean with reversion rates that mimic fast (for example daily) and slow (for example seasonal) random effects. Applying asymptotic expansion method presented by Fouque in 2012, we make an extensive and detailed derivation of the approximation prices for European options. We also present numerical studies on the behavior and accuracy of our first- and second-order asymptotic expansion formulas.     

Common-factor stochastic volatility modelling with observable proxy

Multi-asset modelling is of fundamental importance to financial applications such as risk management and portfolio selection. In this article, we propose a multivariate stochastic volatility modelling framework with a parsimonious and interpretable correlation structure. Building on well-established evidence of common volatility factors among individual assets, we consider a multivariate diffusion process with a common-factor structure in the volatility innovations. Upon substituting an observable market proxy for the common volatility factor, we markedly improve the estimation of several model parameters and latent volatilities. The model is applied to a portfolio of several important constituents of the S&P500 in the financial sector, with the VIX index as the common-factor proxy. We find that the prediction intervals for asset forecasts are comparable to those of more complex dependence models, but that option-pricing uncertainty can be greatly reduced by adopting a common-volatility structure. The Canadian Journal of Statistics 48: 36–61; 2020 © 2020 Statistical Society of Canada     

Computational aspects of pricing foreign exchange options with stochastic volatility and stochastic interest rates

This paper presents an analytic result for the price of a European call option on a foreign exchange currency rate. Market volatility is assumed correlated with the exchange rate and interest rates, domestic and foreign, are assumed to be stochastic. Integrals involving interest rates are derived, characteristic functions are produced, and, with evaluation, the nature of the integrals involved in Fourier inversion is examined. By comparison with FX market data, some of the effects of the nature of stochastic interest rates upon option prices are examined.     

An analytical approximation method for pricing barrier options under the double Heston model

Abstract

The purpose of the paper is to provide an efficient pricing method for single barrier options under the double Heston model. By rewriting the model as a singular and regular perturbed BS model, the double Heston model can separately mimic a fast time-scale and a slow time-scale. With the singular and regular perturbation techniques, we analytically derive the first-order asymptotic expansion of the solution to a barrier option pricing partial differential equation. The convergence and efficiency of the approximate method is verified by Monte Carlo simulation. Numerical results show that the presented asymptotic pricing method is fast and accurate.     

Methodology for stochastic volatility process calibration application to the CAC 40 index

In considering volatility as a stochastic, the aim of this paper is to estimate the four parameters related to a particular stochastic process named P1 and based on a Wiener–Levy process. We present the methodology to estimate its four parameters. We calibrate this theoretical model P1 to the CAC 40 index real data. In the same time, we test the normality of the random variables related to the two Wiener–Levy processes. The calibration is performed using the implemented aforesaid algorithm. We compare the stochastic process P1 with another process named P2 and to the Heston [Closed form solution for options with stochastic volatility with application to bonds and currency options, Rev. Financ. Stud. 6(2) (1993), pp. 327–343] process named H0 and to two other improved Heston processes named H1 and H2. For the empirical study, the same algorithm is used to calibrate the five processes. The calibration is based on a database including the CAC 40 index daily ‘closing fixing’ values for the time period from 3rd January 2005 to 22nd January 2007. The data are divided into 18 classes relative to 18 different contracts of European calls on the CAC 40 index. As a result, we find that, the normality test of the CAC 40 index is rejected which is in accordance with the previous original works dealing with this problem. For the five volatility processes, the normality test is verified almost for the same contracts. We also find that according to the used data, the process P1 and its equivalent H1 are the best for calibration.     

American option pricing under double Heston stochastic volatility model: simulation and strong convergence analysis

In this work, we investigate the double Heston model dynamics which is defined by two independent variance processes with non-Lipschitz diffusions. Next, it is analysed the strong convergence of the volatility processes of the double Heston model. Then, we examine the LSM algorithm to determine the American style option price in the double Heston model. Besides, by performing the antithetic simulation in the original LSM algorithm, we aim to accelerate the algorithm. The accuracy and speed of two algorithms are studied and at last, by some numerical results illustrate the accuracy of the proposed algorithms and examine the effect of the different parameters of the model on the value of the option.     

Estimation of the integrated volatility using noisy high-frequency data with jumps and endogeneity

In this paper, we investigate a new estimator of the integrated volatility of Itô semimartingales in the presence of both market microstructure noise and jumps when sampling times are endogenous. In the first step, our estimation wipes off the effects of the microstructure noise, and in the second step our estimator shrinks the effects of jumps. We provide consistency of the estimator when the jumps have finite variation and infinite variation and establish a central limit theorem for the estimator in a general endogenous time setting when the jumps only have finite variation. Simulation illustrates the performance of the proposed estimator.     

Kalman filter-based modelling and forecasting of stochastic volatility with threshold

We propose a parametric nonlinear time-series model, namely the Autoregressive-Stochastic volatility with threshold (AR-SVT) model with mean equation for forecasting level and volatility. Methodology for estimation of parameters of this model is developed by first obtaining recursive Kalman filter time-update equation and then employing the unrestricted quasi-maximum likelihood method. Furthermore, optimal one-step and two-step-ahead out-of-sample forecasts formulae along with forecast error variances are derived analytically by recursive use of conditional expectation and variance. As an illustration, volatile all-India monthly spices export during the period January 2006 to January 2012 is considered. Entire data analysis is carried out using EViews and matrix laboratory (MATLAB) software packages. The AR-SVT model is fitted and interval forecasts for 10 hold-out data points are obtained. Superiority of this model for describing and forecasting over other competing models for volatility, namely AR-Generalized autoregressive conditional heteroscedastic, AR-Exponential GARCH, AR-Threshold GARCH, and AR-Stochastic volatility models is shown for the data under consideration. Finally, for the AR-SVT model, optimal out-of-sample forecasts along with forecasts of one-step-ahead variances are obtained.     

Non‐parametric Threshold Estimation for Models with Stochastic Diffusion Coefficient and Jumps

Abstract. We consider a stochastic process driven by diffusions and jumps. Given a discrete record of observations, we devise a technique for identifying the times when jumps larger than a suitably defined threshold occurred. This allows us to determine a consistent non‐parametric estimator of the integrated volatility when the infinite activity jump component is Lévy. Jump size estimation and central limit results are proved in the case of finite activity jumps. Some simulations illustrate the applicability of the methodology in finite samples and its superiority on the multipower variations especially when it is not possible to use high frequency data.     

Optimal stopping for Brownian motion with applications to sequential analysis and option pricing

Herman Chernoff made fundamental contributions to analytical and computational methods for solving optimal stopping problems for Brownian motion. He also showed how these optimal stopping problems are closely related to some basic problems in sequential analysis and singular stochastic control. This paper gives a survey of these and related developments and describes some recent applications to option valuation in financial economics.     

Robust portfolio choice for a defined contribution pension plan with stochastic income and interest rate

This paper considers a robust portfolio choice problem for a defined contribution pension plan with stochastic income and stochastic interest rate. The investment objective of the pension plan is to maximize the expected utility of the wealth at the retirement time. We assume that the financial market consists of a stock, a zero-coupon bond and a risk-free asset. And the member of defined contribution pension plan is ambiguity-averse, which means that the member is uncertain about the expected return rate of the bond and stock. Meanwhile, the member's ambiguity-aversion level toward these two financial assets is quite different. The closed-form expressions of the robust optimal investment strategy and the corresponding value function are derived by adopting the stochastic dynamic programming approach. Furthermore, the sensitive analysis of model parameters on the optimal investment strategy are presented. We find that the member's aversion on model ambiguity increases her hedging demand and has remarkable impact on the optimal investment strategy. Moreover, we demonstrate that ignoring model uncertainty will lead to significant utility loss for the ambiguity-averse member, and the model uncertainty about the stock dynamics implies greater effect on the outcome of the investment than the bond.     

Stochastic differential equation systems for an SIS epidemic model with vaccination and immigration

Two Itô stochastic differential equation (SDE) systems are constructed for a Susceptible-Infected-Susceptible epidemic model with temporary vaccination. A constant number of new members enter the population and total size of the population is variable. Some conditions for disease extinction in the stochastic models are established and compared with conditions in deterministic one. It is shown that the two stochastic models are equivalent in the sense that their solutions come from same distribution. In addition, the SDE models are simulated and the equivalence of the two stochastic models is confirmed by numerical examples. The probability distribution for extinction is also obtained numerically, provided there exists a probability for disease persistence whereas the expected duration of epidemic is acquired when extinction occurs with probability 1.