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.     

Option pricing under regime switching

Abstract

This paper develops a family of option pricing models when the underlying stock price dynamic is modelled by a regime switching process in which prices remain in one volatility regime for a random amount of time before switching over into a new regime. Our family includes the regime switching models of Hamilton (Hamilton J 1989 Econometrica 57 357–84), in which volatility influences returns. In addition, our models allow for feedback effects from returns to volatilities. Our family also includes GARCH option models as a special limiting case. Our models are more general than GARCH models in that our variance updating schemes do not only depend on levels of volatility and asset innovations, but also allow for a second factor that is orthogonal to asset innovations. The underlying processes in our family capture the asymmetric response of volatility to good and bad news and thus permit negative (or positive) correlation between returns and volatility. We provide the theory for pricing options under such processes, present an analytical solution for the special case where returns provide no feedback to volatility levels, and develop an efficient algorithm for the computation of American option prices for the general case.     

Implied stopping rules for American basket options from Markovian projection

This work addresses the problem of pricing American basket options in a multivariate setting, which includes among others, the Bachelier and Black–Scholes models. In high dimensions, nonlinear PDE methods for solving the problem become prohibitively costly due to the curse of dimensionality. Instead, this work proposes to use a stopping rule that depends on the dynamics of a low-dimensional Markovian projection of the given basket of assets. From a numerical analysis point of view, we split the given non-smooth high-dimensional problem into two subproblems, namely one dealing with a smooth high-dimensionality integration in the parameter space and the other dealing with a low-dimensional, non-smooth optimal stopping problem in the projected state space. Assuming that we know the density of the forward process and using the Laplace approximation, we first efficiently evaluate the diffusion coefficient corresponding to the low-dimensional Markovian projection of the basket. Then, we approximate the optimal early exercise boundary of the option by solving an HJB PDE in the projected, low-dimensional space. The resulting near-optimal early exercise boundary is used to produce an exercise strategy for the high-dimensional option, thereby providing a lower bound for the price of the American basket option. A corresponding upper bound is also provided. These bounds allow one to assess the accuracy of the proposed pricing method. Indeed, our approximate early exercise strategy provides a straightforward lower bound for the American basket option price. Following a duality argument due to Rogers, we derive a corresponding upper bound solving only the low-dimensional optimal control problem. Numerically, we show the feasibility of the method using baskets with dimensions up to 50. In these examples, the resulting option price relative errors are only of the order of few percent.     

Are regime-shift sources of risk priced in the market?

In this paper we develop a discrete-time pricing model for European options where the log-return of the underlying asset is subject to discontinuous regime shifts in its mean and/or volatility which follow a Markov chain. The model allows for multiple regime shifts whose risk cannot be hedge out and thus must be priced in option market. The paper provides estimates of the price of regime-shift risk coefficients based on a joint estimation procedure of the Markov regime-switching process of the underlying stock and the suggested option pricing model. The results of the paper indicate that bull-to-bear and bear-to-crash regime shifts carry substantial prices of risk. Risk averse investors in the markets price these regime shifts by assigning higher transition (switching) probabilities to them under the risk neutral probability measure than under the physical. Ignoring these sources of risk will lead to substantial option pricing errors. In addition, the paper shows that investors also price reverse regime shifts, like the crash-to-bear and bear-to-bull ones, by assigning smaller transition probabilities under the risk neutral measure than the physical. Finally, the paper evaluates the pricing performance of the model and indicates that it can be successfully employed, in practice, to price European options.     

Valuation of quanto options in a Markovian regime-switching market: A Markov-modulated Gaussian HJM model

We consider the valuation of European quanto call options in an incomplete market where the domestic and foreign forward interest rates are allowed to exhibit regime shifts under the Heath–Jarrow–Morton (HJM) framework, and the foreign price dynamics is exogenously driven by a regime switching jump-diffusion model with Markov-modulated Poisson processes. We derive closed-form solutions for four different types of quanto call options, which include: options struck in a foreign currency, a foreign equity call struck in domestic currency, a foreign equity call option with a guaranteed exchange rate, and an equity-linked foreign exchange-rate call.     

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.     

Principle of equivalent utility and universal variable life insurance pricing

In this paper we study the pricing problem for a class of universal variable life (UVL) insurance products, using the idea of principle of equivalent utility. As the main features of UVL products we allow the (death) benefit to depend on certain indices or assets that are not necessarily tradable (e.g., pension plans), and we also consider the “multiple decrement” cases in which various status of the insured are allowed and the benefit varies in accordance with the status. Following the general theory of indifference pricing, we formulate the pricing problem as stochastic control problems, and derive the corresponding HJB equations for the value functions. In the case of exponential utilities, we show that the prices can be expressed explicitly in terms of the global, bounded solutions of a class of semilinear parabolic PDEs with exponential growth. In the case of general insurance models where multiple decrements and random time benefit payments are all allowed, we show that the price should be determined by the solutions to a system of HJB equations, each component corresponds to the value function of an optimization problem with the particular status of the insurer.     

American option pricing under the double Heston model based on asymptotic expansion

This paper focuses on pricing American put options under the double Heston model proposed by Christoffersen et al. By introducing an explicit exercise rule, we obtain the asymptotic expansion of the solution to the partial differential equation for pricing American put options. We calculate American option price by the sum of the European option price and the early exercise premium. The early exercise premium is calculated by the difference between the American and European option prices based on asymptotic expansions. The European option price is obtained by the efficient COS method. Based on the obtained American option price, the double Heston model is calibrated by minimizing the distance between model and market prices, which yields an optimization problem that is solved by a differential evolution algorithm combined with the Matlab function fmincon.m. Numerical results show that the pricing approach is fast and accurate. Empirical results show that the double Heston model has better performance in pricing short-maturity American put options and capturing the volatility term structure of American put options than the Heston model.     

Risk measures for derivatives with Markov-modulated pure jump processes

We consider a regime-switching HJB approach to evaluate risk measures for derivative securities when the price process of the underlying risky asset is governed by the exponential of a pure jump process with drift and a Markov switching compensator. The pure jump process is flexible enough to incorporate both the infinite, (small), jump activity and the finite, (large), jump activity. The drift and the compensator of the pure jump process switch over time according to the state of a continuous-time hidden Markov chain representing the state of an economy. The market described by our model is incomplete. Hence, there is more than one pricing kernel and there is no perfect hedging strategy for a derivative security. We derive the regime-switching HJB equations for coherent risk measures for the unhedged position of derivative securities, including standard European options and barrier options. For measuring risk inherent in the unhedged option position, we first need to mark the position into the market by valuing the option. We employ a well-known tool in actuarial science, namely, the Esscher transform to select a pricing kernel for valuation of an option and to generate a family of real-world probabilities for risk measurement. We also derive the regime-switching HJB-variational inequalities for coherent risk measures for American-style options.     

American options: the EPV pricing model

Summary We explicitly solve the pricing problem for perpetual American puts and calls, and provide an efficient semi-explicit pricing procedure for options with finite time horizon. Contrary to the standard approach, which uses the price process as a primitive, we model the price process as the expected present value of a stream, which is a monotone function of a Lévy process. Certain processes exhibiting mean-reverting, stochastic volatility and/or switching features can be modeled this way. This specification allows us to consider assets that pay no dividends at all when the level of the underlying stochastic factor is too low, assets that pay dividends at a fixed rate when the underlying stochastic process remains in some range, or capped dividends.The authors are grateful to the anonymous referees for valuable comments and suggestions.     

Pricing Options under Generalized GARCH and Stochastic Volatility Processes

In this paper, we develop an efficient lattice algorithm to price European and American options under discrete time GARCH processes. We show that this algorithm is easily extended to price options under generalized GARCH processes, with many of the existing stochastic volatility bivariate diffusion models appearing as limiting cases. We establish one unifying algorithm that can price options under almost all existing GARCH specifications as well as under a large family of bivariate diffusions in which volatility follows its own, perhaps correlated, process.     

American stochastic volatility call option pricing: A lattice based approach

This study presents a new method of pricing options on assets with stochastic volatility that is lattice based, and can easily accommodate early exercise for American options. Unlike traditional lattice methods, recombination is not a problem in the new model, and it is easily adapted to alternative volatility processes. Approximations are developed for European C.E.V. calls and American stochastic volatility calls. The application of the pricing model to exchange traded calls is also illustrated using a sample of market prices. Modifying the model to price American puts is straightforward, and the approach can easily be extended to other non-recombining lattices.     

Optimal high-frequency trading with limit and market orders

We propose a framework for studying optimal market-making policies in a limit order book (LOB). The bid–ask spread of the LOB is modeled by a tick-valued continuous-time Markov chain. We consider a small agent who continuously submits limit buy/sell orders at best bid/ask quotes, and may also set limit orders at best bid (resp. ask) plus (resp. minus) a tick for obtaining execution order priority, which is a crucial issue in high-frequency trading. The agent faces an execution risk since her limit orders are executed only when they meet counterpart market orders. She is also subject to inventory risk due to price volatility when holding the risky asset. The agent can then also choose to trade with market orders, and therefore obtain immediate execution, but at a less favorable price. The objective of the market maker is to maximize her expected utility from revenue over a short-term horizon by a trade-off between limit and market orders, while controlling her inventory position. This is formulated as a mixed regime switching regular/impulse control problem that we characterize in terms of a quasi-variational system by dynamic programming methods. Calibration procedures are derived for estimating the transition matrix and intensity parameters for the spread and for Cox processes modelling the execution of limit orders. We provide an explicit backward splitting scheme for solving the problem and show how it can be reduced to a system of simple equations involving only the inventory and spread variables. Several computational tests are performed both on simulated and real data, and illustrate the impact and profit when considering execution priority in limit orders and market orders.     

Pricing American options under stochastic volatility and stochastic interest rates

We introduce a new analytical approach to price American options. Using an explicit and intuitive proxy for the exercise rule, we derive tractable pricing formulas using a short-maturity asymptotic expansion. Depending on model parameters, this method can accurately price options with time-to-maturity up to several years. The main advantage of our approach over existing methods lies in its straightforward extension to models with stochastic volatility and stochastic interest rates. We exploit this advantage by providing an analysis of the impact of volatility mean-reversion, volatility of volatility, and correlations on the American put price.     

Machine learning for quantitative finance: fast derivative pricing,hedging and fitting

In this paper, we show how we can deploy machine learning techniques in the context of traditional quant problems. We illustrate that for many classical problems, we can arrive at speed-ups of several orders of magnitude by deploying machine learning techniques based on Gaussian process regression. The price we have to pay for this extra speed is some loss of accuracy. However, we show that this reduced accuracy is often well within reasonable limits and hence very acceptable from a practical point of view. The concrete examples concern fitting and estimation. In the fitting context, we fit sophisticated Greek profiles and summarize implied volatility surfaces. In the estimation context, we reduce computation times for the calculation of vanilla option values under advanced models, the pricing of American options and the pricing of exotic options under models beyond the Black–Scholes setting.     

Optimal discrete hedging of American options using an integrated approach to options with complex embedded decisions

In order to solve the problem of optimal discrete hedging of American options, this paper utilizes an integrated approach in which the writer’s decisions (including hedging decisions) and the holder’s decisions are treated on equal footing. From basic principles expressed in the language of acceptance sets we derive a general pricing and hedging formula and apply it to American options. The result combines the important aspects of the problem into one price. It finds the optimal compromise between risk reduction and transaction costs, i.e. optimally placed rebalancing times. Moreover, it accounts for the interplay between the early exercise and hedging decisions. We then perform a numerical calculation to compare the price of an agent who has exponential preferences and uses our method of optimal hedging against a delta hedger. The results show that the optimal hedging strategy is influenced by the early exercise boundary and that the worst case holder behavior for a sub-optimal hedger significantly deviates from the classical Black–Scholes exercise boundary.     

A new method for generating approximation algorithms for financial mathematics applications

This paper describes a new technique that can be used in financial mathematics for a wide range of situations where the calculation of complicated integrals is required. The numerical schemes proposed here are deterministic in nature but their proof relies on known results from probability theory regarding the weak convergence of probability measures. We adapt those results to unbounded payoffs under certain mild assumptions that are satisfied in finance. Because our approximation schemes avoid repeated simulations and provide computational savings, they can potentially be used when calculating simultaneously the price of several derivatives contingent on the same underlying. We show how to apply the new methods to calculate the price of spread options and American call options on a stock paying a known dividend. The method proves useful for calculations related to the log-Weibull model proposed recently for empirical asset pricing.     

Robustly Hedging Variable Annuities With Guarantees Under Jump and Volatility Risks

Recent variable annuities offer participation in the equity market and attractive protection against downside movements. Accurately quantifying this additional equity market risk and robustly hedging options embedded in the guarantees of variable annuities are new challenges for insurance companies. Due to sensitivities of the benefits to tails of the account value distribution, a simple Black–Scholes model is inadequate in preventing excessive liabilities. A model which realistically describes the real world price dynamics over a long time horizon is essential for the risk management of the variable annuities. In this article, both jump risk and volatility risk are considered for risk management of lookback options embedded in guarantees with a ratchet feature. We evaluate relative performances of delta hedging and dynamic discrete risk minimization hedging strategies. Using the underlying as the hedging instrument, we show that, under a Black–Scholes model, local risk minimization hedging can be significantly better than delta hedging. In addition, we compare risk minimization hedging using the underlying with that of using standard options. We demonstrate that, under a Merton's jump diffusion model, hedging using standard options is superior to hedging using the underlying in terms of the risk reduction. Finally, we consider a market model for volatility risks in which the at‐the‐money implied volatility is a state variable. We compute risk minimization hedging by modeling at‐the‐money Black–Scholes implied volatility explicitly; the hedging effectiveness is evaluated, however, under a joint model for the underlying price and implied volatility. Our computational results suggest that, when implied volatility risk is suitably modeled, risk minimization hedging using standard options, compared to hedging using the underlying, can potentially be more effective in risk reduction under both jump and volatility risks.     

Optimal liquidation under stochastic liquidity

We solve explicitly a two-dimensional singular control problem of finite fuel type for an infinite time horizon. The problem stems from the optimal liquidation of an asset position in a financial market with multiplicative and transient price impact. Liquidity is stochastic in that the volume effect process, which determines the intertemporal resilience of the market in the spirit of Predoiu et al. (SIAM J. Financ. Math. 2:183–212, 2011), is taken to be stochastic, being driven by its own random noise. The optimal control is obtained as the local time of a diffusion process reflected at a non-constant free boundary. To solve the HJB variational inequality and prove optimality, we need a combination of probabilistic arguments and calculus of variations methods, involving Laplace transforms of inverse local times for diffusions reflected at elastic boundaries.     

Pricing vulnerable options with correlated jump-diffusion processes depending on various states of the economy

In this paper, we use a Markov-modulated regime switching approach to model various states of the economy, and study the pricing of vulnerable European options when the dynamics of the underlying asset value and the asset value of the counterparty follow two correlated jump-diffusion processes under regime switching. The correlation is modelled by both the diffusion parts and the pure jump parts which describe the uncertainty of the value of the risky assets. We develop a method to determine an equivalent martingale measure and a parsimonious representation of the risk-neutral density is provided. Based on this, we derive an analytical pricing formula for vulnerable options via two-dimensional Laplace transforms, and implement the formula through numerical Laplace inversion.