Parameter estimation in SDE:s
|
|
- Brice Burke
- 5 years ago
- Views:
Transcription
1 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 of parameters in stochastic differential equations. It contains two main assignments with several sub-assignments. Many of these are labelled Extra. You should solve all the regular problems (you need to solve at least two extra problems if you want to get 2 bonus credits from the computer exercise.) It is preferred if you use Matlab, but you are allowed to use the programming language or package of your choice. If you choose not to use Matlab, please note that you are required to document your code extra carefully. 1 Preparations for the exercise Read chapters 8.1, 11 and 14 in Madsen et. al. [2006 and this instruction. Then you should prepare the computer exercise by writing down the Matlab functions needed for the exercise. Before the computer exercise some of the questions below will be posed. All of the posed question must be answered correctly in order to pass the computer exercise. 2 Catalogue of questions You should be able to answer the following questions before the computer exercise. 1. Suggest suitable moment restrictions for estimating the model using GMM. ε i are i.i.d. normal. y i = β 0 + β 1 x i + ε i, (1) 2. Discretize and find suitable moment restrictions for the CKLS-model. 3. Write down the approximate likelihood function for the CIR model when discretized using the Euler scheme. 4. (If you do exercise 3.4.2) Calculate first, second and third order moments for the Cox-Ingersoll-Ross process using Dynkin s formula (see equation 4). 1
2 3 Computer Exercises 3.1 Generalized Method of Moments Estimation with the Generalized Method of Moments (GMM) is an alternative that should be used mainly when maximum likelihood estimation isn t possible, or when we want to aviod explicit assumptions about the distribution of the innovations. The method has a certain amount of arbitrariness when specifying suitable moment restrictions which sometimes may lead to lower efficiency than Maximum Likelihood 1. On the other hand, for simpler models and with proper moment restrictions, one can show that both methods give identical estimators. If the number of moment restrictions equals the number of parameters, some effort can be spared by solving a system of equations rather than solving an optimization problem. The method is then referred to as the Method of Moments. At a first glance, the method might seem to be a purely theoretical construction since the true parameter values are needed to calculate the weight matrix which in turn are needed to calculate the estimates. It turns out that this is not a problem since it can be shown that the Generalized Method of Moments is a consistent (but not efficient) for all positive definite weight matrices. The pseudo algorithm for the estimation becomes: 1. Estimate the parameters with an arbitrary pos. def. weight matrix W, for instance the identity matrix to get an initial estimate of the parameters. ( 1 ) T ( 1 ) θ = arg min J N (θ) = arg min f(xt, θ) W f(xt, θ) N N 2. Use these parameters to calculate the optimal weight matrix, chosen as the inverse of the covariance matrix for the moment restrictions. 3. Estimate new parameter with the new weight matrix, estimate a new weight matrix from these parameters and so on until convergence Tests Complex models can be compared with LR-tests if we use Maximum Likelihood techniques. There are similar tests for the GMM-method. The test statistic is given by: [ N J N ( θ Restricted ) J N ( θ Unrestricted ) and under the null-hypothesis it has an asymptotic χ 2 (k)-distribution where k is the difference in the number of parameters. 1 This is not always a bad thing, one can for instance choose economically sound restrictions. 2
3 Note: You can not use MLmax to derive confidence intervals for your estimates when using GMM. You should instead use fminunc and compute the Hessian and the Jacobian of the moment restrictions separately. 3.2 Regression As an introduction to GMM we will consider a linear regression model. The data set is the same as in computer exercise 1. Load the data set by writing: > > load regdata.mat Assignment: Use your results from from preparatory question 1 to estimate the parameters. Calculate a 95% confidence interval and compare your results with the results from exercise 1. Assignment: Test the linear model against the quadratic model: y i = β 0 + β 1 x i + β 2 x 2 i + ε i. Conclusions? 3.3 Parameter estimation in stochastic differential equations You will meet two different processes in this part of the computer exercise: the Cox- Ingersoll-Ross (CIR) process and the CKLS-process. the CIR process is often used to model interest rates and is defined as: This can be generalized to the CKLS model given by: dr t = κ(θ r t )dt + σ r t dw t. (2) dr t = κ(θ r t )dt + σr γ t dw t. (3) These models can also be used to model the volatility in stochastic variance models, since they always remain positive. Load the data by writing: > > load cirdata.mat > > load cklsdata.mat By discretizing the stochastic differential equation one can use for instance GMM for parameter estimation. The discretization is usually done with an Euler scheme but of course higher order methods works even better. The advantage of using the discretized process is that we can estimate complex models without much problem. 3
4 Assignment: Estimate the parameters and their covariance matrices in the models above by using GMM. Assignment: Estimate the parameters and their covariance matrices in the models above by using approximate Quasi Maximum Likelihood (by discretizing the model and using the likelihood generated by the approximate model). 3.4 Improved approximate estimators All problems in this Sections are extra assignments Transformation of data for the CIR-model and the Shoji and Ozaki approximated likelihood (Extra) Using the transformation y t = 2 r t we obtain a state independent diffussion term. Itos formula gives [ 1 dy t = κ(θ r t ) rt 2r 3/2 2 σ2 r t dt + 1 rt σdw t rt t [ 2 = κ(θ y2 t y t 4 ) + 4 σ 2 y2 t dt + σdw yt 3 t 8 [ 2 = κθ κ y t y t σ 2 dt + σdw t y t 4 [ ) 2 = (κθ σ2 κ y t dt + σdw t 4 2 y t = µ(y t )dt + σdw t 4
5 We can now use the Shoji and Ozaki [1998 method to approximate the likelihood with a Gaussian likelihood. N log(p ytk y tk 1 (y tk )) k=2 N ( (ytk m k ) 2 1 ) 2v k 2 log(2πv k) k=2 m k = y tk 1 + a k K k + σ2 c k (K k b k k ) b k 2b 2 k v k = σ2 (exp(2b k k ) 1) 2b k k = t k t k 1 K k = exp(b k k ) 1 a k = µ(y tk 1 ) = 2 y tk 1 b k = µ (y tk 1 ) = 2 yt 2 k 1 c k = µ (y tk 1 ) = 4 y 3 t k 1 ) (κθ σ2 κ y t k ) (κθ σ2 κ 4 2 ) (κθ σ2 4 Assignment: Estimate the parameters and their covariance matrices in the CIR model by using the Shoji and Ozaki approximated likelihood Simulated Maximum Likelihood (Extra) Sometimes it is possible to use Maximum Likelihood on a simulated likelihood in stochastic differential equations. Assume that we have observations y n, n = 1,... N from some model. Also assume that the sampling interval is. We now want to find the likelihood N L(θ) = p(y i y i 1, θ). i=1 If this is not available in closed form, one might approximate it using simulation. First we discretize the dynamics of the stochastic differential equation using a scheme that gives a simple transition probability, preferably Euler. Then, divide the interval into M subintervals of length δ = /M. The idea is the to simulate K trajectories on a grid of size δ up to subinterval M 1 starting in y n. Do this for every time n, resulting in K samples at every n. Under the Euler discretization, the transition density from M 1 to M is Gaussian with mean µ n,k and standard deviation σ n,k, both given by the 5
6 model in question. Then we approximate the likelihood by the following: L(θ) = N p(y n+1 y n, θ) n=1 ( N 1 K n=1 ) K ϕ(y n+1, µ n,k θ, σ n,k θ ) k=1 where ϕ(y, µ, σ) is the density of the Normal distribution. The log likelihood is given by ( ( )) N N 1 K l(θ) = log(l(θ)) = log(p(y n+1 y n ), θ) log ϕ(y n+1, µ n,k θ, σ n,k θ ) K n=1 We then use numerical optimization techniques to maximize the loglikelihood. In this exercise we keep the number of intermediate steps fairly small, say M = 2 or 3. If M is too big, the variance of the sample will become large and make the estimation difficult. The remedy for this is to use an importance sampler. n=1 k=1 Assignment: Estimate the parameters and their covariance matrices in the models above by using Simulated Maximum Likelihood. Assignment: model. Do the same thing for the transformed data using the transformed CIR Hint: It is usually a good idea to use the same sequence of random numbers each time the likelihood is evaluated. This is referred to as Common Random Numbers and is a way to avoid the Monte Carlo error in the minimization. Practically this is done by drawing a large number of random numbers, in this case N (M 1) K and use these as input to the likelihood Exact Moments (GMM/EF) (Extra) Sometimes the bias that discretization brings can be avoided by calculating the moments exactly. The moments can be calculated using the so called Dynkin s formula, which is nothing more than an application of the Itô formula. Dynkin s formula is given by: [ t E [f(x t ) X t 1 = f(x 0 ) + E Af(X s )ds X t 1 (4) where the operator A (standard Itô s formula) is defined as, n f Af = µ i + 1 n ( ) σσ T 2 f. x i 2 i,j x i x j i=1 6 i,j=1 0
7 and where the process {X t } is the solution to dx(t) = µ(t, X(t))dt+σ(t, X(t))dW (t). For our purposes {X t } is a one-dimensional process which makes Dynkin s formula a lot easier (n = 1). Hint: To show how Dynkins formula can be used we will study an example with f(x) = x 2. Suppose that we want to calculate the moments of the CIR process. Ax 2 = κ(θ x) x2 x + 1 ) ((σ 2 x) 2 x 2 2 x 2 = κ(θ x) 2x σ2 x 2 = 2κx 2 + (2κθ + σ 2 )x By plugging this into the formula we get: E [ [ t X Xt 2 = x E X ( 2κXs 2 + (2κθ + σ 2 )X s )ds 0 Finally we assume that the expectation and the time integral can be interchanged, followed by taking the derivatives with respect to t on both sides. E X [X 2 t t = 2κE X [ X 2 t + (2κθ + σ 2 )E X [X t. Extra: Estimate the parameters and their covariance matrix for the Cox-Ingersoll- Ross process by using exact moment, that is using f(x) = {x, x 2, x 3 }. 3.5 Exact likelihood for the CIR model (Extra) The CIR process has a known transition density. If we let k = t k t k 1, c = 2κ/(σ 2 (1 exp( κ k ))), and Y = 2cX, then Y tk Y tk 1 is distributed as a noncentral chi-squared with 2κθ/σ 2 degrees of freedom and noncentrality parameter Y tk 1 exp(κ k ). Using this we can write up the density as p Xtk X tk 1 (x tk x tk 1 ) = ce u v ( v u) q/2 Iq (2 uv) (5) where u = cx tk 1 exp( κ k ), v = cx tk, q = 2κθ/σ 2 1 and I q (z) is the modified Bessel function of the first kind of order q (besseli(q,z) in Matlab). Extra: Estimate the parameters and their covariance matrix for the Cox-Ingersoll- Ross model by using the exact likelihood. Compare the result with the previous estimates where you used GMM, Simulated ML, Shoji and Ozaki and Exact moments. 7
8 4 Feedback Comments and ideas relating to the computer exercise are always welcome. Send them to Erik Lindström, or call MATLAB routines fminsearch Simplex based optimisation routine. fminunc Numerical minimization of a multidimensional function. The routine is based on quasi-newton (BFGS) and it can return the minimising parameter value, the minimal function value and the Hessian. The cryptic name comes from function minimization unconditional. MLmax Customized Quasi-Newton based optimization algorithm for maximum likelihood estimation. Maximises the likelihood function by using the score function s quadratic variations to estimate Fisher Information matrix. Needs the log-likelihood returned as a vector. >> [xout,logl,covm=mlmax(@lnl,x0,indata) Referenser Madsen, H., Nielsen, J. N., Lindström, E., Baadsgaard, M. and Holst, J. Statistics in Finance, IMM;DTU; Lyngby och Matematisk Statistik, LTH, Lund Nolsøe, K., Madsen, H., Nielsen, J. N. and Baadsgaard, M. Lecture Notes in estimation functions, IMM;DTU; Lyngby Shoji, I. and Ozaki, T. (1998). Estimation for nonlinear stochastic differential equations by a local linearization method, Stochastic Analysis and Applications, 16,
Computer 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 informationComputer Exercise 2 Simulation
Lund University with Lund Institute of Technology Valuation of Derivative Assets Centre for Mathematical Sciences, Mathematical Statistics Spring 2010 Computer Exercise 2 Simulation This lab deals with
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 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 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 informationIntro to GLM Day 2: GLM and Maximum Likelihood
Intro to GLM Day 2: GLM and Maximum Likelihood Federico Vegetti Central European University ECPR Summer School in Methods and Techniques 1 / 32 Generalized Linear Modeling 3 steps of GLM 1. Specify the
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 informationThe method of Maximum Likelihood.
Maximum Likelihood The method of Maximum Likelihood. In developing the least squares estimator - no mention of probabilities. Minimize the distance between the predicted linear regression and the observed
More informationApplication of MCMC Algorithm in Interest Rate Modeling
Application of MCMC Algorithm in Interest Rate Modeling Xiaoxia Feng and Dejun Xie Abstract Interest rate modeling is a challenging but important problem in financial econometrics. This work is concerned
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 informationEC316a: Advanced Scientific Computation, Fall Discrete time, continuous state dynamic models: solution methods
EC316a: Advanced Scientific Computation, Fall 2003 Notes Section 4 Discrete time, continuous state dynamic models: solution methods We consider now solution methods for discrete time models in which decisions
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 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 informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Mike Giles (Oxford) Monte Carlo methods 2 1 / 24 Lecture outline
More informationM.I.T Fall Practice Problems
M.I.T. 15.450-Fall 2010 Sloan School of Management Professor Leonid Kogan Practice Problems 1. Consider a 3-period model with t = 0, 1, 2, 3. There are a stock and a risk-free asset. The initial stock
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 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 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 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 informationModule 2: Monte Carlo Methods
Module 2: Monte Carlo Methods Prof. Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute MC Lecture 2 p. 1 Greeks In Monte Carlo applications we don t just want to know the expected
More informationCalibration of Interest Rates
WDS'12 Proceedings of Contributed Papers, Part I, 25 30, 2012. ISBN 978-80-7378-224-5 MATFYZPRESS Calibration of Interest Rates J. Černý Charles University, Faculty of Mathematics and Physics, Prague,
More informationLecture 2 - Calibration of interest rate models and optimization
- Calibration of interest rate models and optimization Elisabeth Larsson Uppsala University, Uppsala, Sweden March 2015 E. Larsson, March 2015 (1 : 23) Introduction to financial instruments Introduction
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 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 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 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 information12. Conditional heteroscedastic models (ARCH) MA6622, Ernesto Mordecki, CityU, HK, 2006.
12. Conditional heteroscedastic models (ARCH) MA6622, Ernesto Mordecki, CityU, HK, 2006. References for this Lecture: Robert F. Engle. Autoregressive Conditional Heteroscedasticity with Estimates of Variance
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 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 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 informationFinal exam solutions
EE365 Stochastic Control / MS&E251 Stochastic Decision Models Profs. S. Lall, S. Boyd June 5 6 or June 6 7, 2013 Final exam solutions This is a 24 hour take-home final. Please turn it in to one of the
More informationImplementing an Agent-Based General Equilibrium Model
Implementing an Agent-Based General Equilibrium Model 1 2 3 Pure Exchange General Equilibrium We shall take N dividend processes δ n (t) as exogenous with a distribution which is known to all agents There
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 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 informationLecture IV Portfolio management: Efficient portfolios. Introduction to Finance Mathematics Fall Financial mathematics
Lecture IV Portfolio management: Efficient portfolios. Introduction to Finance Mathematics Fall 2014 Reduce the risk, one asset Let us warm up by doing an exercise. We consider an investment with σ 1 =
More informationModel Estimation. Liuren Wu. Fall, Zicklin School of Business, Baruch College. Liuren Wu Model Estimation Option Pricing, Fall, / 16
Model Estimation Liuren Wu Zicklin School of Business, Baruch College Fall, 2007 Liuren Wu Model Estimation Option Pricing, Fall, 2007 1 / 16 Outline 1 Statistical dynamics 2 Risk-neutral dynamics 3 Joint
More informationdt + ρσ 2 1 ρ2 σ 2 B i (τ) = 1 e κ iτ κ i
A 2D Gaussian model (akin to Brigo & Mercurio Section 4.2) Suppose where dx(t) = ( κ1 0 0 κ 2 ) ( X1 (t) X 2 (t) In this case we find (BLACKBOARD) that r(t) = δ 0 + X 1 (t) + X 2 (t) ) ( σ1 0 dt + ρσ 2
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 informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Abdul-Lateef Haji-Ali Based on slides by: Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Haji-Ali (Oxford)
More informationUser 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 informationRISK-NEUTRAL VALUATION AND STATE SPACE FRAMEWORK. JEL Codes: C51, C61, C63, and G13
RISK-NEUTRAL VALUATION AND STATE SPACE FRAMEWORK JEL Codes: C51, C61, C63, and G13 Dr. Ramaprasad Bhar School of Banking and Finance The University of New South Wales Sydney 2052, AUSTRALIA Fax. +61 2
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 informationMonte Carlo Methods in Financial Engineering
Paul Glassennan Monte Carlo Methods in Financial Engineering With 99 Figures
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 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 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 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 informationdt+ ρσ 2 1 ρ2 σ 2 κ i and that A is a rather lengthy expression that we may or may not need. (Brigo & Mercurio Lemma Thm , p. 135.
A 2D Gaussian model (akin to Brigo & Mercurio Section 4.2) Suppose where ( κ1 0 dx(t) = 0 κ 2 r(t) = δ 0 +X 1 (t)+x 2 (t) )( X1 (t) X 2 (t) ) ( σ1 0 dt+ ρσ 2 1 ρ2 σ 2 )( dw Q 1 (t) dw Q 2 (t) ) In this
More informationLikelihood Methods of Inference. Toss coin 6 times and get Heads twice.
Methods of Inference Toss coin 6 times and get Heads twice. p is probability of getting H. Probability of getting exactly 2 heads is 15p 2 (1 p) 4 This function of p, is likelihood function. Definition:
More informationدرس هفتم یادگیري ماشین. (Machine Learning) دانشگاه فردوسی مشهد دانشکده مهندسی رضا منصفی
یادگیري ماشین توزیع هاي نمونه و تخمین نقطه اي پارامترها Sampling Distributions and Point Estimation of Parameter (Machine Learning) دانشگاه فردوسی مشهد دانشکده مهندسی رضا منصفی درس هفتم 1 Outline Introduction
More informationQuasi-Monte Carlo for Finance
Quasi-Monte Carlo for Finance Peter Kritzer Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy of Sciences Linz, Austria NCTS, Taipei, November 2016 Peter Kritzer
More informationNumerical Simulation of Stochastic Differential Equations: Lecture 1, Part 2. Integration For deterministic h : R R,
Numerical Simulation of Stochastic Differential Equations: Lecture, Part Des Higham Department of Mathematics University of Strathclyde Lecture, part : SDEs Ito stochastic integrals Ito SDEs Examples of
More information"Vibrato" Monte Carlo evaluation of Greeks
"Vibrato" Monte Carlo evaluation of Greeks (Smoking Adjoints: part 3) Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Oxford-Man Institute of Quantitative Finance MCQMC 2008,
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 informationMaximum Likelihood Estimation
Maximum Likelihood Estimation EPSY 905: Fundamentals of Multivariate Modeling Online Lecture #6 EPSY 905: Maximum Likelihood In This Lecture The basics of maximum likelihood estimation Ø The engine that
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 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 information[D7] PROBABILITY DISTRIBUTION OF OUTSTANDING LIABILITY FROM INDIVIDUAL PAYMENTS DATA Contributed by T S Wright
Faculty and Institute of Actuaries Claims Reserving Manual v.2 (09/1997) Section D7 [D7] PROBABILITY DISTRIBUTION OF OUTSTANDING LIABILITY FROM INDIVIDUAL PAYMENTS DATA Contributed by T S Wright 1. Introduction
More informationGov 2001: Section 5. I. A Normal Example II. Uncertainty. Gov Spring 2010
Gov 2001: Section 5 I. A Normal Example II. Uncertainty Gov 2001 Spring 2010 A roadmap We started by introducing the concept of likelihood in the simplest univariate context one observation, one variable.
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 informationEco504 Spring 2010 C. Sims FINAL EXAM. β t 1 2 φτ2 t subject to (1)
Eco54 Spring 21 C. Sims FINAL EXAM There are three questions that will be equally weighted in grading. Since you may find some questions take longer to answer than others, and partial credit will be given
More informationChapter 7: Estimation Sections
1 / 40 Chapter 7: Estimation Sections 7.1 Statistical Inference Bayesian Methods: Chapter 7 7.2 Prior and Posterior Distributions 7.3 Conjugate Prior Distributions 7.4 Bayes Estimators Frequentist Methods:
More information1. You are given the following information about a stationary AR(2) model:
Fall 2003 Society of Actuaries **BEGINNING OF EXAMINATION** 1. You are given the following information about a stationary AR(2) model: (i) ρ 1 = 05. (ii) ρ 2 = 01. Determine φ 2. (A) 0.2 (B) 0.1 (C) 0.4
More informationComputational Finance Least Squares Monte Carlo
Computational Finance Least Squares Monte Carlo School of Mathematics 2019 Monte Carlo and Binomial Methods In the last two lectures we discussed the binomial tree method and convergence problems. One
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 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 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 informationAmath 546/Econ 589 Univariate GARCH Models
Amath 546/Econ 589 Univariate GARCH Models Eric Zivot April 24, 2013 Lecture Outline Conditional vs. Unconditional Risk Measures Empirical regularities of asset returns Engle s ARCH model Testing for ARCH
More informationBrooks, Introductory Econometrics for Finance, 3rd Edition
P1.T2. Quantitative Analysis Brooks, Introductory Econometrics for Finance, 3rd Edition Bionic Turtle FRM Study Notes Sample By David Harper, CFA FRM CIPM and Deepa Raju www.bionicturtle.com Chris Brooks,
More informationStochastic Modelling in Finance
in Finance Department of Mathematics and Statistics University of Strathclyde Glasgow, G1 1XH April 2010 Outline and Probability 1 and Probability 2 Linear modelling Nonlinear modelling 3 The Black Scholes
More informationOPTIMAL PORTFOLIO CONTROL WITH TRADING STRATEGIES OF FINITE
Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 005 Seville, Spain, December 1-15, 005 WeA11.6 OPTIMAL PORTFOLIO CONTROL WITH TRADING STRATEGIES OF
More informationLecture 3. Sergei Fedotov Introduction to Financial Mathematics. Sergei Fedotov (University of Manchester) / 6
Lecture 3 Sergei Fedotov 091 - Introduction to Financial Mathematics Sergei Fedotov (University of Manchester) 091 010 1 / 6 Lecture 3 1 Distribution for lns(t) Solution to Stochastic Differential Equation
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 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 informationAn Efficient Numerical Scheme for Simulation of Mean-reverting Square-root Diffusions
Journal of Numerical Mathematics and Stochastics,1 (1) : 45-55, 2009 http://www.jnmas.org/jnmas1-5.pdf JNM@S Euclidean Press, LLC Online: ISSN 2151-2302 An Efficient Numerical Scheme for Simulation of
More informationA THREE-FACTOR CONVERGENCE MODEL OF INTEREST RATES
Proceedings of ALGORITMY 01 pp. 95 104 A THREE-FACTOR CONVERGENCE MODEL OF INTEREST RATES BEÁTA STEHLÍKOVÁ AND ZUZANA ZÍKOVÁ Abstract. A convergence model of interest rates explains the evolution of the
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 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 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 informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Abdul-Lateef Haji-Ali Based on slides by: Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Haji-Ali (Oxford)
More informationCalibration Lecture 4: LSV and Model Uncertainty
Calibration Lecture 4: LSV and Model Uncertainty March 2017 Recap: Heston model Recall the Heston stochastic volatility model ds t = rs t dt + Y t S t dw 1 t, dy t = κ(θ Y t ) dt + ξ Y t dw 2 t, where
More informationPricing Variance Swaps under Stochastic Volatility Model with Regime Switching - Discrete Observations Case
Pricing Variance Swaps under Stochastic Volatility Model with Regime Switching - Discrete Observations Case Guang-Hua Lian Collaboration with Robert Elliott University of Adelaide Feb. 2, 2011 Robert Elliott,
More informationChapter 7: Estimation Sections
Chapter 7: Estimation Sections 7.1 Statistical Inference Bayesian Methods: 7.2 Prior and Posterior Distributions 7.3 Conjugate Prior Distributions Frequentist Methods: 7.5 Maximum Likelihood Estimators
More informationContents Critique 26. portfolio optimization 32
Contents Preface vii 1 Financial problems and numerical methods 3 1.1 MATLAB environment 4 1.1.1 Why MATLAB? 5 1.2 Fixed-income securities: analysis and portfolio immunization 6 1.2.1 Basic valuation of
More informationWrite legibly. Unreadable answers are worthless.
MMF 2021 Final Exam 1 December 2016. This is a closed-book exam: no books, no notes, no calculators, no phones, no tablets, no computers (of any kind) allowed. Do NOT turn this page over until you are
More informationVaR Estimation under Stochastic Volatility Models
VaR Estimation under Stochastic Volatility Models Chuan-Hsiang Han Dept. of Quantitative Finance Natl. Tsing-Hua University TMS Meeting, Chia-Yi (Joint work with Wei-Han Liu) December 5, 2009 Outline Risk
More 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 informationHomework 3: Asset Pricing
Homework 3: Asset Pricing Mohammad Hossein Rahmati November 1, 2018 1. Consider an economy with a single representative consumer who maximize E β t u(c t ) 0 < β < 1, u(c t ) = ln(c t + α) t= The sole
More informationEmpirical Approach to the Heston Model Parameters on the Exchange Rate USD / COP
Empirical Approach to the Heston Model Parameters on the Exchange Rate USD / COP ICASQF 2016, Cartagena - Colombia C. Alexander Grajales 1 Santiago Medina 2 1 University of Antioquia, Colombia 2 Nacional
More informationCash Accumulation Strategy based on Optimal Replication of Random Claims with Ordinary Integrals
arxiv:1711.1756v1 [q-fin.mf] 6 Nov 217 Cash Accumulation Strategy based on Optimal Replication of Random Claims with Ordinary Integrals Renko Siebols This paper presents a numerical model to solve the
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 informationIntroduction to the Maximum Likelihood Estimation Technique. September 24, 2015
Introduction to the Maximum Likelihood Estimation Technique September 24, 2015 So far our Dependent Variable is Continuous That is, our outcome variable Y is assumed to follow a normal distribution having
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 informationMonte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMSN50)
Monte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMSN50) Magnus Wiktorsson Centre for Mathematical Sciences Lund University, Sweden Lecture 5 Sequential Monte Carlo methods I January
More informationHigh Dimensional Bayesian Optimisation and Bandits via Additive Models
1/20 High Dimensional Bayesian Optimisation and Bandits via Additive Models Kirthevasan Kandasamy, Jeff Schneider, Barnabás Póczos ICML 15 July 8 2015 2/20 Bandits & Optimisation Maximum Likelihood inference
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 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 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 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 informationApplied Stochastic Processes and Control for Jump-Diffusions
Applied Stochastic Processes and Control for Jump-Diffusions Modeling, Analysis, and Computation Floyd B. Hanson University of Illinois at Chicago Chicago, Illinois siam.. Society for Industrial and Applied
More information