From Characteristic Functions and Fourier Transforms to PDFs/CDFs and Option Prices
|
|
- Norma Ball
- 5 years ago
- Views:
Transcription
1 From Characteristic Functions and Fourier Transforms to PDFs/CDFs and Option Prices Liuren Wu Zicklin School of Business, Baruch College Fall, 2007 Liuren Wu Fourier Transforms Option Pricing, Fall, / 21
2 From a dynamic process to a terminal distribution If we specify a dynamic process for X and its initial conditions, we can compute the distribution of X over a certain time horizon, [t, T ]. The reverse is not true. Given a conditional distribution for X over [t, T ], there can be many processes that can generate this distribution There are many routes to Rome. Knowing the dynamic process is important for hedging practices Dynamic hedging is more likely to work if the underlying process is continuous instead of discontinuous. The objective of this note: Given a dynamic process for X, derive the probability density (PDF) of X and price options based on integration of terminal payoffs over the probability densities (instead of based on dynamic hedging arguments and PDEs). Liuren Wu Fourier Transforms Option Pricing, Fall, / 21
3 Formalizing the idea Assume X is a Markov process: Any information about X up to time t is summarized by X t. Let f (X T X t ) denote the conditional distribution of X T conditional on time-t information (filtration F t ) under the risk-neutral measure Q. Let Π(X T ) denote the payoff of a contingent claim (derivative), which we assume is a function of X T. Then, its time-t value is p t = E Q t [e R ] T r t s ds Π(X T ) = e R T r t s ds Π(X T )f (X T X t )dx T Instead of solving PDEs, binomial trees, or simulating the process, we focus on doing the integration based on the risk-neutral densities. In most of the examples, I assume deterministic interest rates: p t = e rτ Π(X T )f (X T X t )dx T X where r now is the time-t continuously compounded spot rate of maturity T t. X Liuren Wu Fourier Transforms Option Pricing, Fall, / 21
4 The road map For most models (processes) discussed in this class on security prices S t, we can derive the characteristic function or Fourier transform of the log security return (ln S T /S t ) (semi-)analytically. Given the characteristic function (CF), we just need one numerical integration to obtain the probability density function (PDF) or cumulative density function (CDF). Given the Fourier transforms (FT), we just need one numerical integration to obtain the value of vanilla options. The integration is one-dimensional in both cases no matter how many dimensions/factors the security price S t is composed of. We can apply fast Fourier inversion (FFT) to make the numerical integration (for both PDF and option prices) very fast. Solving PDEs becomes difficult as the dimension increases. Simulation works (slowly) for high-dimensional cases and is reserved for pricing exotics. Liuren Wu Fourier Transforms Option Pricing, Fall, / 21
5 Characteristic function: definition In probability theory, the characteristic function (CF) of any random variable X completely defines its probability distribution. On the real line it is given by the following formula: φ X (u) E [ e iux ] = e iux f X (x)dx = e iux df X (x), u R where u is a real number, i is the imaginary unit, and E denotes the expected value, f X (x) denotes the probability density function (PDF), and F X (x) denotes the cumulative density function (CDF). CF is well defined on the whole real line (u). For option pricing, we extend the definition to the complex plane, u D C, where D denotes the subset of the complex plane on which the expectation is well defined. φ X (u) under this extended definition is called the generalized Fourier transform. Ω Liuren Wu Fourier Transforms Option Pricing, Fall, / 21
6 Example: The Black-Scholes model Under BSM, the log security return follows, s t ln S t /S 0 = (µ 12 ) σ2 t + σw t The return is normally distributed with mean ( µ 1 2 σ2) t and variance σ 2 t. ( ) The PDF is f s (x) = 1 2πσ 2 t exp (x (µ 1 2 σ2 )t) 2 2σ 2 t. The CF is: φ s (u) = E [ ] e iust = e iumean 1 2 u2variance = e (iuµt iu 1 2 σ2 t 1 2 u2 σ 2 t) = e (iuµt 1 2 σ2 (iu+u 2 )t) Under Q, µ = r q. Liuren Wu Fourier Transforms Option Pricing, Fall, / 21
7 The inversion: From CF to PDF and CDF There is a bijection between CDF and CFs: Two distinct probability distributions never share the same CF. Given a CF φ, it is possible to reconstruct the corresponding CDF: 1 F X (y) F X (x) = lim τ 2π +τ In general this is an improper integral... Another form of the inversion F X (x) = π The inversion formula for PDF: f X (x) = τ e iux e iuy φ X (u)du iu e iux φ X ( u) e iux φ X (u) du. iu e iux φ X (u)du = 1 π 0 e iux φ X (u)du. All the integrals here should be understood as a principal value if there is no separate convergence at the limits. Exercise: Based on the BSM model, try a simple discretization of the integral and see how close you can get to the analytical soln. Can you think of some better ways to improve the accuracy and speed of the numerical integration? Liuren Wu Fourier Transforms Option Pricing, Fall, / 21
8 Proofs Preliminary results: For fixed x, e iux = cos ux + i sin ux, 1 R e iuζ π iu du = π 1 R sin uζ du = sgn(ζ), u 2 R 0 sin uζ R du = sgn(ζ), π u sgn(y x)df (y) = R x df (y) + R x df (y) = 1 2F (x), φ(u)&φ( u) are complex conjugates. CDF inversion: I = Hence, F (x) = π I. PDF inversion: = = = e iux φ X ( u) e iux φ X (u) du 0 iu e iux e iuz e iux e iuz df (z)du 0 iu 2 sin u(x z) df (z)du 0 u Z 2 sin u(x z) dudf (z) = πsgn(x z)df (z) 0 u = π(2f (x) 1) f (x) = F (x) = 1 2π 0 e iux φ ( u) + e iux φ (u) du = 1 π e iux φ (u) du 0 Liuren Wu Fourier Transforms Option Pricing, Fall, / 21
9 Fourier transforms and inversions of European options Take a European call option as an example. We perform the following rescaling and change of variables: c(k) = c(k, t)/s 0 = e rt E Q [ 0 (e s t e k ] )1 st k. We represent the option price as percentage of the spot, and we quote the price in terms of log strike instead of strike. We can derive the Fourier transform of the call option in terms of the Fourier transform (CF) of the log return ln S t /S 0. Hence, if we know the CF of the return, we would know the transform of the option. Then, we can use numerical inversion to obtain option prices directly. There are two ways of doing this. Liuren Wu Fourier Transforms Option Pricing, Fall, / 21
10 I. The CDF analog Treat c (k) = e rt E Q 0 as a CDF. The option transform: χ I c(u) [( e s t e k) 1 st k] = e rt ( e s t e k) 1 st xdf (s) e iuk dc(k) = e rt φ s (u i) iu + 1, u R. Thus, if we know the CF of the return, φ s (u), we know the transform of the option, χ I c(u). The inversion formula is analogous to the inversion of a CDF: c (x) = 1 2 e qt + 1 2π 0 e iux χ I c ( u) e iux χ I c (u) du. iu Note the slight difference from the CDF inversion. Use quadrature methods for the numerical integration. The literature often writes: c (x) = e qt Q1 (x) e rt e x Q 2 (x). Then, we must invert twice. References: Duffie, Pan, Singleton, 2000, Transform Analysis and Asset Pricing for Affine Jump Diffusions, Econometrica, 68(6), Singleton, 2001, Estimation of Affine Asset Pricing Models Using the Empirical Characteristic Function, Journal of Econometrics, 102, Liuren Wu Fourier Transforms Option Pricing, Fall, / 21
11 Proofs: The option transform χ I c (u) e iuk dc (k) = e iuk Z c (k) c (k) iue iuk dk Checking the boundary conditions, we have c ( ) = 0 (when strike is infinity) and c () = e rt F t /S t = e qt when the strike is zero. Hence, e iu c ( ) = 0 and we will carry the other non-convergent limit e iu e qt. χ I c (u) = e iu e qt c (k) iue iuk dk = e iu e qt iue rt " e st e k # 1 s t k df (s) e iuk dk Z " = e iu e qt iue rt e st e k # 1 s t k eiuk dk df (s) Z " Z s = e iu e qt iue rt t e iuk+st e (iu+1)k # dk df (s) = 2 e iu e qt iue rt 4 e s t eiuk iu e 3 (iu+1)k st 5 df (s). iu + 1 We need to check the boundary again. lim k e (iu+1)k = 0 given the real component e. The other boundary is non-convergent e s t e iu, which we pull out and take the expectation to have iue rt e s t e iu df (s) = e iu e rt F t /S t = e iu e qt, iu which cancels out the other nonconvergent term. Z " # χ I e (iu+1)st c (u) = iue rt e(iu+1)s t df (s) = e rt e (iu+1)s t rt e φ (u i) df (s) =. iu iu + 1 iu + 1 iu + 1 Liuren Wu Fourier Transforms Option Pricing, Fall, / 21
12 Proofs: The option transform inversion The proof for the option transform inversion is similar to that for the CDF, except for the boundary conditions: Z x sgn (y x) dc (y) = dc (y) + dc (y) = c (x) + c () + c ( ) c (x) = e qt 2c (x). x In the case of CDF, F ( ) = 1, and F () = 0. In the case of call option prices, c ( ) = 0 (when strike is infinity), and c () = e rt F t /S t = e qt (when strike is zero). Hence, the inversion formula is I e iux χ (u) e iux χ ( u) du 0 iu = (as before) = πsgn (x z) df (z) = π e qt 2c (x). Thus, c (x) = 1 2 e qt + 1 2π I. Liuren Wu Fourier Transforms Option Pricing, Fall, / 21
13 II. The PDF analog Treat c(k) analogous to a PDF. The option transform: χ II c (z) e izk c(k)dk = e rt φ s (z i) (iz) (iz + 1) with z = u iα, α D R + for the option transform to be well defined. The range of α depends on payoff structure and model. The exact value choice of α is a numerical issue. Carr and Madan (1999, Journal of Computational Finance) refer to α as the dampening coefficient. Given the transform on return φ s(u), we know the transform on call. The inversion is analogous to that for a PDF: c(k) = 1 2 iα+ iα e izk χ II c (z)dz = e αk π 0 e iuk χ II c (u iα)du. References: Carr&Wu, Time-Changed Levy Processes and Option Pricing, JFE, 2004, 17(1), Liuren Wu Fourier Transforms Option Pricing, Fall, / 21
14 Proofs χ II (z) e izk c (k) dk = e rt " = e rt " = e rt " Z s t 2 = e rt 4 e s t eizk iz e st e k 1 s t k df (s) # e izk dk e st e k # 1 s t k eizk dk df (s) e izk+st e (iz+1)k # dk df (s) e(iz+1)k iz st 5 df (s) We need to consider the boundary conditions at k =. lim k e (iz+1)k = 0 as long as the real component of izu is greater than 1. lim k e izk = 0 as long as the real component of iz is greater than 0. Hence, taken together, we need the real component of iz to be greater than zero. If we write z = u iα, with both u and α real, we have iz = iu + α. Hence, we need α > 0 for the above boundary condition to converge. Given that u i > 0, we have Z " χ II (z) = e rt e (iz+1)st iz e(iz+1)s t iz + 1 # df (s) = e rt e (iz+1)s t rt e φ (z i) df (s) =. iu (iz + 1) iz (iz + 1) Liuren Wu Fourier Transforms Option Pricing, Fall, / 21
15 Fast Fourier Transform (FFT) FFT is an efficient algorithm for computing discrete Fourier coefficients. The discrete Fourier transform is a mapping of f = (f 0,, f N 1 ) on the vector of Fourier coefficients d = (d 0,, d N 1 ), such that d j = 1 N 1 f k e jk 2π N i, j = 0, 1,, N 1. N k=0 FFT allows the efficient calculation of d if N is an even number, say N = 2 n, n N. The algorithm reduces the number of multiplcations in the required N summations from an order of 2 2n to that of n2 n 1, a very considerable reduction. By a suitable discretization, we can approximate the inversion of a PDF (also option price) in the above form to take advantage of the computational efficiency of FFT. Liuren Wu Fourier Transforms Option Pricing, Fall, / 21
16 Return PDF inversion Compare the PDF inversion with the FFT form: f X (x) = 1 π 0 e iux φ X (u)du. d j = 1 N 1 f k e jk 2π N i N k=0 Discretize the integral: f X (x) 1 π N k=0 e iu k x φ X (u k ) u Set η = u, u k = ηk. Set x j = b + λj with λ = 2π/(ηN) being the return grid and b being a parameter that controls the return range. To center return around zero, set b = λn/2. The PDF becomes f X (x j ) 1 N 1 f k e jk 2π N i, f m = N N π eiu k b φ X (u k )η. (1) k=0 with j = 0, 1,, N 1. The summation has the FFT form and can hence be computed efficiently. Liuren Wu Fourier Transforms Option Pricing, Fall, / 21
17 Call value inversion Compare the call inversion (method II) with the FFT form: c(k) = e νk π 0 e iuk χ II c (u iν)du. Discretize the integral: c(k) e νk π Set η = u, u m = ηm. d j = 1 N 1 f m e jm 2π N i N m=0 N m=0 e iumk χ II c (u m iν) u Set k j = b + λj with λ = 2π/(ηN) being the return grid and b being a parameter that controls the return range. To center return around zero, set b = λn/2. The call value becomes c(k j ) 1 N 1 f m e jm 2π N i, N m=0 f m = N π e νk j +iu mb χ II c (u m )η. with j = 0, 1,, N 1. The summation has the FFT form and can hence be computed efficiently. Liuren Wu Fourier Transforms Option Pricing, Fall, / 21
18 FFT implementation To implement the FFT, we need to fix the following parameters N = 2 n : The number of summation grids. Setting it to be the power of 2 can speed up the FFT calculation. η = u: The discrete summation grid width. The smaller the grid, the better the approximation. However, given N, η also determines the strike grid λ = 2π/(ηN). The finer the summation grid η, the coarser the strike spacing returned from the FFT calculation. There is a trade off: If we want to have more option value calculated at a finer grid of strikes, we would need to use a coarser summation grid and hence less accuracy. The lower and upper bound truncation b = λn/2 is also determined by the summation grid choice. FFT generates option values at N strikes simultaneously. However, if the strike grid is larger, many of the returned strikes are out of the interesting region. Liuren Wu Fourier Transforms Option Pricing, Fall, / 21
19 Fractional FFT Fractional FFT (FRFT) separates the integration grid choice from the strike grids. With appropriate control, it can generate more accurate option values given the same amount of calculation. The method can efficiently compute, N 1 d j = f m e jmζi, j = 0, 1,..., N 1, m=0 for any value of the parameter ζ. The standard FFT can be seen as a special case for ζ = 2π/N. Therefore, we can use the FRFT method to compute, c(k j ) 1 N 1 f m e jmηλi, N m=0 f m = N π e νk j +iu mb χ II c (u m )η. without the trade-off between the summation grid η and the strike spacing λ. We require ηλ = 2π/N under standard FFT. Liuren Wu Fourier Transforms Option Pricing, Fall, / 21
20 Fractional FFT implementation Let d = D(f, ζ) denote the FRFT operation, with D(f) = D(f, 2π/N) being the standard FFT as a special case. An N-point FRFT can be implemented by invoking three 2N-point FFT procedures. Define the following 2N-point vectors: ( ( ) ) N 1 y = f n e iπn2 ζ, n=0 (0)N 1 n=0, (2) ( ( ) N 1 z = e iπn2 ζ, (e ) ) N 1 iπ(n n)2 α. (3) n=0 The FRFT is given by, D k (h, ζ) = (e ) N 1 iπk2 ζ k=0 n=0 D 1 k (D j (y) D j (z)), (4) where D 1 k ( ) denotes the inverse FFT operation and denotes element-by-element vector multiplication. Liuren Wu Fourier Transforms Option Pricing, Fall, / 21
21 Fractional FFT implementation Due to the multiple application of the FFT operations, an N-point FRFT procedure demands a similar number of elementary operations as a 4N-point FFT procedure. Given the free choices on λ and η, FRFT can be applied more efficiently. Using a smaller N with FRFT can achieve the same option pricing accuracy as using a much larger N with FFT. The accuracy improvement is larger when we have a better understanding of the model and model parameters so that we can set the boundaries more tightly. Caveat: The more freedom also asks for more discretion and caution in applying this method to generate robust results in all situations. This concern becomes especially important for model estimation, during which the trial model parameters can vary greatly. Reference: Chourdakis, 2005, Option pricing using fractional FFT, JCF, 8(2). Liuren Wu Fourier Transforms Option Pricing, Fall, / 21
Time-changed Brownian motion and option pricing
Time-changed Brownian motion and option pricing Peter Hieber Chair of Mathematical Finance, TU Munich 6th AMaMeF Warsaw, June 13th 2013 Partially joint with Marcos Escobar (RU Toronto), Matthias Scherer
More informationModel Estimation. Liuren Wu. Fall, Zicklin School of Business, Baruch College. Liuren Wu Model Estimation Option Pricing, Fall, / 16
Model Estimation Liuren Wu Zicklin School of Business, Baruch College Fall, 2007 Liuren Wu Model Estimation Option Pricing, Fall, 2007 1 / 16 Outline 1 Statistical dynamics 2 Risk-neutral dynamics 3 Joint
More informationA Highly Efficient Shannon Wavelet Inverse Fourier Technique for Pricing European Options
A Highly Efficient Shannon Wavelet Inverse Fourier Technique for Pricing European Options Luis Ortiz-Gracia Centre de Recerca Matemàtica (joint work with Cornelis W. Oosterlee, CWI) Models and Numerics
More informationThe Black-Scholes Model
The Black-Scholes Model Liuren Wu Options Markets (Hull chapter: 12, 13, 14) Liuren Wu ( c ) The Black-Scholes Model colorhmoptions Markets 1 / 17 The Black-Scholes-Merton (BSM) model Black and Scholes
More informationThe Black-Scholes Model
The Black-Scholes Model Liuren Wu Options Markets Liuren Wu ( c ) The Black-Merton-Scholes Model colorhmoptions Markets 1 / 18 The Black-Merton-Scholes-Merton (BMS) model Black and Scholes (1973) and Merton
More informationStochastic Volatility (Working Draft I)
Stochastic Volatility (Working Draft I) Paul J. Atzberger General comments or corrections should be sent to: paulatz@cims.nyu.edu 1 Introduction When using the Black-Scholes-Merton model to price derivative
More informationUsing Lévy Processes to Model Return Innovations
Using Lévy Processes to Model Return Innovations Liuren Wu Zicklin School of Business, Baruch College Option Pricing Liuren Wu (Baruch) Lévy Processes Option Pricing 1 / 32 Outline 1 Lévy processes 2 Lévy
More informationLecture 1: Stochastic Volatility and Local Volatility
Lecture 1: Stochastic Volatility and Local Volatility Jim Gatheral, Merrill Lynch Case Studies in Financial Modelling Course Notes, Courant Institute of Mathematical Sciences, Fall Term, 2003 Abstract
More informationAN ANALYTICALLY TRACTABLE UNCERTAIN VOLATILITY MODEL
AN ANALYTICALLY TRACTABLE UNCERTAIN VOLATILITY MODEL FABIO MERCURIO BANCA IMI, MILAN http://www.fabiomercurio.it 1 Stylized facts Traders use the Black-Scholes formula to price plain-vanilla options. An
More informationExact Sampling of Jump-Diffusion Processes
1 Exact Sampling of Jump-Diffusion Processes and Dmitry Smelov Management Science & Engineering Stanford University Exact Sampling of Jump-Diffusion Processes 2 Jump-Diffusion Processes Ubiquitous in finance
More informationMATH3075/3975 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS
MATH307/37 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS School of Mathematics and Statistics Semester, 04 Tutorial problems should be used to test your mathematical skills and understanding of the lecture material.
More informationPricing European Options by Stable Fourier-Cosine Series Expansions
Pricing European Options by Stable Fourier-Cosine Series Expansions arxiv:171.886v2 [q-fin.cp] 8 Jan 217 Chunfa Wang Fiance School of Zhejiang University of Finance and Economics, Hangzhou, China, cfwang@zufe.edu.cn
More informationSlides for Risk Management
Slides for Risk Management Introduction to the modeling of assets Groll Seminar für Finanzökonometrie Prof. Mittnik, PhD Groll (Seminar für Finanzökonometrie) Slides for Risk Management Prof. Mittnik,
More informationAnalytical Option Pricing under an Asymmetrically Displaced Double Gamma Jump-Diffusion Model
Analytical Option Pricing under an Asymmetrically Displaced Double Gamma Jump-Diffusion Model Advances in Computational Economics and Finance Univerity of Zürich, Switzerland Matthias Thul 1 Ally Quan
More informationMSC FINANCIAL ENGINEERING PRICING I, AUTUMN LECTURE 6: EXTENSIONS OF BLACK AND SCHOLES RAYMOND BRUMMELHUIS DEPARTMENT EMS BIRKBECK
MSC FINANCIAL ENGINEERING PRICING I, AUTUMN 2010-2011 LECTURE 6: EXTENSIONS OF BLACK AND SCHOLES RAYMOND BRUMMELHUIS DEPARTMENT EMS BIRKBECK In this section we look at some easy extensions of the Black
More informationBinomial model: numerical algorithm
Binomial model: numerical algorithm S / 0 C \ 0 S0 u / C \ 1,1 S0 d / S u 0 /, S u 3 0 / 3,3 C \ S0 u d /,1 S u 5 0 4 0 / C 5 5,5 max X S0 u,0 S u C \ 4 4,4 C \ 3 S u d / 0 3, C \ S u d 0 S u d 0 / C 4
More informationStochastic modelling of electricity markets Pricing Forwards and Swaps
Stochastic modelling of electricity markets Pricing Forwards and Swaps Jhonny Gonzalez School of Mathematics The University of Manchester Magical books project August 23, 2012 Clip for this slide Pricing
More informationValuation of performance-dependent options in a Black- Scholes framework
Valuation of performance-dependent options in a Black- Scholes framework Thomas Gerstner, Markus Holtz Institut für Numerische Simulation, Universität Bonn, Germany Ralf Korn Fachbereich Mathematik, TU
More informationOption Pricing and Calibration with Time-changed Lévy processes
Option Pricing and Calibration with Time-changed Lévy processes Yan Wang and Kevin Zhang Warwick Business School 12th Feb. 2013 Objectives 1. How to find a perfect model that captures essential features
More informationPricing Variance Swaps under Stochastic Volatility Model with Regime Switching - Discrete Observations Case
Pricing Variance Swaps under Stochastic Volatility Model with Regime Switching - Discrete Observations Case Guang-Hua Lian Collaboration with Robert Elliott University of Adelaide Feb. 2, 2011 Robert Elliott,
More informationA GENERAL FORMULA FOR OPTION PRICES IN A STOCHASTIC VOLATILITY MODEL. Stephen Chin and Daniel Dufresne. Centre for Actuarial Studies
A GENERAL FORMULA FOR OPTION PRICES IN A STOCHASTIC VOLATILITY MODEL Stephen Chin and Daniel Dufresne Centre for Actuarial Studies University of Melbourne Paper: http://mercury.ecom.unimelb.edu.au/site/actwww/wps2009/no181.pdf
More information2.1 Mathematical Basis: Risk-Neutral Pricing
Chapter Monte-Carlo Simulation.1 Mathematical Basis: Risk-Neutral Pricing Suppose that F T is the payoff at T for a European-type derivative f. Then the price at times t before T is given by f t = e r(t
More informationDynamic Relative Valuation
Dynamic Relative Valuation Liuren Wu, Baruch College Joint work with Peter Carr from Morgan Stanley October 15, 2013 Liuren Wu (Baruch) Dynamic Relative Valuation 10/15/2013 1 / 20 The standard approach
More informationA Consistent Pricing Model for Index Options and Volatility Derivatives
A Consistent Pricing Model for Index Options and Volatility Derivatives 6th World Congress of the Bachelier Society Thomas Kokholm Finance Research Group Department of Business Studies Aarhus School of
More informationFourier Space Time-stepping Method for Option Pricing with Lévy Processes
FST method Extensions Indifference pricing Fourier Space Time-stepping Method for Option Pricing with Lévy Processes Vladimir Surkov University of Toronto Computational Methods in Finance Conference University
More informationThe data-driven COS method
The data-driven COS method Á. Leitao, C. W. Oosterlee, L. Ortiz-Gracia and S. M. Bohte Delft University of Technology - Centrum Wiskunde & Informatica Reading group, March 13, 2017 Reading group, March
More information3.4 Copula approach for modeling default dependency. Two aspects of modeling the default times of several obligors
3.4 Copula approach for modeling default dependency Two aspects of modeling the default times of several obligors 1. Default dynamics of a single obligor. 2. Model the dependence structure of defaults
More informationP&L Attribution and Risk Management
P&L Attribution and Risk Management Liuren Wu Options Markets (Hull chapter: 15, Greek letters) Liuren Wu ( c ) P& Attribution and Risk Management Options Markets 1 / 19 Outline 1 P&L attribution via the
More informationFinancial Derivatives Section 5
Financial Derivatives Section 5 The Black and Scholes Model Michail Anthropelos anthropel@unipi.gr http://web.xrh.unipi.gr/faculty/anthropelos/ University of Piraeus Spring 2018 M. Anthropelos (Un. of
More informationM5MF6. Advanced Methods in Derivatives Pricing
Course: Setter: M5MF6 Dr Antoine Jacquier MSc EXAMINATIONS IN MATHEMATICS AND FINANCE DEPARTMENT OF MATHEMATICS April 2016 M5MF6 Advanced Methods in Derivatives Pricing Setter s signature...........................................
More informationThe data-driven COS method
The data-driven COS method Á. Leitao, C. W. Oosterlee, L. Ortiz-Gracia and S. M. Bohte Delft University of Technology - Centrum Wiskunde & Informatica CMMSE 2017, July 6, 2017 Álvaro Leitao (CWI & TUDelft)
More informationLOGNORMAL MIXTURE SMILE CONSISTENT OPTION PRICING
LOGNORMAL MIXTURE SMILE CONSISTENT OPTION PRICING FABIO MERCURIO BANCA IMI, MILAN http://www.fabiomercurio.it Daiwa International Workshop on Financial Engineering, Tokyo, 26-27 August 2004 1 Stylized
More informationLinearity-Generating Processes, Unspanned Stochastic Volatility, and Interest-Rate Option Pricing
Linearity-Generating Processes, Unspanned Stochastic Volatility, and Interest-Rate Option Pricing Liuren Wu, Baruch College Joint work with Peter Carr and Xavier Gabaix at New York University Board of
More informationImplied Lévy Volatility
Joint work with José Manuel Corcuera, Peter Leoni and Wim Schoutens July 15, 2009 - Eurandom 1 2 The Black-Scholes model The Lévy models 3 4 5 6 7 Delta Hedging at versus at Implied Black-Scholes Volatility
More informationQueens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 2017 Instructor: Dr. Sateesh Mane.
Queens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 217 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 217 13 Lecture 13 November 15, 217 Derivation of the Black-Scholes-Merton
More informationValuation of European Call Option via Inverse Fourier Transform
ISSN 2255-9094 (online) ISSN 2255-9086 (print) December 2017, vol. 20, pp. 91 96 doi: 10.1515/itms-2017-0016 https://www.degruyter.com/view/j/itms Valuation of European Call Option via Inverse Fourier
More informationFast and accurate pricing of discretely monitored barrier options by numerical path integration
Comput Econ (27 3:143 151 DOI 1.17/s1614-7-991-5 Fast and accurate pricing of discretely monitored barrier options by numerical path integration Christian Skaug Arvid Naess Received: 23 December 25 / Accepted:
More informationChapter 15: Jump Processes and Incomplete Markets. 1 Jumps as One Explanation of Incomplete Markets
Chapter 5: Jump Processes and Incomplete Markets Jumps as One Explanation of Incomplete Markets It is easy to argue that Brownian motion paths cannot model actual stock price movements properly in reality,
More informationSaddlepoint Approximation Methods for Pricing. Financial Options on Discrete Realized Variance
Saddlepoint Approximation Methods for Pricing Financial Options on Discrete Realized Variance Yue Kuen KWOK Department of Mathematics Hong Kong University of Science and Technology Hong Kong * This is
More informationTEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING
TEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING Semih Yön 1, Cafer Erhan Bozdağ 2 1,2 Department of Industrial Engineering, Istanbul Technical University, Macka Besiktas, 34367 Turkey Abstract.
More informationMATH 476/567 ACTUARIAL RISK THEORY FALL 2016 PROFESSOR WANG
MATH 476/567 ACTUARIAL RISK THEORY FALL 206 PROFESSOR WANG Homework 5 (max. points = 00) Due at the beginning of class on Tuesday, November 8, 206 You are encouraged to work on these problems in groups
More informationA Brief Introduction to Stochastic Volatility Modeling
A Brief Introduction to Stochastic Volatility Modeling Paul J. Atzberger General comments or corrections should be sent to: paulatz@cims.nyu.edu Introduction When using the Black-Scholes-Merton model to
More informationSolving the Black-Scholes Equation
Solving the Black-Scholes Equation An Undergraduate Introduction to Financial Mathematics J. Robert Buchanan 2014 Initial Value Problem for the European Call The main objective of this lesson is solving
More informationValuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility model
Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility model 1(23) Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility
More informationUQ, STAT2201, 2017, Lectures 3 and 4 Unit 3 Probability Distributions.
UQ, STAT2201, 2017, Lectures 3 and 4 Unit 3 Probability Distributions. Random Variables 2 A random variable X is a numerical (integer, real, complex, vector etc.) summary of the outcome of the random experiment.
More informationOption Pricing for a Stochastic-Volatility Jump-Diffusion Model
Option Pricing for a Stochastic-Volatility Jump-Diffusion Model Guoqing Yan and Floyd B. Hanson Department of Mathematics, Statistics, and Computer Science University of Illinois at Chicago Conference
More informationHomework Assignments
Homework Assignments Week 1 (p 57) #4.1, 4., 4.3 Week (pp 58-6) #4.5, 4.6, 4.8(a), 4.13, 4.0, 4.6(b), 4.8, 4.31, 4.34 Week 3 (pp 15-19) #1.9, 1.1, 1.13, 1.15, 1.18 (pp 9-31) #.,.6,.9 Week 4 (pp 36-37)
More informationSimple Robust Hedging with Nearby Contracts
Simple Robust Hedging with Nearby Contracts Liuren Wu and Jingyi Zhu Baruch College and University of Utah October 22, 2 at Worcester Polytechnic Institute Wu & Zhu (Baruch & Utah) Robust Hedging with
More informationVaR Estimation under Stochastic Volatility Models
VaR Estimation under Stochastic Volatility Models Chuan-Hsiang Han Dept. of Quantitative Finance Natl. Tsing-Hua University TMS Meeting, Chia-Yi (Joint work with Wei-Han Liu) December 5, 2009 Outline Risk
More informationPath Dependent British Options
Path Dependent British Options Kristoffer J Glover (Joint work with G. Peskir and F. Samee) School of Finance and Economics University of Technology, Sydney 18th August 2009 (PDE & Mathematical Finance
More informationPricing Barrier Options under Local Volatility
Abstract Pricing Barrier Options under Local Volatility Artur Sepp Mail: artursepp@hotmail.com, Web: www.hot.ee/seppar 16 November 2002 We study pricing under the local volatility. Our research is mainly
More information1.12 Exercises EXERCISES Use integration by parts to compute. ln(x) dx. 2. Compute 1 x ln(x) dx. Hint: Use the substitution u = ln(x).
2 EXERCISES 27 2 Exercises Use integration by parts to compute lnx) dx 2 Compute x lnx) dx Hint: Use the substitution u = lnx) 3 Show that tan x) =/cos x) 2 and conclude that dx = arctanx) + C +x2 Note:
More informationPricing Bermudan options in Lévy process models
1 1 Dept. of Industrial & Enterprise Systems Engineering University of Illinois at Urbana-Champaign Joint with Xiong Lin Bachelier Finance Society 6th Congress 6/24/2010 American options in Lévy models
More informationRisk management. Introduction to the modeling of assets. Christian Groll
Risk management Introduction to the modeling of assets Christian Groll Introduction to the modeling of assets Risk management Christian Groll 1 / 109 Interest rates and returns Interest rates and returns
More informationLarge Deviations and Stochastic Volatility with Jumps: Asymptotic Implied Volatility for Affine Models
Large Deviations and Stochastic Volatility with Jumps: TU Berlin with A. Jaquier and A. Mijatović (Imperial College London) SIAM conference on Financial Mathematics, Minneapolis, MN July 10, 2012 Implied
More informationLocal vs Non-local Forward Equations for Option Pricing
Local vs Non-local Forward Equations for Option Pricing Rama Cont Yu Gu Abstract When the underlying asset is a continuous martingale, call option prices solve the Dupire equation, a forward parabolic
More informationEFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS
Commun. Korean Math. Soc. 23 (2008), No. 2, pp. 285 294 EFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS Kyoung-Sook Moon Reprinted from the Communications of the Korean Mathematical Society
More informationProbability in Options Pricing
Probability in Options Pricing Mark Cohen and Luke Skon Kenyon College cohenmj@kenyon.edu December 14, 2012 Mark Cohen and Luke Skon (Kenyon college) Probability Presentation December 14, 2012 1 / 16 What
More informationBarrier options. In options only come into being if S t reaches B for some 0 t T, at which point they become an ordinary option.
Barrier options A typical barrier option contract changes if the asset hits a specified level, the barrier. Barrier options are therefore path-dependent. Out options expire worthless if S t reaches the
More informationPDE Methods for the Maximum Drawdown
PDE Methods for the Maximum Drawdown Libor Pospisil, Jan Vecer Columbia University, Department of Statistics, New York, NY 127, USA April 1, 28 Abstract Maximum drawdown is a risk measure that plays an
More informationApplying stochastic time changes to Lévy processes
Applying stochastic time changes to Lévy processes Liuren Wu Zicklin School of Business, Baruch College Option Pricing Liuren Wu (Baruch) Stochastic time changes Option Pricing 1 / 38 Outline 1 Stochastic
More informationLeverage Effect, Volatility Feedback, and Self-Exciting Market Disruptions
On-line Appendix for Leverage Effect, Volatility Feedback, and Self-Exciting Market Disruptions Peter Carr 1 and Liuren Wu 2 I. Model Summary We fix a filtered probability space {Ω,F,P,(F t ) t 0 } and
More informationA No-Arbitrage Theorem for Uncertain Stock Model
Fuzzy Optim Decis Making manuscript No (will be inserted by the editor) A No-Arbitrage Theorem for Uncertain Stock Model Kai Yao Received: date / Accepted: date Abstract Stock model is used to describe
More information1.1 Basic Financial Derivatives: Forward Contracts and Options
Chapter 1 Preliminaries 1.1 Basic Financial Derivatives: Forward Contracts and Options A derivative is a financial instrument whose value depends on the values of other, more basic underlying variables
More informationOption Pricing for a Stochastic-Volatility Jump-Diffusion Model with Log-Uniform Jump-Amplitudes
Option Pricing for a Stochastic-Volatility Jump-Diffusion Model with Log-Uniform Jump-Amplitudes Floyd B. Hanson and Guoqing Yan Department of Mathematics, Statistics, and Computer Science University of
More informationUSC Math. Finance April 22, Path-dependent Option Valuation under Jump-diffusion Processes. Alan L. Lewis
USC Math. Finance April 22, 23 Path-dependent Option Valuation under Jump-diffusion Processes Alan L. Lewis These overheads to be posted at www.optioncity.net (Publications) Topics Why jump-diffusion models?
More informationSolving the Black-Scholes Equation
Solving the Black-Scholes Equation An Undergraduate Introduction to Financial Mathematics J. Robert Buchanan 2010 Initial Value Problem for the European Call rf = F t + rsf S + 1 2 σ2 S 2 F SS for (S,
More informationEfficient Static Replication of European Options under Exponential Lévy Models
CIRJE-F-539 Efficient Static Replication of European Options under Exponential Lévy Models Akihiko Takahashi University of Tokyo Akira Yamazaki Mizuho-DL Financial Technology Co., Ltd. January 8 CIRJE
More informationImplementing Models in Quantitative Finance: Methods and Cases
Gianluca Fusai Andrea Roncoroni Implementing Models in Quantitative Finance: Methods and Cases vl Springer Contents Introduction xv Parti Methods 1 Static Monte Carlo 3 1.1 Motivation and Issues 3 1.1.1
More informationTangent Lévy Models. Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford.
Tangent Lévy Models Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford June 24, 2010 6th World Congress of the Bachelier Finance Society Sergey
More informationImplied Volatility Surface
Implied Volatility Surface Liuren Wu Zicklin School of Business, Baruch College Fall, 2007 Liuren Wu Implied Volatility Surface Option Pricing, Fall, 2007 1 / 22 Implied volatility Recall the BSM formula:
More informationEquity correlations implied by index options: estimation and model uncertainty analysis
1/18 : estimation and model analysis, EDHEC Business School (joint work with Rama COT) Modeling and managing financial risks Paris, 10 13 January 2011 2/18 Outline 1 2 of multi-asset models Solution to
More informationIEOR 165 Lecture 1 Probability Review
IEOR 165 Lecture 1 Probability Review 1 Definitions in Probability and Their Consequences 1.1 Defining Probability A probability space (Ω, F, P) consists of three elements: A sample space Ω is the set
More informationFixed-Income Options
Fixed-Income Options Consider a two-year 99 European call on the three-year, 5% Treasury. Assume the Treasury pays annual interest. From p. 852 the three-year Treasury s price minus the $5 interest could
More informationFast Pricing and Calculation of Sensitivities of OTM European Options Under Lévy Processes
Fast Pricing and Calculation of Sensitivities of OTM European Options Under Lévy Processes Sergei Levendorskĭi Jiayao Xie Department of Mathematics University of Leicester Toronto, June 24, 2010 Levendorskĭi
More informationA distributed Laplace transform algorithm for European options
A distributed Laplace transform algorithm for European options 1 1 A. J. Davies, M. E. Honnor, C.-H. Lai, A. K. Parrott & S. Rout 1 Department of Physics, Astronomy and Mathematics, University of Hertfordshire,
More informationHedging under Arbitrage
Hedging under Arbitrage Johannes Ruf Columbia University, Department of Statistics Modeling and Managing Financial Risks January 12, 2011 Motivation Given: a frictionless market of stocks with continuous
More informationValuation of Equity / FX Instruments
Technical Paper: Valuation of Equity / FX Instruments MathConsult GmbH Altenberger Straße 69 A-4040 Linz, Austria 14 th October, 2009 1 Vanilla Equity Option 1.1 Introduction A vanilla equity option is
More informationNumerical Evaluation of American Options Written on Two Underlying Assets using the Fourier Transform Approach
1 / 26 Numerical Evaluation of American Options Written on Two Underlying Assets using the Fourier Transform Approach Jonathan Ziveyi Joint work with Prof. Carl Chiarella School of Finance and Economics,
More informationFINANCIAL OPTION ANALYSIS HANDOUTS
FINANCIAL OPTION ANALYSIS HANDOUTS 1 2 FAIR PRICING There is a market for an object called S. The prevailing price today is S 0 = 100. At this price the object S can be bought or sold by anyone for any
More informationSADDLEPOINT APPROXIMATIONS TO OPTION PRICES 1. By L. C. G. Rogers and O. Zane University of Bath and First Chicago NBD
The Annals of Applied Probability 1999, Vol. 9, No. 2, 493 53 SADDLEPOINT APPROXIMATIONS TO OPTION PRICES 1 By L. C. G. Rogers and O. Zane University of Bath and First Chicago NBD The use of saddlepoint
More informationApplication of an Interval Backward Finite Difference Method for Solving the One-Dimensional Heat Conduction Problem
Application of an Interval Backward Finite Difference Method for Solving the One-Dimensional Heat Conduction Problem Malgorzata A. Jankowska 1, Andrzej Marciniak 2 and Tomasz Hoffmann 2 1 Poznan University
More informationMASM006 UNIVERSITY OF EXETER SCHOOL OF ENGINEERING, COMPUTER SCIENCE AND MATHEMATICS MATHEMATICAL SCIENCES FINANCIAL MATHEMATICS.
MASM006 UNIVERSITY OF EXETER SCHOOL OF ENGINEERING, COMPUTER SCIENCE AND MATHEMATICS MATHEMATICAL SCIENCES FINANCIAL MATHEMATICS May/June 2006 Time allowed: 2 HOURS. Examiner: Dr N.P. Byott This is a CLOSED
More informationUnified Credit-Equity Modeling
Unified Credit-Equity Modeling Rafael Mendoza-Arriaga Based on joint research with: Vadim Linetsky and Peter Carr The University of Texas at Austin McCombs School of Business (IROM) Recent Advancements
More informationMachine Learning for Quantitative Finance
Machine Learning for Quantitative Finance Fast derivative pricing Sofie Reyners Joint work with Jan De Spiegeleer, Dilip Madan and Wim Schoutens Derivative pricing is time-consuming... Vanilla option pricing
More informationAdvanced Topics in Derivative Pricing Models. Topic 4 - Variance products and volatility derivatives
Advanced Topics in Derivative Pricing Models Topic 4 - Variance products and volatility derivatives 4.1 Volatility trading and replication of variance swaps 4.2 Volatility swaps 4.3 Pricing of discrete
More informationOptimal Hedging of Variance Derivatives. John Crosby. Centre for Economic and Financial Studies, Department of Economics, Glasgow University
Optimal Hedging of Variance Derivatives John Crosby Centre for Economic and Financial Studies, Department of Economics, Glasgow University Presentation at Baruch College, in New York, 16th November 2010
More informationEfficiently pricing barrier options in a Markov-switching framework
Efficiently pricing barrier options in a Markov-switching framework Peter Hieber HVB-Institute for Mathematical Finance, Technische Universität München, Parkring 11, 85748 Garching-Hochbrück, Germany,
More informationLévy models in finance
Lévy models in finance Ernesto Mordecki Universidad de la República, Montevideo, Uruguay PASI - Guanajuato - June 2010 Summary General aim: describe jummp modelling in finace through some relevant issues.
More informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Abdul-Lateef Haji-Ali Based on slides by: Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Haji-Ali (Oxford)
More information6. Numerical methods for option pricing
6. Numerical methods for option pricing Binomial model revisited Under the risk neutral measure, ln S t+ t ( ) S t becomes normally distributed with mean r σ2 t and variance σ 2 t, where r is 2 the riskless
More informationAdvanced Numerical Methods
Advanced Numerical Methods Solution to Homework One Course instructor: Prof. Y.K. Kwok. When the asset pays continuous dividend yield at the rate q the expected rate of return of the asset is r q under
More informationFrom Discrete Time to Continuous Time Modeling
From Discrete Time to Continuous Time Modeling Prof. S. Jaimungal, Department of Statistics, University of Toronto 2004 Arrow-Debreu Securities 2004 Prof. S. Jaimungal 2 Consider a simple one-period economy
More informationModeling of Price. Ximing Wu Texas A&M University
Modeling of Price Ximing Wu Texas A&M University As revenue is given by price times yield, farmers income risk comes from risk in yield and output price. Their net profit also depends on input price, but
More informationSYSM 6304: Risk and Decision Analysis Lecture 6: Pricing and Hedging Financial Derivatives
SYSM 6304: Risk and Decision Analysis Lecture 6: Pricing and Hedging Financial Derivatives M. Vidyasagar Cecil & Ida Green Chair The University of Texas at Dallas Email: M.Vidyasagar@utdallas.edu October
More informationHow do Variance Swaps Shape the Smile?
How do Variance Swaps Shape the Smile? A Summary of Arbitrage Restrictions and Smile Asymptotics Vimal Raval Imperial College London & UBS Investment Bank www2.imperial.ac.uk/ vr402 Joint Work with Mark
More informationA Simple Robust Link Between American Puts and Credit Insurance
A Simple Robust Link Between American Puts and Credit Insurance Liuren Wu at Baruch College Joint work with Peter Carr Ziff Brothers Investments, April 2nd, 2010 Liuren Wu (Baruch) DOOM Puts & Credit Insurance
More informationNUMERICAL METHODS OF PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS FOR OPTION PRICE
Trends in Mathematics - New Series Information Center for Mathematical Sciences Volume 13, Number 1, 011, pages 1 5 NUMERICAL METHODS OF PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS FOR OPTION PRICE YONGHOON
More informationThe Forward PDE for American Puts in the Dupire Model
The Forward PDE for American Puts in the Dupire Model Peter Carr Ali Hirsa Courant Institute Morgan Stanley New York University 750 Seventh Avenue 51 Mercer Street New York, NY 10036 1 60-3765 (1) 76-988
More informationSTOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL
STOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL YOUNGGEUN YOO Abstract. Ito s lemma is often used in Ito calculus to find the differentials of a stochastic process that depends on time. This paper will introduce
More information