Lecture on Parameter Estimation for Stochastic Differential Equations. Erik Lindström
|
|
- Randall Hall
- 5 years ago
- Views:
Transcription
1 Lecture on Parameter Estimation for Stochastic Differential Equations Erik Lindström
2 Recap We are interested in the parameters θ in the Stochastic Integral Equations X(t) = X(0) + t 0 µ θ (s, X(s))ds + t 0 σ θ (s, X(s))dW(s) (1)
3 Recap We are interested in the parameters θ in the Stochastic Integral Equations X(t) = X(0) + t Why? Model validation 0 µ θ (s, X(s))ds + t 0 σ θ (s, X(s))dW(s) (1)
4 Recap We are interested in the parameters θ in the Stochastic Integral Equations X(t) = X(0) + t Why? Model validation Risk management 0 µ θ (s, X(s))ds + t 0 σ θ (s, X(s))dW(s) (1)
5 Recap We are interested in the parameters θ in the Stochastic Integral Equations X(t) = X(0) + t 0 µ θ (s, X(s))ds + t 0 σ θ (s, X(s))dW(s) Why? Model validation Risk management Advanced hedging (Greeks and quadratic hedging (P/Q)) (1)
6 Some asymptotics Consider the arithmetic Brownian motion dx(t) = µdt + σdw(t) (2)
7 Some asymptotics Consider the arithmetic Brownian motion dx(t) = µdt + σdw(t) (2) The drift is estimated by computing the mean, and compensating for the sampling δ = t n+1 t n ˆµ = 1 N 1 X(t n+1 ) X(t n ). (3) δn n=0
8 Some asymptotics Consider the arithmetic Brownian motion dx(t) = µdt + σdw(t) (2) The drift is estimated by computing the mean, and compensating for the sampling δ = t n+1 t n ˆµ = 1 N 1 X(t n+1 ) X(t n ). (3) δn n=0 Expanding this expression reveals that the MLE is given by ˆµ = X(t N) X(t 0 ) t N t 0 = µ + σ W(t N) W(t 0 ) t N t 0. (4)
9 Some asymptotics Consider the arithmetic Brownian motion dx(t) = µdt + σdw(t) (2) The drift is estimated by computing the mean, and compensating for the sampling δ = t n+1 t n ˆµ = 1 N 1 X(t n+1 ) X(t n ). (3) δn n=0 Expanding this expression reveals that the MLE is given by ˆµ = X(t N) X(t 0 ) t N t 0 = µ + σ W(t N) W(t 0 ) t N t 0. (4) The MLE for the diffusion (σ) parameter is given by ˆσ 2 = N 1 1 (X(t n+1 ) X(t n ) ˆµδ) 2 d σ 2 χ2 (N 1) δ(n 1) N 1 n=0
10 A simple method Many data sets are sampled at high frequency, making the bias due to discretization of the SDEs some of the schemes in Chapter 12 acceptable.
11 A simple method Many data sets are sampled at high frequency, making the bias due to discretization of the SDEs some of the schemes in Chapter 12 acceptable. The simplest discretization, the Explicit Euler method, would for the stochastic differential equation dx(t) = µ(t, X(t))dt + σ(t, X(t))dW(t) (6) correspond the Discretized Maximum Likelihood (DML) estimator given by N 1 ˆθ DML = arg max log ϕ (X(t n+1 ), X(t n ) + µ(t n, X(t n )), Σ(t n, X(t θ Θ n=1 (7) where ϕ(x, m, P) is the density for a multivariate Normal distribution with argument x, mean m and covariance P and Σ(t, X(t)) = σ(t, X(t))σ(t, X(t)) T. (8)
12 Consistency The DMLE is generally NOT consistent.
13 Consistency The DMLE is generally NOT consistent. Approximate ML estimators (13.5) are, provided enough computational resources are allocated Simulation based estimators Fokker-Planck based estimators Series expansions.
14 Consistency The DMLE is generally NOT consistent. Approximate ML estimators (13.5) are, provided enough computational resources are allocated Simulation based estimators Fokker-Planck based estimators Series expansions. GMM-type estimators (13.6) are consistent if the moments are correctly specified (which is a non-trivial problem!)
15 Simultion based estimators Discretely observed SDEs are Markov processes
16 Simultion based estimators Discretely observed SDEs are Markov processes Then it follows that p θ (x t x s ) = E θ [p θ (x t x τ ) F(s)], t > τ > s (9) This is the Pedersen algorithm.
17 Simultion based estimators Discretely observed SDEs are Markov processes Then it follows that p θ (x t x s ) = E θ [p θ (x t x τ ) F(s)], t > τ > s (9) This is the Pedersen algorithm. Improved by Durham-Gallant (2002) and Lindström (2012)
18 Simultion based estimators Discretely observed SDEs are Markov processes Then it follows that p θ (x t x s ) = E θ [p θ (x t x τ ) F(s)], t > τ > s (9) This is the Pedersen algorithm. Improved by Durham-Gallant (2002) and Lindström (2012) Works very well for Multivariate models! and is easily (...) extended to Levy driven SDEs.
19 Some key points Naive implementation only provides a point wise estimate - use CRNs or importance sampling Variance reduction helps (antithetic variates, control variates) The near optimal importance sampler is a Bridge process, as it reduces variance AND improves the asymptotics. There is a version that is completely bias free, albeit somewhat restrictive in terms of the class of feasible models.
20 Fokker-Planck Consider the expectation E [h(x(t)) F(0)] = h(x(t))p(x(t) x(0))dx(t) (10) and then E [h(x(t)) F(0)] (11) t
21 Fokker-Planck Consider the expectation E [h(x(t)) F(0)] = h(x(t))p(x(t) x(0))dx(t) (10) and then E [h(x(t)) F(0)] (11) t Two possible ways to compute this, direct and using the Itō formula.
22 Fokker-Planck Consider the expectation E [h(x(t)) F(0)] = h(x(t))p(x(t) x(0))dx(t) (10) and then E [h(x(t)) F(0)] (11) t Two possible ways to compute this, direct and using the Itō formula. Equating these yields where p t (x(t) x(0)) = A p(x(t) x(0)) (12) A p(x(t)) = x(t) (µ( )p(x(t)))+1 2 x 2 (t) 2 ( ) σ 2 ( )p(x(t)). (13)
23 Example of the Fokker-Planck equation From (Lindström, 2007) p(s,x s ;t,x t ) Grid Time Figure: Fokker-Planck equation computed for the CKLS process
24 Comments on the PDE approach Generally better than the Monte Carlo method in low dimensional problems Durham Gallant Poulsen Order2 Pade(1,1) Order4 Pade(2,2) 10 2 MAE time Figure: Comparing Monte Carlo, 2nd order and 4th order numerical approximations of the Fokker-Planck equation
25 Discussion Fokker-Planck is the preferred method is the state space is non-trivial (see Pedersen et. al, 2011)
26 Discussion Fokker-Planck is the preferred method is the state space is non-trivial (see Pedersen et. al, 2011) Successfully used in 1-d and 2-d problems
27 Discussion Fokker-Planck is the preferred method is the state space is non-trivial (see Pedersen et. al, 2011) Successfully used in 1-d and 2-d problems but the ``curse of dimensionality'' will eventually make the method infeasable
28 Series expansion The solution to the Fokker-Planck equation when dx(t) = µdt + σdw(t) (14) is p(x(t) x(0)) = N(x(t); x(0) + µt, σ 2 t).
29 Series expansion The solution to the Fokker-Planck equation when dx(t) = µdt + σdw(t) (14) is p(x(t) x(0)) = N(x(t); x(0) + µt, σ 2 t). Hermite polynomials are the orthogonal polynomial basis when using a Gaussian as weight function.
30 Series expansion The solution to the Fokker-Planck equation when dx(t) = µdt + σdw(t) (14) is p(x(t) x(0)) = N(x(t); x(0) + µt, σ 2 t). Hermite polynomials are the orthogonal polynomial basis when using a Gaussian as weight function. This is used in the 'series expansion approach', see e.g. (Ait-Sahalia, 2002, 2008)
31 Key ideas Transform from X Y Z where Z is approximately standard Gaussian. We assume that dx(t) = µ(x(t))dt + σ(x(t))dw(t) (15)
32 Key ideas Transform from X Y Z where Z is approximately standard Gaussian. We assume that dx(t) = µ(x(t))dt + σ(x(t))dw(t) (15) First step. Y(t) = du σ(u) (16)
33 Key ideas Transform from X Y Z where Z is approximately standard Gaussian. We assume that dx(t) = µ(x(t))dt + σ(x(t))dw(t) (15) First step. Y(t) = It then follows that du σ(u) (16) dy(t) = µ Y (Y(t))dt + dw(t) (17)
34 Key ideas Transform from X Y Z where Z is approximately standard Gaussian. We assume that dx(t) = µ(x(t))dt + σ(x(t))dw(t) (15) First step. Y(t) = It then follows that du σ(u) (16) Second step: Transform dy(t) = µ Y (Y(t))dt + dw(t) (17) Z(t k ) = Y(t k) Y(t k 1 ) t k t k 1. (18)
35 Expansion A Hermite expansion for the density p Z at order J is given by p J Z (z y(0), t k t k 1 ) = ϕ(z) J η j (t k t k 1, y 0 )H j (z) (19) j=0 where H j (z) = e z2 /2 dj dz j e z2 /2. (20)
36 Expansion A Hermite expansion for the density p Z at order J is given by p J Z (z y(0), t k t k 1 ) = ϕ(z) where J η j (t k t k 1, y 0 )H j (z) (19) j=0 H j (z) = e z2 /2 dj dz j e z2 /2. (20) The coefficients are computed by projecting the density onto the basis functions H j (z) (recall Hilbert space theory) η j (t, y 0 ) = 1 H j! j (z)p J Z (z y(0), t k t k 1 )dz (21)
37 Practical concerns The series expansion can be extremely accurate. The standard approach is to compute η j by Taylor expansion up to order (t k t k 1 ) K = ( t) K Some restrictions 'so-called reducible diffusion' when using the method for multivariate diffusions.
38 Other alternatives - GMM/EF What about non-likelihood methods? The model is governed by some p-dimensional parameter.
39 Other alternatives - GMM/EF What about non-likelihood methods? The model is governed by some p-dimensional parameter. Suppose some set of features are important, h l (x), l = 1,..., q p
40 Other alternatives - GMM/EF What about non-likelihood methods? The model is governed by some p-dimensional parameter. Suppose some set of features are important, h l (x), l = 1,..., q p Compute h 1 (x(t)) E θ [h 1 (X t )] f(x(t); θ) =... (22) h q (x t ) E θ [h q (X t )] and form J N (θ) = ( 1 N T ( N 1 f(x(n); θ)) W N n=1 ) N f(x(n); θ) n=1 (23)
41 GMM The Generalized Methods of Moments (GMM) estimators is then given by ˆθ = arg min J N (θ) (24)
42 GMM The Generalized Methods of Moments (GMM) estimators is then given by ˆθ = arg min J N (θ) (24) It can be shown that N (ˆθN θ 0 ) N(0, Σ) (25) where Σ = ( Γ T N Ω 1 N Γ N) 1. (26) and Γ N and Ω N are estimates of [( )] f(x, θ) Γ = E θ T, Ω = Var [f(x, θ)]. (27)
43 Quasi Likelihood What if we treat the distribution as Gaussian? This is called quasi likelihood. Asymptotics can be derived from the GMM asymptotics (Höök & Lindström, 2016) showed that this can be extremely efficient from a computational point of view, essentially O(1) for N observations.
44 References Pedersen, A. R. (1995). A new approach to maximum likelihood estimation for stochastic differential equations based on discrete observations. Scandinavian journal of statistics, AïtSahalia, Y. (2002). Maximum Likelihood Estimation of Discretely Sampled Diffusions: A Closedform Approximation Approach. Econometrica, 70(1), Lindström, E. (2007). Estimating parameters in diffusion processes using an approximate maximum likelihood approach. Annals of Operations Research, 151(1),
45 References Durham, G. B. & Gallant, A. R. (2002). Numerical techniques for maximum likelihood estimation of continuous-time diffusion processes. Journal of Business & Economic Statistics, 20(3), Lindström, E. (2012). A regularized bridge sampler for sparsely sampled diffusions. Statistics and Computing, 22(2), Whitaker, G. A., Golightly, A., Boys, R. J., & Sherlock, C. (2016). Improved bridge constructs for stochastic differential equations. Statistics and Computing, Höök, L. J., & Lindström, E. (2016). Efficient computation of the quasi likelihood function for discretely observed diffusion processes. Computational Statistics & Data Analysis.
46
Lecture 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 informationarxiv: v1 [math.pr] 1 Jul 2013
ESTIMATION OF FIRST PASSAGE TIME DENSITIES OF DIFFUSIONS PROCESSES THROUGH TIME-VARYING BOUNDARIES arxiv:307.0336v [math.pr] Jul 03 Imene Allab and Francois Watier Department of Mathematics, Université
More informationDSGE Methods. Estimation of DSGE models: GMM and Indirect Inference. Willi Mutschler, M.Sc.
DSGE Methods Estimation of DSGE models: GMM and Indirect Inference Willi Mutschler, M.Sc. Institute of Econometrics and Economic Statistics University of Münster willi.mutschler@wiwi.uni-muenster.de Summer
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 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 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 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 informationNumerical Simulation of Stochastic Differential Equations: Lecture 2, Part 1. Recap: SDE. Euler Maruyama. Lecture 2, Part 1: Euler Maruyama
Numerical Simulation of Stochastic Differential Equations: Lecture, Part 1 Des Higham Department of Mathematics Universit of Strathclde Lecture, Part 1: Euler Maruama Definition of Euler Maruama Method
More informationDSGE-Models. Limited Information Estimation General Method of Moments and Indirect Inference
DSGE-Models General Method of Moments and Indirect Inference Dr. Andrea Beccarini Willi Mutschler, M.Sc. Institute of Econometrics and Economic Statistics University of Münster willi.mutschler@uni-muenster.de
More informationGillespie s Algorithm and its Approximations. Des Higham Department of Mathematics and Statistics University of Strathclyde
Gillespie s Algorithm and its Approximations Des Higham Department of Mathematics and Statistics University of Strathclyde djh@maths.strath.ac.uk The Three Lectures 1 Gillespie s algorithm and its relation
More informationList of projects. FMS020F NAMS002 Statistical inference for partially observed stochastic processes, 2016
List of projects FMS020F NAMS002 Statistical inference for partially observed stochastic processes, 206 Work in groups of two (if this is absolutely not possible for some reason, please let the lecturers
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 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 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 informationIntroduction to Estimation Methods for Time Series models Lecture 2
Introduction to Estimation Methods for Time Series models Lecture 2 Fulvio Corsi SNS Pisa Fulvio Corsi Introduction to Estimation () Methods for Time Series models Lecture 2 SNS Pisa 1 / 21 Estimators:
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 informationGeneralized Method of Moment
Generalized Method of Moment CHUNG-MING KUAN Department of Finance & CRETA National Taiwan University June 16, 2010 C.-M. Kuan (Finance & CRETA, NTU Generalized Method of Moment June 16, 2010 1 / 32 Lecture
More informationSimulation and Parametric Estimation of SDEs
Simulation and Parametric Estimation of SDEs Jhonny Gonzalez School of Mathematics The University of Manchester Magical books project August 23, 2012 Motivation The problem Simulation of SDEs SDEs driven
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 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 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 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 informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond January 18, 2017 Contents 1 Batch and Recursive Estimation 2 Towards Bayesian Filtering 3 Kalman Filter and Bayesian Filtering and Smoothing
More informationGeometric projection of stochastic differential equations
Geometric projection of stochastic differential equations John Armstrong (King s College London) Damiano Brigo (Imperial) August 9, 2018 Idea: Projection Idea: Projection Projection gives a method of systematically
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 informationSkew Brownian motion with dry friction: The Pugachev Sveshnikov equation approach
Skew Brownian motion with dry friction: The Pugachev Sveshnikov equation approach arxiv:175.198v2 [math.pr] 27 Jul 217 Sergey Berezin and Oleg Zayats Department of Applied Mathematics Peter the Great St.Petersburg
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 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 informationChapter 1: A Brief Review of Maximum Likelihood, GMM, and Numerical Tools. Joan Llull. Microeconometrics IDEA PhD Program
Chapter 1: A Brief Review of Maximum Likelihood, GMM, and Numerical Tools Joan Llull Microeconometrics IDEA PhD Program Maximum Likelihood Chapter 1. A Brief Review of Maximum Likelihood, GMM, and Numerical
More informationIntroduction to the Estimation of Stochastic Differential Equations Based on Discrete Observations
Stochastic Finance 4 Autumn School & International Conference Introduction to the Estimation of Stochastic Differential Equations Based on Discrete Observations João Nicolau Instituto Superior de Economia
More informationADVANCED FINANCIAL ECONOMETRICS PROF. MASSIMO GUIDOLIN
Massimo Guidolin Massimo.Guidolin@unibocconi.it Dept. of Finance ADVANCED FINANCIAL ECONOMETRICS PROF. MASSIMO GUIDOLIN a.a. 14/15 p. 1 LECTURE 3: REVIEW OF BASIC ESTIMATION METHODS: GMM AND OTHER EXTREMUM
More informationA nonparametric method of multi-step ahead forecasting in diffusion processes
A nonparametric method of multi-step ahead forecasting in diffusion processes Mariko Yamamura a, Isao Shoji b a School of Pharmacy, Kitasato University, Minato-ku, Tokyo, 108-8641, Japan. b Graduate School
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 informationInverse Langevin approach to time-series data analysis
Inverse Langevin approach to time-series data analysis Fábio Macêdo Mendes Anníbal Dias Figueiredo Neto Universidade de Brasília Saratoga Springs, MaxEnt 2007 Outline 1 Motivation 2 Outline 1 Motivation
More informationSTA414/2104 Statistical Methods for Machine Learning II
STA414/2104 Statistical Methods for Machine Learning II Murat A. Erdogdu & David Duvenaud Department of Computer Science Department of Statistical Sciences Lecture 3 Slide credits: Russ Salakhutdinov Announcements
More informationXiaohu Wang, Peter C. B. Phillips and Jun Yu. January 2011 COWLES FOUNDATION DISCUSSION PAPER NO. 1778
BIAS IN ESTIMATING MULTIVARIATE AND UNIVARIATE DIFFUSIONS By Xiaohu Wang, Peter C. B. Phillips and Jun Yu January 20 COWLES FOUNDATION DISCUSSION PAPER NO. 778 COWLES FOUNDATION FOR RESEARCH IN ECONOMICS
More informationMonte-Carlo Methods and Stochastic Processes
Monte-Carlo Methods and Stochastic Processes From Linear to Non-Linear EMMANUEL GOBET ECOLE POLYTECHNIQUE - UNIVERSITY PARIS-SACLAY CMAP, PALAISEAU CEDEX, FRANCE CRC Press Taylor & Francis Group 6000 Broken
More informationVerona Course April Lecture 1. Review of probability
Verona Course April 215. Lecture 1. Review of probability Viorel Barbu Al.I. Cuza University of Iaşi and the Romanian Academy A probability space is a triple (Ω, F, P) where Ω is an abstract set, F is
More informationOn the Power of Tests for Regime Switching
On the Power of Tests for Regime Switching joint work with Drew Carter and Ben Hansen Douglas G. Steigerwald UC Santa Barbara May 2015 D. Steigerwald (UCSB) Regime Switching May 2015 1 / 42 Motivating
More informationTMS165/MSA350 Stochastic Calculus, Lecture on Applications
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
More informationEcon 583 Final Exam Fall 2008
Econ 583 Final Exam Fall 2008 Eric Zivot December 11, 2008 Exam is due at 9:00 am in my office on Friday, December 12. 1 Maximum Likelihood Estimation and Asymptotic Theory Let X 1,...,X n be iid random
More informationProblems 5: Continuous Markov process and the diffusion equation
Problems 5: Continuous Markov process and the diffusion equation Roman Belavkin Middlesex University Question Give a definition of Markov stochastic process. What is a continuous Markov process? Answer:
More informationOPTIMAL DESIGN INPUTS FOR EXPERIMENTAL CHAPTER 17. Organization of chapter in ISSO. Background. Linear models
CHAPTER 17 Slides for Introduction to Stochastic Search and Optimization (ISSO)by J. C. Spall OPTIMAL DESIGN FOR EXPERIMENTAL INPUTS Organization of chapter in ISSO Background Motivation Finite sample
More informationSome Terminology and Concepts that We will Use, But Not Emphasize (Section 6.2)
Some Terminology and Concepts that We will Use, But Not Emphasize (Section 6.2) Statistical analysis is based on probability theory. The fundamental object in probability theory is a probability space,
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 informationEstimation for the standard and geometric telegraph process. Stefano M. Iacus. (SAPS VI, Le Mans 21-March-2007)
Estimation for the standard and geometric telegraph process Stefano M. Iacus University of Milan(Italy) (SAPS VI, Le Mans 21-March-2007) 1 1. Telegraph process Consider a particle moving on the real line
More informationObtaining Critical Values for Test of Markov Regime Switching
University of California, Santa Barbara From the SelectedWorks of Douglas G. Steigerwald November 1, 01 Obtaining Critical Values for Test of Markov Regime Switching Douglas G Steigerwald, University of
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 informationNumerical treatment of stochastic delay differential equations
Numerical treatment of stochastic delay differential equations Evelyn Buckwar Department of Mathematics Heriot-Watt University, Edinburgh Evelyn Buckwar Roma, 16.10.2007 Overview An example: machine chatter
More informationBayesian Inference for Discretely Sampled Diffusion Processes: A New MCMC Based Approach to Inference
Bayesian Inference for Discretely Sampled Diffusion Processes: A New MCMC Based Approach to Inference Osnat Stramer 1 and Matthew Bognar 1 Department of Statistics and Actuarial Science, University of
More informationFall 2017 STAT 532 Homework Peter Hoff. 1. Let P be a probability measure on a collection of sets A.
1. Let P be a probability measure on a collection of sets A. (a) For each n N, let H n be a set in A such that H n H n+1. Show that P (H n ) monotonically converges to P ( k=1 H k) as n. (b) For each n
More informationGeneralized Method of Moments (GMM) Estimation
Econometrics 2 Fall 2004 Generalized Method of Moments (GMM) Estimation Heino Bohn Nielsen of29 Outline of the Lecture () Introduction. (2) Moment conditions and methods of moments (MM) estimation. Ordinary
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 informationContinuum Limit of Forward Kolmogorov Equation Friday, March 06, :04 PM
Continuum Limit of Forward Kolmogorov Equation Friday, March 06, 2015 2:04 PM Please note that one of the equations (for ordinary Brownian motion) in Problem 1 was corrected on Wednesday night. And actually
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 informationGaussian, Markov and stationary processes
Gaussian, Markov and stationary processes Gonzalo Mateos Dept. of ECE and Goergen Institute for Data Science University of Rochester gmateosb@ece.rochester.edu http://www.ece.rochester.edu/~gmateosb/ November
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 informationComputer Intensive Methods in Mathematical Statistics
Computer Intensive Methods in Mathematical Statistics Department of mathematics johawes@kth.se Lecture 16 Advanced topics in computational statistics 18 May 2017 Computer Intensive Methods (1) Plan of
More informationMarkov chain Monte Carlo algorithms for SDE parameter estimation
Markov chain Monte Carlo algorithms for SDE parameter estimation Andrew Golightly and Darren J. Wilkinson Abstract This chapter considers stochastic differential equations for Systems Biology models derived
More informationAccurate Maximum Likelihood Estimation for Parametric Population Analysis. Bob Leary UCSD/SDSC and LAPK, USC School of Medicine
Accurate Maximum Likelihood Estimation for Parametric Population Analysis Bob Leary UCSD/SDSC and LAPK, USC School of Medicine Why use parametric maximum likelihood estimators? Consistency: θˆ θ as N ML
More informationThe Wiener Itô Chaos Expansion
1 The Wiener Itô Chaos Expansion The celebrated Wiener Itô chaos expansion is fundamental in stochastic analysis. In particular, it plays a crucial role in the Malliavin calculus as it is presented in
More informationIEOR 165 Lecture 7 1 Bias-Variance Tradeoff
IEOR 165 Lecture 7 Bias-Variance Tradeoff 1 Bias-Variance Tradeoff Consider the case of parametric regression with β R, and suppose we would like to analyze the error of the estimate ˆβ in comparison to
More informationCOMS 4721: Machine Learning for Data Science Lecture 1, 1/17/2017
COMS 4721: Machine Learning for Data Science Lecture 1, 1/17/2017 Prof. John Paisley Department of Electrical Engineering & Data Science Institute Columbia University OVERVIEW This class will cover model-based
More informationHarnack Inequalities and Applications for Stochastic Equations
p. 1/32 Harnack Inequalities and Applications for Stochastic Equations PhD Thesis Defense Shun-Xiang Ouyang Under the Supervision of Prof. Michael Röckner & Prof. Feng-Yu Wang March 6, 29 p. 2/32 Outline
More informationReflected Brownian Motion
Chapter 6 Reflected Brownian Motion Often we encounter Diffusions in regions with boundary. If the process can reach the boundary from the interior in finite time with positive probability we need to decide
More informationentropy Continuous-Discrete Path Integral Filtering Entropy 2009, 11, ; doi: /e ISSN
Entropy 9,, 4-43; doi:.339/e34 Article Continuous-Discrete Path Integral Filtering OPEN ACCESS entropy ISSN 99-43 www.mdpi.com/journal/entropy Bhashyam Balaji Radar Systems Section, Defence Research and
More informationEstimation, Inference, and Hypothesis Testing
Chapter 2 Estimation, Inference, and Hypothesis Testing Note: The primary reference for these notes is Ch. 7 and 8 of Casella & Berger 2. This text may be challenging if new to this topic and Ch. 7 of
More informationProbabilistic Graphical Models
Probabilistic Graphical Models Brown University CSCI 2950-P, Spring 2013 Prof. Erik Sudderth Lecture 12: Gaussian Belief Propagation, State Space Models and Kalman Filters Guest Kalman Filter Lecture by
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 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 informationMarket Risk. MFM Practitioner Module: Quantitiative Risk Management. John Dodson. February 8, Market Risk. John Dodson.
MFM Practitioner Module: Quantitiative Risk Management February 8, 2017 This week s material ties together our discussion going back to the beginning of the fall term about risk measures based on the (single-period)
More informationSupplemental Material for KERNEL-BASED INFERENCE IN TIME-VARYING COEFFICIENT COINTEGRATING REGRESSION. September 2017
Supplemental Material for KERNEL-BASED INFERENCE IN TIME-VARYING COEFFICIENT COINTEGRATING REGRESSION By Degui Li, Peter C. B. Phillips, and Jiti Gao September 017 COWLES FOUNDATION DISCUSSION PAPER NO.
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 informationInferring biological dynamics Iterated filtering (IF)
Inferring biological dynamics 101 3. Iterated filtering (IF) IF originated in 2006 [6]. For plug-and-play likelihood-based inference on POMP models, there are not many alternatives. Directly estimating
More informationIf we want to analyze experimental or simulated data we might encounter the following tasks:
Chapter 1 Introduction If we want to analyze experimental or simulated data we might encounter the following tasks: Characterization of the source of the signal and diagnosis Studying dependencies Prediction
More informationControlled Diffusions and Hamilton-Jacobi Bellman Equations
Controlled Diffusions and Hamilton-Jacobi Bellman Equations Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 2014 Emo Todorov (UW) AMATH/CSE 579, Winter
More informationApproximation of BSDEs using least-squares regression and Malliavin weights
Approximation of BSDEs using least-squares regression and Malliavin weights Plamen Turkedjiev (turkedji@math.hu-berlin.de) 3rd July, 2012 Joint work with Prof. Emmanuel Gobet (E cole Polytechnique) Plamen
More informationHomework # , Spring Due 14 May Convergence of the empirical CDF, uniform samples
Homework #3 36-754, Spring 27 Due 14 May 27 1 Convergence of the empirical CDF, uniform samples In this problem and the next, X i are IID samples on the real line, with cumulative distribution function
More informationEconometrics I, Estimation
Econometrics I, Estimation Department of Economics Stanford University September, 2008 Part I Parameter, Estimator, Estimate A parametric is a feature of the population. An estimator is a function of the
More informationEconomics 583: Econometric Theory I A Primer on Asymptotics
Economics 583: Econometric Theory I A Primer on Asymptotics Eric Zivot January 14, 2013 The two main concepts in asymptotic theory that we will use are Consistency Asymptotic Normality Intuition consistency:
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 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 informationF & B Approaches to a simple model
A6523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 215 http://www.astro.cornell.edu/~cordes/a6523 Lecture 11 Applications: Model comparison Challenges in large-scale surveys
More informationFinal Exam November 24, Problem-1: Consider random walk with drift plus a linear time trend: ( t
Problem-1: Consider random walk with drift plus a linear time trend: y t = c + y t 1 + δ t + ϵ t, (1) where {ϵ t } is white noise with E[ϵ 2 t ] = σ 2 >, and y is a non-stochastic initial value. (a) Show
More informationNon-orthogonal Hilbert Projections in Trend Regression
Non-orthogonal Hilbert Projections in Trend Regression Peter C. B. Phillips Cowles Foundation for Research in Economics Yale University, University of Auckland & University of York and Yixiao Sun Department
More informationSequential Monte Carlo Methods for Bayesian Computation
Sequential Monte Carlo Methods for Bayesian Computation A. Doucet Kyoto Sept. 2012 A. Doucet (MLSS Sept. 2012) Sept. 2012 1 / 136 Motivating Example 1: Generic Bayesian Model Let X be a vector parameter
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 informationModelling Non-linear and Non-stationary Time Series
Modelling Non-linear and Non-stationary Time Series Chapter 2: Non-parametric methods Henrik Madsen Advanced Time Series Analysis September 206 Henrik Madsen (02427 Adv. TS Analysis) Lecture Notes September
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 informationGMM-based inference in the AR(1) panel data model for parameter values where local identi cation fails
GMM-based inference in the AR() panel data model for parameter values where local identi cation fails Edith Madsen entre for Applied Microeconometrics (AM) Department of Economics, University of openhagen,
More information13. Parameter Estimation. ECE 830, Spring 2014
13. Parameter Estimation ECE 830, Spring 2014 1 / 18 Primary Goal General problem statement: We observe X p(x θ), θ Θ and the goal is to determine the θ that produced X. Given a collection of observations
More informationIntroduction to Estimation Methods for Time Series models. Lecture 1
Introduction to Estimation Methods for Time Series models Lecture 1 Fulvio Corsi SNS Pisa Fulvio Corsi Introduction to Estimation () Methods for Time Series models Lecture 1 SNS Pisa 1 / 19 Estimation
More informationComparison of continuous and discrete-time data-based modeling for hypoelliptic systems
Comparison of continuous and discrete- data-based modeling for hypoelliptic systems Fei Lu, Kevin K. Lin, and Alexandre J. Chorin arxiv:1605.02273v2 [math.na] 7 Dec 2016 Abstract We compare two approaches
More informationKolmogorov 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 informationLAN property for sde s with additive fractional noise and continuous time observation
LAN property for sde s with additive fractional noise and continuous time observation Eulalia Nualart (Universitat Pompeu Fabra, Barcelona) joint work with Samy Tindel (Purdue University) Vlad s 6th birthday,
More informationThe moment-generating function of the log-normal distribution using the star probability measure
Noname manuscript No. (will be inserted by the editor) The moment-generating function of the log-normal distribution using the star probability measure Yuri Heymann Received: date / Accepted: date Abstract
More informationON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER
ON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER GERARDO HERNANDEZ-DEL-VALLE arxiv:1209.2411v1 [math.pr] 10 Sep 2012 Abstract. This work deals with first hitting time densities of Ito processes whose
More informationTheoretical Statistics. Lecture 17.
Theoretical Statistics. Lecture 17. Peter Bartlett 1. Asymptotic normality of Z-estimators: classical conditions. 2. Asymptotic equicontinuity. 1 Recall: Delta method Theorem: Supposeφ : R k R m is differentiable
More informationBayesian Regression Linear and Logistic Regression
When we want more than point estimates Bayesian Regression Linear and Logistic Regression Nicole Beckage Ordinary Least Squares Regression and Lasso Regression return only point estimates But what if we
More informationStochastic Modelling in Climate Science
Stochastic Modelling in Climate Science David Kelly Mathematics Department UNC Chapel Hill dtbkelly@gmail.com November 16, 2013 David Kelly (UNC) Stochastic Climate November 16, 2013 1 / 36 Why use stochastic
More information