Free boundary and optimal stopping problems for American Asian options

We give a complete and self-contained proof of the existence of a strong solution to the free boundary and optimal stopping problems for pricing American path-dependent options. The framework is sufficiently general to include geometric Asian options with nonconstant volatility and recent path-dependent volatility models.      

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 form formula for integration. Moreover, these two methods solve the backward dynamic programing 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 the 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 that this model is only bidimensional, the whole history of the process impacts on the price, and how to handle 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.     

Additive and multiplicative duals for American option pricing

We investigate and compare two dual formulations of the American option pricing problem based on two decompositions of supermartingales: the additive dual of Haugh and Kogan (Oper. Res. 52:258–270, 2004) and Rogers (Math. Finance 12:271–286, 2002) and the multiplicative dual of Jamshidian (Minimax optimality of Bermudan and American claims and their Monte- Carlo upper bound approximation. NIB Capital, The Hague, 2003). Both provide upper bounds on American option prices; we show how to improve these bounds iteratively and use this to show that any multiplicative dual can be improved by an additive dual and vice versa. This iterative improvement converges to the optimal value function. We also compare bias and variance under the two dual formulations as the time horizon grows; either method may have smaller bias, but the variance of the multiplicative method typically grows much faster than that of the additive method. We show that in the case of a discrete state space, the additive dual coincides with the dual of the optimal stopping problem in the sense of linear programming duality and the multiplicative method arises through a nonlinear duality.      

Enhanced policy iteration for American options via scenario selection

Kolodko and Schoenmakers (2006 Kolodko, A and Schoenmakers, J. 2006. Iterative construction of the optimal Bermudan stopping time. Finan. Stochast., 10: 27–49. [Crossref], [Web of Science ®] , [Google Scholar]) and Bender and Schoenmakers (2006 Bender, C, Kolodko, A and Schoenmakers, J. 2006. Iterating cancelable snowballs and related exotics in a many-factor Libor model. Risk, September: 126–130.  [Google Scholar]) introduced a policy iteration that allows the achievement of a tight lower approximations of the price for early exercise options via a nested Monte Carlo simulation in a Markovian setting. In this paper we enhance the algorithm by a scenario selection method. It is demonstrated by numerical examples that the scenario selection can significantly reduce the number of inner simulations actually performed, and thus can greatly speed up the method (by up to a factor of 15 in some examples). Moreover, it is shown that the modified algorithm retains the desirable properties of the original, such as the monotone improvement property, termination after a finite number of iteration steps, and numerical stability.     

Primal–dual quasi-Monte Carlo simulation with dimension reduction for pricing American options

The pricing of American options is one of the most challenging problems in financial engineering due to the involved optimal stopping time problem, which can be solved by using dynamic programming (DP). But applying DP is not always practical, especially when the state space is high dimensional. However, the curse of dimensionality can be overcome by Monte Carlo (MC) simulation. We can get lower and upper bounds by MC to ensure that the true price falls into a valid confidence interval. During the recent decades, progress has been made in using MC simulation to obtain both the lower bound by least-squares Monte Carlo method (LSM) and the upper bound by duality approach. However, there are few works on pricing American options using quasi-Monte Carlo (QMC) methods, especially to compute the upper bound. For comparing the sample variances and standard errors in the numerical experiments, randomized QMC (RQMC) methods are usually used. In this paper, we propose to use RQMC to replace MC simulation to compute both the lower bound (by the LSM) and the upper bound (by the duality approach). Moreover, we propose to use dimension reduction techniques, such as the Brownian bridge, principal component analysis, linear transformation and the gradients based principle component analysis. We perform numerical experiments on American–Asian options and American max-call options under the Black–Scholes model and the variance gamma model, in which the options have the path-dependent feature or are written on multiple underlying assets. We find that RQMC in combination with dimension reduction techniques can significantly increase the efficiency in computing both the lower and upper bounds, resulting in better estimates and tighter confidence intervals of the true price than pure MC simulation.