User s Guide for the Matlab Library Implementing Closed Form MLE for Diffusions
|
|
- Simon Short
- 5 years ago
- Views:
Transcription
1 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 Abstract This document explains the use of the attached Matlab code for estimating the parameters of diffusions using closed-form maximum-likelihood. Copyright c Yacine Aït-Sahalia. Redistribution and use in source and binary forms with or without modification are permitted provided that the following conditions are met: i Redistributions of source code must retain the above copyright notice this list of conditions and the following disclaimer; ii Redistributions in binary form must reproduce the above copyright notice this list of conditions and the following disclaimer in the documentation and/or other materials provided with the distribution. DISCLAIMER: THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CON- TRIBUTORS AS IS AND ANY EXPRESS OR IMPLIED WARRANTIES INCLUDING BUT NOT LIMITED TO THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PAR- TICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT INDIRECT INCIDENTAL SPECIAL EXEM- PLARY OR CONSEQUENTIAL DAMAGES INCLUDING BUT NOT LIMITED TO PROCURE- MENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE DATA OR PROFITS; OR BUSI- NESS INTERRUPTION HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY WHETHER IN CONTRACT STRICT LIABILITY OR TORT INCLUDING NEGLIGENCE OR OTHERWISE ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. 1 References This code implements in Matlab the closed-form maximum-likelihood estimation method for diffusions developed in: Aït-Sahalia Y Transition Densities for Interest Rate and Other Nonlinear Diffusions Journal of Finance Aït-Sahalia Y. 00 Maximum-Likelihood Estimation of Discretely-Sampled Diffusions: A Closed- Form Approximation Approach Econometrica Financial support from the NSF under grants SES SES DMS and SES is gratefully acknowledged. Jia Li provided excellent research assistance in developing this Matlab code. Princeton NJ Phone: yacine@princeton.edu. 1
2 Aït-Sahalia Y. 008 Closed-Form Likelihood Expansions for Multivariate Diffusions Annals of Statistics Extensions to models with latent factors can be found in: Aït-Sahalia Y. and Kimmel R. 007 Maximum Likelihood Estimation of Stochastic Volatility Models Journal of Financial Economics Aït-Sahalia Y. and Kimmel R. 010 Estimating Affine Multifactor Term Structure Models Using Closed-Form Likelihood Expansions Journal of Financial Economics Installation This code uses anonymous functions which require Matlab version 7 R14 released in 004 or later..1 Installing the Package The only thing you need to do is to download the file ClosedFormMLE.zip and unzip it somewhere in your computer. You have two options: 1. You can work under the folder ClosedFormMLE\ you may rename this folder this does not matter keeping its subfolders in place and save all your files in that folder.. Or you add this library to your Matlab s search path. We recommend the first option for casual users of Matlab. This package contains the following folders and files: Type Name Usage Folder DERIVESTsuite Numerical differentiation Folder Models A library of closed-form log-likelihood functions Matlab function mymle.m Implementing MLE Matlab function mymle summary.m Post estimation analysis Matlab function example cev.m Demo: CEV model estimation Matlab function logdensityloglik.m A helper function Note: the package DERIVESTsuite was written by John D Errico and downloaded from: It is enclosed here as a convenience for endusers. Please refer to the package for its license. The folder Models contains the formulae for the closed-form approximations of log likelihood functions for various stochastic differential equations SDE that have been implemented in this toolbox. Section 3 below gives the list of models that are curently implemented. Since this is an ongoing project we will keep adding models into the library. Interested users can visit the author s website for updates and request additional models by . The folder DERIVESTsuite contains routines for numerical differentiation necessary to compute Fisher s information so that we can obtain the asymptotic standard errors of the estimates. Neither one of these folders need to be edited or modified. The two functions mymle.m and mymle summary.m need not be modified either. The simplest way to use this code is to first run example cev.m and then modify the file example cev.m as needed: estimate a different model replace the simulated data with real data etc.
3 . An Example: Estimating the CEV Model on Simulated Data Start Matlab and switch your current directory to ClosedFormMLE\. In Matlab open the file example cev.m and run it by pressing the F5 key in Matlab s editor: Debug >> Run Step 1 of the code simulates a time series according to the CEV model i.e. Model U3 see with parameters a = 0.08 b = 0.5 c = 0.7 d = 1.5. The sampling interval is set to be 1/5 which should be thought of as weekly data. There are n = 500 observations saved in the variable x. Note: since the sample size is small the estimation is fast but the estimates are inaccurate. You may also want to experiment with n = to get statistically better results. Step does the estimation with the function mymle.m. The code is output=mymle@modelu3xdelparam0; is the function handle of the log likelihood function x is our data del is our sampling interval = 1/5 and param0 is the initial value of the optimization procedure underlying MLE. The user may need to experiment with different initial values to find the global maximum. In this example we set the initial value as the true value of the parameters that is param0=[ ]. See Section 3 for the ordering of parameters in each model. If we wanted to estimate the other model say Model U1 see the syntax would be output=mymle@modelu1xdelparam0foru1; where the initial value param0foru1 should be 3-dimensional since Model U1 has only 3 parameters. The output is a Matlab structure which packages the results of the estimation. You do not have to worry if you are not familiar with this data structure. We provide the output in step 3. As the code runs it first maximize the log likelihood function and then prints the MLE estimates on the screen. It also tells you whether the maximization procedure converges under the default tolerance. After that it moves onto computing the standard error and then prints the estimates standard error constructed from the inverse of Hessian and the misspecification-robust standard error i.e. the sandwich estimate. Step 3 does the post estimation analysis. The code is mymle summaryoutput; where you take the output structure from step directly as your input in this step. The code does two things. First it reprints the table of estimates standard errors robust standard errors and exit status of the maximization procedure as below: 3
4 Parameter 1 Parameter Parameter Parameter 4 Figure 1: Example: Log-Likelihood Near the Solution ****************************** COEF SE Robust SE ****************************** fminsearch converged to a solution. Secondly it plots the marginal log likelihood for each parameter in a neighborhood of the estimates. See Figure 1. The blue curve is the likelihood function while each red star is our estimate for the corresponding parameter. This gives us a way to visually check whether the maximization is successful of course this says nothing about global optimality!..3 Syntax of the Two Functions for reference only.3.1 mymle Syntax: output = mymlelogdensity x del param0 options; Inputs: logdensity is the function handle of the model @ModelB4... x: is the data. This is a n k matrix where n is the number of observations and k is the dimensionality of your model. For univariate models k = 1 and for bivariate models k =... 4
5 del: is the sampling interval of the time series param0 : initial parameter value of the maximization procedure. The dimensionality and ordering of its entries should be consistent with the underlying model listed in Section 3. options : this input is optional. The input is a Matlab option structure for its optimization toolbox. See Matlab s help document of optimset for more details. Output is structure consisting the following fields: exitmsg: exit status of the optimization procedure param: the MLE estimates see Section 3 for the ordering of parameters. variance: covariance estimates variance robust: robust covariance estimates se: standard errors se robust: robust standard errors objfun: the function handle of the negative log likelihood function used internally in the optimization procedure loglik: the optimized value of log likelihood exitflag: exit status reported by fminsearch..3. mymle summary Syntax: mymle summarymle output The input of this function is the output structure from mymle. 3 List of Models Currently Implemented The models below are in no particular order. The 9 canonical affine models and semi-affine models are implemented separately see below. There is of course some duplication since U4 and U7 are affine for instance. 3.1 Univariate Models Univariate models have the following canonical form: dx t = µ X t ; parameters dt + σ X t ; parameters dw t Model U1 µ x = ax + bx σ x = σx 3/ Parameters: a b σ. 5
6 3.1. Model U µ x = a + bx σ x = dx Parameters: a b d Model U3 µ x = b a x σ x = cx ρ Parameters: a b c ρ. Note: do not let the parameter ρ go below 1. Special cases for ρ = 1 and 3/ Model U4 µ x = κ α x σ x = σx 1/ Parameters: κ α σ Model U5 µ x = θ 0 + θ 1 x + θ x + θ 3 x 3 σ x = γx ρ Parameters: γ ρ θ 0 θ 1 θ θ 3. Note: do not let the parameter ρ go below 1. Special cases for ρ = 1 and 3/ Model U6 µ x = a + bx + cx + dx 3 σ x = f Parameters: a b c d f Model U7 µ x = κ α x σ x = σ Parameters: κ α σ. 6
7 3.1.8 Model U8 µ x = a 1 x + a 0 + a 1 x + a x σ x = σx ρ Parameters: a 1 a 0 a 1 a σ ρ. Note: do not let the parameter ρ go below 1. Special cases for ρ = 1 and 3/ Model U9 µ x = a 1 x + a 0 + a 1 x + a x σ x = b 0 + b 1 x + b x b3 1/ Parameters: a 1 a 0 a 1 a b 0 b 1 b b Model U10 Note: do not let the parameter ρ go below 1. µ x = a 1 x + a 0 + a 1 x + a x σ x = b 0 + b 1 x + b x b3 Parameters: a 1 a 0 a 1 a b 0 b 1 b b Model U11 µ x = a + bx σ x = f + dx Parameters: a b f d Model U1 µ x = β x αx3 σ x = γx 1/ Parameters: β α γ Model U13 µ x = a 0 + a 1 x + a 1x + a x + a 3 x 3 σ x = σx ρ Parameters: a 0 a 1 a 1 a a 3 σ ρ. Note: do not let the parameter ρ go below 1. Special cases for ρ = 1 and 3/. 7
8 Model U14 µ x = h + ax + bx σ x = vx 3/ Parameters: h a b v. 3. Bivariate Models Consider a -dimensional process X t = X 1t X t determined by or written componentwise dx1t dx t = µ1 X 1t X t dt + µ X 1t X t dx t = µ X t dt + Σ X t dw t Σ 11 X 1t X t Σ 1 X 1t X t Σ 1 X 1t X t Σ X 1t X t dw1t. The Brownian motions are independent. Any correlation between the components is introduced by the off-diagonal terms in the Σ matrix. dw t 3..1 Model B1 a + bx µ x 1 x = c + dx Σ x 1 x = Parameters: r h a b c d. r x 0 h 1 r x 3.. Model B a0 + a 1 x 1 + a x µ x 1 x = b 0 + b 1 x 1 + b x Σ x 1 x = c 0 + c 1 x 1 + c x 0 0 d 0 + d 1 x 1 + d x Parameters: a 0 a 1 a b 0 b 1 b c 0 c 1 c d 0 d 1 d Model B3 µ x / µ x 1 x = α + βx x 0 Σ x 1 x = σρx γ σ 1 ρ x γ Parameters: µ α β σ ρ γ. 8
9 3..4 Model B4 a 0 + a 1 x µ x 1 x = b a x + λgx β a + f x a 1 r a + f x a r a + f x a Σ x 1 x = Parameters: a 0 a 1 a b λ g β f r. 0 gx β 3..5 Model B5 bx1 µ x 1 x = c dx hx 1 x 0 Σ x 1 x = gr x g 1 r x Parameters: b c d h g r Model B6 This is the Heston model with x 1 being the log of the asset price and x its stochastic volatility. m x / µ x 1 x = a bx x 0 Σ x 1 x = sr x s 1 r x Parameters: m a b s r. or Note that this model is observationally equivalent to the ones with 1 r x r x Σ x 1 x = 0 s x Σ x 1 x = r x 1 r x s x Model B7 0 µ x 1 x = a 1 a x x1 γ Σ x 1 x = ηx 1 x γ x 0 Parameters: a 1 a γ η. 9
10 3..8 Model B8 a + bx1 µ x 1 x = cx dx γ 1 Σ x 1 x = ex 0 0 f Parameters: a b c d γ f Model B9 θ5 1 µ x 1 x = ex θ θ 1 x e x 0 Σ x 1 x = θ 3 θ 4 1 θ 4 θ 3 Parameters: θ 1 θ θ 3 θ 4 θ Model B10 b1 a 1 x 1 µ x 1 x = b a x g 1 0 Σ x 1 x = 0 g x Parameters: a 1 a b 1 b g 1 g Model B11 k1 + k x µ x 1 x = κ θ x 1 ρ x ρ x Σ x 1 x = 0 σx Parameters: k 1 k ρ κ θ σ Model B1 ax1 µ x 1 x = bx cx 1 e x 0 Σ x 1 x = dr d 1 r Parameters: a b c d r. 10
11 3..13 Model B13 b11 a 1 x 1 + b 1 a x µ x 1 x = b 1 a 1 x 1 + b a x σ 11 σ 1 Σ x 1 x = σ 1 σ Parameters: a 1 a b 11 b 1 b 1 b σ 11 σ 1 σ 1 σ Model B14 k1 x x 1 µ x 1 x = k θ x σ 1 x1 0 Σ x 1 x = 0 σ x Parameters: k 1 k θ σ 1 σ Model B15 a + bx1 µ x 1 x = fx 1 + dx x1 0 Σ x 1 x = h 1 + gx1 Parameters: g h a b f d Model B16 a + bx1 + gx µ x 1 x = d + ηx 1 + fx x1 0 Σ x 1 x = 0 x Parameters: a b g d η f Model B17 a00 a 1 + a x / + n 0 1 g 1 + nu 1 g 1 a 1 + a x b+d µ x 1 x = a 01 + a 11 x + nu 1 g 11 a 1 + a x b+d 1 g Σ x 1 x = 1 a1 + a x g 1 a1 + a x 0 g 11 a 1 + a x b Parameters: a 00 a 01 a 11 A 1 A n 0 nu 1 g 1 g 11 b d. 11
12 3..18 Model B Model B19 b1 x 1 µ x 1 x = a + b x g 11 e x1 0 Σ x 1 x = g r g 1 r Parameters: b 1 a b g 11 g r. b1 x 1 µ x 1 x = a + b x e x 0 Σ x 1 x = g r g 1 r Parameters: b 1 a b g r Model B0 a1 + b 1 x 1 µ x 1 x = a + b x x 0 Σ x 1 x = gr x g 1 r x P arameters : a 1 b 1 a b g r Model B1 a 1 b 1 x 1 µ x 1 x = a 1 b 1 x 1 + a b x x1 0 Σ x 1 x = g 1 x1 g x1 Parameters: a 1 b 1 a 1 a b g 1 g. 3.. Model B k1 + k x µ x 1 x = k 3 a x 1 r x r x Σ x 1 x = 0 s x b Parameters: k 1 k k 3 a r s b. 3.3 Trivariate Models Again the Brownian motions are independent. Any correlation between the components is introduced by the off-diagonal terms in the Σ matrix. Some models are specified by the local variance matrix ΣΣ directly ΣΣ is what matters not Σ. 1
13 3.3.1 Model T1 µ x 1 x x 3 = Σ x 1 x x 3 = k 11 m 1 x 1 + k 1 m x k 1 m 1 x 1 + k m x k 33 m 3 x 3 a 1 + b 1 x 1 1/ a + b x 1/ a 3 + b 3 x 3 Parameters: m 1 m m 3 k 11 k 1 k 1 k k 33 a 1 b 1 a b a 3 b Model T µ x 1 x x 3 = Σ x 1 x x 3 Σ x 1 x x 3 = a 1 x + x 3 x 1 a + b x a 3 + b 3 x 3 s 1x 1 s 1 s r 1 x 1 s 1 s 3 r 13 x 1 s s s 3 r 3 s 3 Parameters: a 1 a b a 3 b 3 s 1 s s 3 r 1 r 13 r Model T3 µ x 1 x x 3 = Σ x 1 x x 3 Σ x 1 x x 3 = k 11 g 1 x 1 + k 1 g x + k 13 g 3 x 3 k g x k 33 g 3 x 3 s 11 x 1 s 1 x s 13 x 3 0 s x s 33 x 3 Parameters: g 1 g g 3 k 11 k 1 k 13 k k 33 s 11 s 1 s 13 s s Model T4 b 1 x 1 µ x 1 x x 3 = b x b 3 x 3 Σ x 1 x x 3 = d 11 x d 1 x x 1 d x 0 d 31 x 3 x 1 d 3 x 3 x d 33 x 3 Parameters: b 1 b b 3 d 11 d d 33 d 1 d 31 d 3 13
14 3.3.5 Model T5 µ x 1 x x 3 = x 3 x / v 1 v x r 1 r x 3 x 0 0 ρ Σ x 1 x x 3 = 1 v 3 x 1 ρ 1 v 3 x 0 ρ r 3 ρ 13 r 3 ρ 13ρ 1 3 r 1 ρ 3 1 ρ 13 ρ 3 ρ 13ρ ρ 1 Parameters: v 1 v v 3 r 1 r r 3 ρ 1 ρ 13 ρ Model T6 µ x 1 x x 3 = Σ x 1 x x 3 = a 1 b 1 x 1 a 1 b 1 x 1 + a b x a 31 b 1 x 1 + a 3 b x + a 3 b 3 x 3 x1 0 0 g 1 x1 g x1 0 g 31 x1 g 3 x1 g 33 x1 Parameters: a 1 b 1 a 1 a b a 31 a 3 a 3 b 3 g 1 g g 31 g 3 g Canonical Affine Models Dimensions 1 to Model AY01 µ x 1 = κ 11 x 1 σ x 1 = 1 Pamameters: κ Model AY11 µ x 1 = κ 11 θ 1 x 1 σ x 1 = x 1/ 1 Parameters: θ 1 κ Model AY0 κ 11 x 1 µ x 1 x = κ 1 x 1 κ x 1 0 Σ x 1 x = 0 1 Parameters: κ 11 κ 1 κ. 14
15 3.4.4 Model AY1 κ 11 θ 1 x 1 µ x 1 x = κ 1 θ 1 x 1 κ x x 1/ 1 0 Σ x 1 x = β 1 x 1 1/ Parameters: θ 1 κ 11 κ 1 κ β Model AY κ11 θ 1 x 1 + κ 1 θ x µ x 1 x = κ 1 θ 1 x 1 + κ θ x x 1/ 1 0 Σ x 1 x = 0 x 1/ Parameters: θ 1 θ κ 11 κ 1 κ 1 κ Model AY03 κ 11 x 1 µ x 1 x x 3 = κ 1 x 1 κ x κ 31 x 1 κ 3 x κ 33 x Σ x 1 x x 3 = Parameters: κ 11 κ 1 κ κ 31 κ 3 κ Model AY13 µ x 1 x x 3 = Σ x 1 x x 3 = κ 11 θ 1 x 1 κ 1 θ 1 x 1 κ x κ 3 x 3 κ 31 θ 1 x 1 κ 3 x κ 33 x 3 x 1/ β 1 x 1 1/ β 31 x 1 1/ Parameters: θ 1 κ 11 κ 1 κ κ 3 κ 31 κ 3 κ 33 β 1 β
16 3.4.8 Model AY3 µ x 1 x x 3 = Σ x 1 x x 3 = κ 11 θ 1 x 1 + κ 1 θ x κ 1 θ 1 x 1 + κ θ x κ 31 θ 1 x 1 + κ 3 θ x κ 33 x 3 x 1/ x 1/ β 31 x 1 + β 3 x 1/ Parameters: θ 1 θ κ 11 κ 1 κ 1 κ κ 31 κ 3 κ 33 β 31 β Model AY33 µ x 1 x x 3 = Σ x 1 x x 3 = κ 11 θ 1 x 1 + κ 1 θ x + κ 13 θ 3 x 3 κ 1 θ 1 x 1 + κ θ x + κ 3 θ 3 x 3 κ 31 θ 1 x 1 + κ 3 θ x + κ 33 θ 3 x 3 x 1/ x 1/ x 1/ 3 Parameters: θ 1 θ θ 3 κ 11 κ 1 κ 13 κ 1 κ κ 3 κ 31 κ 3 κ Semi-Affine Models Dimensions 1 to Model SA11 µ x = a + bx + d x σ x = x Parameters: a b d 3.5. Model SA1 a + b 11 x 1 + d 1 x1 µ x 1 x = b 1 x 1 + b x + d 1 + fx1 x1 0 Σ x 1 x = fx1 Parameters: a b 11 d 1 b 1 b d f Model SA a 1 + b 11 x 1 + b 1 x + d 1 x1 µ x 1 x = a + b 1 x 1 + b x + d x x1 0 Σ x 1 x = 0 x Parameters: a 1 b 11 b 1 d 1 a b 1 b d 16
17 3.5.4 Model SA13 a + b 11 x 1 + d 1 x1 µ x 1 x x 3 = b 1 x 1 + b x + b 3 x 3 + d 1 + e1 x 1 b 31 x 1 + b 3 x + b 33 x 3 + d e x 1 x1 0 0 Σ x 1 x x 3 = e1 x e x Parameters: a b 11 d 1 b 1 b b 3 d e 1 b 31 b 3 b 33 d 3 e Model SA3 µ x 1 x x 3 = a 1 + b 11 x 1 + b 1 x + d 1 x1 a + b 1 x 1 + b x + d x b 31 x 1 + b 3 x + b 33 x 3 + d e1 x 1 + e x x1 0 0 Σ x 1 x x 3 = 0 x e1 x 1 + e x Parameters: a 1 b 11 b 1 d 1 a b 1 b d b 31 b 3 b 33 d 3 e 1 e Model SA33 µ x 1 x x 3 = a 1 + b 11 x 1 + b 1 x + b 13 x 3 + d 1 x1 a + b 1 x 1 + b x + b 3 x 3 + d x a 3 + b 31 x 1 + b 3 x + b 33 x 3 + d 3 x3 x1 0 0 Σ x 1 x x 3 = 0 x x3 Parameters: a 1 b 11 b 1 b 13 d 1 a b 1 b b 3 d a 3 b 31 b 3 b 33 d 3. 17
MLEMVD: 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 informationCalibration and Simulation of Interest Rate Models in MATLAB Kevin Shea, CFA Principal Software Engineer MathWorks
Calibration and Simulation of Interest Rate Models in MATLAB Kevin Shea, CFA Principal Software Engineer MathWorks 2014 The MathWorks, Inc. 1 Outline Calibration to Market Data Calibration to Historical
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 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 informationSupplementary Appendix to The Risk Premia Embedded in Index Options
Supplementary Appendix to The Risk Premia Embedded in Index Options Torben G. Andersen Nicola Fusari Viktor Todorov December 214 Contents A The Non-Linear Factor Structure of Option Surfaces 2 B Additional
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 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 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 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 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 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 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 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 information1 Residual life for gamma and Weibull distributions
Supplement to Tail Estimation for Window Censored Processes Residual life for gamma and Weibull distributions. Gamma distribution Let Γ(k, x = x yk e y dy be the upper incomplete gamma function, and let
More informationInt. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session CPS001) p approach
Int. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session CPS001) p.5901 What drives short rate dynamics? approach A functional gradient descent Audrino, Francesco University
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 informationSimulating more interesting stochastic processes
Chapter 7 Simulating more interesting stochastic processes 7. Generating correlated random variables The lectures contained a lot of motivation and pictures. We'll boil everything down to pure algebra
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 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 information1 Roy model: Chiswick (1978) and Borjas (1987)
14.662, Spring 2015: Problem Set 3 Due Wednesday 22 April (before class) Heidi L. Williams TA: Peter Hull 1 Roy model: Chiswick (1978) and Borjas (1987) Chiswick (1978) is interested in estimating regressions
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 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 informationEmpirical Distribution Testing of Economic Scenario Generators
1/27 Empirical Distribution Testing of Economic Scenario Generators Gary Venter University of New South Wales 2/27 STATISTICAL CONCEPTUAL BACKGROUND "All models are wrong but some are useful"; George Box
More informationTwo-dimensional COS method
Two-dimensional COS method Marjon Ruijter Winterschool Lunteren 22 January 2013 1/29 Introduction PhD student since October 2010 Prof.dr.ir. C.W. Oosterlee). CWI national research center for 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 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 informationDuangporn Jearkpaporn, Connie M. Borror Douglas C. Montgomery and George C. Runger Arizona State University Tempe, AZ
Process Monitoring for Correlated Gamma Distributed Data Using Generalized Linear Model Based Control Charts Duangporn Jearkpaporn, Connie M. Borror Douglas C. Montgomery and George C. Runger Arizona State
More informationTEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING
TEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING Semih Yön 1, Cafer Erhan Bozdağ 2 1,2 Department of Industrial Engineering, Istanbul Technical University, Macka Besiktas, 34367 Turkey Abstract.
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 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 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 information13.3 A Stochastic Production Planning Model
13.3. A Stochastic Production Planning Model 347 From (13.9), we can formally write (dx t ) = f (dt) + G (dz t ) + fgdz t dt, (13.3) dx t dt = f(dt) + Gdz t dt. (13.33) The exact meaning of these expressions
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 informationCorresponding author: Gregory C Chow,
Co-movements of Shanghai and New York stock prices by time-varying regressions Gregory C Chow a, Changjiang Liu b, Linlin Niu b,c a Department of Economics, Fisher Hall Princeton University, Princeton,
More informationOnline Appendix to Financing Asset Sales and Business Cycles
Online Appendix to Financing Asset Sales usiness Cycles Marc Arnold Dirk Hackbarth Tatjana Xenia Puhan August 31, 2015 University of St. allen, Rosenbergstrasse 52, 9000 St. allen, Switzerl. Telephone:
More informationRohini Kumar. Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque)
Small time asymptotics for fast mean-reverting stochastic volatility models Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque) March 11, 2011 Frontier Probability Days,
More 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 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 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 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 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 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 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 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 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 information5. Itô Calculus. Partial derivative are abstractions. Usually they are called multipliers or marginal effects (cf. the Greeks in option theory).
5. Itô Calculus Types of derivatives Consider a function F (S t,t) depending on two variables S t (say, price) time t, where variable S t itself varies with time t. In stard calculus there are three types
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 informationMarkowitz portfolio theory
Markowitz portfolio theory Farhad Amu, Marcus Millegård February 9, 2009 1 Introduction Optimizing a portfolio is a major area in nance. The objective is to maximize the yield and simultaneously minimize
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 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 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 informationThe Two Sample T-test with One Variance Unknown
The Two Sample T-test with One Variance Unknown Arnab Maity Department of Statistics, Texas A&M University, College Station TX 77843-343, U.S.A. amaity@stat.tamu.edu Michael Sherman Department of Statistics,
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 informationEconomathematics. Problem Sheet 1. Zbigniew Palmowski. Ws 2 dw s = 1 t
Economathematics Problem Sheet 1 Zbigniew Palmowski 1. Calculate Ee X where X is a gaussian random variable with mean µ and volatility σ >.. Verify that where W is a Wiener process. Ws dw s = 1 3 W t 3
More 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 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 informationKeywords: China; Globalization; Rate of Return; Stock Markets; Time-varying parameter regression.
Co-movements of Shanghai and New York Stock prices by time-varying regressions Gregory C Chow a, Changjiang Liu b, Linlin Niu b,c a Department of Economics, Fisher Hall Princeton University, Princeton,
More informationLecture 4. Finite difference and finite element methods
Finite difference and finite element methods Lecture 4 Outline Black-Scholes equation From expectation to PDE Goal: compute the value of European option with payoff g which is the conditional expectation
More 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 informationPremia 14 HESTON MODEL CALIBRATION USING VARIANCE SWAPS PRICES
Premia 14 HESTON MODEL CALIBRATION USING VARIANCE SWAPS PRICES VADIM ZHERDER Premia Team INRIA E-mail: vzherder@mailru 1 Heston model Let the asset price process S t follows the Heston stochastic volatility
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 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 informationSimulation Wrap-up, Statistics COS 323
Simulation Wrap-up, Statistics COS 323 Today Simulation Re-cap Statistics Variance and confidence intervals for simulations Simulation wrap-up FYI: No class or office hours Thursday Simulation wrap-up
More informationBACHELIER FINANCE SOCIETY. 4 th World Congress Tokyo, Investments and forward utilities. Thaleia Zariphopoulou The University of Texas at Austin
BACHELIER FINANCE SOCIETY 4 th World Congress Tokyo, 26 Investments and forward utilities Thaleia Zariphopoulou The University of Texas at Austin 1 Topics Utility-based measurement of performance Utilities
More informationSPDE and portfolio choice (joint work with M. Musiela) Princeton University. Thaleia Zariphopoulou The University of Texas at Austin
SPDE and portfolio choice (joint work with M. Musiela) Princeton University November 2007 Thaleia Zariphopoulou The University of Texas at Austin 1 Performance measurement of investment strategies 2 Market
More 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 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 informationEmpirical Analysis of the US Swap Curve Gough, O., Juneja, J.A., Nowman, K.B. and Van Dellen, S.
WestminsterResearch http://www.westminster.ac.uk/westminsterresearch Empirical Analysis of the US Swap Curve Gough, O., Juneja, J.A., Nowman, K.B. and Van Dellen, S. This is a copy of the final version
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 informationWITH SKETCH ANSWERS. Postgraduate Certificate in Finance Postgraduate Certificate in Economics and Finance
WITH SKETCH ANSWERS BIRKBECK COLLEGE (University of London) BIRKBECK COLLEGE (University of London) Postgraduate Certificate in Finance Postgraduate Certificate in Economics and Finance SCHOOL OF ECONOMICS,
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 informationLecture Note 9 of Bus 41914, Spring Multivariate Volatility Models ChicagoBooth
Lecture Note 9 of Bus 41914, Spring 2017. Multivariate Volatility Models ChicagoBooth Reference: Chapter 7 of the textbook Estimation: use the MTS package with commands: EWMAvol, marchtest, BEKK11, dccpre,
More informationHomework Assignments for BusAdm 713: Business Forecasting Methods. Assignment 1: Introduction to forecasting, Review of regression
Homework Assignments for BusAdm 713: Business Forecasting Methods Note: Problem points are in parentheses. Assignment 1: Introduction to forecasting, Review of regression 1. (3) Complete the exercises
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 informationEnergy Systems under Uncertainty: Modeling and Computations
Energy Systems under Uncertainty: Modeling and Computations W. Römisch Humboldt-University Berlin Department of Mathematics www.math.hu-berlin.de/~romisch Systems Analysis 2015, November 11 13, IIASA (Laxenburg,
More informationUsing Monte Carlo Integration and Control Variates to Estimate π
Using Monte Carlo Integration and Control Variates to Estimate π N. Cannady, P. Faciane, D. Miksa LSU July 9, 2009 Abstract We will demonstrate the utility of Monte Carlo integration by using this algorithm
More informationPricing Barrier Options under Local Volatility
Abstract Pricing Barrier Options under Local Volatility Artur Sepp Mail: artursepp@hotmail.com, Web: www.hot.ee/seppar 16 November 2002 We study pricing under the local volatility. Our research is mainly
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 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 informationAnalyzing Oil Futures with a Dynamic Nelson-Siegel Model
Analyzing Oil Futures with a Dynamic Nelson-Siegel Model NIELS STRANGE HANSEN & ASGER LUNDE DEPARTMENT OF ECONOMICS AND BUSINESS, BUSINESS AND SOCIAL SCIENCES, AARHUS UNIVERSITY AND CENTER FOR RESEARCH
More informationProxy Function Fitting: Some Implementation Topics
OCTOBER 2013 ENTERPRISE RISK SOLUTIONS RESEARCH OCTOBER 2013 Proxy Function Fitting: Some Implementation Topics Gavin Conn FFA Moody's Analytics Research Contact Us Americas +1.212.553.1658 clientservices@moodys.com
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 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 informationWeb-based Supplementary Materials for. A space-time conditional intensity model. for invasive meningococcal disease occurence
Web-based Supplementary Materials for A space-time conditional intensity model for invasive meningococcal disease occurence by Sebastian Meyer 1,2, Johannes Elias 3, and Michael Höhle 4,2 1 Department
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 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 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 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 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 informationOptimal asset allocation under forward performance criteria Oberwolfach, February 2007
Optimal asset allocation under forward performance criteria Oberwolfach, February 2007 Thaleia Zariphopoulou The University of Texas at Austin 1 References Indifference valuation in binomial models (with
More informationAn Implementation of Markov Regime Switching GARCH Models in Matlab
An Implementation of Markov Regime Switching GARCH Models in Matlab Thomas Chuffart Aix-Marseille University (Aix-Marseille School of Economics), CNRS & EHESS Abstract MSGtool is a MATLAB toolbox which
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 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 informationA No-Arbitrage Theorem for Uncertain Stock Model
Fuzzy Optim Decis Making manuscript No (will be inserted by the editor) A No-Arbitrage Theorem for Uncertain Stock Model Kai Yao Received: date / Accepted: date Abstract Stock model is used to describe
More 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 informationOptimal investments under dynamic performance critria. Lecture IV
Optimal investments under dynamic performance critria Lecture IV 1 Utility-based measurement of performance 2 Deterministic environment Utility traits u(x, t) : x wealth and t time Monotonicity u x (x,
More informationA discretionary stopping problem with applications to the optimal timing of investment decisions.
A discretionary stopping problem with applications to the optimal timing of investment decisions. Timothy Johnson Department of Mathematics King s College London The Strand London WC2R 2LS, UK Tuesday,
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 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 informationWindow Width Selection for L 2 Adjusted Quantile Regression
Window Width Selection for L 2 Adjusted Quantile Regression Yoonsuh Jung, The Ohio State University Steven N. MacEachern, The Ohio State University Yoonkyung Lee, The Ohio State University Technical Report
More information