Showing papers in "arXiv: Computational Finance in 2012"

Online Portfolio Selection: A Survey

Abstract: Online portfolio selection is a fundamental problem in computational finance, which has been extensively studied across several research communities, including finance, statistics, artificial intelligence, machine learning, and data mining, etc. This article aims to provide a comprehensive survey and a structural understanding of published online portfolio selection techniques. From an online machine learning perspective, we first formulate online portfolio selection as a sequential decision problem, and then survey a variety of state-of-the-art approaches, which are grouped into several major categories, including benchmarks, "Follow-the-Winner" approaches, "Follow-the-Loser" approaches, "Pattern-Matching" based approaches, and "Meta-Learning Algorithms". In addition to the problem formulation and related algorithms, we also discuss the relationship of these algorithms with the Capital Growth theory in order to better understand the similarities and differences of their underlying trading ideas. This article aims to provide a timely and comprehensive survey for both machine learning and data mining researchers in academia and quantitative portfolio managers in the financial industry to help them understand the state-of-the-art and facilitate their research and practical applications. We also discuss some open issues and evaluate some emerging new trends for future research directions.

A new approach to unbiased estimation for SDE's

Abstract: In this paper, we introduce a new approach to constructing unbiased estimators when computing expectations of path functionals associated with stochastic differential equations (SDEs). Our randomization idea is closely related to multi-level Monte Carlo and provides a simple mechanism for constructing a finite variance unbiased estimator with "square root convergence rate" whenever one has available a scheme that produces strong error of order greater than 1/2 for the path functional under consideration.

Antithetic multilevel Monte Carlo estimation for multi-dimensional SDEs without L\'{e}vy area simulation

Abstract: In this paper we introduce a new multilevel Monte Carlo (MLMC) estimator for multi-dimensional SDEs driven by Brownian motions Giles has previously shown that if we combine a numerical approximation with strong order of convergence $O(\Delta t)$ with MLMC we can reduce the computational complexity to estimate expected values of functionals of SDE solutions with a root-mean-square error of $\epsilon$ from $O(\epsilon^{-3})$ to $O(\epsilon^{-2})$ However, in general, to obtain a rate of strong convergence higher than $O(\Delta t^{1/2})$ requires simulation, or approximation, of L\'{e}vy areas In this paper, through the construction of a suitable antithetic multilevel correction estimator, we are able to avoid the simulation of L\'{e}vy areas and still achieve an $O(\Delta t^2)$ multilevel correction variance for smooth payoffs, and almost an $O(\Delta t^{3/2})$ variance for piecewise smooth payoffs, even though there is only $O(\Delta t^{1/2})$ strong convergence This results in an $O(\epsilon^{-2})$ complexity for estimating the value of European and Asian put and call options

Multilevel Monte Carlo methods for applications in finance

Abstract: Since Giles introduced the multilevel Monte Carlo path simulation method [18], there has been rapid development of the technique for a variety of applications in computational finance. This paper surveys the progress so far, highlights the key features in achieving a high rate of multilevel variance convergence, and suggests directions for future research.

First order strong approximations of scalar SDEs with values in a domain

Abstract: We are interested in strong approximations of one-dimensional SDEs which have non-Lipschitz coefficients and which take values in a domain. Under a set of general assumptions we derive an implicit scheme that preserves the domain of the SDEs and is strongly convergent with rate one. Moreover, we show that this general result can be applied to many SDEs we encounter in mathematical finance and bio-mathematics. We will demonstrate flexibility of our approach by analysing classical examples of SDEs with sublinear coefficients (CIR, CEV models and Wright-Fisher diffusion) and also with superlinear coefficients (3/2-volatility, Ait-Sahalia model). Our goal is to justify an efficient Multi-Level Monte Carlo (MLMC) method for a rich family of SDEs, which relies on good strong convergence properties.

Minimax Option Pricing Meets Black-Scholes in the Limit

Abstract: Option contracts are a type of financial derivative that allow investors to hedge risk and speculate on the variation of an asset's future market price. In short, an option has a particular payout that is based on the market price for an asset on a given date in the future. In 1973, Black and Scholes proposed a valuation model for options that essentially estimates the tail risk of the asset price under the assumption that the price will fluctuate according to geometric Brownian motion. More recently, DeMarzo et al., among others, have proposed more robust valuation schemes, where we can even assume an adversary chooses the price fluctuations. This framework can be considered as a sequential two-player zero-sum game between the investor and Nature. We analyze the value of this game in the limit, where the investor can trade at smaller and smaller time intervals. Under weak assumptions on the actions of Nature (an adversary), we show that the minimax option price asymptotically approaches exactly the Black-Scholes valuation. The key piece of our analysis is showing that Nature's minimax optimal dual strategy converges to geometric Brownian motion in the limit.

Error estimates for binomial approximations of game put options

Abstract: We construct algorithms via binomial approximations for computation of prices of game put options and obtain estimates of approximation errors.

From characteristic functions to implied volatility expansions

Abstract: For any strictly positive martingale $S = \exp(X)$ for which $X$ has a characteristic function, we provide an expansion for the implied volatility. This expansion is explicit in the sense that it involves no integrals, but only polynomials in the log strike. We illustrate the versatility of our expansion by computing the approximate implied volatility smile in three well-known martingale models: one finite activity exponential L\'evy model (Merton), one infinite activity exponential L\'evy model (Variance Gamma), and one stochastic volatility model (Heston). Finally, we illustrate how our expansion can be used to perform a model-free calibration of the empirically observed implied volatility surface.

Perturbative Expansion Technique for Non-linear FBSDEs with Interacting Particle Method

Abstract: In this paper, we propose an efficient Monte Carlo implementation of non-linear FBSDEs as a system of interacting particles inspired by the ideas of branching diffusion method. It will be particularly useful to investigate large and complex systems, and hence it is a good complement of our previous work presenting an analytical perturbation procedure for generic non-linear FBSDEs. There appear multiple species of particles, where the first one follows the diffusion of the original underlying state, and the others the Malliavin derivatives with a grading structure. The number of branching points are capped by the order of perturbation, which is expected to make the scheme less numerically intensive. The proposed method can be applied to semi-linear problems, such as American and Bermudan options, Credit Value Adjustment (CVA), and even fully non-linear issues, such as the optimal portfolio problems in incomplete and/or constrained markets, feedbacks from large investors, and also the analysis of various risk measures.

Second Order Multiscale Stochastic Volatility Asymptotics: Stochastic Terminal Layer Analysis & Calibration

Abstract: Multiscale stochastic volatility models have been developed as an efficient way to capture the principle effects on derivative pricing and portfolio optimization of randomly varying volatility. The recent book Fouque, Papanicolaou, Sircar and Solna (2011, CUP) analyzes models in which the volatility of the underlying is driven by two diffusions -- one fast mean-reverting and one slow-varying, and provides a first order approximation for European option prices and for the implied volatility surface, which is calibrated to market data. Here, we present the full second order asymptotics, which are considerably more complicated due to a terminal layer near the option expiration time. We find that, to second order, the implied volatility approximation depends quadratically on log-moneyness, capturing the convexity of the implied volatility curve seen in data. We introduce a new probabilistic approach to the terminal layer analysis needed for the derivation of the second order singular perturbation term, and calibrate to S&P 500 options data.

Counterparty Risk Valuation: A Marked Branching Diffusion Approach

Abstract: The purpose of this paper is to design an algorithm for the computation of the counterparty risk which is competitive in regards of a brute force "Monte-Carlo of Monte-Carlo" method (with nested simulations). This is achieved using marked branching diffusions describing a Galton-Watson random tree. Such an algorithm leads at the same time to a computation of the (bilateral) counterparty risk when we use the default-risky or counterparty-riskless option values as mark-to-market. Our method is illustrated by various numerical examples.

Maximum likelihood approach for several stochastic volatility models

Abstract: Volatility measures the amplitude of price fluctuations. Despite it is one of the most important quantities in finance, volatility is not directly observable. Here we apply a maximum likelihood method which assumes that price and volatility follow a two-dimensional diffusion process where volatility is the stochastic diffusion coefficient of the log-price dynamics. We apply this method to the simplest versions of the expOU, the OU and the Heston stochastic volatility models and we study their performance in terms of the log-price probability, the volatility probability, and its Mean First-Passage Time. The approach has some predictive power on the future returns amplitude by only knowing current volatility. The assumed models do not consider long-range volatility auto-correlation and the asymmetric return-volatility cross-correlation but the method still arises very naturally these two important stylized facts. We apply the method to different market indexes and with a good performance in all cases.

Fast computation of vanilla prices in time-changed models and implied volatilities using rational approximations

Abstract: We present a new numerical method to price vanilla options quickly in time-changed Brownian motion models. The method is based on rational function approximations of the Black-Scholes formula. Detailed numerical results are given for a number of widely used models. In particular, we use the variance-gamma model, the CGMY model and the Heston model without correlation to illustrate our results. Comparison to the standard fast Fourier transform method with respect to accuracy and speed appears to favour the newly developed method in the cases considered. We present error estimates for the option prices. Additionally, we use this method to derive a procedure to compute, for a given set of arbitrage-free European call option prices, the corresponding Black-Scholes implied volatility surface. To achieve this, rational function approximations of the inverse of the Black-Scholes formula are used. We are thus able to work out implied volatilities more efficiently than one can by the use of other common methods. Error estimates are presented for a wide range of parameters.

An FBSDE Approach to American Option Pricing with an Interacting Particle Method

Abstract: In the paper, we propose a new calculation scheme for American options in the framework of a forward backward stochastic differential equation (FBSDE). The well-known decomposition of an American option price with that of a European option of the same maturity and the remaining early exercise premium can be cast into the form of a decoupled non-linear FBSDE. We numerically solve the FBSDE by applying an interacting particle method recently proposed by Fujii and Takahashi (2012d), which allows one to perform a Monte Carlo simulation in a fully forward-looking manner. We perform the fourth-order analysis for the Black-Scholes (BS) model and the third-order analysis for the Heston model. The comparison to those obtained from existing tree algorithms shows the effectiveness of the particle method.

Quasi-Monte Carlo methods for the Heston model

Abstract: In this paper, we discuss the application of quasi-Monte Carlo methods to the Heston model. We base our algorithms on the Broadie-Kaya algorithm, an exact simulation scheme for the Heston model. As the joint transition densities are not available in closed-form, the Linear Transformation method due to Imai and Tan, a popular and widely applicable method to improve the effectiveness of quasi-Monte Carlo methods, cannot be employed in the context of path-dependent options when the underlying price process follows the Heston model. Consequently, we tailor quasi-Monte Carlo methods directly to the Heston model. The contributions of the paper are threefold: We firstly show how to apply quasi-Monte Carlo methods in the context of the Heston model and the SVJ model, secondly that quasi-Monte Carlo methods improve on Monte Carlo methods, and thirdly how to improve the effectiveness of quasi-Monte Carlo methods by using bridge constructions tailored to the Heston and SVJ models. Finally, we provide some extensions for computing greeks, barrier options, multidimensional and multi-asset pricing, and the 3/2 model.

An introduction to particle integration methods: with applications to risk and insurance

Abstract: Interacting particle methods are increasingly used to sample from complex and high-dimensional distributions. These stochastic particle integration techniques can be interpreted as an universal acceptance-rejection sequential particle sampler equipped with adaptive and interacting recycling mechanisms. Practically, the particles evolve randomly around the space independently and to each particle is associated a positive potential function. Periodically, particles with high potentials duplicate at the expense of low potential particle which die. This natural genetic type selection scheme appears in numerous applications in applied probability, physics, Bayesian statistics, signal processing, biology, and information engineering. It is the intention of this paper to introduce them to risk modeling. From a purely mathematical point of view, these stochastic samplers can be interpreted as Feynman-Kac particle integration methods. These functional models are natural mathematical extensions of the traditional change of probability measures, commonly used to design an importance sampling strategy. In this article, we provide a brief introduction to the stochastic modeling and the theoretical analysis of these particle algorithms. Then we conclude with an illustration of a subset of such methods to resolve important risk measure and capital estimation in risk and insurance modelling.

The solution of discretionary stopping problems with applications to the optimal timing of investment decisions

Abstract: We present a methodology for obtaining explicit solutions to infinite time horizon optimal stopping problems involving general, one-dimensional, Ito diffusions, payoff functions that need not be smooth and state-dependent discounting. This is done within a framework based on dynamic programming techniques employing variational inequalities and links to the probabilistic approaches employing $r$-excessive functions and martingale theory. The aim of this paper is to facilitate the the solution of a wide variety of problems, particularly in finance or economics.

Strong Convergence for Euler-Maruyama and Milstein Schemes with Asymptotic Method

Abstract: Motivated by weak convergence results in the paper of Takahashi and Yoshida (2005), we show strong convergence for an accelerated Euler-Maruyama scheme applied to perturbed stochastic differential equations. The Milstein scheme with the same acceleration is also discussed as an extended result. The theoretical results can be applied to analyzing the multi-level Monte Carlo method originally developed by M.B. Giles. Several numerical experiments for the SABR stochastic volatility model are presented in order to confirm the efficiency of the schemes.

Numerical Valuation of Derivatives in High-Dimensional Settings via PDE Expansions

Abstract: In this article, we propose a new numerical approach to high-dimensional partial differential equations (PDEs) arising in the valuation of exotic derivative securities. The proposed method is extended from Reisinger and Wittum (2007) and uses principal component analysis (PCA) of the underlying process in combination with a Taylor expansion of the value function into solutions to low-dimensional PDEs. The approximation is related to anchored analysis of variance (ANOVA) decompositions and is expected to be accurate whenever the covariance matrix has one or few dominating eigenvalues. A main purpose of the present article is to give a careful analysis of the numerical accuracy and computational complexity compared to state-of-the-art Monte Carlo methods on the example of Bermudan swaptions and Ratchet floors, which are considered difficult benchmark problems. We are able to demonstrate that for problems with medium to high dimensionality and moderate time horizons the presented PDE method delivers results comparable in accuracy to the MC methods considered here in similar or (often significantly) faster runtime.

Quantile Mechanics 3: Series Representations and Approximation of some Quantile Functions appearing in Finance

Abstract: It has long been agreed by academics that the inversion method is the method of choice for generating random variates, given the availability of the quantile function However for several probability distributions arising in practice a satisfactory method of approximating these functions is not available The main focus of this paper will be to develop Taylor and asymptotic series expansions for the quantile functions belonging to the following probability distributions; Variance Gamma, Generalized Inverse Gaussian, Hyperbolic and alpha-Stable As a secondary matter, based on these analytic expressions we briefly investigate the problem of approximating the quantile function

Estimation of the shape parameter of a generalized Pareto distribution based on a transformation to Pareto distributed variables

Abstract: Random variables of the generalized Pareto distribution, can be transformed to that of the Pareto distribution. Explicit expressions exist for the maximum likelihood estimators of the parameters of the Pareto distribution. The performance of the estimation of the shape parameter of generalized Pareto distributed using transformed observations, based on the probability weighted method is tested. It was found to improve the performance of the probability weighted estimator and performs good with respect to bias and MSE.

Conditional sampling for barrier option pricing under the Heston model

Abstract: We propose a quasi-Monte Carlo algorithm for pricing knock-out and knock-in barrier options under the Heston (1993) stochastic volatility model. This is done by modifying the LT method from Imai and Tan (2006) for the Heston model such that the first uniform variable does not influence the stochastic volatility path and then conditionally modifying its marginals to fulfill the barrier condition(s). We show this method is unbiased and never does worse than the unconditional algorithm. Additionally the conditioning is combined with a root finding method to also force positive payouts. The effectiveness of this method is shown by extensive numerical results.

The Smile of certain L\'evy-type Models

Abstract: We consider a class of assets whose risk-neutral pricing dynamics are described by an exponential L\'evy-type process subject to default. The class of processes we consider features locally-dependent drift, diffusion and default-intensity as well as a locally-dependent L\'evy measure. Using techniques from regular perturbation theory and Fourier analysis, we derive a series expansion for the price of a European-style option. We also provide precise conditions under which this series expansion converges to the exact price. Additionally, for a certain subclass of assets in our modeling framework, we derive an expansion for the implied volatility induced by our option pricing formula. The implied volatility expansion is exact within its radius of convergence. As an example of our framework, we propose a class of CEV-like L\'evy-type models. Within this class, approximate option prices can be computed by a single Fourier integral and approximate implied volatilities are explicit (i.e., no integration is required). Furthermore, the class of CEV-like L\'evy-type models is shown to provide a tight fit to the implied volatility surface of S{&}P500 index options.

On a Symmetrization of Diffusion Processes

Abstract: The latter author, together with collaborators, proposed a numerical scheme to calculate the price of barrier options. The scheme is based on a symmetrization of diffusion process. The present paper aims to give a mathematical credit to the use of the numerical scheme for Heston or SABR type stochastic volatility models. This will be done by showing a fairly general result on the symmetrization (in multi-dimension/multi-reflections). Further applications (to time-inhomogeneous diffusions/ to time dependent boundaries/to curved boundaries) are also discussed.

Perturbative Expansion of FBSDE in an Incomplete Market with Stochastic Volatility

Abstract: In this work, we apply our newly proposed perturbative expansion technique to a quadratic growth FBSDE appearing in an incomplete market with stochastic volatility that is not perfectly hedgeable. By combining standard asymptotic expansion technique for the underlying volatility process, we derive explicit expression for the solution of the FBSDE up to the third order of volatility-of-volatility, which can be directly translated into the optimal investment strategy. We compare our approximation with the exact solution, which is known to be derived by the Cole-Hopf transformation in this popular setup. The result is very encouraging and shows good accuracy of the approximation up to quite long maturities. Since our new methodology can be extended straightforwardly to multi-dimensional setups, we expect it will open real possibilities to obtain explicit optimal portfolios or hedging strategies under realistic assumptions.

Approximating stochastic volatility by recombinant trees

Abstract: A general method to construct recombinant tree approximations for stochastic volatility models is developed and applied to the Heston model for stock price dynamics. In this application, the resulting approximation is a four tuple Markov process. The first two components are related to the stock and volatility processes and take values in a two-dimensional binomial tree. The other two components of the Markov process are the increments of random walks with simple values in $\{-1,+1\}$. The resulting efficient option pricing equations are numerically implemented for general American and European options including the standard put and calls, barrier, lookback and Asian-type pay-offs. The weak and extended weak convergences are also proved.

A No-Arbitrage Model of Liquidity in Financial Markets involving Brownian Sheets

Abstract: We consider a dynamic market model where buyers and sellers submit limit orders. If at a given moment in time, the buyer is unable to complete his entire order due to the shortage of sell orders at the required limit price, the unmatched part of the order is recorded in the order book. Subsequently these buy unmatched orders may be matched with new incoming sell orders. The resulting demand curve constitutes the sole input to our model. The clearing price is then mechanically calculated using the market clearing condition. We use a Brownian sheet to model the demand curve, and provide some theoretical assumptions under which such a model is justified. Our main result is the proof that if there exists a unique equivalent martingale measure for the clearing price, then under some mild assumptions there is no arbitrage. We use the Ito- Wentzell formula to obtain that result, and also to characterize the dynamics of the demand curve and of the clearing price in the equivalent measure. We find that the volatility of the clearing price is (up to a stochastic factor) inversely proportional to the sum of buy and sell order flow density (evaluated at the clearing price), which confirms the intuition that volatility is inversely proportional to volume. We also demonstrate that our approach is implementable. We use real order book data and simulate option prices under a particularly simple parameterization of our model. The no-arbitrage conditions we obtain are applicable to a wide class of models, in the same way that the Heath-Jarrow-Morton conditions apply to a wide class of interest rate models.

Sparsifying Defaults: Optimal Bailout Policies for Financial Networks in Distress

Abstract: The events of the last few years revealed an acute need for tools to systematically model and analyze large financial networks. Many applications of such tools include the forecasting of systemic failures and analyzing probable effects of economic policy decisions. We consider optimizing the amount and structure of a bailout in a borrower-lender network: Given a fixed amount of cash to be injected into the system, how should it be distributed among the nodes in order to achieve the smallest overall amount of unpaid liabilities or the smallest number of nodes in default? We develop an exact algorithm for the problem of minimizing the amount of unpaid liabilities, by showing that it is equivalent to a linear program. For the problem of minimizing the number of defaults, we develop an approximate algorithm using a reweighted l1 minimization approach. We illustrate this algorithm using an example with synthetic data for which the optimal solution can be calculated exactly, and show through numerical simulation that the solutions calculated by our algorithm are close to optimal.

The Exact Smile of some Local Volatility Models

Abstract: We introduce a new class of local volatility models. Within this framework, we obtain expressions for both (i) the price of any European option and (ii) the induced implied volatility smile. As an illustration of our framework, we perform specific pricing and implied volatility computations for a CEV-like example. Numerical examples are provided.

Dynamic Hedging Using Generated Genetic Programming Implied Volatility Models

Abstract: The purpose of this paper is to improve the accuracy of dynamic hedging using implied volatilities generated by genetic programming. Using real data from S&P500 index options, the genetic programming's ability to forecast Black and Scholes implied volatility is compared between static and dynamic training-subset selection methods. The performance of the best generated GP implied volatilities is tested in dynamic hedging and compared with Black-Scholes model. Based on MSE total, the dynamic training of GP yields better results than those obtained from static training with fixed samples. According to hedging errors, the GP model is more accurate almost in all hedging strategies than the BS model, particularly for in-the-money call options and at-the-money put options.