New robust inference for predictive regressions
|
|
- Barbra Riley
- 5 years ago
- Views:
Transcription
1 New robust inference for predictive regressions Anton Skrobotov Russian Academy of National Economy and Public Administration and Innopolis University based on joint work with Rustam Ibragimov and Jihyun Kim Research supported by a grant from Russian Science Foundaton (Project No ) The IAAE grant is gratefully acknowledged Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
2 Motivation Motivation Endogeneity Consider the following predictive regression for t = 1,..., T, where x t is some covariate. y t = α + x t 1 + u t, (1) x t = ρx t 1 + v t, (2) Our purpose to test the null hypothesis of no predictability of y t (e.g., stock returns). In other words, we want to test H 0 : = 0. If the process x t be stationary then we estimate by OLS and construct t-ratio, t( ˆ OLS ), which converges to standard Normal distribution. Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
3 Motivation Motivation Endogeneity Consider the following predictive regression for t = 1,..., T, where x t is some covariate. y t = α + x t 1 + u t, (3) x t = ρx t 1 + v t, (4) If, howewer, x t is (near) non-stationary, ( x t = 1 c ) x t 1 + v t, T we can use standard Normal inference only if the u t and v t are uncorrelated. Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
4 Motivation Motivation Endogeneity Consider the following predictive regression y t = α + x t 1 + u t, (5) x t = ρx t 1 + v t, (6) More presicely, let ξ i = (u i, v i ), and (appropriately normalizing) u t converges to a bivariate Wiener process ([ [rt ]] 1 [rt ] i=1 ξ i B(s) = (U(r), V (r)) ) with covariance matrix ( ) ω 2 Ω = y ω xy ω xy It can be shown that t( ˆ OLS ) P + Q, where Q is normal and P is function of Ornstein-Uhlenbeck process (near-unit root distributions). If ω xy = 0, then P vanishes from the limiting distribution, so that t( ˆ OLS ) N(0, 1). Otherwise, there are serious size distortions. ω 2 x Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
5 Motivation Motivation Non-stationary volatility Consider the following predictive regression y t = α + x t 1 + u t, (7) x t = ρx t 1 + v t, (8) u t = σ t 1 ε t, where ε t is MDS w.r.t. ltration F t 1, s.t. E(ε 2 t F t 1 ) = 1. Therefore, E(u y,2 t F t 1 ) = σ 2 t 1. The volatility process σ t may be (nearly) non-stationary: σ t = ω(z t ), where z t is (near) unit root process, or σ t = ω(t/t ) (xed or random), or σ t = ω(z t / T ) (see Choi et al., 2016). Then the limiting distribution of t( ˆ OLS ) is non-normal, even when there is no endogeneity (no dependence between u t and v t ). Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
6 Choi et al., 2016 The robust test for no predictability First assume α = 0 in predictive regression. Also assume MDS assumption of u t (no non-stationary volatility). Choi et al., 2016 proposed so called Cauchy estimator of, ( T ) 1 T ˆ T = x t 1 sgn(x t 1 )y t, (9) t=1 where sgn( ) is a sign function such that sign(x) = 1 for x 0 and sign(x) = 1 for x < 0. Extention 1: α 0 recursive de-meaning (the limiting distribution is the same). Extention 2: u t follows non-stationary volatility assumption Time Change in continuous time framework (the limiting distribution is the same). t=1 Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
7 Choi et al., 2016 The robust test for no predictability It can be shown that under the null hypothesis of = 0, 1 T T t=1 u t 1 0 σ(s)dw (s) = d MN(0, Q), where Q is some (random) variance depending on the volatility process (Q = 1 0 σ(s)2 ds), and therefore 1 T T t=1 sgn(x t 1 )y t = 1 T T t=1 sgn(x t 1 )u t Therefore, the numerator of the Cauchy estimator, ˆδ T := 1 0 σ(s)dw (s) = d MN(0, Q). T sgn(x t 1 )y t (10) t=1 has approximately mixed normal distribution with mean zero and variance Q T. Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
8 Choi et al., 2016 Brief review of t-statistic based robust inference Bakirov and Szekely (2005), Ibragimov and Mueller (2010, 2016): Usual small sample t-test of level α 5%: conservative for independent heterogeneous Gaussian observations (not α = 10%) X j N(µ, σ 2 j ), j = 1,..., q: H 0 : µ = 0 vs. H 1 : µ 0 t-statistic t = q X s X X = q 1 q j=1 X j, s 2 X = (q 1) 1 q j=1 (X j X) 2 cv q (α) = critical value of T q 1 : P ( T q 1 > cv q (α)) = α P ( t > cv(α) H 0 ) P ( t > cv(α) H 0, σ 2 1 = = σ 2 q) = α Holds under heavy tails, mixtures of normals (stable, Student-t) Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
9 Choi et al., 2016 Brief review of t-statistic based robust inference Regression: Assume data can be classied in a nite number q of groups that allow asymptotically independent normal inference about the (scalar) parameter of interest, so that ˆ j idn(, vj 2 ) for j = 1,..., q. Time series example: Divide data into q = 4 consecutive blocks, and estimate the model 4 times. Treat ˆ j as observations for the usual t-statistic, and reject a 5% level test if t-statistic is larger than usual critical value for q 1 degrees of freedom. Results in valid inference by small sample result. Exploits information ˆ j idn(, vj 2 ) in an ecient way. Does not rely on single asymptotic model of sampling variability for estimated standard deviation Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
10 Monte Carlo Results Same design as in Andrews (1991): Linear Regression, 5 regressors, 4 nonconstant regressors are independent draws from stationary Gaussian AR(1), as are the disturbances, + heteroskedasticity. T = 128, 5% level test about coefficient of one nonconstant regressor. t-statistic (q) ˆω 2 QA ˆω2 PW ˆω 2 BT (b) ρ Size ρ Size Adjusted Power
11 Choi et al., 2016 The robust test for no predictability It can be shown that under the null hypothesis of = 0, 1 T T t=1 u t 1 0 σ(s)dw (s) = d MN(0, Q), where Q is some (random) variance depending on the volatility process (Q = 1 0 σ(s)2 ds), and therefore 1 T T t=1 sgn(x t 1 )y t = 1 T T t=1 sgn(x t 1 )u t Therefore, the numerator of the Cauchy estimator, ˆδ T := 1 0 σ(s)dw (s) = d MN(0, Q). T sgn(x t 1 )y t (11) t=1 has approximately mixed normal distribution with mean zero and variance Q T. Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
12 Choi et al., 2016 The robust test for no predictability Consider partition of the data into q 2 approximately equal groups G j = {s : (j 1)T/q < s jt/q}, j = 1,..., q. Then we have q subsamples for which we calculate q numerators of Cauchy estimator of (11), ˆδ T,j, j = 1,..., q. Therefore, we have (ˆδ T,1, ˆδ T,2,..., ˆδ T,q ) N (0, diag(t Q 1,..., T Q q )) (12) under assumption about asymptotically independence of δ T,j. Following Ibragimov and Muller (2010), asymptotic Gaussianity then allow us to construct asymptotically valid (conservative) test of level α 83 of H 0 : = 0 against H 1 : 0 by rejecting H 0 when t IV, exceed (1 α/2) percentile of Student t-distribution with q 1 degrees of freedom, where t IV, is constructed as t IV, = q ˆδ/sˆδ (13) with ˆδ = q 1 q j=1 ˆδ T,j and sˆδ = (q 1) 1 q j=1 (ˆδT,j ˆδ) 2. Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
13 Choi et al., 2016 The robust test for no predictability Consider partition of the data into q 2 approximately equal groups G j = {s : (j 1)T/q < s jt/q}, j = 1,..., q. Then we have q subsamples for which we calculate q numerators of Cauchy estimator of (11), ˆδ T,j, j = 1,..., q. Therefore, we have (ˆδ T,1, ˆδ T,2,..., ˆδ T,q ) N (0, diag(t Q 1,..., T Q q )) (12) under assumption about asymptotically independence of δ T,j. Following Ibragimov and Muller (2010), asymptotic Gaussianity then allow us to construct asymptotically valid (conservative) test of level α 83 of H 0 : = 0 against H 1 : 0 by rejecting H 0 when t IV, exceed (1 α/2) percentile of Student t-distribution with q 1 degrees of freedom, where t IV, is constructed as t IV, = q ˆδ/sˆδ (13) with ˆδ = q 1 q j=1 ˆδ T,j and sˆδ = (q 1) 1 q j=1 (ˆδT,j ˆδ) 2. Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
14 We use continuous time model as in Choi et al. (2016) dy t = T X tdt + du t, (14) dx t = κ T X tdt + σ t dv t, (15) ) du t = σ t (dw t + xλ(dt, dx), (16) where V t and W t are Brownian motion with correlation coecient and α = 0. The continuous data are generated to b eobserved at δ-intervals over T years, so that there are δt daily obserations. See Choi et al. (2016). Compare with Bonferroni Q-test of Campbell and Yogo, 2006 (BQ) and restricted likelihood ratio test of Chen and Deo, 2009 (RLRT) R Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
15 Volatility process: Model CNST. Constant volatility : σ 2 t = σ 2 0, σ 0 = 1. Model SB. Structural break in volatility : σ 0 + (σ 1 σ 0 )I(t 4T/5) with σ 0 = 1 and σ 1 = 4. Model GBM. Geometric Brownian motion: dσt 2 = 1 ω 2 2 T σ2 t dt + ω2 T σt 2 dz t, Z t is Brownian motion which correlated with W t with correlation coecient - and ω is set to be 9. Model RS. Regime switching : σ t = σ 0 (1 s t ) + σ 1 s t, where s t be a homogeneous Markov process indicating the current state of the world which is independent on both Y t and X t with transition matrix with the state space {0, 1} P t = ( ) + ( ) ( exp λ ) T t with λ = 60, σ 0 = 1 and σ 1 = 4. s t is initialized by its invariant distribution. Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
16 Finite Sample Property Constant volatility: σ t = σ Table: Sizes of tests κ = 0 κ = 5 κ = CNST OLS BQ RLRT Cauchy RT q= q= q= q= Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
17 Single break in vlatility: σ t = 1 for t [0, 4T/5] and σ t = 4 for t [4T/5, T ] Table: Sizes of tests κ = 0 κ = 5 κ = SB OLS BQ RLRT Cauchy RT q= q= q= q= Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
18 Regime switching model: σ t = σ 0 (1 s t ) + σ 1 s t, s t is homogeneous Markov process independent of both Y and X Table: Sizes of tests κ = 0 κ = 5 κ = RS OLS BQ RLRT Cauchy RT q= q= q= q= Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
19 Geometric Brownian motion: dσ 2 t = ω2 2 σ2 t dt + ωσ 2 t dz t, Z t BM, correlated with W (with -), ω = 9/ T. Table: Sizes of tests κ = 0 κ = 5 κ = GBM OLS BQ RLRT Cauchy RT q= q= q= q= Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
20 OLS:,Bonf.Q:,RLRT:,CauchyRT:,q=4:,q=8:,q=12:,q=16: (a) CNST, κ = 0, T = 5 (b) CNST, κ = 5, T = (c) CNST, κ = 20, T = 5 Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
21 OLS:,Bonf.Q:,RLRT:,CauchyRT:,q=4:,q=8:,q=12:,q=16: (a) CNST, κ = 0, T = 50 (b) CNST, κ = 5, T = (c) CNST, κ = 20, T = 50 Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
22 OLS:,Bonf.Q:,RLRT:,CauchyRT:,q=4:,q=8:,q=12:,q=16: (a) SB, κ = 0, T = 5 (b) SB, κ = 5, T = (c) SB, κ = 20, T = 5 Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
23 OLS:,Bonf.Q:,RLRT:,CauchyRT:,q=4:,q=8:,q=12:,q=16: (a) SB, κ = 0, T = 50 (b) SB, κ = 5, T = (c) SB, κ = 20, T = 50 Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
24 OLS:,Bonf.Q:,RLRT:,CauchyRT:,q=4:,q=8:,q=12:,q=16: (a) GBM, κ = 0, T = 5 (b) GBM, κ = 5, T = (c) GBM, κ = 20, T = 5 Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
25 OLS:,Bonf.Q:,RLRT:,CauchyRT:,q=4:,q=8:,q=12:,q=16: (a) GBM, κ = 0, T = 50 (b) GBM, κ = 5, T = (c) GBM, κ = 20, T = 50 Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
26 OLS:,Bonf.Q:,RLRT:,CauchyRT:,q=4:,q=8:,q=12:,q=16: (a) RS, κ = 0, T = 5 (b) RS, κ = 5, T = (c) RS, κ = 20, T = 5 Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
27 OLS:,Bonf.Q:,RLRT:,CauchyRT:,q=4:,q=8:,q=12:,q=16: (a) RS, κ = 0, T = 50 (b) RS, κ = 5, T = (c) RS, κ = 20, T = 50 Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
28 Concluding remarks New approach to robust inference in predictive regression Based on instrumental variable estimator (Cauchy estimator) which allows the endogeneity between variables The inference then based on splitting the sample and obtaining robust Student t-test The obtaining approach is robust to a wide class of errors: dependence, heterockedasticity, nonstationary volatility and heavy tails Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
29 Thank you for attention! Anton Skrobotov (RANEPA, Innopolis) New robust inference for predictive regressions June 28, / 27
The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2009, Mr. Ruey S. Tsay. Solutions to Final Exam
The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2009, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (42 pts) Answer briefly the following questions. 1. Questions
More informationThe stochastic calculus
Gdansk A schedule of the lecture Stochastic differential equations Ito calculus, Ito process Ornstein - Uhlenbeck (OU) process Heston model Stopping time for OU process Stochastic differential equations
More informationAssicurazioni Generali: An Option Pricing Case with NAGARCH
Assicurazioni Generali: An Option Pricing Case with NAGARCH Assicurazioni Generali: Business Snapshot Find our latest analyses and trade ideas on bsic.it Assicurazioni Generali SpA is an Italy-based insurance
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 informationGraduate School of Business, University of Chicago Business 41202, Spring Quarter 2007, Mr. Ruey S. Tsay. Solutions to Final Exam
Graduate School of Business, University of Chicago Business 41202, Spring Quarter 2007, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (30 pts) Answer briefly the following questions. 1. Suppose that
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 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 informationThe University of Chicago, Booth School of Business Business 41202, Spring Quarter 2017, Mr. Ruey S. Tsay. Solutions to Final Exam
The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2017, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (40 points) Answer briefly the following questions. 1. Describe
More informationLimit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies
Limit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies George Tauchen Duke University Viktor Todorov Northwestern University 2013 Motivation
More informationM.I.T Fall Practice Problems
M.I.T. 15.450-Fall 2010 Sloan School of Management Professor Leonid Kogan Practice Problems 1. Consider a 3-period model with t = 0, 1, 2, 3. There are a stock and a risk-free asset. The initial stock
More informationThe University of Chicago, Booth School of Business Business 41202, Spring Quarter 2010, Mr. Ruey S. Tsay Solutions to Final Exam
The University of Chicago, Booth School of Business Business 410, Spring Quarter 010, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (4 pts) Answer briefly the following questions. 1. Questions 1
More information12. Conditional heteroscedastic models (ARCH) MA6622, Ernesto Mordecki, CityU, HK, 2006.
12. Conditional heteroscedastic models (ARCH) MA6622, Ernesto Mordecki, CityU, HK, 2006. References for this Lecture: Robert F. Engle. Autoregressive Conditional Heteroscedasticity with Estimates of Variance
More informationExperimental Design and Statistics - AGA47A
Experimental Design and Statistics - AGA47A Czech University of Life Sciences in Prague Department of Genetics and Breeding Fall/Winter 2014/2015 Matúš Maciak (@ A 211) Office Hours: M 14:00 15:30 W 15:30
More informationMultiname and Multiscale Default Modeling
Multiname and Multiscale Default Modeling Jean-Pierre Fouque University of California Santa Barbara Joint work with R. Sircar (Princeton) and K. Sølna (UC Irvine) Special Semester on Stochastics with Emphasis
More informationRobust Econometric Inference for Stock Return Predictability
Robust Econometric Inference for Stock Return Predictability Alex Kostakis (MBS), Tassos Magdalinos (Southampton) and Michalis Stamatogiannis (Bath) Alex Kostakis, MBS Marie Curie, Konstanz (Alex Kostakis,
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 informationChapter 5 Univariate time-series analysis. () Chapter 5 Univariate time-series analysis 1 / 29
Chapter 5 Univariate time-series analysis () Chapter 5 Univariate time-series analysis 1 / 29 Time-Series Time-series is a sequence fx 1, x 2,..., x T g or fx t g, t = 1,..., T, where t is an index denoting
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Generating Random Variables and Stochastic Processes Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More 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 informationA1. Relating Level and Slope to Expected Inflation and Output Dynamics
Appendix 1 A1. Relating Level and Slope to Expected Inflation and Output Dynamics This section provides a simple illustrative example to show how the level and slope factors incorporate expectations regarding
More informationARCH and GARCH models
ARCH and GARCH models Fulvio Corsi SNS Pisa 5 Dic 2011 Fulvio Corsi ARCH and () GARCH models SNS Pisa 5 Dic 2011 1 / 21 Asset prices S&P 500 index from 1982 to 2009 1600 1400 1200 1000 800 600 400 200
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 informationQueens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 2017 Instructor: Dr. Sateesh Mane.
Queens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 2017 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 2017 14 Lecture 14 November 15, 2017 Derivation of the
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 informationSTOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL
STOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL YOUNGGEUN YOO Abstract. Ito s lemma is often used in Ito calculus to find the differentials of a stochastic process that depends on time. This paper will introduce
More informationBooth School of Business, University of Chicago Business 41202, Spring Quarter 2014, Mr. Ruey S. Tsay. Solutions to Midterm
Booth School of Business, University of Chicago Business 41202, Spring Quarter 2014, Mr. Ruey S. Tsay Solutions to Midterm Problem A: (30 pts) Answer briefly the following questions. Each question has
More informationPractical example of an Economic Scenario Generator
Practical example of an Economic Scenario Generator Martin Schenk Actuarial & Insurance Solutions SAV 7 March 2014 Agenda Introduction Deterministic vs. stochastic approach Mathematical model Application
More informationBROWNIAN MOTION Antonella Basso, Martina Nardon
BROWNIAN MOTION Antonella Basso, Martina Nardon basso@unive.it, mnardon@unive.it Department of Applied Mathematics University Ca Foscari Venice Brownian motion p. 1 Brownian motion Brownian motion plays
More informationRisk Neutral Valuation
copyright 2012 Christian Fries 1 / 51 Risk Neutral Valuation Christian Fries Version 2.2 http://www.christian-fries.de/finmath April 19-20, 2012 copyright 2012 Christian Fries 2 / 51 Outline Notation Differential
More informationS t d with probability (1 p), where
Stochastic Calculus Week 3 Topics: Towards Black-Scholes Stochastic Processes Brownian Motion Conditional Expectations Continuous-time Martingales Towards Black Scholes Suppose again that S t+δt equals
More informationFinancial Risk Management
Financial Risk Management Professor: Thierry Roncalli Evry University Assistant: Enareta Kurtbegu Evry University Tutorial exercices #4 1 Correlation and copulas 1. The bivariate Gaussian copula is given
More informationDerivation Of The Capital Asset Pricing Model Part I - A Single Source Of Uncertainty
Derivation Of The Capital Asset Pricing Model Part I - A Single Source Of Uncertainty Gary Schurman MB, CFA August, 2012 The Capital Asset Pricing Model CAPM is used to estimate the required rate of return
More informationContinous time models and realized variance: Simulations
Continous time models and realized variance: Simulations Asger Lunde Professor Department of Economics and Business Aarhus University September 26, 2016 Continuous-time Stochastic Process: SDEs Building
More informationInternet Appendix for Asymmetry in Stock Comovements: An Entropy Approach
Internet Appendix for Asymmetry in Stock Comovements: An Entropy Approach Lei Jiang Tsinghua University Ke Wu Renmin University of China Guofu Zhou Washington University in St. Louis August 2017 Jiang,
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 informationVolatility. Roberto Renò. 2 March 2010 / Scuola Normale Superiore. Dipartimento di Economia Politica Università di Siena
Dipartimento di Economia Politica Università di Siena 2 March 2010 / Scuola Normale Superiore What is? The definition of volatility may vary wildly around the idea of the standard deviation of price movements
More informationConditional Density Method in the Computation of the Delta with Application to Power Market
Conditional Density Method in the Computation of the Delta with Application to Power Market Asma Khedher Centre of Mathematics for Applications Department of Mathematics University of Oslo A joint work
More informationLecture 3. Sergei Fedotov Introduction to Financial Mathematics. Sergei Fedotov (University of Manchester) / 6
Lecture 3 Sergei Fedotov 091 - Introduction to Financial Mathematics Sergei Fedotov (University of Manchester) 091 010 1 / 6 Lecture 3 1 Distribution for lns(t) Solution to Stochastic Differential Equation
More 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 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 informationRMSC 4005 Stochastic Calculus for Finance and Risk. 1 Exercises. (c) Let X = {X n } n=0 be a {F n }-supermartingale. Show that.
1. EXERCISES RMSC 45 Stochastic Calculus for Finance and Risk Exercises 1 Exercises 1. (a) Let X = {X n } n= be a {F n }-martingale. Show that E(X n ) = E(X ) n N (b) Let X = {X n } n= be a {F n }-submartingale.
More informationTesting for a Unit Root with Near-Integrated Volatility
Testing for a Unit Root with Near-Integrated Volatility H. Peter Boswijk Department of Quantitative Economics, University of Amsterdam y January Abstract This paper considers tests for a unit root when
More informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Mike Giles (Oxford) Monte Carlo methods 2 1 / 24 Lecture outline
More informationHelp Session 2. David Sovich. Washington University in St. Louis
Help Session 2 David Sovich Washington University in St. Louis TODAY S AGENDA Today we will cover the Change of Numeraire toolkit We will go over the Fundamental Theorem of Asset Pricing as well EXISTENCE
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 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 informationIMPA Commodities Course : Forward Price Models
IMPA Commodities Course : Forward Price Models Sebastian Jaimungal sebastian.jaimungal@utoronto.ca Department of Statistics and Mathematical Finance Program, University of Toronto, Toronto, Canada http://www.utstat.utoronto.ca/sjaimung
More informationModeling via Stochastic Processes in Finance
Modeling via Stochastic Processes in Finance Dimbinirina Ramarimbahoaka Department of Mathematics and Statistics University of Calgary AMAT 621 - Fall 2012 October 15, 2012 Question: What are appropriate
More informationOption pricing in the stochastic volatility model of Barndorff-Nielsen and Shephard
Option pricing in the stochastic volatility model of Barndorff-Nielsen and Shephard Indifference pricing and the minimal entropy martingale measure Fred Espen Benth Centre of Mathematics for Applications
More informationSTOCHASTIC INTEGRALS
Stat 391/FinMath 346 Lecture 8 STOCHASTIC INTEGRALS X t = CONTINUOUS PROCESS θ t = PORTFOLIO: #X t HELD AT t { St : STOCK PRICE M t : MG W t : BROWNIAN MOTION DISCRETE TIME: = t < t 1
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 informationSupplementary online material to Information tradeoffs in dynamic financial markets
Supplementary online material to Information tradeoffs in dynamic financial markets Efstathios Avdis University of Alberta, Canada 1. The value of information in continuous time In this document I address
More informationPortability, salary and asset price risk: a continuous-time expected utility comparison of DB and DC pension plans
Portability, salary and asset price risk: a continuous-time expected utility comparison of DB and DC pension plans An Chen University of Ulm joint with Filip Uzelac (University of Bonn) Seminar at SWUFE,
More informationLecture 17. The model is parametrized by the time period, δt, and three fixed constant parameters, v, σ and the riskless rate r.
Lecture 7 Overture to continuous models Before rigorously deriving the acclaimed Black-Scholes pricing formula for the value of a European option, we developed a substantial body of material, in continuous
More informationThe rth moment of a real-valued random variable X with density f(x) is. x r f(x) dx
1 Cumulants 1.1 Definition The rth moment of a real-valued random variable X with density f(x) is µ r = E(X r ) = x r f(x) dx for integer r = 0, 1,.... The value is assumed to be finite. Provided that
More informationθ(t ) = T f(0, T ) + σ2 T
1 Derivatives Pricing and Financial Modelling Andrew Cairns: room M3.08 E-mail: A.Cairns@ma.hw.ac.uk Tutorial 10 1. (Ho-Lee) Let X(T ) = T 0 W t dt. (a) What is the distribution of X(T )? (b) Find E[exp(
More informationStochastic Differential Equations in Finance and Monte Carlo Simulations
Stochastic Differential Equations in Finance and Department of Statistics and Modelling Science University of Strathclyde Glasgow, G1 1XH China 2009 Outline Stochastic Modelling in Asset Prices 1 Stochastic
More informationA model of stock price movements
... A model of stock price movements Johan Gudmundsson Thesis submitted for the degree of Master of Science 60 ECTS Master Thesis Supervised by Sven Åberg. Department of Physics Division of Mathematical
More information1. You are given the following information about a stationary AR(2) model:
Fall 2003 Society of Actuaries **BEGINNING OF EXAMINATION** 1. You are given the following information about a stationary AR(2) model: (i) ρ 1 = 05. (ii) ρ 2 = 01. Determine φ 2. (A) 0.2 (B) 0.1 (C) 0.4
More informationOrnstein-Uhlenbeck Theory
Beatrice Byukusenge Department of Technomathematics Lappeenranta University of technology January 31, 2012 Definition of a stochastic process Let (Ω,F,P) be a probability space. A stochastic process is
More informationDependence Structure and Extreme Comovements in International Equity and Bond Markets
Dependence Structure and Extreme Comovements in International Equity and Bond Markets René Garcia Edhec Business School, Université de Montréal, CIRANO and CIREQ Georges Tsafack Suffolk University Measuring
More informationStatistical methods for financial models driven by Lévy processes
Statistical methods for financial models driven by Lévy processes José Enrique Figueroa-López Department of Statistics, Purdue University PASI Centro de Investigación en Matemátics (CIMAT) Guanajuato,
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 informationStochastic modelling of electricity markets Pricing Forwards and Swaps
Stochastic modelling of electricity markets Pricing Forwards and Swaps Jhonny Gonzalez School of Mathematics The University of Manchester Magical books project August 23, 2012 Clip for this slide Pricing
More informationGeneralized Multi-Factor Commodity Spot Price Modeling through Dynamic Cournot Resource Extraction Models
Generalized Multi-Factor Commodity Spot Price Modeling through Dynamic Cournot Resource Extraction Models Bilkan Erkmen (joint work with Michael Coulon) Workshop on Stochastic Games, Equilibrium, and Applications
More informationWeak Convergence to Stochastic Integrals
Weak Convergence to Stochastic Integrals Zhengyan Lin Zhejiang University Join work with Hanchao Wang Outline 1 Introduction 2 Convergence to Stochastic Integral Driven by Brownian Motion 3 Convergence
More informationValue at Risk Ch.12. PAK Study Manual
Value at Risk Ch.12 Related Learning Objectives 3a) Apply and construct risk metrics to quantify major types of risk exposure such as market risk, credit risk, liquidity risk, regulatory risk etc., and
More informationLecture 11: Stochastic Volatility Models Cont.
E4718 Spring 008: Derman: Lecture 11:Stochastic Volatility Models Cont. Page 1 of 8 Lecture 11: Stochastic Volatility Models Cont. E4718 Spring 008: Derman: Lecture 11:Stochastic Volatility Models Cont.
More information25857 Interest Rate Modelling
25857 Interest Rate Modelling UTS Business School University of Technology Sydney Chapter 19. Allowing for Stochastic Interest Rates in the Black-Scholes Model May 15, 2014 1/33 Chapter 19. Allowing for
More informationRobust Econometric Inference for Stock Return Predictability
Robust Econometric Inference for Stock Return Predictability Alex Kostakis (MBS), Tassos Magdalinos (Southampton) and Michalis Stamatogiannis (Bath) Alex Kostakis, MBS 2nd ISNPS, Cadiz (Alex Kostakis,
More information( ) since this is the benefit of buying the asset at the strike price rather
Review of some financial models for MAT 483 Parity and Other Option Relationships The basic parity relationship for European options with the same strike price and the same time to expiration is: C( KT
More informationBooth School of Business, University of Chicago Business 41202, Spring Quarter 2010, Mr. Ruey S. Tsay. Solutions to Midterm
Booth School of Business, University of Chicago Business 41202, Spring Quarter 2010, Mr. Ruey S. Tsay Solutions to Midterm Problem A: (30 pts) Answer briefly the following questions. Each question has
More informationLecture 5a: ARCH Models
Lecture 5a: ARCH Models 1 2 Big Picture 1. We use ARMA model for the conditional mean 2. We use ARCH model for the conditional variance 3. ARMA and ARCH model can be used together to describe both conditional
More informationLecture 5. Predictability. Traditional Views of Market Efficiency ( )
Lecture 5 Predictability Traditional Views of Market Efficiency (1960-1970) CAPM is a good measure of risk Returns are close to unpredictable (a) Stock, bond and foreign exchange changes are not predictable
More informationOptimal Placement of a Small Order Under a Diffusive Limit Order Book (LOB) Model
Optimal Placement of a Small Order Under a Diffusive Limit Order Book (LOB) Model José E. Figueroa-López Department of Mathematics Washington University in St. Louis INFORMS National Meeting Houston, TX
More informationOn modelling of electricity spot price
, Rüdiger Kiesel and Fred Espen Benth Institute of Energy Trading and Financial Services University of Duisburg-Essen Centre of Mathematics for Applications, University of Oslo 25. August 2010 Introduction
More informationStochastic Calculus, Application of Real Analysis in Finance
, Application of Real Analysis in Finance Workshop for Young Mathematicians in Korea Seungkyu Lee Pohang University of Science and Technology August 4th, 2010 Contents 1 BINOMIAL ASSET PRICING MODEL Contents
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 informationCredit Risk : Firm Value Model
Credit Risk : Firm Value Model Prof. Dr. Svetlozar Rachev Institute for Statistics and Mathematical Economics University of Karlsruhe and Karlsruhe Institute of Technology (KIT) Prof. Dr. Svetlozar Rachev
More informationA Robust Test for Normality
A Robust Test for Normality Liangjun Su Guanghua School of Management, Peking University Ye Chen Guanghua School of Management, Peking University Halbert White Department of Economics, UCSD March 11, 2006
More informationBeyond the Black-Scholes-Merton model
Econophysics Lecture Leiden, November 5, 2009 Overview 1 Limitations of the Black-Scholes model 2 3 4 Limitations of the Black-Scholes model Black-Scholes model Good news: it is a nice, well-behaved model
More informationSlides 4. Matthieu Gomez Fall 2017
Slides 4 Matthieu Gomez Fall 2017 How to Compute Stationary Distribution of a Diffusion? Kolmogorov Forward Take a diffusion process dx t = µ(x t )dt + σ(x t )dz t How does the density of x t evolves?
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 informationChapter 14. The Multi-Underlying Black-Scholes Model and Correlation
Chapter 4 The Multi-Underlying Black-Scholes Model and Correlation So far we have discussed single asset options, the payoff function depended only on one underlying. Now we want to allow multiple underlyings.
More informationLocal Volatility Dynamic Models
René Carmona Bendheim Center for Finance Department of Operations Research & Financial Engineering Princeton University Columbia November 9, 27 Contents Joint work with Sergey Nadtochyi Motivation 1 Understanding
More informationAdvanced topics in continuous time finance
Based on readings of Prof. Kerry E. Back on the IAS in Vienna, October 21. Advanced topics in continuous time finance Mag. Martin Vonwald (martin@voni.at) November 21 Contents 1 Introduction 4 1.1 Martingale.....................................
More information1 The continuous time limit
Derivative Securities, Courant Institute, Fall 2008 http://www.math.nyu.edu/faculty/goodman/teaching/derivsec08/index.html Jonathan Goodman and Keith Lewis Supplementary notes and comments, Section 3 1
More informationSaddlepoint Approximation Methods for Pricing. Financial Options on Discrete Realized Variance
Saddlepoint Approximation Methods for Pricing Financial Options on Discrete Realized Variance Yue Kuen KWOK Department of Mathematics Hong Kong University of Science and Technology Hong Kong * This is
More informationStochastic Volatility
Stochastic Volatility A Gentle Introduction Fredrik Armerin Department of Mathematics Royal Institute of Technology, Stockholm, Sweden Contents 1 Introduction 2 1.1 Volatility................................
More informationIlliquidity, Credit risk and Merton s model
Illiquidity, Credit risk and Merton s model (joint work with J. Dong and L. Korobenko) A. Deniz Sezer University of Calgary April 28, 2016 Merton s model of corporate debt A corporate bond is a contingent
More informationAn Analytical Approximation for Pricing VWAP Options
.... An Analytical Approximation for Pricing VWAP Options Hideharu Funahashi and Masaaki Kijima Graduate School of Social Sciences, Tokyo Metropolitan University September 4, 215 Kijima (TMU Pricing of
More informationStochastic Dynamical Systems and SDE s. An Informal Introduction
Stochastic Dynamical Systems and SDE s An Informal Introduction Olav Kallenberg Graduate Student Seminar, April 18, 2012 1 / 33 2 / 33 Simple recursion: Deterministic system, discrete time x n+1 = f (x
More informationdt+ ρσ 2 1 ρ2 σ 2 κ i and that A is a rather lengthy expression that we may or may not need. (Brigo & Mercurio Lemma Thm , p. 135.
A 2D Gaussian model (akin to Brigo & Mercurio Section 4.2) Suppose where ( κ1 0 dx(t) = 0 κ 2 r(t) = δ 0 +X 1 (t)+x 2 (t) )( X1 (t) X 2 (t) ) ( σ1 0 dt+ ρσ 2 1 ρ2 σ 2 )( dw Q 1 (t) dw Q 2 (t) ) In this
More informationA potentially useful approach to model nonlinearities in time series is to assume different behavior (structural break) in different subsamples
1.3 Regime switching models A potentially useful approach to model nonlinearities in time series is to assume different behavior (structural break) in different subsamples (or regimes). If the dates, the
More informationContinuous Time Finance. Tomas Björk
Continuous Time Finance Tomas Björk 1 II Stochastic Calculus Tomas Björk 2 Typical Setup Take as given the market price process, S(t), of some underlying asset. S(t) = price, at t, per unit of underlying
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 informationValuation of derivative assets Lecture 8
Valuation of derivative assets Lecture 8 Magnus Wiktorsson September 27, 2018 Magnus Wiktorsson L8 September 27, 2018 1 / 14 The risk neutral valuation formula Let X be contingent claim with maturity T.
More informationResults for option pricing
Results for option pricing [o,v,b]=optimal(rand(1,100000 Estimators = 0.4619 0.4617 0.4618 0.4613 0.4619 o = 0.46151 % best linear combination (true value=0.46150 v = 1.1183e-005 %variance per uniform
More informationOptimal robust bounds for variance options and asymptotically extreme models
Optimal robust bounds for variance options and asymptotically extreme models Alexander Cox 1 Jiajie Wang 2 1 University of Bath 2 Università di Roma La Sapienza Advances in Financial Mathematics, 9th January,
More informationThe University of Chicago, Booth School of Business Business 41202, Spring Quarter 2012, Mr. Ruey S. Tsay. Solutions to Final Exam
The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2012, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (40 points) Answer briefly the following questions. 1. Consider
More information