TMS165/MSA350 Stochastic Calculus, Lecture on Applications
|
|
- Aubrey McCoy
- 5 years ago
- Views:
Transcription
1 TMS165/MSA35 Stochastic Calculus, Lecture on Applications In this lecture we demonstrate how statistical methods such as the maximum likelihood method likelihood ratio estimation can be applied to the Ornstein-Uhlenbeck (OU) process. 1. Elements of diffusion theory Diffusion processes SDE A time homogeneous diffusion process is the solution X = {X(t) t to an SDE of the form dx(t) = µ(x(t))dt+σ(x(t))db(t) for t, (1) where the drift coefficient µ : the diffusion coefficient σ : are sufficiently nice functions. Here B = {B(t) t denotes a Brownian motion as usual. By definition, the solution X to (1) satisfies X(t) = X()+ Markov property t µ(x(r))dr + t σ(x(r))db(r) for t. () The solution X to the SDE (1) is a Markov process, which is to say that P { X(t) F X s = P{X(t) X(s) for s t. Here Fs X = σ{x(r):r s for s is the filtration generated by the process X itself. While a detailed rigorous proof of the Markov property is extremly complicated, the Markov property is easy to underst from a more heursitic point of view: Using the representation () for both X(t) X(s) we get X(t) = X(s)+ t s = X(s)+lim µ(x(r))dr + t s σ(x(r)) db(r) n µ(x(t i 1 ))(t i t i 1 )+lim n σ(x(t i 1 ))(B(t i ) B(t i 1 )), where s = t < t 1 <... < t n = t is a partion of the interval [s,t] that becomes finer finer in the limit. From this we see the that the only thing from the past Fs X that affects the future value X(t) is the value X(s)=X(t ) of X at time s=t. 1
2 Transition densities The transition density function p(t,x,y) = d P{X(t+s) y X(s)=x for t> dy of the diffusion process X given by (1) or () satisfies the Kolmogorov backward PDE σ(x) p(t,x,y) = t x p(t,x,y)+µ(x) x p(t,x,y). Conversely, under general conditions, a solution to this PDE is the transition density function of the diffusion process given by (1) or () if it is a probability density function as a function of y for each choice of x t>, that is, p(t,x,y) p(t,x,y)dy = 1, in addition satisfies p(t,x,y) as t for x y. In general it is not an easy task or even possible to find an explicit expression for the transition density function. Argubly, the most common systematic way to try to solve the Kolmogorov backward PDE is to consider the Laplace transform ˆp(λ,x,y) = which must satisfy the ODE (for x y), as λ ˆp(λ,x,y) = σ(x) e λt p(t,x,y)dt for λ>, x ˆp(λ,x,y)+µ(x) x ˆp(λ,x,y) e λt t p(t,x,y)dt = λ e λt p(t,x,y)dt (for x y) when p(,x,y)=. The conditions that p is a density function translates to ˆp(λ,x,y)dy = e λt p(t,x,y)dtdy = e λt p(t,x,y)dydt = e λt dt = 1 λ. As the ODE for ˆp(λ,x,y) usually has a unique solution that integrates to 1/λ, in the above fashion, this determines the Laplace transform ˆp(λ, x, y) of p(t, x, y), after which p(t, x, y) is found by inverse Laplace tranformation. Although the above indicated way to find transition densities for diffusions does indeed work for many important equations, the details of the solution are in general too difficult to be attempted on undergraduate level, it is usually more rewarding to search the literature (web) for solutions than to try to derive them oneself.
3 Stationary distribution A stationary density function for X is a probability density function π that satisfies d π(y) = dy P{X(t+s) y X(s)=xπ(x)dx = p(t,x,y)π(x)dx for t>. By the Chapman-Kolmogorov equation below Theorem 5.6 in Klebaner s book, this means that if the process has the stationary density function (distribution) at a certain time, then it has the stationary density function (distribution) at all later times. It can be shown (see Task of Home exercise session 5) that if X is started according to the stationary distribution at time t =, then X is a stationary process. Also, regardless of how X is started, X will converge (in a way we select not to specify) to a stationary process with the stationary marginal distribution as t. As the transition density satisfies the Kolmogorov forward PDE ( ) σ(y) p(t,x,y) = t y p(t,x,y) ( ) µ(y)p(t,x,y), y we see from the above integral equation that ( σ(y) y ( )π(y) y µ(y) = t σ(y) y + )π(y) y µ(y) ( = t σ(y) y + )p(t,x,y)π(x)dx y µ(y) =. If it exists, then the stationary density function is given by π(y) = C { y σ(y) exp µ(z) y σ(z) dz. Here C > is a normalizing constant selected to make π a probability density function, that is, π(y)dy = 1, while y is any constant. The existence of the stationary density function is the same thing as that this normalizing procedure is possible to carry out. Note that it is easy to see that π satisfies the ODE (cf. above) ( σ(y) y )π(y). y µ(y) Finite dimensional distributions The joint probability density function of (X(t ),X(t 1 ),...,X(t n )) for = t < t 1 <... <t n is given by f X(t ),X(t 1 ),...,X(t n)(x,x 1,...,x n ) = π(x ) 3 n p(t i t i 1,x i 1,x i )
4 when the process is started with the stationary density function π (provided that it exists) at time, while n f X(t1 ),...,X(t n)(x 1,...,x n ) = p(t i t i 1,x i 1,x i ) when the process is started at a fixed value X()=x. Euler method We may simulate a weak solution to the SDE at a time grid =t <t 1 <... <t n =T by means of the Euler method, as X(t i ) X(t i 1 )+µ(x(t i 1 ))(t i t i 1 )+σ(x(t i 1 ))ξ i for i = 1,...,n, where {ξ i n are independent rom variables such that ξ i is zero-mean normal distributed with starddeviation t i t i 1, where X() is a rom variable that is indepedent of {ξ i n has the stationary distribution if X is started according to the stationary distribution, while X()=x if X is started at a fixed value x. Likelihood ratios If the diffusion process X satisfies the equation for a P 1 -Brownian motion B dx(t) = µ 1 (X(t))dt+σ(X(t))dB(t) dx(t) = µ (X(t))dt+σ(X(t))dW(t) for a P -Brownian motion W, where P 1 P are two different probability measures, then the likelihood ratio between P P 1 based on an observation {X(t) t [,T] of the process X in the time interval [,T] is given by { dp T µ (X(t)) µ 1 (X(t)) = exp dp 1 σ(x(t)) dx(t) 1 T µ (X(t)) µ 1 (X(t)) σ(x(t)) dt. This likelihood ratio can be used to judge which is the most likely of the above two SDE-models for X based on the observations: If dp /dp 1 is bigger than 1 then the model with the drift µ is the most likely, while a value of dp /dp 1 that is smaller than 1 indicates that the drift µ 1 is the most likely. The likelihood ratio can also be used to estimate parameters of a parametric SDE, as we will see later. 4
5 . The OU process An OU process is the solution X to the SDE dx(t) = µx(t)dt+σdb(t), where µ> σ> are parameters. In other words, the drift is µ(x) = µx the difffusion coefficient σ(x) = σ. The OU process has stationary distribution π(y) = 1 y /( { { σ exp µz σ dz 1 σ exp y ) µz σ dz dy = µ πσ exp that is, a zero-mean normal distribution with starddeviation σ/ µ. { µy σ, As we have shown in Exercise Session 5, the OU process has transition density p(t,x,y) = exp (y e µt x) πσ 1 e µt σ (1 e µt, )/µ thatis,anormaldistributionwithmeane µt xstarddeviationσ 1 e µt / µ. It is also a routine matter to differentiate to check that this function p satisfies the Kolmogorov backward PDE. If X is an OU process dx(t) = µ 1 X(t)dt+σdB(t) for a P µ1 -Brownian motion B, an OU process dx(t) = µ X(t)dt+σdW(t) for a P µ -Brownian motion W, then the corresponding likelihood ratio is given by { dp µ = exp µ T T µ 1 dp µ1 σ X(t)dX(t) µ µ 1 σ X(t) dt. In particular, we can find which is the most likely of the models dx(t) = db(t) dx(t) = X(t)dt+dW(t) by computing the likelihood ratio for µ 1 =, µ =1 σ=1 { dp T 1 = exp X(t)dX(t) 1 T X(t) dt, dp 5
6 then check whether dp 1 /dp > 1, indicating that µ = 1 is the most appropriate model, or dp 1 /dp < 1, indicating that µ= is the most appropriate model. We can also estimate the parameter µ for the equation dx(t) = µx(t)dt+dw(t) by means of maximizing the likelihood { dp µ = exp µ dp T which by differentiation gives the estimate µ = T 3. Application to the OU process X(t)dX(t) µ T /( T ) X(t) dx(t) X(t) dt. X(t) dt, We used the Euler method to simulate an OU process {X(t) t [,1] started according to the stationary distribution, an OU process {Y(t) t [,1] started at zero. In both cases the drift was µ(x) = µ x the diffusion coefficient σ(x) = σ, where µ = σ = 1. We use distance 1 1 between the time points of the simulation grid, so that =t < t 1 <... <t 1 =1, where t i t i 1 = 1 1 for i = 1,...,1. The simulations were carried out by means of the following Mathematica programs. <<Statistics ContinuousDistributions ; dt=1/1; T=1; {mu,sigma={1,1; For[i=; X={om[NormalDistribution[,sigma/Sqrt[*mu]]], i<=t/dt, i++, AppendTo[X, X[[i-1]] - mu*x[[i-1]]*dt + om[normaldistribution[,sigma*sqrt[dt]]]]] For[i=; Y=, i<=t/dt, i++, AppendTo[Y, Y[[i-1]] - mu*y[[i-1]]*dt + om[normaldistribution[,sigma*sqrt[dt]]]]] The results of the simulations are depicted in the following two figures 6
7 Y t OU process dy t Y t dt db t, Y t The joint density functions of (X(t ),...,X(t 1 )) (Y(t ),...,Y(t 1 )) are given by f X(t ),...,X(t 1 )(x,...,x 1 ) 1 = π(x ) p(t i t i 1,x i 1,x i ) = µ πσ exp { 1 µx σ f Y(t1 ),...,Y(t 1 )(y 1,...,y 1 ) = = 1 1 exp (x i e µ/1 x i 1 ) πσ 1 e µ/1 σ (1 e µ/1 )/µ p(t i t i 1,y i 1,y i ) exp (y i e µ/1 y i 1 ) πσ 1 e µ/1 σ (1 e µ/1, )/µ respectively, where y =. Hence we may use the maximum likelihood method to estimate µ σ from our simulated data (pretending that they are unknown), by 7
8 means of maximizing the likelihood µ πσ exp { µx(t ) 1 σ 1 exp (X(t i) e µ/1 X(t i 1 )) πσ 1 e µ/1 σ (1 e µ/1 )/µ exp (Y(t i) e µ/1 Y(t i 1 )) πσ 1 e µ/1 σ (1 e µ/1, )/µ respectively. These maximum likelihood estimates were carried out by means of the following Mathematcia code (with the densities logged to not get numerical underflows). foustationary[mu,sigma,x ] := Sqrt[mu]*Exp[-mu*x^/sigma^]/(Sqrt[Pi]*sigma); pou[mu,sigma,x,y,t ] := Exp[-(y-x*Exp[-mu*t])^/(*sigma^*(1-Exp[-*mu*t])/(*mu))] /(Sqrt[*Pi]*sigma*Sqrt[1-Exp[-*mu*t]]/Sqrt[*mu]); MLStationary[mu,sigma,dt,Data ] := Log[fOUStationary[mu,sigma,Data[[1]]]] + Sum[Log[pOU[mu,sigma,Data[[i-1]],Data[[i]],dt]], {i,,length[data]] MLNonStationary[mu,sigma,dt,Data ] := Sum[Log[pOU[mu,sigma,Data[[i-1]],Data[[i]],dt]], {i,,length[data]] NMaximize[MLStationary[mu,sigma,dt,X],mu>,sigma>, mu,sigma] Out[]:= {89.45, {mu -> , sigma -> NMaximize[MLNonStationary[mu,sigma,dt,Y],mu>,sigma>, mu,sigma] Out[]:= {886.86, {mu ->.99798, sigma -> 1.88 Note how well this fits with the correct values µ=1 σ=1 for the parameters. We calculate the likelihood ratios { dp 1 (X) = exp dp 1 X(t)dX(t) 1 1 X(t) dt { dp 1 (Y) = exp dp 1 Y(t)dY(t) 1 1 Y(t) dt, 8
9 for our simulated processes, in order to fine whether dx(t) = db(t) or dx(t) = X(t)dt+dW(t) dy(t) = db(t) or dy(t) = Y(t)dt+dW(t), respectively, are the most likely models for the data. As both ratios where significantly larger than 1 (see the enclosed Mathematica code), the model with µ=1 was the most likely for both data sets. OUatioTest[Data ] := Exp[Sum[-Data[[i-1]]*(Data[[i]]-Data[[i-1]]), {i,, Length[Data]] - Sum[Data[[i]]^*dt, {i,1,length[data]]/]; {OUatioTest[X], OUatioTest[Y] Out[]:= {1.134, We may estimate the parameter µ by means of maximizing the likelihood ratios (dp µ /dp )(X) (dp µ /dp )(Y), respectively, which gives the µ estimates µ = 1 /( 1 ) X(t) dx(t) X(t) dt µ = 1 /( 1 ) Y(t)dY(t) Y(t) dt, respectively. Both results were very close to the correct µ=1, as the following Mathematica code illustrates: OUatioEst[Data ] := -Sum[Data[[i-1]]*(Data[[i]]-Data[[i-1]]), {i,,length[data]] /Sum[Data[[i]]^*dt, {i,1,length[data]]; {OUatioEst[X], OUatioEst[Y] { ,
Kolmogorov Equations and Markov Processes
Kolmogorov Equations and Markov Processes May 3, 013 1 Transition measures and functions Consider a stochastic process {X(t)} t 0 whose state space is a product of intervals contained in R n. We define
More informationLecture 21: Stochastic Differential Equations
: Stochastic Differential Equations In this lecture, we study stochastic differential equations. See Chapter 9 of [3] for a thorough treatment of the materials in this section. 1. Stochastic differential
More informationp 1 ( Y p dp) 1/p ( X p dp) 1 1 p
Doob s inequality Let X(t) be a right continuous submartingale with respect to F(t), t 1 P(sup s t X(s) λ) 1 λ {sup s t X(s) λ} X + (t)dp 2 For 1 < p
More informationLecture 12: Diffusion Processes and Stochastic Differential Equations
Lecture 12: Diffusion Processes and Stochastic Differential Equations 1. Diffusion Processes 1.1 Definition of a diffusion process 1.2 Examples 2. Stochastic Differential Equations SDE) 2.1 Stochastic
More informationIntroduction to Diffusion Processes.
Introduction to Diffusion Processes. Arka P. Ghosh Department of Statistics Iowa State University Ames, IA 511-121 apghosh@iastate.edu (515) 294-7851. February 1, 21 Abstract In this section we describe
More informationI forgot to mention last time: in the Ito formula for two standard processes, putting
I forgot to mention last time: in the Ito formula for two standard processes, putting dx t = a t dt + b t db t dy t = α t dt + β t db t, and taking f(x, y = xy, one has f x = y, f y = x, and f xx = f yy
More informationHomogenization with stochastic differential equations
Homogenization with stochastic differential equations Scott Hottovy shottovy@math.arizona.edu University of Arizona Program in Applied Mathematics October 12, 2011 Modeling with SDE Use SDE to model system
More informationStochastic Calculus and Black-Scholes Theory MTH772P Exercises Sheet 1
Stochastic Calculus and Black-Scholes Theory MTH772P Exercises Sheet. For ξ, ξ 2, i.i.d. with P(ξ i = ± = /2 define the discrete-time random walk W =, W n = ξ +... + ξ n. (i Formulate and prove the property
More informationIntroduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0
Introduction to multiscale modeling and simulation Lecture 5 Numerical methods for ODEs, SDEs and PDEs The need for multiscale methods Two generic frameworks for multiscale computation Explicit methods
More informationSloppy derivations of Ito s formula and the Fokker-Planck equations
Sloppy derivations of Ito s formula and the Fokker-Planck equations P. G. Harrison Department of Computing, Imperial College London South Kensington Campus, London SW7 AZ, UK email: pgh@doc.ic.ac.uk April
More informationSTATISTICS 385: STOCHASTIC CALCULUS HOMEWORK ASSIGNMENT 4 DUE NOVEMBER 23, = (2n 1)(2n 3) 3 1.
STATISTICS 385: STOCHASTIC CALCULUS HOMEWORK ASSIGNMENT 4 DUE NOVEMBER 23, 26 Problem Normal Moments (A) Use the Itô formula and Brownian scaling to check that the even moments of the normal distribution
More informationLecture 6: Bayesian Inference in SDE Models
Lecture 6: Bayesian Inference in SDE Models Bayesian Filtering and Smoothing Point of View Simo Särkkä Aalto University Simo Särkkä (Aalto) Lecture 6: Bayesian Inference in SDEs 1 / 45 Contents 1 SDEs
More informationIntroduction to numerical simulations for Stochastic ODEs
Introduction to numerical simulations for Stochastic ODEs Xingye Kan Illinois Institute of Technology Department of Applied Mathematics Chicago, IL 60616 August 9, 2010 Outline 1 Preliminaries 2 Numerical
More informationThe multidimensional Ito Integral and the multidimensional Ito Formula. Eric Mu ller June 1, 2015 Seminar on Stochastic Geometry and its applications
The multidimensional Ito Integral and the multidimensional Ito Formula Eric Mu ller June 1, 215 Seminar on Stochastic Geometry and its applications page 2 Seminar on Stochastic Geometry and its applications
More informationLecture 1: Pragmatic Introduction to Stochastic Differential Equations
Lecture 1: Pragmatic Introduction to Stochastic Differential Equations Simo Särkkä Aalto University, Finland (visiting at Oxford University, UK) November 13, 2013 Simo Särkkä (Aalto) Lecture 1: Pragmatic
More informationStochastic differential equation models in biology Susanne Ditlevsen
Stochastic differential equation models in biology Susanne Ditlevsen Introduction This chapter is concerned with continuous time processes, which are often modeled as a system of ordinary differential
More information{σ x >t}p x. (σ x >t)=e at.
3.11. EXERCISES 121 3.11 Exercises Exercise 3.1 Consider the Ornstein Uhlenbeck process in example 3.1.7(B). Show that the defined process is a Markov process which converges in distribution to an N(0,σ
More informationStochastic Calculus February 11, / 33
Martingale Transform M n martingale with respect to F n, n =, 1, 2,... σ n F n (σ M) n = n 1 i= σ i(m i+1 M i ) is a Martingale E[(σ M) n F n 1 ] n 1 = E[ σ i (M i+1 M i ) F n 1 ] i= n 2 = σ i (M i+1 M
More informationLAN property for ergodic jump-diffusion processes with discrete observations
LAN property for ergodic jump-diffusion processes with discrete observations Eulalia Nualart (Universitat Pompeu Fabra, Barcelona) joint work with Arturo Kohatsu-Higa (Ritsumeikan University, Japan) &
More informationConvoluted Brownian motions: a class of remarkable Gaussian processes
Convoluted Brownian motions: a class of remarkable Gaussian processes Sylvie Roelly Random models with applications in the natural sciences Bogotá, December 11-15, 217 S. Roelly (Universität Potsdam) 1
More informationResearch Article A Necessary Characteristic Equation of Diffusion Processes Having Gaussian Marginals
Abstract and Applied Analysis Volume 01, Article ID 598590, 9 pages doi:10.1155/01/598590 Research Article A Necessary Characteristic Equation of Diffusion Processes Having Gaussian Marginals Syeda Rabab
More informationThe concentration of a drug in blood. Exponential decay. Different realizations. Exponential decay with noise. dc(t) dt.
The concentration of a drug in blood Exponential decay C12 concentration 2 4 6 8 1 C12 concentration 2 4 6 8 1 dc(t) dt = µc(t) C(t) = C()e µt 2 4 6 8 1 12 time in minutes 2 4 6 8 1 12 time in minutes
More informationModel Specification Testing in Nonparametric and Semiparametric Time Series Econometrics. Jiti Gao
Model Specification Testing in Nonparametric and Semiparametric Time Series Econometrics Jiti Gao Department of Statistics School of Mathematics and Statistics The University of Western Australia Crawley
More informationHandbook of Stochastic Methods
C. W. Gardiner Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences Third Edition With 30 Figures Springer Contents 1. A Historical Introduction 1 1.1 Motivation I 1.2 Some Historical
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 7 9/25/2013
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.65/15.070J Fall 013 Lecture 7 9/5/013 The Reflection Principle. The Distribution of the Maximum. Brownian motion with drift Content. 1. Quick intro to stopping times.
More informationMathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca November 26th, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
More informationSDE Coefficients. March 4, 2008
SDE Coefficients March 4, 2008 The following is a summary of GARD sections 3.3 and 6., mainly as an overview of the two main approaches to creating a SDE model. Stochastic Differential Equations (SDE)
More informationLecture on Parameter Estimation for Stochastic Differential Equations. Erik Lindström
Lecture on Parameter Estimation for Stochastic Differential Equations Erik Lindström Recap We are interested in the parameters θ in the Stochastic Integral Equations X(t) = X(0) + t 0 µ θ (s, X(s))ds +
More informationLecture 4: Introduction to stochastic processes and stochastic calculus
Lecture 4: Introduction to stochastic processes and stochastic calculus Cédric Archambeau Centre for Computational Statistics and Machine Learning Department of Computer Science University College London
More informationBernardo D Auria Stochastic Processes /10. Notes. Abril 13 th, 2010
1 Stochastic Calculus Notes Abril 13 th, 1 As we have seen in previous lessons, the stochastic integral with respect to the Brownian motion shows a behavior different from the classical Riemann-Stieltjes
More informationLecture 4: Ito s Stochastic Calculus and SDE. Seung Yeal Ha Dept of Mathematical Sciences Seoul National University
Lecture 4: Ito s Stochastic Calculus and SDE Seung Yeal Ha Dept of Mathematical Sciences Seoul National University 1 Preliminaries What is Calculus? Integral, Differentiation. Differentiation 2 Integral
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationApproximating diffusions by piecewise constant parameters
Approximating diffusions by piecewise constant parameters Lothar Breuer Institute of Mathematics Statistics, University of Kent, Canterbury CT2 7NF, UK Abstract We approximate the resolvent of a one-dimensional
More information6. Brownian Motion. Q(A) = P [ ω : x(, ω) A )
6. Brownian Motion. stochastic process can be thought of in one of many equivalent ways. We can begin with an underlying probability space (Ω, Σ, P) and a real valued stochastic process can be defined
More informationExercises. T 2T. e ita φ(t)dt.
Exercises. Set #. Construct an example of a sequence of probability measures P n on R which converge weakly to a probability measure P but so that the first moments m,n = xdp n do not converge to m = xdp.
More informationThe Smoluchowski-Kramers Approximation: What model describes a Brownian particle?
The Smoluchowski-Kramers Approximation: What model describes a Brownian particle? Scott Hottovy shottovy@math.arizona.edu University of Arizona Applied Mathematics October 7, 2011 Brown observes a particle
More informationBackward Stochastic Differential Equations with Infinite Time Horizon
Backward Stochastic Differential Equations with Infinite Time Horizon Holger Metzler PhD advisor: Prof. G. Tessitore Università di Milano-Bicocca Spring School Stochastic Control in Finance Roscoff, March
More informationStochastic Integration.
Chapter Stochastic Integration..1 Brownian Motion as a Martingale P is the Wiener measure on (Ω, B) where Ω = C, T B is the Borel σ-field on Ω. In addition we denote by B t the σ-field generated by x(s)
More informationPoisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation
Poisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation Jingyi Zhu Department of Mathematics University of Utah zhu@math.utah.edu Collaborator: Marco Avellaneda (Courant
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationStochastic Differential Equations.
Chapter 3 Stochastic Differential Equations. 3.1 Existence and Uniqueness. One of the ways of constructing a Diffusion process is to solve the stochastic differential equation dx(t) = σ(t, x(t)) dβ(t)
More informationLecture 12: Detailed balance and Eigenfunction methods
Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility), 4.2-4.4 (explicit examples of eigenfunction methods) Gardiner
More informationA Short Introduction to Diffusion Processes and Ito Calculus
A Short Introduction to Diffusion Processes and Ito Calculus Cédric Archambeau University College, London Center for Computational Statistics and Machine Learning c.archambeau@cs.ucl.ac.uk January 24,
More informationM4A42 APPLIED STOCHASTIC PROCESSES
M4A42 APPLIED STOCHASTIC PROCESSES G.A. Pavliotis Department of Mathematics Imperial College London, UK LECTURE 1 12/10/2009 Lectures: Mondays 09:00-11:00, Huxley 139, Tuesdays 09:00-10:00, Huxley 144.
More informationTheoretical Tutorial Session 2
1 / 36 Theoretical Tutorial Session 2 Xiaoming Song Department of Mathematics Drexel University July 27, 216 Outline 2 / 36 Itô s formula Martingale representation theorem Stochastic differential equations
More informationBernardo D Auria Stochastic Processes /12. Notes. March 29 th, 2012
1 Stochastic Calculus Notes March 9 th, 1 In 19, Bachelier proposed for the Paris stock exchange a model for the fluctuations affecting the price X(t) of an asset that was given by the Brownian motion.
More informationVäxjö University. Some recent simulation techniques for diffusion bridges. School of Mathematics and System Engineering
School of Mathematics and System Engineering Reports from MSI - Rapporter från MSI Växjö University Some recent simulation techniques for diffusion bridges Yadigar Sekerci May 2009 MSI Växjö University
More informationIntroduction. Stochastic Processes. Will Penny. Stochastic Differential Equations. Stochastic Chain Rule. Expectations.
19th May 2011 Chain Introduction We will Show the relation between stochastic differential equations, Gaussian processes and methods This gives us a formal way of deriving equations for the activity of
More informationAffine Processes. Econometric specifications. Eduardo Rossi. University of Pavia. March 17, 2009
Affine Processes Econometric specifications Eduardo Rossi University of Pavia March 17, 2009 Eduardo Rossi (University of Pavia) Affine Processes March 17, 2009 1 / 40 Outline 1 Affine Processes 2 Affine
More informationStochastic Differential Equations
Chapter 5 Stochastic Differential Equations We would like to introduce stochastic ODE s without going first through the machinery of stochastic integrals. 5.1 Itô Integrals and Itô Differential Equations
More informationEstimation of arrival and service rates for M/M/c queue system
Estimation of arrival and service rates for M/M/c queue system Katarína Starinská starinskak@gmail.com Charles University Faculty of Mathematics and Physics Department of Probability and Mathematical Statistics
More informationBrownian Motion and Stochastic Calculus
ETHZ, Spring 17 D-MATH Prof Dr Martin Larsson Coordinator A Sepúlveda Brownian Motion and Stochastic Calculus Exercise sheet 6 Please hand in your solutions during exercise class or in your assistant s
More information2 Quick Review of Continuous Random Variables
CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems CDS b R. M. Murray Stochastic Systems 8 January 27 Reading: This set of lectures provides a brief introduction to stochastic systems. Friedland,
More informationIntroduction to Probability and Stocastic Processes - Part I
Introduction to Probability and Stocastic Processes - Part I Lecture 2 Henrik Vie Christensen vie@control.auc.dk Department of Control Engineering Institute of Electronic Systems Aalborg University Denmark
More informationSolution for Problem 7.1. We argue by contradiction. If the limit were not infinite, then since τ M (ω) is nondecreasing we would have
362 Problem Hints and Solutions sup g n (ω, t) g(ω, t) sup g(ω, s) g(ω, t) µ n (ω). t T s,t: s t 1/n By the uniform continuity of t g(ω, t) on [, T], one has for each ω that µ n (ω) as n. Two applications
More informationWeek 9 Generators, duality, change of measure
Week 9 Generators, duality, change of measure Jonathan Goodman November 18, 013 1 Generators This section describes a common abstract way to describe many of the differential equations related to Markov
More informationBrownian Motion. An Undergraduate Introduction to Financial Mathematics. J. Robert Buchanan. J. Robert Buchanan Brownian Motion
Brownian Motion An Undergraduate Introduction to Financial Mathematics J. Robert Buchanan 2010 Background We have already seen that the limiting behavior of a discrete random walk yields a derivation of
More informationThe Cameron-Martin-Girsanov (CMG) Theorem
The Cameron-Martin-Girsanov (CMG) Theorem There are many versions of the CMG Theorem. In some sense, there are many CMG Theorems. The first version appeared in ] in 944. Here we present a standard version,
More informationInterest Rate Models:
1/17 Interest Rate Models: from Parametric Statistics to Infinite Dimensional Stochastic Analysis René Carmona Bendheim Center for Finance ORFE & PACM, Princeton University email: rcarmna@princeton.edu
More informationShort-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility
Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility José Enrique Figueroa-López 1 1 Department of Statistics Purdue University Statistics, Jump Processes,
More informationStochastic optimal control with rough paths
Stochastic optimal control with rough paths Paul Gassiat TU Berlin Stochastic processes and their statistics in Finance, Okinawa, October 28, 2013 Joint work with Joscha Diehl and Peter Friz Introduction
More informationThis is a Gaussian probability centered around m = 0 (the most probable and mean position is the origin) and the mean square displacement m 2 = n,or
Physics 7b: Statistical Mechanics Brownian Motion Brownian motion is the motion of a particle due to the buffeting by the molecules in a gas or liquid. The particle must be small enough that the effects
More informationQuantization of stochastic processes with applications on Euler-Maruyama schemes
UPTEC F15 063 Examensarbete 30 hp Oktober 2015 Quantization of stochastic processes with applications on Euler-Maruyama schemes Viktor Edward Abstract Quantization of stochastic processes with applications
More informationA new approach for investment performance measurement. 3rd WCMF, Santa Barbara November 2009
A new approach for investment performance measurement 3rd WCMF, Santa Barbara November 2009 Thaleia Zariphopoulou University of Oxford, Oxford-Man Institute and The University of Texas at Austin 1 Performance
More informationMA8109 Stochastic Processes in Systems Theory Autumn 2013
Norwegian University of Science and Technology Department of Mathematical Sciences MA819 Stochastic Processes in Systems Theory Autumn 213 1 MA819 Exam 23, problem 3b This is a linear equation of the form
More informationMultivariate Generalized Ornstein-Uhlenbeck Processes
Multivariate Generalized Ornstein-Uhlenbeck Processes Anita Behme TU München Alexander Lindner TU Braunschweig 7th International Conference on Lévy Processes: Theory and Applications Wroclaw, July 15 19,
More informationStochastic process. X, a series of random variables indexed by t
Stochastic process X, a series of random variables indexed by t X={X(t), t 0} is a continuous time stochastic process X={X(t), t=0,1, } is a discrete time stochastic process X(t) is the state at time t,
More information1 Simulating normal (Gaussian) rvs with applications to simulating Brownian motion and geometric Brownian motion in one and two dimensions
Copyright c 2007 by Karl Sigman 1 Simulating normal Gaussian rvs with applications to simulating Brownian motion and geometric Brownian motion in one and two dimensions Fundamental to many applications
More informationLecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations
Lecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations Simo Särkkä Aalto University, Finland November 18, 2014 Simo Särkkä (Aalto) Lecture 4: Numerical Solution of SDEs November
More informationLecture 6: Multiple Model Filtering, Particle Filtering and Other Approximations
Lecture 6: Multiple Model Filtering, Particle Filtering and Other Approximations Department of Biomedical Engineering and Computational Science Aalto University April 28, 2010 Contents 1 Multiple Model
More informationAn adaptive numerical scheme for fractional differential equations with explosions
An adaptive numerical scheme for fractional differential equations with explosions Johanna Garzón Departamento de Matemáticas, Universidad Nacional de Colombia Seminario de procesos estocásticos Jointly
More informationHandbook of Stochastic Methods
Springer Series in Synergetics 13 Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences von Crispin W Gardiner Neuausgabe Handbook of Stochastic Methods Gardiner schnell und portofrei
More informationFIRST-ORDER ORDINARY DIFFERENTIAL EQUATIONS I: Introduction and Analytic Methods David Levermore Department of Mathematics University of Maryland
FIRST-ORDER ORDINARY DIFFERENTIAL EQUATIONS I: Introduction and Analytic Methods David Levermore Department of Mathematics University of Maryland 4 September 2012 Because the presentation of this material
More informationLeast Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises
Least Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises Hongwei Long* Department of Mathematical Sciences, Florida Atlantic University, Boca Raton Florida 33431-991,
More informationStationary distributions of non Gaussian Ornstein Uhlenbeck processes for beam halos
N Cufaro Petroni: CYCLOTRONS 2007 Giardini Naxos, 1 5 October, 2007 1 Stationary distributions of non Gaussian Ornstein Uhlenbeck processes for beam halos CYCLOTRONS 2007 Giardini Naxos, 1 5 October Nicola
More informationNonparametric Drift Estimation for Stochastic Differential Equations
Nonparametric Drift Estimation for Stochastic Differential Equations Gareth Roberts 1 Department of Statistics University of Warwick Brazilian Bayesian meeting, March 2010 Joint work with O. Papaspiliopoulos,
More informationBirth and Death Processes. Birth and Death Processes. Linear Growth with Immigration. Limiting Behaviour for Birth and Death Processes
DTU Informatics 247 Stochastic Processes 6, October 27 Today: Limiting behaviour of birth and death processes Birth and death processes with absorbing states Finite state continuous time Markov chains
More informationGaussian processes for inference in stochastic differential equations
Gaussian processes for inference in stochastic differential equations Manfred Opper, AI group, TU Berlin November 6, 2017 Manfred Opper, AI group, TU Berlin (TU Berlin) inference in SDE November 6, 2017
More informationLocal vs. Nonlocal Diffusions A Tale of Two Laplacians
Local vs. Nonlocal Diffusions A Tale of Two Laplacians Jinqiao Duan Dept of Applied Mathematics Illinois Institute of Technology Chicago duan@iit.edu Outline 1 Einstein & Wiener: The Local diffusion 2
More informationInfinite-dimensional methods for path-dependent equations
Infinite-dimensional methods for path-dependent equations (Università di Pisa) 7th General AMaMeF and Swissquote Conference EPFL, Lausanne, 8 September 215 Based on Flandoli F., Zanco G. - An infinite-dimensional
More informationLecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Process Approximations
Lecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Process Approximations Simo Särkkä Aalto University Tampere University of Technology Lappeenranta University of Technology Finland November
More informationDISCRETE-TIME STOCHASTIC MODELS, SDEs, AND NUMERICAL METHODS. Ed Allen. NIMBioS Tutorial: Stochastic Models With Biological Applications
DISCRETE-TIME STOCHASTIC MODELS, SDEs, AND NUMERICAL METHODS Ed Allen NIMBioS Tutorial: Stochastic Models With Biological Applications University of Tennessee, Knoxville March, 2011 ACKNOWLEDGEMENT I thank
More informationBrownian Motion. Chapter Stochastic Process
Chapter 1 Brownian Motion 1.1 Stochastic Process A stochastic process can be thought of in one of many equivalent ways. We can begin with an underlying probability space (Ω, Σ,P and a real valued stochastic
More information6 Continuous-Time Birth and Death Chains
6 Continuous-Time Birth and Death Chains Angela Peace Biomathematics II MATH 5355 Spring 2017 Lecture notes follow: Allen, Linda JS. An introduction to stochastic processes with applications to biology.
More informationELEMENTS OF PROBABILITY THEORY
ELEMENTS OF PROBABILITY THEORY Elements of Probability Theory A collection of subsets of a set Ω is called a σ algebra if it contains Ω and is closed under the operations of taking complements and countable
More informationExercise Sheet 7 - Solutions
Exercise Sheet 7 - Solutions Alessandro Gnoatto June 10, 015 1 Exercise 1 1. Let Y t = log X t, then dy t = 1 dx t 1 X t Xt Xt σ dt by Ito s Lemma. Hence This implies that dy t = (µ σ )dt + σdw t. log
More informationWeak solutions of mean-field stochastic differential equations
Weak solutions of mean-field stochastic differential equations Juan Li School of Mathematics and Statistics, Shandong University (Weihai), Weihai 26429, China. Email: juanli@sdu.edu.cn Based on joint works
More informationLikelihood Functions for Stochastic Signals in White Noise* TYRONE E. DUNCAN
INFORMATION AND CONTROL 16, 303-310 (1970) Likelihood Functions for Stochastic Signals in White Noise* TYRONE E. DUNCAN Computer, Information and Control Engineering, The University of Nlichigan, Ann Arbor,
More informationRough paths methods 1: Introduction
Rough paths methods 1: Introduction Samy Tindel Purdue University University of Aarhus 216 Samy T. (Purdue) Rough Paths 1 Aarhus 216 1 / 16 Outline 1 Motivations for rough paths techniques 2 Summary of
More informationExact Simulation of Diffusions and Jump Diffusions
Exact Simulation of Diffusions and Jump Diffusions A work by: Prof. Gareth O. Roberts Dr. Alexandros Beskos Dr. Omiros Papaspiliopoulos Dr. Bruno Casella 28 th May, 2008 Content 1 Exact Algorithm Construction
More informationLecture 12: Detailed balance and Eigenfunction methods
Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 2015 Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility),
More informationContinuous-Time Markov Chain
Continuous-Time Markov Chain Consider the process {X(t),t 0} with state space {0, 1, 2,...}. The process {X(t),t 0} is a continuous-time Markov chain if for all s, t 0 and nonnegative integers i, j, x(u),
More information1. Stochastic Process
HETERGENEITY IN QUANTITATIVE MACROECONOMICS @ TSE OCTOBER 17, 216 STOCHASTIC CALCULUS BASICS SANG YOON (TIM) LEE Very simple notes (need to add references). It is NOT meant to be a substitute for a real
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
More informationOn the Rate of Convergence to Equilibrium for Two-sided Reflected Brownian Motion and for the Ornstein-Uhlenbeck Process
Submitted for Publication On the Rate of Convergence to Equilibrium for Two-sided Reflected Brownian Motion and for the Ornstein-Uhlenbeck Process Peter W. Glynn Rob J. Wang August 14, 018 Abstract This
More informationChain Rule. MATH 311, Calculus III. J. Robert Buchanan. Spring Department of Mathematics
3.33pt Chain Rule MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Spring 2019 Single Variable Chain Rule Suppose y = g(x) and z = f (y) then dz dx = d (f (g(x))) dx = f (g(x))g (x)
More informationStochastic Modelling Unit 1: Markov chain models
Stochastic Modelling Unit 1: Markov chain models Russell Gerrard and Douglas Wright Cass Business School, City University, London June 2004 Contents of Unit 1 1 Stochastic Processes 2 Markov Chains 3 Poisson
More informationDeterministic Dynamic Programming
Deterministic Dynamic Programming 1 Value Function Consider the following optimal control problem in Mayer s form: V (t 0, x 0 ) = inf u U J(t 1, x(t 1 )) (1) subject to ẋ(t) = f(t, x(t), u(t)), x(t 0
More informationAsymptotic properties of maximum likelihood estimator for the growth rate for a jump-type CIR process
Asymptotic properties of maximum likelihood estimator for the growth rate for a jump-type CIR process Mátyás Barczy University of Debrecen, Hungary The 3 rd Workshop on Branching Processes and Related
More informationPath integrals for classical Markov processes
Path integrals for classical Markov processes Hugo Touchette National Institute for Theoretical Physics (NITheP) Stellenbosch, South Africa Chris Engelbrecht Summer School on Non-Linear Phenomena in Field
More information