MLEMVD: A R Package for Maximum Likelihood Estimation of Multivariate Diffusion Models
|
|
- Kristopher Harvey
- 6 years ago
- Views:
Transcription
1 MLEMVD: A R Package for Maximum Likelihood Estimation of Multivariate Diffusion Models Matthew Dixon and Tao Wu 1 Illinois Institute of Technology May 19th id=
2 Time Homogeneous Markovian Diffusions (Sahalia2002) consider the multivariate time-homogenous Markovian diffusion of the form dx t = µ(x t )dt + Σ(X t )dw t, µ(x) : R m R m, Σ(x) : R m R m m W t R m are independent Wiener processes. Key idea: The log of the transition function f X (x x 0, ) can be given in closed form under mild restrictions on the form of µ and Σ.
3 Time Homogeneous Markovian Diffusions We observe X t at times t 0, t 1,..., t n, where denotes the difference between observation times and is assumed independent. Under this finite data, the log-likelihood takes the form: l n (p) := n l X (, x ti x ti 1 ; p), l X := ln f X. i=1 Under a Hermite expansion of l X and application of a number of transformation we get the following compact closed form expression with K terms: l (K) X (x x 0) = m 2 ln(2π ) D v (x) + C 1 X (x x0) + K k=0 C (k) X (x x 0) k k!, D v := 1 2 ln(det[v(x)]).
4 Example: Geometric Brownian Motion Given an initial value X 0 = x 0, X t evolves with drift parameter µ and volatility parameter σ dx t = µx t dt + σx t dw t. The transition function takes the form f X (x t x 0, t) = 1 exp (lnx t lnx 0 (µ σ 2 /2)t) 2 2πσt 2σ 2 t The exact log likelihood function, evaluated over a uniform time series of n observations of the state variable with spacing is given by l n (p) = n i=1 l X (, x ti x ti 1 ; p) = 1 2 n i=1 (ln(2π σ2 x 2 ti ) + (ln[x ti /x ti 1 ] (µ σ2 /2) )) 2 /(σ 2 ).
5 Example: Heston Model The Heston model describes the evolution of the log of the stock price s t = ln S t ds t = (a + by t )dt + Y t dw Q 1 (t), (1) dy t = κ (θ Y t )dt + σ Y t dw Q 2 (t), (2) where a = r d, b = 1 2.
6 Heston Model A key characteristic of the model is that the Wiener processes are correlated dw Q 1 dw Q 2 = ρdt => leverage effect. There are five parameters in the model κ: mean-reversion rate θ: long-term variance σ: volatility of variance ρ: instantaneous correlation between dw Q 1 and dw Q 2 y 0 : initial variance
7 Likelihood function estimation: full observance State vector X t := [s t, Y t ] has a transition density function for the conditional density of X t+ = x given X t = x 0 denoted by f X (, x x 0 ; p). The log likelihood function for observations at times t 0, t 1,..., t n is given by l n (p) = 1 n n l X (, x ti x ti 1 ; p), i=1 l X (, x x 0 ; p) := ln f X (, x x 0 ; p) and is given in closed form.
8 Likelihood function estimation: partial observance X t := [s t ; Y t ] is partially observed and hence l n can not be directly estimated from time series data. Approximate the likelihood of the observed state vector G t := [s t ; C t ], where C t is the ATM constant maturity option price. The transition density function for the conditional density of G t+ = g given G t = g 0 is now denoted by f G (, g g 0 ; p) and the log likelihood function is given by l n (p) := 1 n n l G (, g(t i ) g(t i 1 ); p). i=1
9 Likelihood function estimation: partial observance G t = f (X t ; p) represents the stock and option prices Under a change of variables from G t to X t, the log of the transition density f G can be expressed in terms of f X through a Jacobian J t to give: l G (, g g 0 ; p) := ln J t (, g g 0 ; p)+l X (, f 1 (g; p) f 1 (g 0 ; p); p) Goal: Find a numerical approximation to the inverse map f 1 (G; p)
10 Pricing The mapping between X t and G t is given by [ s f (X t ; p) = t C(S t, Y t, K, ; p) C( ) is the Heston model option price defined above and = T t is the constant time to maturity of the option. ],
11 Calibration: overview Since Y t, as previously mentioned, is unobserved, we seek the inverse [ ] [ ] f 1 s (G t ; p) = t st C 1 = (S t, C t, K, ; p) Use the map from G t to X t to imply the volatility y t from the observed state vector {g t } n i=1 of underlying prices and corresponding constant maturity, ATM option prices. Find the maximum likelihood estimate p. Y t
12 Calibration: Step 1 Fix p and find the corresponding implied volatility y p t i y p t i = arg min c ti C(S ti, y, K, ; p). y This error function is convex and so has a unique solution independent of the initial choice of y in a simple root finding method.
13 Calibration: Step 2 Using (y p t i, p) we compute the Jacobian of G t w.r.t. X t and the log of transition density function l X ( ) l G (, g ti g ti 1 ; p) = lnj ti + l X (, x p t i x p t i p), x p t i := [s ti ; y p t i ] The Jacobian is the vega of the option.
14 Calibration: Step 3 Steps 1 and 2 are now repeated for all remaining times t i and the log likelihood function is evaluated for the combination (y p, p). The value of yt p i 1 is used to initialize the solver for the least absolute error problem.
15 Calibration: Step 4 A new value of p is generated and Steps 1 and 3 are repeated until the likelihood function has been maximized by a numerical solver p arg min l n (p). p F
16 Geometric Brownian Motion log likelihood log likelihood param 1 param 2 Figure: The marginal log likelihood function with respect to each parameter of the geometric brownian diffusion model. The mean is represented by Parameter 1 and the volatility by Parameter 2.
17 Geometric Brownian Motion µ σ Actual Estimated Est. Std. Error Est. Huber Sandwich Error Error in Est. Std. Error Error in Est. Huber Sandwich Error Table: The correct parameters, the estimates and the standard error estimates using 500 simulated stock prices.
18 Geometric Brownian Motion Actual maximum log likelihood Est. maximum log likelihood L2 Norm of Score Error L2 Norm of Hessian Error L2 Norm of Information matrix Table: This table lists the numerical error in a selection of estimated values.
19 Figure: The marginal log likelihood function with respect to each parameter of the Heston model applied to the simulated underlying and prices. Parameters 1-4 represents ρ, κ, θ and σ. Heston Model log likelihood log likelihood param 1 param 2 log likelihood log likelihood param 3 param 4
20 Heston Model ρ κ θ σ Actual Estimated Std. Error Huber Sandwich Error Table: The correct parameters, the estimates and the lower bounds on the standard error using 50 simulated observations of the stock and the ATM option price.
21 Heston Model abs price.model price.sim [ / ] Last Jan 01 Feb Apr 02 May Jul 02 Aug Oct 01 Nov 05 Dec Figure: The absolute error in the calibrated model prices versus the simulated prices.
22 Heston Model abs vol.sim vol.model [ / ] Last Jan 01 Feb Apr 02 May Jul 02 Aug Oct 01 Nov 05 Dec Figure: The absolute error in the implied volatilities versus the simulated volatilities.
23 MLEMVD Package Installation library("devtools") install_github("mfrdixon/mlemvd")
24 MLEVMD Univariate Models Model µ σ constraints U1 x(a + bx) σx 3/2 U2 a + bx dx U3 b(a x) cx d U4 κ(α x) σx 1/2 U5 3 i=0 θ ix i γx ρ ρ 1 U6 a + bx + cx 2 + dx 3 f U7 κ(α x) σ a U8 1 x + a 0 + a 1 x + a 2 x 2 σx p ρ 1 a U9 1 x + a 0 + a 1 x + a 2 x 2 (b 0 + b 1 x + b 2 x b3 ) 1/2 a U10 1 x + a 0 + a 1 x + a 2 x 2 b 0 + b 1 x + b 2 x b3 ρ 1 U11 a + bx f + dx β U12 x αx 3 γx 1/2 a U13 1 x + a 0 + a 1 x + a 2 x 2 + a 3 x 3 σx ρ ρ 1 Table: The specification of various univariate diffusion models currently supported by the package.
25 MLEMVD Bivariate Models Model µ(x1, x2) Σ(x1, x2) constraints ( ) ( ) a + bx2 ρ (x2) B1 0 c + dx2 ( ) ( h (1 ρ 2 )x2 ) a0 + a1x1 + a2x2 c0 + c1x1 + c2x2 0 B2 b0 + b1x1 + b2x2 0 d0 + d1x1 + d2x2 ( ) ( ) µ x2/2 x2 0 B3 α + βx2 σρx γ 2 σ 1 ρ 2 x γ 2 ( ) ( a0 + a1x2 B4 b(a x2) + λgx β 1 ρ 2 a + f (x2 a) ρ ) a + f (x2 a) 2 a + f (x2 a) 0 gx β 2 ( ) ( ) bx1 hx1 x2 0 B5 c dx2 gρ x2 g 1 ρ 2 x2 ( ) ( ) m x2/2 x2 0 B6 a bx2 σ 1 ρ 2 x2 σρ 2a > σ 2 x2 ( ) ( ) 0 2x1 B7 γ 2ηx1 x2 γ a1 a2x2 2 x2 0 ( ) ( a + bx1 dx γ ) B8 1 ex2 0 ( cx2 ) ( 0 f a + bx1 dx γ ) B9 1 ex2 0 ( cx2 ) ( 0 f ) B10 b1(a1 x1) g1 0 b2(a2 x2) 0 g2 x2 ( ) ( k1 + k2x2 1 ρ 2 x2 ρ ) x2 B11 κ(θ x2) 0 σx2 ( ) ( ) ax1 cx1e x2 0 B12 bx2 dr d ( ) ( 1 ) r 2 B13 b11(a1 x1) + b12(a2 x2) σ11 σ12 b21(a1( x1) + b22(a2 ) x2) ( σ21 σ22 ) B14 k1(x2 x1) σ x1 0 ( k2(θ x2) ) ( 0 σ2 x2 ) a + bx1 x1 0 B15 ( fx1 + dx2 ) ( h 1 + gx1 ) a + bx1 + gx2 x1 0 B16 d + ηx1 + fx2 h x2
User s Guide for the Matlab Library Implementing Closed Form MLE for Diffusions
User s Guide for the Matlab Library Implementing Closed Form MLE for Diffusions Yacine Aït-Sahalia Department of Economics and Bendheim Center for Finance Princeton University and NBER This Version: July
More informationCounterparty Credit Risk Simulation
Counterparty Credit Risk Simulation Alex Yang FinPricing http://www.finpricing.com Summary Counterparty Credit Risk Definition Counterparty Credit Risk Measures Monte Carlo Simulation Interest Rate Curve
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 informationSensitivity Analysis on Long-term Cash flows
Sensitivity Analysis on Long-term Cash flows Hyungbin Park Worcester Polytechnic Institute 19 March 2016 Eastern Conference on Mathematical Finance Worcester Polytechnic Institute, Worceseter, MA 1 / 49
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 informationApplication of Moment Expansion Method to Option Square Root Model
Application of Moment Expansion Method to Option Square Root Model Yun Zhou Advisor: Professor Steve Heston University of Maryland May 5, 2009 1 / 19 Motivation Black-Scholes Model successfully explain
More informationConvergence Analysis of Monte Carlo Calibration of Financial Market Models
Analysis of Monte Carlo Calibration of Financial Market Models Christoph Käbe Universität Trier Workshop on PDE Constrained Optimization of Certain and Uncertain Processes June 03, 2009 Monte Carlo Calibration
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 informationEstimation of dynamic term structure models
Estimation of dynamic term structure models Greg Duffee Haas School of Business, UC-Berkeley Joint with Richard Stanton, Haas School Presentation at IMA Workshop, May 2004 (full paper at http://faculty.haas.berkeley.edu/duffee)
More informationChapter 3: Black-Scholes Equation and Its Numerical Evaluation
Chapter 3: Black-Scholes Equation and Its Numerical Evaluation 3.1 Itô Integral 3.1.1 Convergence in the Mean and Stieltjes Integral Definition 3.1 (Convergence in the Mean) A sequence {X n } n ln of random
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 informationHeston Model Version 1.0.9
Heston Model Version 1.0.9 1 Introduction This plug-in implements the Heston model. Once installed the plug-in offers the possibility of using two new processes, the Heston process and the Heston time
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 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 informationHeston Stochastic Local Volatility Model
Heston Stochastic Local Volatility Model Klaus Spanderen 1 R/Finance 2016 University of Illinois, Chicago May 20-21, 2016 1 Joint work with Johannes Göttker-Schnetmann Klaus Spanderen Heston Stochastic
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 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 informationStatistical and Computational Inverse Problems with Applications Part 5B: Electrical impedance tomography
Statistical and Computational Inverse Problems with Applications Part 5B: Electrical impedance tomography Aku Seppänen Inverse Problems Group Department of Applied Physics University of Eastern Finland
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 informationTheoretical Problems in Credit Portfolio Modeling 2
Theoretical Problems in Credit Portfolio Modeling 2 David X. Li Shanghai Advanced Institute of Finance (SAIF) Shanghai Jiaotong University(SJTU) November 3, 2017 Presented at the University of South California
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 informationSlides 4. Matthieu Gomez Fall 2017
Slides 4 Matthieu Gomez Fall 2017 How to Compute Stationary Distribution of a Diffusion? Kolmogorov Forward Take a diffusion process dx t = µ(x t )dt + σ(x t )dz t How does the density of x t evolves?
More informationPolynomial Models in Finance
Polynomial Models in Finance Martin Larsson Department of Mathematics, ETH Zürich based on joint work with Damir Filipović, Anders Trolle, Tony Ware Risk Day Zurich, 11 September 2015 Flexibility Tractability
More informationPricing and hedging with rough-heston models
Pricing and hedging with rough-heston models Omar El Euch, Mathieu Rosenbaum Ecole Polytechnique 1 January 216 El Euch, Rosenbaum Pricing and hedging with rough-heston models 1 Table of contents Introduction
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 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 informationRobust Pricing and Hedging of Options on Variance
Robust Pricing and Hedging of Options on Variance Alexander Cox Jiajie Wang University of Bath Bachelier 21, Toronto Financial Setting Option priced on an underlying asset S t Dynamics of S t unspecified,
More informationOrder Making Fiscal Year 2018 Annual Adjustments to Transaction Fee Rates
This document is scheduled to be published in the Federal Register on 04/20/2018 and available online at https://federalregister.gov/d/2018-08339, and on FDsys.gov 8011-01p SECURITIES AND EXCHANGE COMMISSION
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 informationExtended Libor Models and Their Calibration
Extended Libor Models and Their Calibration Denis Belomestny Weierstraß Institute Berlin Haindorf, 7 Februar 2008 Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar
More informationLikelihood Estimation of Jump-Diffusions
Likelihood Estimation of Jump-Diffusions Extensions from Diffusions to Jump-Diffusions, Implementation with Automatic Differentiation, and Applications Berent Ånund Strømnes Lunde DEPARTMENT OF MATHEMATICS
More informationA Closed-form Solution for Outperfomance Options with Stochastic Correlation and Stochastic Volatility
A Closed-form Solution for Outperfomance Options with Stochastic Correlation and Stochastic Volatility Jacinto Marabel Romo Email: jacinto.marabel@grupobbva.com November 2011 Abstract This article introduces
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 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 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 informationFinancial Risk Management
Financial Risk Management Professor: Thierry Roncalli Evry University Assistant: Enareta Kurtbegu Evry University Tutorial exercices #3 1 Maximum likelihood of the exponential distribution 1. We assume
More informationStochastic Dynamical Systems and SDE s. An Informal Introduction
Stochastic Dynamical Systems and SDE s An Informal Introduction Olav Kallenberg Graduate Student Seminar, April 18, 2012 1 / 33 2 / 33 Simple recursion: Deterministic system, discrete time x n+1 = f (x
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 informationIntroduction to Affine Processes. Applications to Mathematical Finance
and Its Applications to Mathematical Finance Department of Mathematical Science, KAIST Workshop for Young Mathematicians in Korea, 2010 Outline Motivation 1 Motivation 2 Preliminary : Stochastic Calculus
More informationNumerical Methods for Pricing Energy Derivatives, including Swing Options, in the Presence of Jumps
Numerical Methods for Pricing Energy Derivatives, including Swing Options, in the Presence of Jumps, Senior Quantitative Analyst Motivation: Swing Options An electricity or gas SUPPLIER needs to be capable,
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 informationEconophysics V: Credit Risk
Fakultät für Physik Econophysics V: Credit Risk Thomas Guhr XXVIII Heidelberg Physics Graduate Days, Heidelberg 2012 Outline Introduction What is credit risk? Structural model and loss distribution Numerical
More informationIntegrating rational functions (Sect. 8.4)
Integrating rational functions (Sect. 8.4) Integrating rational functions, p m(x) q n (x). Polynomial division: p m(x) The method of partial fractions. p (x) (x r )(x r 2 ) p (n )(x). (Repeated roots).
More informationParameter estimation in SDE:s
Lund University Faculty of Engineering Statistics in Finance Centre for Mathematical Sciences, Mathematical Statistics HT 2011 Parameter estimation in SDE:s This computer exercise concerns some estimation
More informationAD in Monte Carlo for finance
AD in Monte Carlo for finance Mike Giles giles@comlab.ox.ac.uk Oxford University Computing Laboratory AD & Monte Carlo p. 1/30 Overview overview of computational finance stochastic o.d.e. s Monte Carlo
More informationTwo hours. To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER
Two hours MATH20802 To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER STATISTICAL METHODS Answer any FOUR of the SIX questions.
More informationEconomics 883: The Basic Diffusive Model, Jumps, Variance Measures. George Tauchen. Economics 883FS Spring 2015
Economics 883: The Basic Diffusive Model, Jumps, Variance Measures George Tauchen Economics 883FS Spring 2015 Main Points 1. The Continuous Time Model, Theory and Simulation 2. Observed Data, Plotting
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 informationMath 416/516: Stochastic Simulation
Math 416/516: Stochastic Simulation Haijun Li lih@math.wsu.edu Department of Mathematics Washington State University Week 13 Haijun Li Math 416/516: Stochastic Simulation Week 13 1 / 28 Outline 1 Simulation
More informationRough volatility models: When population processes become a new tool for trading and risk management
Rough volatility models: When population processes become a new tool for trading and risk management Omar El Euch and Mathieu Rosenbaum École Polytechnique 4 October 2017 Omar El Euch and Mathieu Rosenbaum
More informationCalibration Lecture 1: Background and Parametric Models
Calibration Lecture 1: Background and Parametric Models March 2016 Motivation What is calibration? Derivative pricing models depend on parameters: Black-Scholes σ, interest rate r, Heston reversion speed
More information1 Introduction. 2 Old Methodology BOARD OF GOVERNORS OF THE FEDERAL RESERVE SYSTEM DIVISION OF RESEARCH AND STATISTICS
BOARD OF GOVERNORS OF THE FEDERAL RESERVE SYSTEM DIVISION OF RESEARCH AND STATISTICS Date: October 6, 3 To: From: Distribution Hao Zhou and Matthew Chesnes Subject: VIX Index Becomes Model Free and Based
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 informationEuropean option pricing under parameter uncertainty
European option pricing under parameter uncertainty Martin Jönsson (joint work with Samuel Cohen) University of Oxford Workshop on BSDEs, SPDEs and their Applications July 4, 2017 Introduction 2/29 Introduction
More information25857 Interest Rate Modelling
25857 Interest Rate Modelling UTS Business School University of Technology Sydney Chapter 19. Allowing for Stochastic Interest Rates in the Black-Scholes Model May 15, 2014 1/33 Chapter 19. Allowing for
More informationEnhancing Insurer Value Via Reinsurance Optimization
Enhancing Insurer Value Via Reinsurance Optimization Actuarial Research Symposium 2004 @UNSW Yuriy Krvavych and Michael Sherris University of New South Wales Sydney, AUSTRALIA Actuarial Research Symposium
More informationSample Path Large Deviations and Optimal Importance Sampling for Stochastic Volatility Models
Sample Path Large Deviations and Optimal Importance Sampling for Stochastic Volatility Models Scott Robertson Carnegie Mellon University scottrob@andrew.cmu.edu http://www.math.cmu.edu/users/scottrob June
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 informationVolatility Smiles and Yield Frowns
Volatility Smiles and Yield Frowns Peter Carr NYU CBOE Conference on Derivatives and Volatility, Chicago, Nov. 10, 2017 Peter Carr (NYU) Volatility Smiles and Yield Frowns 11/10/2017 1 / 33 Interest Rates
More informationESTIMATION OF UTILITY FUNCTIONS: MARKET VS. REPRESENTATIVE AGENT THEORY
ESTIMATION OF UTILITY FUNCTIONS: MARKET VS. REPRESENTATIVE AGENT THEORY Kai Detlefsen Wolfgang K. Härdle Rouslan A. Moro, Deutsches Institut für Wirtschaftsforschung (DIW) Center for Applied Statistics
More informationExploring Volatility Derivatives: New Advances in Modelling. Bruno Dupire Bloomberg L.P. NY
Exploring Volatility Derivatives: New Advances in Modelling Bruno Dupire Bloomberg L.P. NY bdupire@bloomberg.net Global Derivatives 2005, Paris May 25, 2005 1. Volatility Products Historical Volatility
More informationLeverage Effect, Volatility Feedback, and Self-Exciting MarketAFA, Disruptions 1/7/ / 14
Leverage Effect, Volatility Feedback, and Self-Exciting Market Disruptions Liuren Wu, Baruch College Joint work with Peter Carr, New York University The American Finance Association meetings January 7,
More informationMulti-dimensional Term Structure Models
Multi-dimensional Term Structure Models We will focus on the affine class. But first some motivation. A generic one-dimensional model for zero-coupon yields, y(t; τ), looks like this dy(t; τ) =... dt +
More informationA New Class of Non-linear Term Structure Models. Discussion
A New Class of Non-linear Term Structure Models by Eraker, Wang and Wu Discussion Pietro Veronesi The University of Chicago Booth School of Business Main Contribution and Outline of Discussion Main contribution
More informationInterest rate models and Solvency II
www.nr.no Outline Desired properties of interest rate models in a Solvency II setting. A review of three well-known interest rate models A real example from a Norwegian insurance company 2 Interest rate
More informationDr. Maddah ENMG 625 Financial Eng g II 10/16/06
Dr. Maddah ENMG 65 Financial Eng g II 10/16/06 Chapter 11 Models of Asset Dynamics () Random Walk A random process, z, is an additive process defined over times t 0, t 1,, t k, t k+1,, such that z( t )
More information"Pricing Exotic Options using Strong Convergence Properties
Fourth Oxford / Princeton Workshop on Financial Mathematics "Pricing Exotic Options using Strong Convergence Properties Klaus E. Schmitz Abe schmitz@maths.ox.ac.uk www.maths.ox.ac.uk/~schmitz Prof. Mike
More informationWeierstrass Institute for Applied Analysis and Stochastics Maximum likelihood estimation for jump diffusions
Weierstrass Institute for Applied Analysis and Stochastics Maximum likelihood estimation for jump diffusions Hilmar Mai Mohrenstrasse 39 1117 Berlin Germany Tel. +49 3 2372 www.wias-berlin.de Haindorf
More informationQuadratic hedging in affine stochastic volatility models
Quadratic hedging in affine stochastic volatility models Jan Kallsen TU München Pittsburgh, February 20, 2006 (based on joint work with F. Hubalek, L. Krawczyk, A. Pauwels) 1 Hedging problem S t = S 0
More informationThe Lognormal Interest Rate Model and Eurodollar Futures
GLOBAL RESEARCH ANALYTICS The Lognormal Interest Rate Model and Eurodollar Futures CITICORP SECURITIES,INC. 399 Park Avenue New York, NY 143 Keith Weintraub Director, Analytics 1-559-97 Michael Hogan Ex
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 informationLecture 7: Computation of Greeks
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
More informationImportance sampling and Monte Carlo-based calibration for time-changed Lévy processes
Importance sampling and Monte Carlo-based calibration for time-changed Lévy processes Stefan Kassberger Thomas Liebmann BFS 2010 1 Motivation 2 Time-changed Lévy-models and Esscher transforms 3 Applications
More informationSTOR Lecture 15. Jointly distributed Random Variables - III
STOR 435.001 Lecture 15 Jointly distributed Random Variables - III Jan Hannig UNC Chapel Hill 1 / 17 Before we dive in Contents of this lecture 1. Conditional pmf/pdf: definition and simple properties.
More informationVolatility Smiles and Yield Frowns
Volatility Smiles and Yield Frowns Peter Carr NYU IFS, Chengdu, China, July 30, 2018 Peter Carr (NYU) Volatility Smiles and Yield Frowns 7/30/2018 1 / 35 Interest Rates and Volatility Practitioners and
More informationGeneralized Multi-Factor Commodity Spot Price Modeling through Dynamic Cournot Resource Extraction Models
Generalized Multi-Factor Commodity Spot Price Modeling through Dynamic Cournot Resource Extraction Models Bilkan Erkmen (joint work with Michael Coulon) Workshop on Stochastic Games, Equilibrium, and Applications
More informationManaging Systematic Mortality Risk in Life Annuities: An Application of Longevity Derivatives
Managing Systematic Mortality Risk in Life Annuities: An Application of Longevity Derivatives Simon Man Chung Fung, Katja Ignatieva and Michael Sherris School of Risk & Actuarial Studies University of
More informationEstimation of Stochastic Volatility Models with Implied. Volatility Indices and Pricing of Straddle Option
Estimation of Stochastic Volatility Models with Implied Volatility Indices and Pricing of Straddle Option Yue Peng Steven C. J. Simon June 14, 29 Abstract Recent market turmoil has made it clear that modelling
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 informationTesting the significance of the RV coefficient
1 / 19 Testing the significance of the RV coefficient Application to napping data Julie Josse, François Husson and Jérôme Pagès Applied Mathematics Department Agrocampus Rennes, IRMAR CNRS UMR 6625 Agrostat
More informationInternet Appendix to Idiosyncratic Cash Flows and Systematic Risk
Internet Appendix to Idiosyncratic Cash Flows and Systematic Risk ILONA BABENKO, OLIVER BOGUTH, and YURI TSERLUKEVICH This Internet Appendix supplements the analysis in the main text by extending the model
More informationRisk, Return, and Ross Recovery
Risk, Return, and Ross Recovery Peter Carr and Jiming Yu Courant Institute, New York University September 13, 2012 Carr/Yu (NYU Courant) Risk, Return, and Ross Recovery September 13, 2012 1 / 30 P, Q,
More informationNBER WORKING PAPER SERIES MAXIMUM LIKELIHOOD ESTIMATION OF STOCHASTIC VOLATILITY MODELS. Yacine Ait-Sahalia Robert Kimmel
NBER WORKING PAPER SERIES MAIMUM LIKELIHOOD ESTIMATION OF STOCHASTIC VOLATILITY MODELS Yacine Ait-Sahalia Robert Kimmel Working Paper 10579 http://www.nber.org/papers/w10579 NATIONAL BUREAU OF ECONOMIC
More informationModelling Credit Spread Behaviour. FIRST Credit, Insurance and Risk. Angelo Arvanitis, Jon Gregory, Jean-Paul Laurent
Modelling Credit Spread Behaviour Insurance and Angelo Arvanitis, Jon Gregory, Jean-Paul Laurent ICBI Counterparty & Default Forum 29 September 1999, Paris Overview Part I Need for Credit Models Part II
More informationStochastic Processes and Brownian Motion
A stochastic process Stochastic Processes X = { X(t) } Stochastic Processes and Brownian Motion is a time series of random variables. X(t) (or X t ) is a random variable for each time t and is usually
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 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 informationCredit Risk using Time Changed Brownian Motions
Credit Risk using Time Changed Brownian Motions Tom Hurd Mathematics and Statistics McMaster University Joint work with Alexey Kuznetsov (New Brunswick) and Zhuowei Zhou (Mac) 2nd Princeton Credit Conference
More informationLattice (Binomial Trees) Version 1.2
Lattice (Binomial Trees) Version 1. 1 Introduction This plug-in implements different binomial trees approximations for pricing contingent claims and allows Fairmat to use some of the most popular binomial
More informationInterest Rate Curves Calibration with Monte-Carlo Simulatio
Interest Rate Curves Calibration with Monte-Carlo Simulation 24 june 2008 Participants A. Baena (UCM) Y. Borhani (Univ. of Oxford) E. Leoncini (Univ. of Florence) R. Minguez (UCM) J.M. Nkhaso (UCM) A.
More information(A note) on co-integration in commodity markets
(A note) on co-integration in commodity markets Fred Espen Benth Centre of Mathematics for Applications (CMA) University of Oslo, Norway In collaboration with Steen Koekebakker (Agder) Energy & Finance
More informationA new approach for scenario generation in risk management
A new approach for scenario generation in risk management Josef Teichmann TU Wien Vienna, March 2009 Scenario generators Scenarios of risk factors are needed for the daily risk analysis (1D and 10D ahead)
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 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 informationAll Investors are Risk-averse Expected Utility Maximizers. Carole Bernard (UW), Jit Seng Chen (GGY) and Steven Vanduffel (Vrije Universiteit Brussel)
All Investors are Risk-averse Expected Utility Maximizers Carole Bernard (UW), Jit Seng Chen (GGY) and Steven Vanduffel (Vrije Universiteit Brussel) First Name: Waterloo, April 2013. Last Name: UW ID #:
More informationSmart Beta: Managing Diversification of Minimum Variance Portfolios
Smart Beta: Managing Diversification of Minimum Variance Portfolios Jean-Charles Richard and Thierry Roncalli Lyxor Asset Management 1, France University of Évry, France Risk Based and Factor Investing
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 informationParameter estimation of diffusion models from discrete observations
221 Parameter estimation of diffusion models from discrete observations Miljenko Huzak Abstract. A short review of diffusion parameter estimations methods from discrete observations is presented. The applicability
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 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 information