Fast narrow bounds on the value of Asian options
|
|
- Sheila Phillips
- 5 years ago
- Views:
Transcription
1 Fast narrow bounds on the value of Asian options G. W. P. Thompson Centre for Financial Research, Judge Institute of Management, University of Cambridge Abstract We consider the problem of finding bounds on the value of fixed-strike and floatingstrike Asian options. A good lower bound for both types was derived in Rogers & Shi 995). We provide an alternative derivation, which leads to a simpler expression for the bound, and also to the bound given by Curran 992) for fixed-strike options; we derive an analogous bound for floating-strike options. Combining these results with a new upper bound allows the accurate valuation of fixed-strike and floating-strike Asian options for typical parameter values. Keywords: Brownian motion, Asian option, fixed-strike, floating-strike. Introduction A fixed-strike Asian call) option with strike K > and maturity T on an asset with price process {S t } is a contract with value T T S t dt K) + at time T. We assume that T = year and that the asset price follows a geometric Brownian motion S t = S expσb t + ct) where c is a constant, S and σ are positive constants and B t is a Brownian motion. The arbitragefree time- value of this option is then equal to exp ρ)e[ S expσb t + αt) dt K) + ], where ρ is the risk-free interest rate assumed constant) and α = ρ 2 σ2 see Baxter & Rennie 996) for an introduction to arbitrage-free valuation). Thus the problem of valuing the option boils down to calculating E[S expσb t +αt) dt K) + ]. Another type of Asian option, the floating-strike option, pays out S t dt S ) +, and leads to the consideration of E[ expσb t + αt) dt expσb + α)) + ]. Several approaches to the problem of valuing Asian options have been tried. Carverhill & Clewlow 99) use a convolution method to compute the distribution of S t dt, and work by Yor 992) and Geman & Yor 993) has lead to a formula for the price as a triple integral. Approximate formulae and bounds in the form of single and double integrals
2 by Levy 992), Levy & Turnbull 992), Curran 992) and Rogers & Shi 995) seem to be faster to evaluate however, and the methods of this paper fall into this category. Using intuition and simple optimization we derive bounds on the value of fixed-strike and floating-strike Asian options which can be computed quickly, and which are accurate for typical parameter values. The method generalizes to other options on sums of lognormal assets: discretely monitored Asians, currency basket options, and swaptions in Gaussian HJM models, for example. The plan of this paper is as follows: in Section 2 we state some useful facts about covariances related to Brownian motion, and in Section 3 present the impressive lower bound of Rogers & Shi 995) for fixed-strike and floating-strike options, and the approximation to the fixed-strike bound given by Curran 992). This approximation is notable since it is very close to the bound of Rogers & Shi 995) and is much easier to compute, involving only one-dimensional integrals rather than a troublesome two-dimensional integral. We give an alternative derivation of the bounds of Rogers & Shi 995), leading to an expression involving only one-dimensional integrals; we also provide a generalization Curran s approximation to floating-strike options. In Section 4 we derive upper bounds to complement the lower bounds, and in Section 5 we present a numerical comparison of the various bounds discussed here, using the parameter values from Curran 992). 2 Useful covariances Like most of the approximation formulae in the literature, we will exploit the high correlation between expσb t + αt) dt and B t dt for the values of σ and α met in practice, and the fact that the second integral is Gaussian. We will need the covariance matrix of the bivariate Gaussian random variable B t, B s ds): ) T E B t, B s ds B t, B s ds) ) ) t t t/2) =. t t/2) 3 Thus the conditional distribution of B t given B s ds = z is normal with mean 3t t/2)z and variance t 3t 2 t/2) 2, and the conditional distribution of B s ds given B t = x is normal with mean t/2)x and variance 3 t t/2)2. We will also need the covariance matrix of B t, B s ds B ): E B t, B s ds B ) T B t, B s ds B ) ) = t t 2 /2 t 2 /2 3 ), 2
3 so, given B s ds B = x, the conditional distribution of B t is normal with mean 3t 2 x/2 and variance t 3t 4 /4. 3 Lower bounds In this section we present two derivations of the bound of Rogers & Shi 995): their own, and an alternative, which yields a simpler expression for the bound. It also leads to the bound of Curran 992) for fixed-strike options. We show how a similar bound for floating-strike options may be derived. The derivation of Rogers & Shi 995) exploits the inequality: E [ A +] = E [ E A + C )] E [ E A C)) +] which holds for any random variables A and C. For fixed-strike Asian option they choose A = S t dt K and C = B t dt. Since the inner expectation is E S exp σb t + αt) ) B s ds dt K and, conditional B t dt = z, B t is normal with mean 3t t/2)z and variance t 3t 2 t/2) 2, we have the bound ) [ V fixed e ρ ] + 3φ 3z Se 3σt t/2)z+αt+ 2 σ2 t 3t 2 t/2) 2 ) dt K dz. 3.) To bound a floating-strike option they use A = S t dt S, C = B t dt B, and get ) [ V floating e ρ 3φ 3z S exp 3σt 2 z/2 + αt + 2 σ2 t 3t 4 /4 )) dt S exp 3σz/2 + α + σ 2 /8 ))] + dz. Both of these formulae are slightly tricky to evaluate since the outer integration has a non-smooth integrand. An alternative approach is to approximate the event that the option eventually ends in-the-money with something more tractable. Let A = {ω : S t dt > K}, and note that [ ) + ] E S exp σb t + αt) dt K = E [S exp σb t + αt) K) IA)] dt. 3.2) If we replace A by some other event A, we no longer have equality in 3.2); the right hand side is now a lower bound. We will use A = { B t dt > γ}. This Gaussian form allows the expectation to be written as a Black-Scholes type formula once γ has been determined, and just leaves us with a one-dimensional integral of a smooth integrand, which should be very fast. 3
4 To determine the optimal value of γ, let N t = σb t + αt + log S and note that for any random variable X with density f X x) γ E expn t ) K; X > γ) dt = Thus the optimal value of γ, γ satisfies E exp N t ) K X = γ) f X γ)) dt. With our choice of X = B t dt, we conclude that EexpN t ) X = γ ) dt = K. 3.3) S exp 3γ σt t/2) + αt + 2 σ2 t 3t 2 t/2) 2)) dt = K, 3.4) which determines γ uniquely. We now have the bound [ Se V fixed e ρ σb E t+αt K ) )] I B s ds > γ dt, and it remains to calculate the expectation. Fix t, ) and let N = σb t + αt + log S and N 2 = B s ds γ. Write µ i = EN i ), σi 2 = VarN i) and c = CovN, N 2 ), then using E [ e N K ) I N 2 > ) ] = e µ + 2 σ2 Φ µ2 + c σ 2 ) KΦ where Φ is the normal distribution function, and substituting µ = αt + log S, µ 2 = γ, σ 2 = σ2 t, σ2 2 = 3, c = σt t/2), we have [ V fixed e ρ Se αt+ γ 2 σ2t ) + σt t/2) Φ / 3 µ2 σ 2 ), )] γ dt KΦ /. 3 Integrating this numerically is significantly easier than integrating 3.). To see that this bound gives the same answers as that of Rogers & Shi 995), let Y = S t dt, Z = B t dt and note that E[EY K Z)) + ] = E[EY K Z))IEY K Z) > )], and since from 3.3) and 3.4), EY K Z) is strictly increasing in Z, we have: EY K Z) > if and only if Z > γ. Thus γ satisfies E[EY K Z)) + ] = E[EY K Z))IZ > γ )] which is just E[Y K)IZ > γ )]. The bound of Curran 992) arises from solving 3.4) approximately, using the following method: let fγ) = E S expσb t + αt) dt B s ds = γ), and note that a reasonable approximation to f is fγ) := S expγσ + α/2), obtained by interchanging the orders of integration and exponentiation. Recall that we seek γ = f K) and observe that if f f 4
5 then f x) f 2x f f x)), the approximation being exact if f f is constant. Thus γ f 2K f f K)), giving 2K γ σ [log S ) ] e 3logK/S) α/2)t t/2)+αt+ 2 σ2 t 3t 2 t/2) 2) dt α/2, which is the continuous limit of the bound given by Curran 992). For the floating-strike option, let A = {ω : S t dt > S } and use an approximation to A of the form A = { B t dt B > γ}. With this choice, γ, the optimal value of γ, satisfies [ ] E S exp σb t + αt) S exp σb + α)) dt B s ds B = γ = giving exp 3γ σ/2 α σ 2 /8 ) exp 3γ σt 2 /2 + αt + 2 σ2 t 3t 4 /4 )) dt =, 3.5) which has a unique solution. Our lower bound for the floating-strike case is thus [ S V floating e ρ E exp σbt + αt) S exp σb + α) ) )] I B s ds B > γ dt, which reduces to [ V floating e ρ Se αt+ γ 2 σ2t σt 2 ) /2 Φ / dt Se α+ γ )] /2 2 σ2 Φ 3 /. 3 Again this bound gives the same answers as that of Rogers & Shi 995). To see this, let Y = S t dt, Z = B t dt, and note first that γ E Y S Z B = γ) = 3σt 2 /2) expat, γ)) dt + 3σ/2) expa, γ)), where At, γ) = log S 3γσt 2 /2 + αt + 2 σ2 t 3t 4 /4 ). Since γ solves EY S Z B = γ ) =, we have expat, γ )) dt = expa, γ )), so γ E Y S Z B = γ) = γ=γ 3σ t 2 )/2) expat, γ )) dt >. Thus E Y S Z B = γ) > if and only if γ > γ. A similar argument to the fixedstrike case completes the proof. We can now generalize Curran s formula to the case of floating-strike options by solving 3.5) approximately. Let fγ) = e 3γσ/2 α σ2 /8 E eσbt+αt B t dt B = γ) and let 5
6 fγ) = expσγ α/2) be the approximation to f obtained by interchanging the orders of integration and exponentiation. Then the solution to 3.5) is given approximately by γ f 2 f f )), giving γ σ [α/2 + log 2 e α/4 σ2 /8 4 Upper bounds exp αt 3αt 2 /4 + 2 σ2 t 3t 4 /4 )) )] dt. In this section we derive a new upper bound on the value of fixed-strike and floating-strike Asian options in the form of a double integral. Rogers & Shi 995) obtained an upper bound by considering the error made by their lower bound see Section 3). As an indication of their relative accuracy, if σ =.3, ρ =.9, S =, and T = year and k =, the lower bounds of Section 3 are both and the upper bound of this section is By comparison, the upper bound given by Rogers & Shi 995) is The inequality underlying the new bound is the following: let X be a random variable, and let f t ω) be a random function with f tω) dt = for all ω. Then [ ) + ] [ ) + ] E S expσb t + αt) dt X = E S expσb t + αt) Xf t ) dt [ ] E S expσb t + αt) Xf t ) + dt = E [ S exp σb t + αt) Xf t ) +] dt. 4.) For both the fixed-strike and floating-strike cases we will use f t = µ t + σb t B s ds), where µ t is a deterministic function satisfying µ t dt = ; and derive an expression for µ t which is approximately optimal in each case. As the bounds have a very similar derivation we will concentrate on the fixed-strike option and give the appropriate modifications for the floating-strike case at the end of the section. Take X = K in 4.) and first consider the choice f t = µ t. To choose µ t we will minimize the right hand side of 4.) over the set of deterministic functions f t such that f t dt =. Let Lλ, {f t }) = E[ S expσb t + αt) Kf t ) + dt λ f t dt )] be the Lagrangian, and consider stationarity with respect to {f t } for the unconstrained problem. This gives the condition KP S exp σb t + αt) Kf t ) λ) ɛ t dt =, where {ɛ t } is some small deterministic perturbation. Thus we see that PS expσb t + αt) Kf t ) must be independent of t. Equivalently we have logkf t /S) αt = γσ t for some 6
7 constant γ. Thus the optimal choice for f t is f t = S/K) exp σγ t + αt ), where the constant γ is chosen so that f t dt =. Since f t dt is monotone increasing in γ, the correct value for γ is easy to estimate numerically. If instead f t = µ t + σb t B s ds), the condition for stationarity with respect to small deterministic perturbations is [ P S expσb t + αt) K µ t + σ B t ))] B s ds = λ, for all t, 4.2) but this cannot easily be re-arranged to give the dependence of µ t on λ. Instead we will use the approximation expσb t ) + σb t which should be reasonable for small σ. This leads to the condition PS expαt) + σb t ) Kf t ) = λ for all t. Letting N t = S expαt) + S expαt)σ Kσ)B t + Kσ B s ds, we conclude that PN t Kµ t ) must be independent of t. Using the facts about the joint distribution of B t, B s ds) given in Section 2, we deduce that where Imposing µ t dt = gives µ t = K S exp αt) + γ v t ), 4.3) v t = VarN t ) = c 2 t t + 2Kσ)c t t t/2) + Kσ) 2 /3, 4.4) c t = S expαt)σ Kσ. 4.5) γ = ) / K Se α )/α vt dt. 4.6) We estimate the integral vt dt numerically. We now know the constant γ and hence the function µ t for our upper bound: [ V fixed e ρ E S exp σb t + αt) K µ t + σ B t ))) + ] B s ds dt. Conditioning on B t = x, this becomes ) V fixed e ρ x φ E [ at, x) + bt, x)n) +] dx dt, 4.7) t t 7
8 where N has a N, ) distribution, and the functions a and b are given by at, x) = S exp σx + αt) Kµ t + σx) + Kσ t/2)x, 4.8) bt, x) = Kσ 3 t t/2)2. 4.9) The calculation of E[a + bn) + ] is straightforward and gives aφa/b) + bφa/b). In the form of 4.7) the integrand is badly behaved near, ) so we perform the change of variables v = t, w = x/ t, giving V fixed e ρ [ 2vφw) at, x)φ ) at, x) + bt, x)φ bt, x) )] at, x) dw dv. bt, x) This expression, combined with 4.3), 4.4), 4.5), 4.6), 4.8) and 4.9), constitutes the upper bound in the fixed-strike case. For the case of a floating-strike option, we now take X = S. Setting Z = B t dt, the condition for stationarity with respect to small deterministic perturbations analogous to 4.2) is E [ S expσb + α); S expσb t + αt) µ t + σ B t Z))S expσb + α))] = λ, t. We approximate this by P[S expσb t + αt) µ t + σb t Z))S expσb + α)] = λ, t, and further, using the approximation expσb t B )) + σb t B ), by P[expαt )) + σb t B ) µ t + σb t Z))] = λ, t. This implies that for some γ, µ t = exp α t )) + γ v t where v t = Var[expαt )) + σb t B )) σb t Z)]. Since µ t dt = we must have γ = exp α))/α)/ vt dt. Our bound is thus [ Se V floating e ρ σb E t+αt µ t + σ B t Z)) Se σb +α ) ] + dt [ ) ] = e ρ E e Nt) N 2 t)e N + 3t) dt, where N t) = σb t + αt + log S, N 2 t) = µ t + σb t B s ds) and N 3 t) = σb + α + log S. If we condition on N 2 t) = x we can perform the remaining expectation analytically. Let 8
9 µ i t) = EN i t)), σ ij t) = CovN i t), N j t)) and denote by tildes the conditional distributions given N 2 t) = x: µ i t, x) = µ i + x µ 2 )σ i2 /σ 22 and σ ij = σ ij σ i2 σ j2 /σ 22. Finally let v 2 = VarN t) N 3 t) N 2 t) = x) = σ 2 σ 3 + σ 33. Our upper bound on the price of a floating-strike Asian option is then V floating e ρ where we take logx) = for x. x )[e µ ) µ2 φ + 2 σ µ µ 3 logx) + σ σ 3 Φ σ22 σ22 v )] xe µ 3+ 2 σ 33 µ µ 3 logx) + σ 3 σ 33 ) Φ dx dt, v 5 Numerical Results In Table we consider fixed-strike options and show the upper bound of Rogers & Shi 995), the upper and lower bounds derived in Sections 3 and 4, the approximation of Curran 992) and the Monte-Carlo results of Levy & Turnbull 992). All calculations assume ρ =.9, an initial stock price of S = and an expiry time of year. For the lower bound and the new upper bound, the approximate time taken on an HP 9/73) is parenthesized; for the Monte-Carlo studies, the estimated standard error is bracketed beneath. In Table 2 we show how the upper and lower bounds of Rogers & Shi 995) compare to the upper bound of Section 4 and the generalization to floating-strike options of Curran s lower bound, described in Section 3. The approximate time taken is parenthesized. References Baxter, M. W. & Rennie, A. 996), Financial Calculus, Cambridge University Press. Carverhill, A. & Clewlow, L. 99), Flexible convolution, Risk 34), Curran, M. 992), Beyond average intelligence, Risk 5), 6. Geman, H. & Yor, M. 993), Bessel processes, Asian options and Perpetuities, Mathematical Finance 34), Levy, E. 992), Pricing European average rate currency options, J. Intnl. Money and Finance, Levy, E. & Turnbull, S. 992), Average Intelligence, Risk 52),
10 Rogers, L. C. G. & Shi, Z. 995), The value of an Asian option, J. Appl. Prob. 32, Yor, M. 992), Some aspects of Brownian motion, Part : Some Special Functionals, Birkäuser.
11 Volatility σ Strike K Curran lower R S lower M C result Upper bound R S upper ).9) [.].3) ).) [.].9) ).) [.].9) ).8) [.].9) ).7) [.].2) ).8) [.].2) ).9) [.].24) ).9) [.].24) ).8) [.].28) ).9) [.3].28) ).8) [.3].63) ).8) [.3].64) Table Comparison of various bounds on fixed-strike Asian option prices for S =, ρ =.9, and an expiry time of year. Parenthesized numbers are computation times in seconds, bracketed numbers are estimates of standard errors from Curran 992)).
12 Volatility σ Interest Rate ρ Generalized R S lower Upper bound R S upper ).36).27) ).54).25) ).78).24) ).24).3) ).32).33) ).44).28) ).8).28) ).22).33) ).32).33) Table 2 Comparison of bounds on floating-strike Asian option prices for S = with an expiry time of year. Parenthesized numbers are computation times in seconds. The parameter values are those used in Rogers & Shi 995). 2
Probability 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 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 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 informationAnalytical formulas for local volatility model with stochastic. Mohammed Miri
Analytical formulas for local volatility model with stochastic rates Mohammed Miri Joint work with Eric Benhamou (Pricing Partners) and Emmanuel Gobet (Ecole Polytechnique Modeling and Managing Financial
More informationStochastic Processes and Stochastic Calculus - 9 Complete and Incomplete Market Models
Stochastic Processes and Stochastic Calculus - 9 Complete and Incomplete Market Models Eni Musta Università degli studi di Pisa San Miniato - 16 September 2016 Overview 1 Self-financing portfolio 2 Complete
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 informationThe stochastic calculus
Gdansk A schedule of the lecture Stochastic differential equations Ito calculus, Ito process Ornstein - Uhlenbeck (OU) process Heston model Stopping time for OU process Stochastic differential equations
More informationNEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 MAS3904. Stochastic Financial Modelling. Time allowed: 2 hours
NEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 Stochastic Financial Modelling Time allowed: 2 hours Candidates should attempt all questions. Marks for each question
More informationMultiname and Multiscale Default Modeling
Multiname and Multiscale Default Modeling Jean-Pierre Fouque University of California Santa Barbara Joint work with R. Sircar (Princeton) and K. Sølna (UC Irvine) Special Semester on Stochastics with Emphasis
More informationAMH4 - ADVANCED OPTION PRICING. Contents
AMH4 - ADVANCED OPTION PRICING ANDREW TULLOCH Contents 1. Theory of Option Pricing 2 2. Black-Scholes PDE Method 4 3. Martingale method 4 4. Monte Carlo methods 5 4.1. Method of antithetic variances 5
More informationEstimating the Greeks
IEOR E4703: Monte-Carlo Simulation Columbia University Estimating the Greeks c 207 by Martin Haugh In these lecture notes we discuss the use of Monte-Carlo simulation for the estimation of sensitivities
More informationRisk Neutral Measures
CHPTER 4 Risk Neutral Measures Our aim in this section is to show how risk neutral measures can be used to price derivative securities. The key advantage is that under a risk neutral measure the discounted
More informationPractical example of an Economic Scenario Generator
Practical example of an Economic Scenario Generator Martin Schenk Actuarial & Insurance Solutions SAV 7 March 2014 Agenda Introduction Deterministic vs. stochastic approach Mathematical model Application
More informationThe British Russian Option
The British Russian Option Kristoffer J Glover (Joint work with G. Peskir and F. Samee) School of Finance and Economics University of Technology, Sydney 25th June 2010 (6th World Congress of the BFS, Toronto)
More information2 Control variates. λe λti λe e λt i where R(t) = t Y 1 Y N(t) is the time from the last event to t. L t = e λr(t) e e λt(t) Exercises
96 ChapterVI. Variance Reduction Methods stochastic volatility ISExSoren5.9 Example.5 (compound poisson processes) Let X(t) = Y + + Y N(t) where {N(t)},Y, Y,... are independent, {N(t)} is Poisson(λ) with
More informationOption Pricing Models for European Options
Chapter 2 Option Pricing Models for European Options 2.1 Continuous-time Model: Black-Scholes Model 2.1.1 Black-Scholes Assumptions We list the assumptions that we make for most of this notes. 1. The underlying
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 informationComputer Exercise 2 Simulation
Lund University with Lund Institute of Technology Valuation of Derivative Assets Centre for Mathematical Sciences, Mathematical Statistics Fall 2017 Computer Exercise 2 Simulation This lab deals with pricing
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Simulating Stochastic Differential Equations Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationFinancial Risk Management
Financial Risk Management Professor: Thierry Roncalli Evry University Assistant: Enareta Kurtbegu Evry University Tutorial exercices #4 1 Correlation and copulas 1. The bivariate Gaussian copula is given
More informationLIBOR models, multi-curve extensions, and the pricing of callable structured derivatives
Weierstrass Institute for Applied Analysis and Stochastics LIBOR models, multi-curve extensions, and the pricing of callable structured derivatives John Schoenmakers 9th Summer School in Mathematical Finance
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 informationS t d with probability (1 p), where
Stochastic Calculus Week 3 Topics: Towards Black-Scholes Stochastic Processes Brownian Motion Conditional Expectations Continuous-time Martingales Towards Black Scholes Suppose again that S t+δt equals
More informationAn Analytical Approximation for Pricing VWAP Options
.... An Analytical Approximation for Pricing VWAP Options Hideharu Funahashi and Masaaki Kijima Graduate School of Social Sciences, Tokyo Metropolitan University September 4, 215 Kijima (TMU Pricing of
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 informationThe value of foresight
Philip Ernst Department of Statistics, Rice University Support from NSF-DMS-1811936 (co-pi F. Viens) and ONR-N00014-18-1-2192 gratefully acknowledged. IMA Financial and Economic Applications June 11, 2018
More informationLecture 4. Finite difference and finite element methods
Finite difference and finite element methods Lecture 4 Outline Black-Scholes equation From expectation to PDE Goal: compute the value of European option with payoff g which is the conditional expectation
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 informationThe University of Chicago, Booth School of Business Business 41202, Spring Quarter 2011, Mr. Ruey S. Tsay. Solutions to Final Exam.
The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2011, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (32 pts) Answer briefly the following questions. 1. Suppose
More informationLocal Volatility Dynamic Models
René Carmona Bendheim Center for Finance Department of Operations Research & Financial Engineering Princeton University Columbia November 9, 27 Contents Joint work with Sergey Nadtochyi Motivation 1 Understanding
More information1 Implied Volatility from Local Volatility
Abstract We try to understand the Berestycki, Busca, and Florent () (BBF) result in the context of the work presented in Lectures and. Implied Volatility from Local Volatility. Current Plan as of March
More informationComputational Finance
Path Dependent Options Computational Finance School of Mathematics 2018 The Random Walk One of the main assumption of the Black-Scholes framework is that the underlying stock price follows a random walk
More informationAsymptotic methods in risk management. Advances in Financial Mathematics
Asymptotic methods in risk management Peter Tankov Based on joint work with A. Gulisashvili Advances in Financial Mathematics Paris, January 7 10, 2014 Peter Tankov (Université Paris Diderot) Asymptotic
More informationlast problem outlines how the Black Scholes PDE (and its derivation) may be modified to account for the payment of stock dividends.
224 10 Arbitrage and SDEs last problem outlines how the Black Scholes PDE (and its derivation) may be modified to account for the payment of stock dividends. 10.1 (Calculation of Delta First and Finest
More informationEconomathematics. Problem Sheet 1. Zbigniew Palmowski. Ws 2 dw s = 1 t
Economathematics Problem Sheet 1 Zbigniew Palmowski 1. Calculate Ee X where X is a gaussian random variable with mean µ and volatility σ >.. Verify that where W is a Wiener process. Ws dw s = 1 3 W t 3
More informationKØBENHAVNS UNIVERSITET (Blok 2, 2011/2012) Naturvidenskabelig kandidateksamen Continuous time finance (FinKont) TIME ALLOWED : 3 hours
This question paper consists of 3 printed pages FinKont KØBENHAVNS UNIVERSITET (Blok 2, 211/212) Naturvidenskabelig kandidateksamen Continuous time finance (FinKont) TIME ALLOWED : 3 hours This exam paper
More informationApproximation Methods in Derivatives Pricing
Approximation Methods in Derivatives Pricing Minqiang Li Bloomberg LP September 24, 2013 1 / 27 Outline of the talk A brief overview of approximation methods Timer option price approximation Perpetual
More informationBlack-Scholes Option Pricing
Black-Scholes Option Pricing The pricing kernel furnishes an alternate derivation of the Black-Scholes formula for the price of a call option. Arbitrage is again the foundation for the theory. 1 Risk-Free
More informationLecture Notes for Chapter 6. 1 Prototype model: a one-step binomial tree
Lecture Notes for Chapter 6 This is the chapter that brings together the mathematical tools (Brownian motion, Itô calculus) and the financial justifications (no-arbitrage pricing) to produce the derivative
More informationDefinition Pricing Risk management Second generation barrier options. Barrier Options. Arfima Financial Solutions
Arfima Financial Solutions Contents Definition 1 Definition 2 3 4 Contenido Definition 1 Definition 2 3 4 Definition Definition: A barrier option is an option on the underlying asset that is activated
More informationExtended Libor Models and Their Calibration
Extended Libor Models and Their Calibration Denis Belomestny Weierstraß Institute Berlin Vienna, 16 November 2007 Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Vienna, 16 November
More informationOptimal trading strategies under arbitrage
Optimal trading strategies under arbitrage Johannes Ruf Columbia University, Department of Statistics The Third Western Conference in Mathematical Finance November 14, 2009 How should an investor trade
More informationThe Use of Importance Sampling to Speed Up Stochastic Volatility Simulations
The Use of Importance Sampling to Speed Up Stochastic Volatility Simulations Stan Stilger June 6, 1 Fouque and Tullie use importance sampling for variance reduction in stochastic volatility simulations.
More informationStock Loan Valuation Under Brownian-Motion Based and Markov Chain Stock Models
Stock Loan Valuation Under Brownian-Motion Based and Markov Chain Stock Models David Prager 1 1 Associate Professor of Mathematics Anderson University (SC) Based on joint work with Professor Qing Zhang,
More informationChapter 14. The Multi-Underlying Black-Scholes Model and Correlation
Chapter 4 The Multi-Underlying Black-Scholes Model and Correlation So far we have discussed single asset options, the payoff function depended only on one underlying. Now we want to allow multiple underlyings.
More informationRohini Kumar. Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque)
Small time asymptotics for fast mean-reverting stochastic volatility models Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque) March 11, 2011 Frontier Probability Days,
More informationOption Pricing. Chapter Discrete Time
Chapter 7 Option Pricing 7.1 Discrete Time In the next section we will discuss the Black Scholes formula. To prepare for that, we will consider the much simpler problem of pricing options when there are
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 informationOptimal stopping problems for a Brownian motion with a disorder on a finite interval
Optimal stopping problems for a Brownian motion with a disorder on a finite interval A. N. Shiryaev M. V. Zhitlukhin arxiv:1212.379v1 [math.st] 15 Dec 212 December 18, 212 Abstract We consider optimal
More informationMARTINGALES AND LOCAL MARTINGALES
MARINGALES AND LOCAL MARINGALES If S t is a (discounted) securtity, the discounted P/L V t = need not be a martingale. t θ u ds u Can V t be a valid P/L? When? Winter 25 1 Per A. Mykland ARBIRAGE WIH SOCHASIC
More informationApplication of Stochastic Calculus to Price a Quanto Spread
Application of Stochastic Calculus to Price a Quanto Spread Christopher Ting http://www.mysmu.edu/faculty/christophert/ Algorithmic Quantitative Finance July 15, 2017 Christopher Ting July 15, 2017 1/33
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 informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Generating Random Variables and Stochastic Processes Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationSimulating Stochastic Differential Equations
IEOR E4603: Monte-Carlo Simulation c 2017 by Martin Haugh Columbia University Simulating Stochastic Differential Equations In these lecture notes we discuss the simulation of stochastic differential equations
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 informationThe Black-Scholes Equation using Heat Equation
The Black-Scholes Equation using Heat Equation Peter Cassar May 0, 05 Assumptions of the Black-Scholes Model We have a risk free asset given by the price process, dbt = rbt The asset price follows a geometric
More informationCourse information FN3142 Quantitative finance
Course information 015 16 FN314 Quantitative finance This course is aimed at students interested in obtaining a thorough grounding in market finance and related empirical methods. Prerequisite If taken
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 informationMarkovian Projection, Heston Model and Pricing of European Basket Optio
Markovian Projection, Heston Model and Pricing of European Basket Options with Smile July 7, 2009 European Options on geometric baskets European Options on arithmetic baskets Pricing via Moment Matching
More informationLecture Note 8 of Bus 41202, Spring 2017: Stochastic Diffusion Equation & Option Pricing
Lecture Note 8 of Bus 41202, Spring 2017: Stochastic Diffusion Equation & Option Pricing We shall go over this note quickly due to time constraints. Key concept: Ito s lemma Stock Options: A contract giving
More informationIMPA Commodities Course : Forward Price Models
IMPA Commodities Course : Forward Price Models Sebastian Jaimungal sebastian.jaimungal@utoronto.ca Department of Statistics and Mathematical Finance Program, University of Toronto, Toronto, Canada http://www.utstat.utoronto.ca/sjaimung
More information2.1 Mean-variance Analysis: Single-period Model
Chapter Portfolio Selection The theory of option pricing is a theory of deterministic returns: we hedge our option with the underlying to eliminate risk, and our resulting risk-free portfolio then earns
More informationAsset-based Estimates for Default Probabilities for Commercial Banks
Asset-based Estimates for Default Probabilities for Commercial Banks Statistical Laboratory, University of Cambridge September 2005 Outline Structural Models Structural Models Model Inputs and Outputs
More informationA new PDE approach for pricing arithmetic average Asian options
A new PDE approach for pricing arithmetic average Asian options Jan Večeř Department of Mathematical Sciences, Carnegie Mellon University, Pittsburgh, PA 15213. Email: vecer@andrew.cmu.edu. May 15, 21
More informationAmerican Option Pricing Formula for Uncertain Financial Market
American Option Pricing Formula for Uncertain Financial Market Xiaowei Chen Uncertainty Theory Laboratory, Department of Mathematical Sciences Tsinghua University, Beijing 184, China chenxw7@mailstsinghuaeducn
More informationNumerical schemes for SDEs
Lecture 5 Numerical schemes for SDEs Lecture Notes by Jan Palczewski Computational Finance p. 1 A Stochastic Differential Equation (SDE) is an object of the following type dx t = a(t,x t )dt + b(t,x t
More informationEMPIRICAL EVIDENCE ON ARBITRAGE BY CHANGING THE STOCK EXCHANGE
Advances and Applications in Statistics Volume, Number, This paper is available online at http://www.pphmj.com 9 Pushpa Publishing House EMPIRICAL EVIDENCE ON ARBITRAGE BY CHANGING THE STOCK EXCHANGE JOSÉ
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 informationHedging under arbitrage
Hedging under arbitrage Johannes Ruf Columbia University, Department of Statistics AnStAp10 August 12, 2010 Motivation Usually, there are several trading strategies at one s disposal to obtain a given
More informationDeriving the Black-Scholes Equation and Basic Mathematical Finance
Deriving the Black-Scholes Equation and Basic Mathematical Finance Nikita Filippov June, 7 Introduction In the 97 s Fischer Black and Myron Scholes published a model which would attempt to tackle the issue
More informationSPDE and portfolio choice (joint work with M. Musiela) Princeton University. Thaleia Zariphopoulou The University of Texas at Austin
SPDE and portfolio choice (joint work with M. Musiela) Princeton University November 2007 Thaleia Zariphopoulou The University of Texas at Austin 1 Performance measurement of investment strategies 2 Market
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 informationarxiv: v2 [q-fin.pr] 23 Nov 2017
VALUATION OF EQUITY WARRANTS FOR UNCERTAIN FINANCIAL MARKET FOAD SHOKROLLAHI arxiv:17118356v2 [q-finpr] 23 Nov 217 Department of Mathematics and Statistics, University of Vaasa, PO Box 7, FIN-6511 Vaasa,
More informationKing s College London
King s College London University Of London This paper is part of an examination of the College counting towards the award of a degree. Examinations are governed by the College Regulations under the authority
More informationBROWNIAN MOTION Antonella Basso, Martina Nardon
BROWNIAN MOTION Antonella Basso, Martina Nardon basso@unive.it, mnardon@unive.it Department of Applied Mathematics University Ca Foscari Venice Brownian motion p. 1 Brownian motion Brownian motion plays
More information1 The continuous time limit
Derivative Securities, Courant Institute, Fall 2008 http://www.math.nyu.edu/faculty/goodman/teaching/derivsec08/index.html Jonathan Goodman and Keith Lewis Supplementary notes and comments, Section 3 1
More informationReplication and Absence of Arbitrage in Non-Semimartingale Models
Replication and Absence of Arbitrage in Non-Semimartingale Models Matematiikan päivät, Tampere, 4-5. January 2006 Tommi Sottinen University of Helsinki 4.1.2006 Outline 1. The classical pricing model:
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 informationOption Pricing Formula for Fuzzy Financial Market
Journal of Uncertain Systems Vol.2, No., pp.7-2, 28 Online at: www.jus.org.uk Option Pricing Formula for Fuzzy Financial Market Zhongfeng Qin, Xiang Li Department of Mathematical Sciences Tsinghua University,
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 informationCharacterization of the Optimum
ECO 317 Economics of Uncertainty Fall Term 2009 Notes for lectures 5. Portfolio Allocation with One Riskless, One Risky Asset Characterization of the Optimum Consider a risk-averse, expected-utility-maximizing
More informationCalculating Implied Volatility
Statistical Laboratory University of Cambridge University of Cambridge Mathematics and Big Data Showcase 20 April 2016 How much is an option worth? A call option is the right, but not the obligation, to
More informationKing s College London
King s College London University Of London This paper is part of an examination of the College counting towards the award of a degree. Examinations are governed by the College Regulations under the authority
More informationCrashcourse Interest Rate Models
Crashcourse Interest Rate Models Stefan Gerhold August 30, 2006 Interest Rate Models Model the evolution of the yield curve Can be used for forecasting the future yield curve or for pricing interest rate
More informationArbitrages and pricing of stock options
Arbitrages and pricing of stock options Gonzalo Mateos Dept. of ECE and Goergen Institute for Data Science University of Rochester gmateosb@ece.rochester.edu http://www.ece.rochester.edu/~gmateosb/ November
More informationContinuous Time Finance. Tomas Björk
Continuous Time Finance Tomas Björk 1 II Stochastic Calculus Tomas Björk 2 Typical Setup Take as given the market price process, S(t), of some underlying asset. S(t) = price, at t, per unit of underlying
More informationInterest Rate Volatility
Interest Rate Volatility III. Working with SABR Andrew Lesniewski Baruch College and Posnania Inc First Baruch Volatility Workshop New York June 16-18, 2015 Outline Arbitrage free SABR 1 Arbitrage free
More information(Appeared in: Quarterly Journal of the Institute of Actuaries of Australia 2: 2-17, 1996.) FROM COMPOUND INTEREST TO ASIAN OPTIONS
(Appeared in: Quarterly Journal of the Institute of Actuaries of Australia 2: 2-17, 1996.) FROM COMPOUND INTEREST TO ASIAN OPTIONS by Daniel Dufresne, University of Montreal Abstract This paper describes
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 informationFinancial Engineering with FRONT ARENA
Introduction The course A typical lecture Concluding remarks Problems and solutions Dmitrii Silvestrov Anatoliy Malyarenko Department of Mathematics and Physics Mälardalen University December 10, 2004/Front
More informationBasic Arbitrage Theory KTH Tomas Björk
Basic Arbitrage Theory KTH 2010 Tomas Björk Tomas Björk, 2010 Contents 1. Mathematics recap. (Ch 10-12) 2. Recap of the martingale approach. (Ch 10-12) 3. Change of numeraire. (Ch 26) Björk,T. Arbitrage
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 informationHedging of swaptions in a Lévy driven Heath-Jarrow-Morton framework
Hedging of swaptions in a Lévy driven Heath-Jarrow-Morton framework Kathrin Glau, Nele Vandaele, Michèle Vanmaele Bachelier Finance Society World Congress 2010 June 22-26, 2010 Nele Vandaele Hedging of
More informationMath489/889 Stochastic Processes and Advanced Mathematical Finance Solutions to Practice Problems
Math489/889 Stochastic Processes and Advanced Mathematical Finance Solutions to Practice Problems Steve Dunbar No Due Date: Practice Only. Find the mode (the value of the independent variable with the
More informationA Moment Matching Approach To The Valuation Of A Volume Weighted Average Price Option
A Moment Matching Approach To The Valuation Of A Volume Weighted Average Price Option Antony Stace Department of Mathematics and MASCOS University of Queensland 15th October 2004 AUSTRALIAN RESEARCH COUNCIL
More informationJDEP 384H: Numerical Methods in Business
Chapter 4: Numerical Integration: Deterministic and Monte Carlo Methods Chapter 8: Option Pricing by Monte Carlo Methods JDEP 384H: Numerical Methods in Business Instructor: Thomas Shores Department of
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 19 11/20/2013. Applications of Ito calculus to finance
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.7J Fall 213 Lecture 19 11/2/213 Applications of Ito calculus to finance Content. 1. Trading strategies 2. Black-Scholes option pricing formula 1 Security
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 informationIlliquidity, Credit risk and Merton s model
Illiquidity, Credit risk and Merton s model (joint work with J. Dong and L. Korobenko) A. Deniz Sezer University of Calgary April 28, 2016 Merton s model of corporate debt A corporate bond is a contingent
More informationOption Approach to Risk-shifting Incentive Problem with Mutually Correlated Projects
Option Approach to Risk-shifting Incentive Problem with Mutually Correlated Projects Hiroshi Inoue 1, Zhanwei Yang 1, Masatoshi Miyake 1 School of Management, T okyo University of Science, Kuki-shi Saitama
More information