A data and digital-contracts driven method for pricing complex derivatives

Abstract

In this paper, we propose a data and digital-contracts driven (DDCD) method for pricing various complex options. The DDCD method is a combination of nonparametric and parametric methods. In general, nonparametric data driven methods use observed data as training data of a learning network directly. Different from these, in the proposed DDCD method, some European-style digital contracts (DCs) of the underlying assets are added as auxiliary information to guide the learning process of the pricing formula. The DCs can be obtained by using the observed data according to parametric methods. Thus, the DCs are actually used as the hints of the pricing formula, and then the DDCD method has superior pricing accuracy to the common data driven method in practical applications. Some Monte Carlo simulation experiments are performed and the results demonstrate that the proposed method not only has the advantages of generalization and superior accuracy, as the non-parametric method has, but also has the property of robustness to financial data with noise, as the parametric method has.     

Simulation-based pricing of convertible bonds

We propose and empirically investigate a pricing model for convertible bonds based on Monte Carlo simulation. The method uses parametric representations of the early exercise decisions and consists of two stages. Pricing convertible bonds with the proposed Monte Carlo approach allows us to better capture both the dynamics of the underlying state variables and the rich set of real-world convertible bond specifications. Furthermore, using the simulation model proposed, we present an empirical pricing study of the US market, using 32 convertible bonds and 69 months of daily market prices. Our results do not confirm the evidence of previous studies that market prices of convertible bonds are on average lower than prices generated by a theoretical model. Similarly, our study is not supportive of a strong positive relationship between moneyness and mean pricing error, as argued in the literature.     

美式期权的Monte Carlo模拟法定价理论及应用

近年来随着计算机技术的飞速发展,美式期权的Monte Carlo模拟法定价取得了实质性的突破。本文分析介绍了美式期权的Monte Carlo模拟法定价理论及在此基础上推导出的线性回归MonteCarlo模拟法定价公式及其在实际的应用。     

Efficient and accurate quadratic approximation methods for pricing Asian strike options

This study is on valuing Asian strike options and presents efficient and accurate quadratic approximation methods that work extremely well, both with regard to the volatility of a wide range of underlying assets, and longer average time windows. We demonstrate that most of the well-known quadratic approximation methods used in the literature for pricing Asian strike options are special cases of our model, with the numerical results demonstrating that our method significantly outperforms the other quadratic approximation methods examined here. Using our method for the calculation of hundreds of Asian strike options, the pricing errors (in terms of the root mean square errors) are reasonably small. Compared with the Monte Carlo benchmark method, our method is shown to be rapid and accurate. We further extend our method to the valuing of quanto forward-starting Asian strike options, with the pricing accuracy of these options being largely the same as the pricing of plain vanilla Asian strike options.     

Smoothing the payoff for efficient computation of Basket option prices

We consider the problem of pricing basket options in a multivariate Black–Scholes or Variance-Gamma model. From a numerical point of view, pricing such options corresponds to moderate and high-dimensional numerical integration problems with non-smooth integrands. Due to this lack of regularity, higher order numerical integration techniques may not be directly available, requiring the use of methods like Monte Carlo specifically designed to work for non-regular problems. We propose to use the inherent smoothing property of the density of the underlying in the above models to mollify the payoff function by means of an exact conditional expectation. The resulting conditional expectation is unbiased and yields a smooth integrand, which is amenable to the efficient use of adaptive sparse-grid cubature. Numerical examples indicate that the high-order method may perform orders of magnitude faster than Monte Carlo or Quasi Monte Carlo methods in dimensions up to 35.     

Efficient willow tree method for European-style and American-style moving average barrier options pricing

Moving average options are widely traded in financial markets, but exiting methods for pricing this type of option are too slow. This paper proposes two efficient willow tree methods for pricing European-style and American-style moving average barrier options (MABOs). We first solve the finite-dimensional partial differential equation model for discretely monitored MABOs by willow tree methods, and then compute the value of continuously monitored MABOs by Richardson’s two-point extrapolation. Our new willow tree method employs the interpolation error minimization technique to reduce complexity. The corresponding convergence rate and error bounds are also analyzed. It shows that our proposed methods can provide the same accuracy as the binomial tree approach and Monte Carlo simulation, but require much less computing time. The numerical experiments support our claims.     

Option pricing for GARCH-type models with generalized hyperbolic innovations

In this paper, we provide a new dynamic asset pricing model for plain vanilla options and we discuss its ability to produce minimum mispricing errors on equity option books. Given the historical measure, the dynamics of assets being modeled by Garch-type models with generalized hyperbolic innovations and the pricing kernel is an exponential affine function of the state variables, we show that the risk-neutral distribution is unique and again implies a generalized hyperbolic dynamics with changed parameters. We provide an empirical test for our pricing methodology on two data sets of options, respectively written on the French CAC 40 and the American SP 500. Then, using our theoretical result associated with Monte Carlo simulations, we compare this approach with natural competitors in order to test its efficiency. More generally, our empirical investigations analyse the ability of specific parametric innovations to reproduce market prices in the context of an exponential affine specification of the stochastic discount factor.     

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.     

On an efficient multiple time step Monte Carlo simulation of the SABR model

In this paper, we will present a multiple time step Monte Carlo simulation technique for pricing options under the Stochastic Alpha Beta Rho model. The proposed method is an extension of the one time step Monte Carlo method that we proposed in an accompanying paper Leitao et al. [Appl. Math. Comput. 2017, 293, 461–479], for pricing European options in the context of the model calibration. A highly efficient method results, with many very interesting and nontrivial components, like Fourier inversion for the sum of log-normals, stochastic collocation, Gumbel copula, correlation approximation, that are not yet seen in combination within a Monte Carlo simulation. The present multiple time step Monte Carlo method is especially useful for long-term options and for exotic options.     

Pricing exotic options using MSL-MC

Today, better numerical approximations are required for multi-dimensional SDEs to improve on the poor performance of the standard Monte Carlo pricing method. With this aim in mind, this paper presents a method (MSL-MC) to price exotic options using multi-dimensional SDEs (e.g. stochastic volatility models). Usually, it is the weak convergence property of numerical discretizations that is most important, because, in financial applications, one is mostly concerned with the accurate estimation of expected payoffs. However, in the recently developed Multilevel Monte Carlo path simulation method (ML-MC), the strong convergence property plays a crucial role. We present a modification to the ML-MC algorithm that can be used to achieve better savings. To illustrate these, various examples of exotic options are given using a wide variety of payoffs, stochastic volatility models and the new Multischeme Multilevel Monte Carlo method (MSL-MC). For standard payoffs, both European and Digital options are presented. Examples are also given for complex payoffs, such as combinations of European options (Butterfly Spread, Strip and Strap options). Finally, for path-dependent payoffs, both Asian and Variance Swap options are demonstrated. This research shows how the use of stochastic volatility models and the θ scheme can improve the convergence of the MSL-MC so that the computational cost to achieve an accuracy of O(ε) is reduced from O(ε^?3) to O(ε^?2) for a payoff under global and non-global Lipschitz conditions.     

Dual pricing of multi-exercise options under volume constraints

In this paper, we study the pricing problem of multi-exercise options under volume constraints. The volume constraint is modelled by an adapted process with values in the positive integers, which describes the maximal number of rights to be exercised at a given time. We derive a representation of the marginal value of an additional nth right as a standard single stopping problem with a modified cash-flow process. This representation then leads to a dual pricing formula, which generalizes a result by Meinshausen and Hambly (Math. Finance 14:557–583, 2004) from the standard multi-exercise option (with at most one right per time step) to general constraints. We also state an explicit Monte Carlo algorithm for computing confidence intervals for the price of multi-exercise options under volume constraints and present numerical results for the pricing of a swing contract in an electricity market.     

A generalization of reset options pricing formulae with stochastic interest rates

A generalization of reset call options with predetermined dates is derived in the case of time-dependent volatility and time-dependent interest rate by applying martingale method and change of nume?aire or change of probability measure. An analytical pricing formula for the reset call option is also obtained when the interest rate follows an extended Vasicek’s model. Numerical results show that the correlated coefficient between the stock price and interest rate is almost unacted on the price of reset call option with short maturity and Monte Carlo method is inefficient. Monte Carlo method should be only used if there is no closed-formed solution for option pricing.     

Nonparametric Estimation of State-Price Densities Implicit in Financial Asset Prices

Implicit in the prices of traded financial assets are Arrow–Debreu prices or, with continuous states, the state-price density (SPD). We construct a nonparametric estimator for the SPD implicit in option prices and we derive its asymptotic sampling theory. This estimator provides an arbitrage-free method of pricing new, complex, or illiquid securities while capturing those features of the data that are most relevant from an asset-pricing perspective, for example, negative skewness and excess kurtosis for asset returns, and volatility "smiles" for option prices. We perform Monte Carlo experiments and extract the SPD from actual S&P 500 option prices.     

Option Pricing Under Autoregressive Random Variance Models

Abstract

The autoregressive random variance (ARV) model introduced by Taylor (1980, 1982, 1986) is a popular version of stochastic volatility (SV) models and a discrete-time simplification of the continuous-time diffusion SV models. This paper introduces a valuation model for options under a discrete-time ARV model with general stock and volatility innovations. It employs the discretetime version of the Esscher transform to determine an equivalent martingale measure under an incomplete market. Various parametric cases of the ARV models, are considered, namely, the log-normal ARV models, the jump-type Poisson ARV models, and the gamma ARV models, and more explicit pricing formulas of a European call option under these parametric cases are provided. A Monte Carlo experiment for some parametric cases is also conducted.     

Conditional beta pricing models: A nonparametric approach

We propose a two-stage procedure to estimate conditional beta pricing models that allows for flexibility in the dynamics of asset betas and market prices of risk (MPR). First, conditional betas are estimated nonparametrically for each asset and period using the time-series of previous data. Then, time-varying MPR are estimated from the cross-section of returns and betas. We prove the consistency and asymptotic normality of the estimators. We also perform Monte Carlo simulations for the conditional version of the three-factor model of Fama and French (1993) and show that nonparametrically estimated betas outperform rolling betas under different specifications of beta dynamics. Using return data on the 25 size and book-to-market sorted portfolios, we find that the nonparametric procedure produces a better fit of the three-factor model to the data, less biased estimates of MPR and lower pricing errors than the Fama–MacBeth procedure with betas estimated under several alternative parametric specifications.     

Regression-based algorithms for life insurance contracts with surrender guarantees

We present a general framework for pricing life insurance contracts embedding a surrender option. The model allows for several sources of risk, such as uncertainty in mortality, interest rates and other financial factors. We describe and compare two numerical schemes based on the Least Squares Monte Carlo method, emphasizing underlying modeling assumptions and computational issues.     

Multilevel Monte Carlo for exponential Lévy models

We apply the multilevel Monte Carlo method for option pricing problems using exponential Lévy models with a uniform timestep discretisation. For lookback and barrier options, we derive estimates of the convergence rate of the error introduced by the discrete monitoring of the running supremum of a broad class of Lévy processes. We then use these to obtain upper bounds on the multilevel Monte Carlo variance convergence rate for the variance gamma, NIG and \(\alpha\)-stable processes. We also provide an analysis of a trapezoidal approximation for Asian options. Our method is illustrated by numerical experiments.     

Efficient willow tree method for variable annuities valuation and risk management☆

Variable annuities (VAs) with various guarantees are popular retirement products in the past decades. However, due to the sophistication of the embedded guarantees, most existing methods only focus on the one of embedded guarantees underlying one specified stochastic model. The method to evaluate VAs with all guarantees and manage its risk is very limited, except for the Monte Carlo method. In this paper, we propose an efficient willow tree method to evaluate VAs embedded with all popular guarantees on the market underlying various stochastic models. Moreover, our tree structure is also applicable to compute dollar delta, value at risk (VaR) and conditional tail expectation (CTE) in hedging and risk-based capital calculation. Numerical experiments demonstrate the accuracy and efficiency of our method in pricing and managing the risk of VAs.     

Multilevel dual approach for pricing American style derivatives

In this article we propose a novel approach to reduce the computational complexity of the dual method for pricing American options. We consider a sequence of martingales that converges to a given target martingale and decompose the original dual representation into a sum of representations that correspond to different levels of approximation to the target martingale. By next replacing in each representation true conditional expectations with their Monte Carlo estimates, we arrive at what one may call a multilevel dual Monte Carlo algorithm. The analysis of this algorithm reveals that the computational complexity of getting the corresponding target upper bound, due to the target martingale, can be significantly reduced. In particular, it turns out that using our new approach, we may construct a multilevel version of the well-known nested Monte Carlo algorithm of Andersen and Broadie (Manag. Sci. 50:1222–1234, 2004) that is, regarding complexity, virtually equivalent to a non-nested algorithm. The performance of this multilevel algorithm is illustrated by a numerical example.     

Rapid and accurate development of prices and Greeks for nth to default credit swaps in the Li model

New techniques are introduced for pricing nth to default credit swaps in the Li model. We demonstrate the use of importance sampling to greatly increase the rate of convergence of Monte Carlo simulations for pricing. This technique is combined with the likelihood ratio and pathwise methods for computing the sensitivities of these products to changes in the hazard rates of the underlying obligors. In particular the extension of the pathwise method has wider significance in that it is shown that the method can be used even when the pay-off is discontinuous.