Moving average options: Machine Learning and Gauss-Hermite quadrature for a double non-Markovian problem

Evaluating moving average options is a tough computational challenge for the energy and commodity market as the payoff of the option depends on the prices of a certain underlying observed on a moving window so, when a long window is considered, the pricing problem becomes high dimensional. We present an efficient method for pricing Bermudan style moving average options, based on Gaussian Process Regression and Gauss-Hermite quadrature, thus named GPR-GHQ. Specifically, the proposed algorithm proceeds backward in time and, at each time-step, the continuation value is computed only in a few points by using Gauss-Hermite quadrature, and then it is learned through Gaussian Process Regression. We test the proposed approach in the Black-Scholes model, where the GPR-GHQ method is made even more efficient by exploiting the positive homogeneity of the continuation value, which allows one to reduce the problem size. Positive homogeneity is also exploited to develop a binomial Markov chain, which is able to deal efficiently with medium-long windows. Secondly, we test GPR-GHQ in the Clewlow-Strickland model, the reference framework for modeling prices of energy commodities. Finally, we consider a challenging problem which involves double non-Markovian feature, that is the rough-Bergomi model. In this case, the pricing problem is even harder since the whole history of the volatility process impacts the future distribution of the process. The manuscript includes a numerical investigation, which displays that GPR-GHQ is very accurate and it is able to handle options with a very long window, thus overcoming the problem of high dimensionality.

Machine Learning for Pricing American Options in High-Dimensional Markovian and non-Markovian models

In this paper we propose two efficient techniques which allow one to compute the price of American basket options. In particular, we consider a basket of assets that follow a multi-dimensional Black-Scholes dynamics. The proposed techniques, called GPR Tree (GRP-Tree) and GPR Exact Integration (GPR-EI), are both based on Machine Learning, exploited together with binomial trees or with a closed formula for integration. Moreover, these two methods solve the backward dynamic programming problem considering a Bermudan approximation of the American option. On the exercise dates, the value of the option is first computed as the maximum between the exercise value and the continuation value and then approximated by means of Gaussian Process Regression. The two methods mainly differ in the approach used to compute the continuation value: a single step of binomial tree or integration according to the probability density of the process. Numerical results show that these two methods are accurate and reliable in handling American options on very large baskets of assets. Moreover we also consider the rough Bergomi model, which provides stochastic volatility with memory. Despite this model is only bidimensional, the whole history of the process impacts on the price, and handling all this information is not obvious at all. To this aim, we present how to adapt the GPR-Tree and GPR-EI methods and we focus on pricing American options in this non-Markovian framework.

Variance Reduction Applied to Machine Learning for Pricing Bermudan/American Options in High Dimension

In this paper we propose an efficient method to compute the price of multi-asset American options, based on Machine Learning, Monte Carlo simulations and variance reduction technique. Specifically, the options we consider are written on a basket of assets, each of them following a Black-Scholes dynamics. In the wake of Ludkovskis approach (2018), we implement here a backward dynamic programming algorithm which considers a finite number of uniformly distributed exercise dates. On these dates, the option value is computed as the maximum between the exercise value and the continuation value, which is obtained by means of Gaussian process regression technique and Monte Carlo simulations. Such a method performs well for low dimension baskets but it is not accurate for very high dimension baskets. In order to improve the dimension range, we employ the European option price as a control variate, which allows us to treat very large baskets and moreover to reduce the variance of price estimators. Numerical tests show that the proposed algorithm is fast and reliable, and it can handle also American options on very large baskets of assets, overcoming the problem of the curse of dimensionality.

Gaussian Process Regression for Pricing Variable Annuities with Stochastic Volatility and Interest Rate

In this paper we investigate price and Greeks computation of a Guaranteed Minimum Withdrawal Benefit (GMWB) Variable Annuity (VA) when both stochastic volatility and stochastic interest rate are considered together in the Heston Hull-White model. We consider a numerical method the solves the dynamic control problem due to the computing of the optimal withdrawal. Moreover, in order to speed up the computation, we employ Gaussian Process Regression (GPR). Starting from observed prices previously computed for some known combinations of model parameters, it is possible to approximate the whole price function on a defined domain. The regression algorithm consists of algorithm training and evaluation. The first step is the most time demanding, but it needs to be performed only once, while the latter is very fast and it requires to be performed only when predicting the target function. The developed method, as well as for the calculation of prices and Greeks, can also be employed to compute the no-arbitrage fee, which is a common practice in the Variable Annuities sector. Numerical experiments show that the accuracy of the values estimated by GPR is high with very low computational cost. Finally, we stress out that the analysis is carried out for a GMWB annuity but it could be generalized to other insurance products.

Taxation of a GMWB Variable Annuity in a Stochastic Interest Rate Model

Modeling taxation of Variable Annuities has been frequently neglected but accounting for it can significantly improve the explanation of the withdrawal dynamics and lead to a better modeling of the financial cost of these insurance products. The importance of including a model for taxation has first been observed by Moenig and Bauer (2016) while considering a GMWB Variable Annuity. In particular, they consider the simple Black-Scholes dynamics to describe the underlying security. Nevertheless, GMWB are long term products and thus accounting for stochastic interest rate has relevant effects on both the financial evaluation and the policy holder behavior, as observed by Gouden`ege et al. (2018). In this paper we investigate the outcomes of these two elements together on GMWB evaluation. To this aim, we develop a numerical framework which allows one to efficiently compute the fair value of a policy. Numerical results show that accounting for both taxation and stochastic interest rate has a determinant impact on the withdrawal strategy and on the cost of GMWB contracts. In addition, it can explain why these products are so popular with people looking for a protected form of investment for retirement.

Computing Credit Valuation Adjustment solving coupled PIDEs in the Bates model

Credit value adjustment (CVA) is the charge applied by financial institutions to the counterparty to cover the risk of losses on a counterpart default event. In this paper we estimate such a premium under the Bates stochastic model (Bates [4]), which considers an underlying affected by both stochastic volatility and random jumps. We propose an efficient method which improves the finite-difference Monte Carlo (FDMC) approach introduced by de Graaf et al. [11]. In particular, the method we propose consists in replacing the Monte Carlo step of the FDMC approach with a finite difference step and the whole method relies on the efficient solution of two coupled partial integro-differential equations (PIDE) which is done by employing the Hybrid Tree-Finite Difference method developed by Briani et al. [6, 7, 8]. Moreover, the direct application of the hybrid techniques in the original FDMC approach is also considered for comparison purposes. Several numerical tests prove the effectiveness and the reliability of the proposed approach when both European and American options are considered.

Pricing and Hedging GMWB in the Heston and in the Black-Scholes with Stochastic Interest Rate Models

Valuing Guaranteed Minimum Withdrawal Benefit (GMWB) has attracted significant attention from both the academic field and real world financial markets. As remarked by Yang and Dai, the Black and Scholes framework seems to be inappropriate for such a long maturity products. Also Chen Vetzal and Forsyth in showed that the price of these products is very sensitive to interest rate and volatility parameters. We propose here to use a stochastic volatility model (Heston model) and a Black Scholes model with stochastic interest rate (Hull White model). For this purpose we present four numerical methods for pricing GMWB variables annuities: a hybrid tree-finite difference method and a Hybrid Monte Carlo method, an ADI finite difference scheme, and a Standard Monte Carlo method. These methods are used to determine the no-arbitrage fee for the most popul