Importance sampling and Monte Carlo-based calibration for time-changed Lévy processes
|
|
- Zoe Stafford
- 5 years ago
- Views:
Transcription
1 Importance sampling and Monte Carlo-based calibration for time-changed Lévy processes Stefan Kassberger Thomas Liebmann BFS 2010
2 1 Motivation 2 Time-changed Lévy-models and Esscher transforms 3 Applications and examples Importance sampling Sensitivities w.r.t. Esscher transform parameters Simulation of a Normal Tempered Stable process Calibration Example 4 Summary
3 Motivation Variety of tractable Lévy-models can be represented as time-changed Brownian motion Esscher-transform well-established for Lévy-models Kassberger and Liebmann (2009) apply independent Esscher transforms to Brownian motion and subordinator in the TCL-context Exploit above idea in the context of Monte Carlo simulation: Variance reduction through importance sampling Calculating sensitivities by likelihood-ratio methods Extending sampling algorithms to more general classes of distributions Monte-Carlo based model calibration
4 Setup n-dimensional Brownian motion (B t ), subordinator (T (t)). Fix τ and denote by κ(u) = loge Q [exp(ut (τ))] the cumulant generating function of T (τ). Define equivalent measure S η,γ via ds η,γ dq = exp ( η B T (τ) 1 ) 2 η 2 T (τ) + γt (τ) κ(γ). This transform is composed of two Esscher transforms: 1 One with parameter η R n on the Brownian motion B. Shifts the drift of B from 0 under Q to η under S η,γ. 2 One with parameter γ on the subordinator. Its cgf under S η,γ is given by κ S η,γ (u) = κ(u + γ) κ(γ).
5 Outline 1 Motivation 2 Time-changed Lévy-models and Esscher transforms 3 Applications and examples Importance sampling Sensitivities w.r.t. Esscher transform parameters Simulation of a Normal Tempered Stable process Calibration Example 4 Summary
6 Importance sampling Set X(t) = B T (t) + mt (t). Take expectations under the transformed measure: [ E Q f((x(t))t [0,τ] ) ] ( = E Sη,γ [exp η B T (τ) + 1 ) ] 2 η 2 T (τ) γt (τ) + κ(γ) f((x(t)) t [0,τ] ). If the subordinator cannot be efficiently sampled under Q but under S η,γ for certain values of γ, we can choose η for variance reduction via importance sampling (compare example for Normal Tempered Stable process).
7 Importance sampling Goal: Estimate ( E Sη,γ [exp η B T (τ) + 1 ) ] 2 η 2 T (τ) γt (τ) + κ(γ) f((x(t)) t [0,τ] ) (1) by simulation. Simulate appropriately discretized version of (X(t)) : (X(t i )) i=0,...,n with 0 = t 0 < < t N = τ. Proceed as follows: Simulate (T (t i )) i=0,...,n and iid N n (0,diag(1 n )) rvs W 1,...,W N. Set Y (t i ) = Y (t i 1 ) + W i T (ti ) T (t i 1 ) + η(t i t i 1 ). Set X(t i ) = Y (t i ) + mt (t i ). Repeat M times to arrive at set of sample paths indexed by k. Estimator for (1) is ( M 1 exp η Y k (τ) + 1 ) M k=1 2 η 2 T k (τ) γt k (τ) + κ(γ) f((x k (t i )) i=0,...,n ).
8 Importance sampling: Variance of the Estimator Let f(x) ce β x for some c 0 and β R n. Bound for variance of the estimator [ ( M 1 Var Sη,γ exp η Y k (τ) + 1 ) ] M k=1 2 η 2 T k (τ) γt k (τ) + κ(γ) f(x k (τ)) c2 M exp( κ( β 2 + β η 2 γ + 2β m) + κ(γ) ) κ is increasing and convex. For given β, c and γ, the bound is minimal if η = β. If γ can be chosen freely, minimum is attained for γ = β 2 /2 + β m. If further m = 0, this holds for γ = β 2 /2 = η 2 /2, i.e. an Esscher transform on B T (τ) with parameter β.
9 Outline 1 Motivation 2 Time-changed Lévy-models and Esscher transforms 3 Applications and examples Importance sampling Sensitivities w.r.t. Esscher transform parameters Simulation of a Normal Tempered Stable process Calibration Example 4 Summary
10 Sensitivities with the LR-method Sensitivities w.r.t. Esscher transform parameters η and γ can be estimated by the likelihood ratio method If the expectations are finite in some neighborhood of γ 0, d dγ E [ S η0,γ f((x(t))t [0,τ] ) ] γ=γ0 = E Sη,γ [ dsη0,γ 0 dq dq ( T (τ) κ (γ 0 ) ) ] f((x(t)) t [0,τ] ) ds η,γ If the expectations are finite in some neighborhood of η 0, d dη E S η,γ0 [ f((x(t))t [0,τ] ) ] η=η0 = E Sη,γ [ dsη0,γ 0 dq dq ds η,γ ( BT (τ) η 0 T (τ) ) f((x(t)) t [0,τ] ) ]
11 Outline 1 Motivation 2 Time-changed Lévy-models and Esscher transforms 3 Applications and examples Importance sampling Sensitivities w.r.t. Esscher transform parameters Simulation of a Normal Tempered Stable process Calibration Example 4 Summary
12 Simulation of a Normal Tempered Stable process Let T be a Tempered Stable TS(κ,a,(2γ) κ ) random variable with κ (0,1), a,γ > 0 and cumulant generating function κ κ,a,γ (u) = a(2γ) κ a(2γ 2u) κ. Algorithm that allows direct sampling is unknown Idea: Simulate stable random variables and perform Esscher transform Algorithm: 1 Sample n iid stable TS(κ,a,0) rvs T k 2 Set importance weights to w k = exp( γt k + κ κ,a,γ (γ)) = exp( γt k + a(2γ) κ ) In Normal Tempered Stable model, η still available for IS
13 Outline 1 Motivation 2 Time-changed Lévy-models and Esscher transforms 3 Applications and examples Importance sampling Sensitivities w.r.t. Esscher transform parameters Simulation of a Normal Tempered Stable process Calibration Example 4 Summary
14 MC-based model calibration MC-based model calibration time-consuming and unstable Change of parameters often requires re-simulation Noise due to finite sample size often precludes gradient-based methods and slows down convergence of optimizer Idea: Use of two Esscher transforms often allows calibration based on a single set of paths or at least helps significantly reduce number of simulation runs Simulate BM and subordinator, keep paths in memory, and then apply Esscher transforms in conjunction with basic transformation (shifting & scaling) Facilitates use of gradient-based optimization algorithms (more exact & stable gradients)
15 Outline 1 Motivation 2 Time-changed Lévy-models and Esscher transforms 3 Applications and examples Importance sampling Sensitivities w.r.t. Esscher transform parameters Simulation of a Normal Tempered Stable process Calibration Example 4 Summary
16 Example: MC-calibration of NIG model Objective: Calibrate model based on NIG(α,β,δ,µ) paths, i.e. find α,β,δ,µ such that observed derivatives prices are replicated Proceed as follows (easily generalizes to appropriate time-scaling): 1 Fix c. Simulate IG(1,c) paths (IG t ) and standard Brownian paths (B t ) 2 Fix α,β,δ,µ 3 Calculate b = δ α 2 β 2 and choose γ such that c 2 2γ = b. 4 Calculate Esscher weights corresponding to γ. 5 Under S 0,γ, X = βδ 2 IG 1 + δ B IG1 + µ is NIG(α,β,δ,µ)-distributed. 6 Calculate option prices / evaluate objective function under S 0,γ, 7 Repeat steps 2-7 if necessary η (unused above) can be employed for importance sampling.
17 Example: Pricing a binary down-and-out call via MC Let the log-return process (X t ) be a Lévy-process with X 1 NIG(α,0,δ,µ) with µ = µ(α,δ) chosen such that (S 0 exp(x t )) is a martingale (assume riskless interest rate r = 0 and S 0 = 100). Price a binary down-and-out call (BDOC) in this model via Esscher-based MC. Payoff of BDOC maturing in 1 year: I(min(S t,t [0,1]) 80) I(S 1 100) Price BDOC via MC: Discretize using 250 time-steps. Simulate Brownian and IG paths. Use Esscher-based MC: Smooth dependence of BDOC-price on parameters (α, δ). Fast calculation of option price surface as only one simulation run is needed.
18 BDOC price surface BDOC price δ α Figure: BDOC price as function of α and δ
19 Example: Calculating sensitivities w.r.t. model parameters Sensitivities based on plain MC notoriously noisy Esscher-based MC makes sensitivity estimates more stable and provides faster convergence In above setup, calculate sensitivity of BDOC-price w.r.t. model parameter δ for α 0 = 2 and δ 0 = 1 First method: Finite differences with resampling (same seed for random number generator in both runs) Second method: Finite differences with Esscher (only one sample plus Esscher transform)
20 Convergence of parameter sensitivities Finite differences (Esscher) Finite differences (resampling) 0 1.0E+3 2.0E+3 4.1E+3 8.2E+3 1.6E+4 3.3E+4 6.6E+4 1.3E+5 2.6E+5 5.2E+5 1.0E Figure: Sensitivities w.r.t. δ as function of sample size
21 Summary Apply independent Esscher transforms to subordinator and BM Gives higher degree of flexibility than Esscher transform applied to TCBM itself Subclass of structure preserving transforms (Kassberger and Liebmann (2009)) Variety of applications: 1 IS with upper bounds for variance 2 IS with a mixture of importance distributions 3 sensitivities w.r.t. Esscher parameters via LR-methods 4 development of sampling algorithms 5 model calibration Approach also works for subordinated stable processes
22 Thanks for your attention!
Equity 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 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 informationRisk Neutral Valuation
copyright 2012 Christian Fries 1 / 51 Risk Neutral Valuation Christian Fries Version 2.2 http://www.christian-fries.de/finmath April 19-20, 2012 copyright 2012 Christian Fries 2 / 51 Outline Notation Differential
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 informationEffectiveness of CPPI Strategies under Discrete Time Trading
Effectiveness of CPPI Strategies under Discrete Time Trading S. Balder, M. Brandl 1, Antje Mahayni 2 1 Department of Banking and Finance, University of Bonn 2 Department of Accounting and Finance, Mercator
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 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 informationCS 774 Project: Fall 2009 Version: November 27, 2009
CS 774 Project: Fall 2009 Version: November 27, 2009 Instructors: Peter Forsyth, paforsyt@uwaterloo.ca Office Hours: Tues: 4:00-5:00; Thurs: 11:00-12:00 Lectures:MWF 3:30-4:20 MC2036 Office: DC3631 CS
More informationNormal Inverse Gaussian (NIG) Process
With Applications in Mathematical Finance The Mathematical and Computational Finance Laboratory - Lunch at the Lab March 26, 2009 1 Limitations of Gaussian Driven Processes Background and Definition IG
More informationLimit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies
Limit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies George Tauchen Duke University Viktor Todorov Northwestern University 2013 Motivation
More informationValue at Risk Ch.12. PAK Study Manual
Value at Risk Ch.12 Related Learning Objectives 3a) Apply and construct risk metrics to quantify major types of risk exposure such as market risk, credit risk, liquidity risk, regulatory risk etc., and
More informationOption Pricing for Discrete Hedging and Non-Gaussian Processes
Option Pricing for Discrete Hedging and Non-Gaussian Processes Kellogg College University of Oxford A thesis submitted in partial fulfillment of the requirements for the MSc in Mathematical Finance November
More informationPricing of some exotic options with N IG-Lévy input
Pricing of some exotic options with N IG-Lévy input Sebastian Rasmus, Søren Asmussen 2 and Magnus Wiktorsson Center for Mathematical Sciences, University of Lund, Box 8, 22 00 Lund, Sweden {rasmus,magnusw}@maths.lth.se
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 informationMINIMAL PARTIAL PROXY SIMULATION SCHEMES FOR GENERIC AND ROBUST MONTE-CARLO GREEKS
MINIMAL PARTIAL PROXY SIMULATION SCHEMES FOR GENERIC AND ROBUST MONTE-CARLO GREEKS JIUN HONG CHAN AND MARK JOSHI Abstract. In this paper, we present a generic framework known as the minimal partial proxy
More informationMLEMVD: A R Package for Maximum Likelihood Estimation of Multivariate Diffusion Models
MLEMVD: A R Package for Maximum Likelihood Estimation of Multivariate Diffusion Models Matthew Dixon and Tao Wu 1 Illinois Institute of Technology May 19th 2017 1 https://papers.ssrn.com/sol3/papers.cfm?abstract
More informationOptimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models
Optimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models José E. Figueroa-López 1 1 Department of Statistics Purdue University University of Missouri-Kansas City Department of Mathematics
More informationNumerical Methods in Option Pricing (Part III)
Numerical Methods in Option Pricing (Part III) E. Explicit Finite Differences. Use of the Forward, Central, and Symmetric Central a. In order to obtain an explicit solution for the price of the derivative,
More informationOptimal Trade Execution: Mean Variance or Mean Quadratic Variation?
Optimal Trade Execution: Mean Variance or Mean Quadratic Variation? Peter Forsyth 1 S. Tse 2 H. Windcliff 2 S. Kennedy 2 1 Cheriton School of Computer Science University of Waterloo 2 Morgan Stanley New
More informationForward Monte-Carlo Scheme for PDEs: Multi-Type Marked Branching Diffusions
Forward Monte-Carlo Scheme for PDEs: Multi-Type Marked Branching Diffusions Pierre Henry-Labordère 1 1 Global markets Quantitative Research, SOCIÉTÉ GÉNÉRALE Outline 1 Introduction 2 Semi-linear PDEs 3
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 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 informationStochastic modelling of electricity markets Pricing Forwards and Swaps
Stochastic modelling of electricity markets Pricing Forwards and Swaps Jhonny Gonzalez School of Mathematics The University of Manchester Magical books project August 23, 2012 Clip for this slide Pricing
More informationOption Pricing and Calibration with Time-changed Lévy processes
Option Pricing and Calibration with Time-changed Lévy processes Yan Wang and Kevin Zhang Warwick Business School 12th Feb. 2013 Objectives 1. How to find a perfect model that captures essential features
More informationUnified Credit-Equity Modeling
Unified Credit-Equity Modeling Rafael Mendoza-Arriaga Based on joint research with: Vadim Linetsky and Peter Carr The University of Texas at Austin McCombs School of Business (IROM) Recent Advancements
More 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 informationA new approach to LIBOR modeling
A new approach to LIBOR modeling Antonis Papapantoleon FAM TU Vienna Based on joint work with Martin Keller-Ressel and Josef Teichmann Istanbul Workshop on Mathematical Finance Istanbul, Turkey, 18 May
More informationGamma. The finite-difference formula for gamma is
Gamma The finite-difference formula for gamma is [ P (S + ɛ) 2 P (S) + P (S ɛ) e rτ E ɛ 2 ]. For a correlation option with multiple underlying assets, the finite-difference formula for the cross gammas
More information- 1 - **** d(lns) = (µ (1/2)σ 2 )dt + σdw t
- 1 - **** These answers indicate the solutions to the 2014 exam questions. Obviously you should plot graphs where I have simply described the key features. It is important when plotting graphs to label
More information1 Rare event simulation and importance sampling
Copyright c 2007 by Karl Sigman 1 Rare event simulation and importance sampling Suppose we wish to use Monte Carlo simulation to estimate a probability p = P (A) when the event A is rare (e.g., when p
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 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 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 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 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 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 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 informationA Numerical Approach to the Estimation of Search Effort in a Search for a Moving Object
Proceedings of the 1. Conference on Applied Mathematics and Computation Dubrovnik, Croatia, September 13 18, 1999 pp. 129 136 A Numerical Approach to the Estimation of Search Effort in a Search for a Moving
More informationOptimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing
Optimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing Prof. Chuan-Ju Wang Department of Computer Science University of Taipei Joint work with Prof. Ming-Yang Kao March 28, 2014
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Other Miscellaneous Topics and Applications of Monte-Carlo Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationDRAFT. 1 exercise in state (S, t), π(s, t) = 0 do not exercise in state (S, t) Review of the Risk Neutral Stock Dynamics
Chapter 12 American Put Option Recall that the American option has strike K and maturity T and gives the holder the right to exercise at any time in [0, T ]. The American option is not straightforward
More informationShort-time asymptotics for ATM option prices under tempered stable processes
Short-time asymptotics for ATM option prices under tempered stable processes José E. Figueroa-López 1 1 Department of Statistics Purdue University Probability Seminar Purdue University Oct. 30, 2012 Joint
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 informationOn modelling of electricity spot price
, Rüdiger Kiesel and Fred Espen Benth Institute of Energy Trading and Financial Services University of Duisburg-Essen Centre of Mathematics for Applications, University of Oslo 25. August 2010 Introduction
More informationStrategies for Improving the Efficiency of Monte-Carlo Methods
Strategies for Improving the Efficiency of Monte-Carlo Methods Paul J. Atzberger General comments or corrections should be sent to: paulatz@cims.nyu.edu Introduction The Monte-Carlo method is a useful
More informationA start of Variational Methods for ERGM Ranran Wang, UW
A start of Variational Methods for ERGM Ranran Wang, UW MURI-UCI April 24, 2009 Outline A start of Variational Methods for ERGM [1] Introduction to ERGM Current methods of parameter estimation: MCMCMLE:
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 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 informationinduced by the Solvency II project
Asset Les normes allocation IFRS : new en constraints assurance induced by the Solvency II project 36 th International ASTIN Colloquium Zürich September 005 Frédéric PLANCHET Pierre THÉROND ISFA Université
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 informationHedging Under Jump Diffusions with Transaction Costs. Peter Forsyth, Shannon Kennedy, Ken Vetzal University of Waterloo
Hedging Under Jump Diffusions with Transaction Costs Peter Forsyth, Shannon Kennedy, Ken Vetzal University of Waterloo Computational Finance Workshop, Shanghai, July 4, 2008 Overview Overview Single factor
More informationA Macro-Finance Model of the Term Structure: the Case for a Quadratic Yield Model
Title page Outline A Macro-Finance Model of the Term Structure: the Case for a 21, June Czech National Bank Structure of the presentation Title page Outline Structure of the presentation: Model Formulation
More informationMultilevel quasi-monte Carlo path simulation
Multilevel quasi-monte Carlo path simulation Michael B. Giles and Ben J. Waterhouse Lluís Antoni Jiménez Rugama January 22, 2014 Index 1 Introduction to MLMC Stochastic model Multilevel Monte Carlo Milstein
More informationCPSC 540: Machine Learning
CPSC 540: Machine Learning Monte Carlo Methods Mark Schmidt University of British Columbia Winter 2019 Last Time: Markov Chains We can use Markov chains for density estimation, d p(x) = p(x 1 ) p(x }{{}
More informationOptimum Thresholding for Semimartingales with Lévy Jumps under the mean-square error
Optimum Thresholding for Semimartingales with Lévy Jumps under the mean-square error José E. Figueroa-López Department of Mathematics Washington University in St. Louis Spring Central Sectional Meeting
More informationVolatility. Roberto Renò. 2 March 2010 / Scuola Normale Superiore. Dipartimento di Economia Politica Università di Siena
Dipartimento di Economia Politica Università di Siena 2 March 2010 / Scuola Normale Superiore What is? The definition of volatility may vary wildly around the idea of the standard deviation of price movements
More informationExtend the ideas of Kan and Zhou paper on Optimal Portfolio Construction under parameter uncertainty
Extend the ideas of Kan and Zhou paper on Optimal Portfolio Construction under parameter uncertainty George Photiou Lincoln College University of Oxford A dissertation submitted in partial fulfilment for
More informationChapter 15: Jump Processes and Incomplete Markets. 1 Jumps as One Explanation of Incomplete Markets
Chapter 5: Jump Processes and Incomplete Markets Jumps as One Explanation of Incomplete Markets It is easy to argue that Brownian motion paths cannot model actual stock price movements properly in reality,
More informationAlgorithmic and High-Frequency Trading
LOBSTER June 2 nd 2016 Algorithmic and High-Frequency Trading Julia Schmidt Overview Introduction Market Making Grossman-Miller Market Making Model Trading Costs Measuring Liquidity Market Making using
More informationMarket interest-rate models
Market interest-rate models Marco Marchioro www.marchioro.org November 24 th, 2012 Market interest-rate models 1 Lecture Summary No-arbitrage models Detailed example: Hull-White Monte Carlo simulations
More informationIntroduction Credit risk
A structural credit risk model with a reduced-form default trigger Applications to finance and insurance Mathieu Boudreault, M.Sc.,., F.S.A. Ph.D. Candidate, HEC Montréal Montréal, Québec Introduction
More informationConditional Density Method in the Computation of the Delta with Application to Power Market
Conditional Density Method in the Computation of the Delta with Application to Power Market Asma Khedher Centre of Mathematics for Applications Department of Mathematics University of Oslo A joint work
More informationLecture 17: More on Markov Decision Processes. Reinforcement learning
Lecture 17: More on Markov Decision Processes. Reinforcement learning Learning a model: maximum likelihood Learning a value function directly Monte Carlo Temporal-difference (TD) learning COMP-424, Lecture
More informationFair Valuation of Insurance Contracts under Lévy Process Specifications Preliminary Version
Fair Valuation of Insurance Contracts under Lévy Process Specifications Preliminary Version Rüdiger Kiesel, Thomas Liebmann, Stefan Kassberger University of Ulm and LSE June 8, 2005 Abstract The valuation
More informationAsymptotic Methods in Financial Mathematics
Asymptotic Methods in Financial Mathematics José E. Figueroa-López 1 1 Department of Mathematics Washington University in St. Louis Statistics Seminar Washington University in St. Louis February 17, 2017
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 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 informationOptimized Least-squares Monte Carlo (OLSM) for Measuring Counterparty Credit Exposure of American-style Options
Optimized Least-squares Monte Carlo (OLSM) for Measuring Counterparty Credit Exposure of American-style Options Kin Hung (Felix) Kan 1 Greg Frank 3 Victor Mozgin 3 Mark Reesor 2 1 Department of Applied
More informationMonte Carlo Simulations
Monte Carlo Simulations Lecture 1 December 7, 2014 Outline Monte Carlo Methods Monte Carlo methods simulate the random behavior underlying the financial models Remember: When pricing you must simulate
More informationOn CAT Options and Bonds
Advanced Modeling in Finance and Insurance Linz, 24th of September 2008 1 Introduction Why are Catastrophes Dangerous? CAT Products CAT Products from the Investor s Point of View 2 Models Models for the
More informationNear-expiration behavior of implied volatility for exponential Lévy models
Near-expiration behavior of implied volatility for exponential Lévy models José E. Figueroa-López 1 1 Department of Statistics Purdue University Financial Mathematics Seminar The Stevanovich Center for
More informationDYNAMIC CDO TERM STRUCTURE MODELLING
DYNAMIC CDO TERM STRUCTURE MODELLING Damir Filipović (joint with Ludger Overbeck and Thorsten Schmidt) Vienna Institute of Finance www.vif.ac.at PRisMa 2008 Workshop on Portfolio Risk Management TU Vienna,
More informationComputational Finance. Computational Finance p. 1
Computational Finance Computational Finance p. 1 Outline Binomial model: option pricing and optimal investment Monte Carlo techniques for pricing of options pricing of non-standard options improving accuracy
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 informationOptimal robust bounds for variance options and asymptotically extreme models
Optimal robust bounds for variance options and asymptotically extreme models Alexander Cox 1 Jiajie Wang 2 1 University of Bath 2 Università di Roma La Sapienza Advances in Financial Mathematics, 9th January,
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 informationFast Convergence of Regress-later Series Estimators
Fast Convergence of Regress-later Series Estimators New Thinking in Finance, London Eric Beutner, Antoon Pelsser, Janina Schweizer Maastricht University & Kleynen Consultants 12 February 2014 Beutner Pelsser
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 informationELEMENTS OF MONTE CARLO SIMULATION
APPENDIX B ELEMENTS OF MONTE CARLO SIMULATION B. GENERAL CONCEPT The basic idea of Monte Carlo simulation is to create a series of experimental samples using a random number sequence. According to the
More informationParametric Inference and Dynamic State Recovery from Option Panels. Torben G. Andersen
Parametric Inference and Dynamic State Recovery from Option Panels Torben G. Andersen Joint work with Nicola Fusari and Viktor Todorov The Third International Conference High-Frequency Data Analysis in
More informationArbitrage Bounds for Volatility Derivatives as Free Boundary Problem. Bruno Dupire Bloomberg L.P. NY
Arbitrage Bounds for Volatility Derivatives as Free Boundary Problem Bruno Dupire Bloomberg L.P. NY bdupire@bloomberg.net PDE and Mathematical Finance, KTH, Stockholm August 16, 25 Variance Swaps Vanilla
More informationHow to hedge Asian options in fractional Black-Scholes model
How to hedge Asian options in fractional Black-Scholes model Heikki ikanmäki Jena, March 29, 211 Fractional Lévy processes 1/36 Outline of the talk 1. Introduction 2. Main results 3. Methodology 4. Conclusions
More informationRapid computation of prices and deltas of nth to default swaps in the Li Model
Rapid computation of prices and deltas of nth to default swaps in the Li Model Mark Joshi, Dherminder Kainth QUARC RBS Group Risk Management Summary Basic description of an nth to default swap Introduction
More informationMonte Carlo Methods in Structuring and Derivatives Pricing
Monte Carlo Methods in Structuring and Derivatives Pricing Prof. Manuela Pedio (guest) 20263 Advanced Tools for Risk Management and Pricing Spring 2017 Outline and objectives The basic Monte Carlo algorithm
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 informationConstruction and behavior of Multinomial Markov random field models
Graduate Theses and Dissertations Iowa State University Capstones, Theses and Dissertations 2010 Construction and behavior of Multinomial Markov random field models Kim Mueller Iowa State University Follow
More informationFinancial Engineering and Structured Products
550.448 Financial Engineering and Structured Products Week of March 31, 014 Structured Securitization Liability-Side Cash Flow Analysis & Structured ransactions Assignment Reading (this week, March 31
More informationRisk Measurement in Credit Portfolio Models
9 th DGVFM Scientific Day 30 April 2010 1 Risk Measurement in Credit Portfolio Models 9 th DGVFM Scientific Day 30 April 2010 9 th DGVFM Scientific Day 30 April 2010 2 Quantitative Risk Management Profit
More informationOption pricing in the stochastic volatility model of Barndorff-Nielsen and Shephard
Option pricing in the stochastic volatility model of Barndorff-Nielsen and Shephard Indifference pricing and the minimal entropy martingale measure Fred Espen Benth Centre of Mathematics for Applications
More informationEstimating Mixed Logit Models with Large Choice Sets. Roger H. von Haefen, NC State & NBER Adam Domanski, NOAA July 2013
Estimating Mixed Logit Models with Large Choice Sets Roger H. von Haefen, NC State & NBER Adam Domanski, NOAA July 2013 Motivation Bayer et al. (JPE, 2007) Sorting modeling / housing choice 250,000 individuals
More informationLarge Deviations and Stochastic Volatility with Jumps: Asymptotic Implied Volatility for Affine Models
Large Deviations and Stochastic Volatility with Jumps: TU Berlin with A. Jaquier and A. Mijatović (Imperial College London) SIAM conference on Financial Mathematics, Minneapolis, MN July 10, 2012 Implied
More 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 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 informationA Robust Option Pricing Problem
IMA 2003 Workshop, March 12-19, 2003 A Robust Option Pricing Problem Laurent El Ghaoui Department of EECS, UC Berkeley 3 Robust optimization standard form: min x sup u U f 0 (x, u) : u U, f i (x, u) 0,
More informationRegression estimation in continuous time with a view towards pricing Bermudan options
with a view towards pricing Bermudan options Tagung des SFB 649 Ökonomisches Risiko in Motzen 04.-06.06.2009 Financial engineering in times of financial crisis Derivate... süßes Gift für die Spekulanten
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 informationShort & Long Run impact of volatility on the effect monetary shocks
Short & Long Run impact of volatility on the effect monetary shocks Fernando Alvarez University of Chicago & NBER Inflation: Drivers & Dynamics Conference 218 Cleveland Fed Alvarez Volatility & Monetary
More informationCPSC 540: Machine Learning
CPSC 540: Machine Learning Monte Carlo Methods Mark Schmidt University of British Columbia Winter 2018 Last Time: Markov Chains We can use Markov chains for density estimation, p(x) = p(x 1 ) }{{} d p(x
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 informationROM SIMULATION Exact Moment Simulation using Random Orthogonal Matrices
ROM SIMULATION Exact Moment Simulation using Random Orthogonal Matrices Bachelier Finance Society Meeting Toronto 2010 Henley Business School at Reading Contact Author : d.ledermann@icmacentre.ac.uk Alexander
More information