CH Sean Han QF, NTHU, Taiwan BFS2010. (Joint work with T.-Y. Chen and W.-H. Liu)

Similar documents
State-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter

Lecture 10 Estimating Nonlinear Regression Models

OBJECTIVES OF TIME SERIES ANALYSIS

Regression with Time Series Data

Unit Root Time Series. Univariate random walk

Linear Combinations of Volatility Forecasts for the WIG20 and Polish Exchange Rates

Econ Autocorrelation. Sanjaya DeSilva

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t

Introduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate.

GMM - Generalized Method of Moments

An introduction to the theory of SDDP algorithm

Linear Gaussian State Space Models

(a) Set up the least squares estimation procedure for this problem, which will consist in minimizing the sum of squared residuals. 2 t.

Time series Decomposition method

Chapter 5. Heterocedastic Models. Introduction to time series (2008) 1

PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD

Vehicle Arrival Models : Headway

Financial Econometrics Jeffrey R. Russell Midterm Winter 2009 SOLUTIONS

What Ties Return Volatilities to Price Valuations and Fundamentals? On-Line Appendix

An Introduction to Malliavin calculus and its applications

Distribution of Estimates

Volatility. Many economic series, and most financial series, display conditional volatility

Smoothing. Backward smoother: At any give T, replace the observation yt by a combination of observations at & before T

Introduction to Probability and Statistics Slides 4 Chapter 4

Modeling and Forecasting Volatility Autoregressive Conditional Heteroskedasticity Models. Economic Forecasting Anthony Tay Slide 1

Recursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H.

Bias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé

Wisconsin Unemployment Rate Forecast Revisited

Lecture 33: November 29

Option pricing and implied volatilities in a 2-hypergeometric stochastic volatility model

Vectorautoregressive Model and Cointegration Analysis. Time Series Analysis Dr. Sevtap Kestel 1

Distribution of Least Squares

ESTIMATION OF DYNAMIC PANEL DATA MODELS WHEN REGRESSION COEFFICIENTS AND INDIVIDUAL EFFECTS ARE TIME-VARYING

L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms

5. NONLINEAR MODELS [1] Nonlinear (NL) Regression Models

Chapter 2. Models, Censoring, and Likelihood for Failure-Time Data

Financial Econometrics Introduction to Realized Variance

EMS SCM joint meeting. On stochastic partial differential equations of parabolic type

REVIEW OF MAXIMUM LIKELIHOOD ESTIMATION

Licenciatura de ADE y Licenciatura conjunta Derecho y ADE. Hoja de ejercicios 2 PARTE A

Optimal Investment under Dynamic Risk Constraints and Partial Information

Simulation of BSDEs and. Wiener Chaos Expansions

1. Consider a pure-exchange economy with stochastic endowments. The state of the economy

Two Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017

Simulation of BSDEs and. Wiener Chaos Expansions

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds

AMartingaleApproachforFractionalBrownian Motions and Related Path Dependent PDEs

Estimation for Parameters in Partially Observed Linear Stochastic System

Institute for Mathematical Methods in Economics. University of Technology Vienna. Singapore, May Manfred Deistler

Utility maximization in incomplete markets

Testing the Random Walk Model. i.i.d. ( ) r

(10) (a) Derive and plot the spectrum of y. Discuss how the seasonality in the process is evident in spectrum.

Vanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law

Properties of Autocorrelated Processes Economics 30331

Sample Autocorrelations for Financial Time Series Models. Richard A. Davis Colorado State University Thomas Mikosch University of Copenhagen

Testing for a Single Factor Model in the Multivariate State Space Framework

Optimal Investment Strategy Insurance Company

Robust estimation based on the first- and third-moment restrictions of the power transformation model

HPCFinance research project 8

4.1 Other Interpretations of Ridge Regression

Exponential Smoothing

The General Linear Test in the Ridge Regression

Probabilistic Robotics

Outline. lse-logo. Outline. Outline. 1 Wald Test. 2 The Likelihood Ratio Test. 3 Lagrange Multiplier Tests

Generalized Least Squares

2.160 System Identification, Estimation, and Learning. Lecture Notes No. 8. March 6, 2006

J. Martin van Zyl Department of Mathematical Statistics and Actuarial Science, University of the Free State, PO Box 339, Bloemfontein, South Africa

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)

Pure Jump Lévy Processes and Self-decomposability in Financial Modelling

12: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME. Σ j =

Reliability of Technical Systems

A Bayesian Approach to Spectral Analysis

Lecture 5. Time series: ECM. Bernardina Algieri Department Economics, Statistics and Finance

Optimal Portfolio under Fractional Stochastic Environment

References are appeared in the last slide. Last update: (1393/08/19)

Understanding the asymptotic behaviour of empirical Bayes methods

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing

Essential Microeconomics : OPTIMAL CONTROL 1. Consider the following class of optimization problems

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.

Exponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits

Supplementary Document

Stochastic Model for Cancer Cell Growth through Single Forward Mutation

Structural Break Detection in Time Series Models

Answers to Exercises in Chapter 7 - Correlation Functions

On a Discrete-In-Time Order Level Inventory Model for Items with Random Deterioration

ON THE NUMBER OF FAMILIES OF BRANCHING PROCESSES WITH IMMIGRATION WITH FAMILY SIZES WITHIN RANDOM INTERVAL

Chapter 14 Wiener Processes and Itô s Lemma. Options, Futures, and Other Derivatives, 9th Edition, Copyright John C. Hull

Cointegration and Implications for Forecasting

On a Fractional Stochastic Landau-Ginzburg Equation

Types of Exponential Smoothing Methods. Simple Exponential Smoothing. Simple Exponential Smoothing

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

Basic notions of probability theory (Part 2)

Sequential Importance Resampling (SIR) Particle Filter

3.1 More on model selection

Maintenance Models. Prof. Robert C. Leachman IEOR 130, Methods of Manufacturing Improvement Spring, 2011

arxiv: v1 [math.ca] 15 Nov 2016

6. COMPUTATION OF CENTILES AND Z-SCORES FOR VELOCITIES BASED ON WEIGHT, LENGTH AND HEAD CIRCUMFERENCE

DEPARTMENT OF ECONOMICS AND FINANCE COLLEGE OF BUSINESS AND ECONOMICS UNIVERSITY OF CANTERBURY CHRISTCHURCH, NEW ZEALAND

Transcription:

CH Sean Han QF, NTHU, Taiwan BFS2010 (Join work wih T.-Y. Chen and W.-H. Liu)

Risk Managemen in Pracice: Value a Risk (VaR) / Condiional Value a Risk (CVaR) Volailiy Esimaion: Correced Fourier Transform Mehod Esimae Exreme Probabiliy by Efficien Imporance Sampling Backesing for VaR Esimaion

Le r( ) be an asse reurn a ime. Is α 100% VaR, denoed by VaR α, is defined by he 1 α 100% percenile of r( ). Tha is, ( ) ( ( ) ) 1 P r VaR α = α Tha is a risk conroller has a α 100% confidence ha he asse price will no drop below VaR α a ime.

Mahemaically, i is no a coheren risk measure* because i doesn' saisfy he risk diversificaion principal. Insead, CVaR does!. Pracically, VaR is commonly required by financial regulaions (Basel Ⅱ Accord). * Arzner P., F. Delbaen, J.-M. Eber, and D. Heah, Coheren Measures of Risk, Mahemaical Finance, 9 (1999): 203-28.

Riskmerics: normal assumpion under EWMA model. Hisorical Simulaion: generae scenarios Model Dependen Approach: Discree-Time vs. Coninuous-Time Models

Assume a diffusion process du d dw ( ) = µ ( ) + σ ( ), Task: Given reurn ime series he volailiy σ ( ). u ( ), esimae * Malliavin and Mancino(2002, 2009)

Compue he Fourier coefficiens of du by Then, 1 2π a ( ) ( ) 0 du = du, 2π 0 1 2π a ( ) cos ( ) ( ) k du = k du, π 0 1 2π b ( ) sin ( ) ( ) k du = k du. π 0 ( ) ( ) bk du ak du u ( ) = a0 + cos( k ) + sin ( k ). k = 1 k k

Fourier coefficiens of variance σ N k ( 2 ) π *( ) * ( ) *( ) * a lim ( ) k σ = as du as+ k du + bs du bs k du +, N 2N + 1 N k ( 2 ) π *( ) * ( ) *( ) * b lim ( ) k σ = as du bs+ k du bs du as k du +, N 2N + 1 where n 0 s= N s= N ( ) is any posiive ineger so ha N 2 2 2 N k k k = 0 ( ) = a ( ) ( k ) + a ( ) ( k ) σ σ cos σ sin. 2,

2 Reconsruc he ime series variance σ ( ). Finally, σ 2 ( ) 2 is an approximaion of σ ( ) as N N approaches infiniy, which can be given by classical Fourier-Fejer inversion formula. σ = lim σ in prob. ( ) ( ) 2 2 N N

We add a funcion ino he final compuaion of ime series variance in order o smooh i. N 2 2 2 k k N k 0 2 ( ) = ( k ) a ( ) ( k ) + b ( ) ( k ) σ lim ϕ δ σ cos σ sin, = ϕ x = ( ) sin x ( ) where 2 is a funcion in order o x smooh he rajecory and δ is a smoohing parameer. Reno (2008) alers he boundary effec in he Fourier ransform mehod.

Idea: (Nonlinear) Leas Squares Mehod for firs-order correcion r σ δε (( a by ˆ ) ) ex p + 2 δ ε. Then by MLE o regress ou a and 2 r = a + byˆ + ε 2 ln ln. δ b

Assuming ha he driving volailiy process is governed by he Ornsein-Uhlenbeck process, dy = α m Y d + βdw ( ). (1) We use he correced esimaor a + byˆ o furher esimae model parameers αβ of Yˆ by means of maximum likelihood mehod.,,m ( )

For a given se of observaions he likelihood funcion is Y1, Y2,..., YN ( ) N 1 1 ( ) ( ( ) ) 2 L αβ,, m = exp Y 2 2 + 1 αm + 1 α Y, = 1 2πβ 2β where denoes he lengh of discreized ime inerval.

By maximizing he righ hand side over he parameers ( αβ,,m), we obain he following maximum likelihood esimaors N N 1 N 1 Y Y ( N 1) YY + 1 1 = 2 = 1 = 1 ˆ α = 1, N 1 2 N 1 2 Y ( N 1) Y = 1 = 1 N 1 1 ˆ= β Y+ 1 ( αm + ( 1 α ) Y), N = 1 N N 1 N 1 N 1 2 Y Y Y YY + 1 1 = 2 = 1 = 1 = 1 ˆm=, 1 2 ˆ α N N 1 2 Y ( N 1) Y = 1 = 1 2 (3) (4) (5)

Le he sochasic volailiy model ds ( ) = µ Sd + exp Y 2 SdW1, dy ( - ) = α m Y d + βdw2. To empirically es our price correcion scheme, we se model parameers as follows: µ = 0.01, S = 50, Y = 2, m= 2, α=5, β = 1, 0 0 and wih he discreizaion lengh so as o generae volailiy series and asse price series. S =1/5000 σ = exp( Y 2)

Two crieria are used for performance comparison: Mean squared errors (MSE) and Maximum absolue errors (MAE). Comparison resuls are shown below: Fourier mehod Correced Fourier mehod Mean squared error 0.0324 0.0025 Maximum absolue error 0.3504 0.1563

Given a Markovian dynamic model of an asse price, is reurn process is r = ln S S. S ( ) Given a loss hreshold D, he exreme probabiliy is defined by P( 0, S ) ( ) 0; D = E I rt D S0. Noe: solve VaR α from P ( 0, S ) 0; VaR α = 1 α. =. (Expeced Shorfall) CVaR E rt rt VaR α T T 0

Given he Black-Scholes Model under measure P, choose dp Q ( h) exp hw saisfing Then = T = T dp E [ S ] exp 0 ( ). T = S D µ D h =, σ σt he exreme probabiliy P( 0, S ) ( ) ( ) 0 = E I rt D QT h S0. 2 h T becomes (6) The unbiased imporance sampling esimaor N 1 ( i) i of P( 0, S ) is I r D Q h. (7) 0 N i = 1 ( ) ( ) T T ( ) 2

Theorem: Under he Black-Scholes model, he proposed imporance sampling esimaor is asympoically opimal or efficien under some scaling scenarios in ime and space. Proof: The variance rae of he proposed imporance sampling scheme approaches zero.

Sochasic volailiy model: ds = µ Sd + σsdw1 σ exp( 2) = Y dy = α( m -Y ) d + βdw 2 Ergodic propery of he averaged variance process 1 T T 0.. 2 as 2 ε f ( Y ) d σ,for ε 0 where ε denoes a small ime scale and Y ε denoes a fas mean-revering process.

E.g. 1 2ν α =, β =. ε ε So ha he imporance sampling as foremenioned can be applied. CVaR esimaion can be easily solved.

c ρ = VaR 99% Given model parameers of sochasic volailiy: m 5, α 5, β 1, S0 50, = = = = = Y 0 3, µ =0 ρ c = VaR 99% CVaR N. Approx. IS 0.8-0.0339-0.0347-0.0386 (7.5073E-05) 0.4-0.0335-0.0343-0.0378 (7.3833E-05) 0-0.0323-0.0331-0.0367 (7.2498E-05) -0.4-0.0317-0.0325-0.0366 (7.2739E-05) -0.8-0.0310-0.0319-0.0351 (6.9643E-05)

Daa sample period: 1998.01.05-2009.07.24

Daa sample period: 2005.01.03-2009.07.24

Daa sample period: 2005.01.03-2009.07.24 RiskMerics Significance 1% Significance 5% LRuc Rejec VaR Model LRuc Rejec VaR Model LRind Rejec VaR Model LRind Don' Rejec VaR Model LRcc Rejec VaR Model LRcc Rejec VaR Model Hisorical Simulaion Significance 1% Significance 5% LRuc Rejec VaR Model LRuc Rejec VaR Model LRind Don' Rejec VaR Model LRind Don' Rejec VaR Model LRcc Rejec VaR Model LRcc Rejec VaR Model SV Significance 1% Significance 5% LRuc Don' Rejec VaR Model LRuc Rejec VaR Model LRind Don' Rejec VaR Model LRind Don' Rejec VaR Model LRcc Don' Rejec VaR Model LRcc Rejec VaR Model GARCH(1,1) Significance 1% Significance 5% LRuc Rejec VaR Model LRuc Rejec VaR Model LRind Don' Rejec VaR Model LRind Rejec VaR Model LRcc Rejec VaR Model LRcc Rejec VaR Model

Remove boundary effec of Fourier ransform mehod for volailiy esimaion. (efficien) imporance sampling mehods are invesigaed. VaR backesing resuls for FX and equiy daa.

Thank You