Lecture 7: Computation of Greeks

Similar documents
Estimating the Greeks

2 f. f t S 2. Delta measures the sensitivityof the portfolio value to changes in the price of the underlying

Importance Sampling for Option Pricing. Steven R. Dunbar. Put Options. Monte Carlo Method. Importance. Sampling. Examples.

1.1 Basic Financial Derivatives: Forward Contracts and Options

IEOR E4703: Monte-Carlo Simulation

Queens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Spring 2018 Instructor: Dr. Sateesh Mane.

Monte Carlo Methods in Option Pricing. UiO-STK4510 Autumn 2015

M5MF6. Advanced Methods in Derivatives Pricing

In this lecture we will solve the final-value problem derived in the previous lecture 4, V (1) + rs = rv (t < T )

Mathematics of Finance Final Preparation December 19. To be thoroughly prepared for the final exam, you should

Numerical schemes for SDEs

Definition Pricing Risk management Second generation barrier options. Barrier Options. Arfima Financial Solutions

MAS3904/MAS8904 Stochastic Financial Modelling

Conditional Density Method in the Computation of the Delta with Application to Power Market

Lecture Note 8 of Bus 41202, Spring 2017: Stochastic Diffusion Equation & Option Pricing

Chapter 3 Common Families of Distributions. Definition 3.4.1: A family of pmfs or pdfs is called exponential family if it can be expressed as

"Vibrato" Monte Carlo evaluation of Greeks

1 Geometric Brownian motion

Queens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 2017 Instructor: Dr. Sateesh Mane.

The Use of Importance Sampling to Speed Up Stochastic Volatility Simulations

Barrier 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.

Financial Risk Management

Lévy models in finance

Monte Carlo Based Numerical Pricing of Multiple Strike-Reset Options

Change of Measure (Cameron-Martin-Girsanov Theorem)

Rapid computation of prices and deltas of nth to default swaps in the Li Model

NEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 MAS3904. Stochastic Financial Modelling. Time allowed: 2 hours

STOCHASTIC INTEGRALS

SYSM 6304: Risk and Decision Analysis Lecture 6: Pricing and Hedging Financial Derivatives

MATH3075/3975 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS

Brownian Motion. Richard Lockhart. Simon Fraser University. STAT 870 Summer 2011

Stochastic Differential Equations in Finance and Monte Carlo Simulations

Capital Allocation Principles

Monte Carlo Simulations

Tangent Lévy Models. Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford.

Time-changed Brownian motion and option pricing

All Investors are Risk-averse Expected Utility Maximizers

King s College London

All Investors are Risk-averse Expected Utility Maximizers. Carole Bernard (UW), Jit Seng Chen (GGY) and Steven Vanduffel (Vrije Universiteit Brussel)

FIN FINANCIAL INSTRUMENTS SPRING 2008

FINANCIAL OPTION ANALYSIS HANDOUTS

Financial Risk Management

Module 2: Monte Carlo Methods

Computing Greeks with Multilevel Monte Carlo Methods using Importance Sampling

Logarithmic derivatives of densities for jump processes

Analytical formulas for local volatility model with stochastic. Mohammed Miri

Module 10:Application of stochastic processes in areas like finance Lecture 36:Black-Scholes Model. Stochastic Differential Equation.

STOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL

Lecture 17. The model is parametrized by the time period, δt, and three fixed constant parameters, v, σ and the riskless rate r.

King s College London

Valuation of performance-dependent options in a Black- Scholes framework

Math 181 Lecture 15 Hedging and the Greeks (Chap. 14, Hull)

AN ANALYTICALLY TRACTABLE UNCERTAIN VOLATILITY MODEL

A Continuity Correction under Jump-Diffusion Models with Applications in Finance

MASM006 UNIVERSITY OF EXETER SCHOOL OF ENGINEERING, COMPUTER SCIENCE AND MATHEMATICS MATHEMATICAL SCIENCES FINANCIAL MATHEMATICS.

SYSM 6304 Risk and Decision Analysis Lecture 2: Fitting Distributions to Data

Generating Random Variables and Stochastic Processes

Local Volatility Dynamic Models

Chapter 3: Black-Scholes Equation and Its Numerical Evaluation

KØBENHAVNS UNIVERSITET (Blok 2, 2011/2012) Naturvidenskabelig kandidateksamen Continuous time finance (FinKont) TIME ALLOWED : 3 hours

Analysis of the sensitivity to discrete dividends : A new approach for pricing vanillas

Short-time-to-expiry expansion for a digital European put option under the CEV model. November 1, 2017

A GENERAL FORMULA FOR OPTION PRICES IN A STOCHASTIC VOLATILITY MODEL. Stephen Chin and Daniel Dufresne. Centre for Actuarial Studies

PDE Methods for the Maximum Drawdown

Asset-or-nothing digitals

AMH4 - ADVANCED OPTION PRICING. Contents

Queens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 2017 Instructor: Dr. Sateesh Mane.

Advanced Topics in Derivative Pricing Models. Topic 4 - Variance products and volatility derivatives

The Lognormal Interest Rate Model and Eurodollar Futures

Computational Finance

Stochastic Modelling in Finance

Exact replication under portfolio constraints: a viability approach

MLEMVD: A R Package for Maximum Likelihood Estimation of Multivariate Diffusion Models

Calculating Implied Volatility

Replication and Absence of Arbitrage in Non-Semimartingale Models

Option Pricing Models for European Options

Lecture Quantitative Finance Spring Term 2015

STOCHASTIC VOLATILITY AND OPTION PRICING

Bluff Your Way Through Black-Scholes

Pricing Barrier Options under Local Volatility

American Option Pricing Formula for Uncertain Financial Market

Chapter 15: Jump Processes and Incomplete Markets. 1 Jumps as One Explanation of Incomplete Markets

2.1 Mathematical Basis: Risk-Neutral Pricing

Probability in Options Pricing

AD in Monte Carlo for finance

Calibration Lecture 1: Background and Parametric Models

Arbitrages and pricing of stock options

Lecture 18. More on option pricing. Lecture 18 1 / 21

Math 416/516: Stochastic Simulation

The stochastic calculus

Stochastic Processes and Stochastic Calculus - 9 Complete and Incomplete Market Models

Math Computational Finance Option pricing using Brownian bridge and Stratified samlping

Exercise. Show the corrected sample variance is an unbiased estimator of population variance. S 2 = n i=1 (X i X ) 2 n 1. Exercise Estimation

Lecture 6: Option Pricing Using a One-step Binomial Tree. Thursday, September 12, 13

IEOR E4602: Quantitative Risk Management

Path Dependent British Options

Derivatives Pricing and Stochastic Calculus

Posterior Inference. , where should we start? Consider the following computational procedure: 1. draw samples. 2. convert. 3. compute properties

S t d with probability (1 p), where

The value of foresight

Transcription:

Lecture 7: Computation of Greeks Ahmed Kebaier kebaier@math.univ-paris13.fr HEC, Paris

Outline 1 The log-likelihood approach

Motivation The pathwise method requires some restrictive regularity assumptions on the payoff function of the option price, at least one time differentiable. This is not the case of some particular payoffs such as the digital options with a payoff function h T = h(x x T ) where h(x) = 1 {x K}. The Greek of such option cannot be evaluated using the pathwise method. Note also that for the same arguments we cannot use the pathwise method to evaluate the Greek payoff Γ of a given vanilla Call option.

General idea Let us assume that the family of random variables (X (θ)) indexed by a parameter θ Θ (Θ an open set of R), admits positive density function p(θ, y).

General idea Let us assume that the family of random variables (X (θ)) indexed by a parameter θ Θ (Θ an open set of R), admits positive density function p(θ, y). f (θ) = Eφ(X (θ)) = φ(y)p(θ, y)dy. R d

General idea Let us assume that the family of random variables (X (θ)) indexed by a parameter θ Θ (Θ an open set of R), admits positive density function p(θ, y). f (θ) = Eφ(X (θ)) = φ(y)p(θ, y)dy. R d Theorem 1 Assume that the density function p(θ, y) taking values in r Θ R d satisfies i) θ p(θ, y) is differentiable on Θ almost every where

General idea Let us assume that the family of random variables (X (θ)) indexed by a parameter θ Θ (Θ an open set of R), admits positive density function p(θ, y). f (θ) = Eφ(X (θ)) = φ(y)p(θ, y)dy. R d Theorem 1 Assume that the density function p(θ, y) taking values in r Θ R d satisfies i) θ p(θ, y) is differentiable on Θ almost every where ii) g : R d R a measurable function such that R d φ(y)g(y)dy < and θ Θ θ p(θ, y) g(y),

General idea Let us assume that the family of random variables (X (θ)) indexed by a parameter θ Θ (Θ an open set of R), admits positive density function p(θ, y). f (θ) = Eφ(X (θ)) = φ(y)p(θ, y)dy. R d Theorem 1 Assume that the density function p(θ, y) taking values in r Θ R d satisfies i) θ p(θ, y) is differentiable on Θ almost every where ii) g : R d R a measurable function such that R φ(y)g(y)dy < and θ Θ d θ p(θ, y) g(y),then ( θ Θ, f (θ) = E φ(x (θ)) log p ) (θ, X (θ)) θ

Proof. According to assumptions i) and ii) we have f (θ) = θ R Eφ(X (θ)) = φ(y) p(θ, y)dy. d θ Rewriting the above expression we get ( f θ p(θ, y) (θ) = φ(y) p(θ, y)dy = E φ(x (θ)) R d p(θ, y) which completes the proof. ) p(θ, X (θ)), p(θ, X (θ)) θ

Black Scholes Delta Recall that the Black-Scholes model with parameters r, σ, T is given by the explicit solution ( S T = S 0 exp (r σ 2 /2)T + σ ) T G, where G N (0, 1). Consequently, for any given measurable function we get ( ( Ef (S T ) = Ef S 0 exp (r σ 2 /2)T + σ )) T G ( ( = f S 0 exp (r σ 2 /2)T + σ )) T y R 1 exp ( y 2 ) dy. 2π 2

Black Scholes Delta Using the following change of variable ( u = S 0 exp (r σ 2 /2)T + σ ) T y with du = σ ( T S 0 exp (r σ 2 /2)T + σ ) T y dy = uσ T dy and y = ζ(u) where ζ(u) = ( log(u/s 0 ) (r σ 2 /2)T ) /σ T

Black Scholes Delta Using the following change of variable ( u = S 0 exp (r σ 2 /2)T + σ ) T y with du = σ ( T S 0 exp (r σ 2 /2)T + σ ) T y dy = uσ T dy and y = ζ(u) where ζ(u) = ( log(u/s 0 ) (r σ 2 /2)T ) /σ T we get Ef (S T ) = R f (y) 1 uσ T 2π exp ( ζ(u)2 2 ) du.

Black Scholes Delta Therefore we deduce that the density of the random variable S T is given by ) 1 g(x) = ( xσ T 2π exp ζ(x)2. 2

Black Scholes Delta Therefore we deduce that the density of the random variable S T is given by ) 1 g(x) = ( xσ T 2π exp ζ(x)2. 2 Then log g(x) S 0 = g S 0 (x) = log(x/s 0) (r σ 2 /2)T g(x) S 0 σ 2. T

Black Scholes Delta Hence, according to the above theorem = e rt E(S T K) + S 0 ( = e rt log(s T /S 0 ) (r σ 2 ) /2)T E (S T K) + S 0 σ 2 T ( ) = e rt G E (S T K) + S 0 σ. T

Example: Path-dependent deltas Let us consider the case of an Asian option. Since the associated payoff involves the vector (S t1, S t2,, S tn )

Example: Path-dependent deltas Let us consider the case of an Asian option. Since the associated payoff involves the vector (S t1, S t2,, S tn ) we need to make explicit the joint density of this vector if we aim to apply the log-likelihood ratio method for the computation of sensitivities associated to the Asian option price.

Example: Path-dependent deltas Let us consider the case of an Asian option. Since the associated payoff involves the vector (S t1, S t2,, S tn ) we need to make explicit the joint density of this vector if we aim to apply the log-likelihood ratio method for the computation of sensitivities associated to the Asian option price. Using the Markovian property of the Brownian motion we can rewrite the density associated to the above vector as follows g(x 1,, x n ) = g 1 (x 1 S 0 )g 2 (x 2 x 1 ) g n (x n x n 1 ), with g j (x j x j 1 ) = 1 x i σ t j t j 1 ϕ(ζ j (x j x j 1 )), ϕ(u) = 1 2π exp ) ( u2 2 where

Example: Path-dependent deltas ζ j (x j x j 1 ) = log(x j/x j 1 ) (r σ 2 /2)(t j t j 1 ) σ t j t j 1.

Example: Path-dependent deltas ζ j (x j x j 1 ) = log(x j/x j 1 ) (r σ 2 /2)(t j t j 1 ) σ t j t j 1. Consequently, log g S 0 (S t1,, S tn ) = log g 1 S 0 (S t1 S 0 ) = G 1 S 0 σ t 1, where G 1 is the same Gaussian used for the simulation of S t1.

Example: Path-dependent deltas ζ j (x j x j 1 ) = log(x j/x j 1 ) (r σ 2 /2)(t j t j 1 ) σ t j t j 1. Consequently, log g S 0 (S t1,, S tn ) = log g 1 S 0 (S t1 S 0 ) = G 1 S 0 σ t 1, where G 1 is the same Gaussian used for the simulation of S t1. Finally, the value of the delta Asian call option using the Log-likelihood ratio method is given by ( ) = e rt E ( 1 n G 1 S ti K) + n S 0 σ. t 1 i=1