V Time Varying Covariance and Correlation

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

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

System of Linear Differential Equations

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

OBJECTIVES OF TIME SERIES ANALYSIS

Financial Econometrics Jeffrey R. Russell Midterm Winter 2009 SOLUTIONS

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

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles

Generalized Least Squares

Robert Kollmann. 6 September 2017

Exponential Smoothing

Notes on Kalman Filtering

10. State Space Methods

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

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

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

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

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

Comparing Means: t-tests for One Sample & Two Related Samples

Summer Term Albert-Ludwigs-Universität Freiburg Empirische Forschung und Okonometrie. Time Series Analysis

Concourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations

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

Affine term structure models

I. Return Calculations (20 pts, 4 points each)

Linear Gaussian State Space Models

GMM - Generalized Method of Moments

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

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

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

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

Math 10B: Mock Mid II. April 13, 2016

Two Coupled Oscillators / Normal Modes

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds

Math 334 Fall 2011 Homework 11 Solutions

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

DEPARTMENT OF STATISTICS

Let us start with a two dimensional case. We consider a vector ( x,

THE UNIVERSITY OF TEXAS AT AUSTIN McCombs School of Business

Vehicle Arrival Models : Headway

1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx.

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

20. Applications of the Genetic-Drift Model

Unit Root Time Series. Univariate random walk

Christos Papadimitriou & Luca Trevisan November 22, 2016

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

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

ECON 482 / WH Hong Time Series Data Analysis 1. The Nature of Time Series Data. Example of time series data (inflation and unemployment rates)

Asymmetry and Leverage in Conditional Volatility Models*

Laplace transfom: t-translation rule , Haynes Miller and Jeremy Orloff

Vector autoregression VAR. Case 1

13.3 Term structure models

Regression with Time Series Data

Chapter 6. Systems of First Order Linear Differential Equations

Announcements: Warm-up Exercise:

EE363 homework 1 solutions

US AND LATIN AMERICAN STOCK MARKET LINKAGES

Forward guidance. Fed funds target during /15/2017

Granger Causality Among Pre-Crisis East Asian Exchange Rates. (Running Title: Granger Causality Among Pre-Crisis East Asian Exchange Rates)

Distribution of Estimates

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

Estimation Uncertainty

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION

Problem Set 5. Graduate Macro II, Spring 2017 The University of Notre Dame Professor Sims

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

Distribution of Least Squares

ACE 562 Fall Lecture 8: The Simple Linear Regression Model: R 2, Reporting the Results and Prediction. by Professor Scott H.

Maximum Likelihood Estimation of Time-Varying Loadings in High-Dimensional Factor Models. Jakob Guldbæk Mikkelsen, Eric Hillebrand and Giovanni Urga

Matlab and Python programming: how to get started

Model Reduction for Dynamical Systems Lecture 6

14 Autoregressive Moving Average Models

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

Exponentially Weighted Moving Average (EWMA) Chart Based on Six Delta Initiatives

CHAPTER 12 DIRECT CURRENT CIRCUITS

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

Georey E. Hinton. University oftoronto. Technical Report CRG-TR February 22, Abstract

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

t + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that

( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively:

Advanced time-series analysis (University of Lund, Economic History Department)

Time series Decomposition method

Modeling the Volatility of Shanghai Composite Index

Explaining Total Factor Productivity. Ulrich Kohli University of Geneva December 2015

Some Basic Information about M-S-D Systems

Dynamic models for largedimensional. Yields on U.S. Treasury securities (3 months to 10 years) y t

Lecture 3: Exponential Smoothing

Lecture #31, 32: The Ornstein-Uhlenbeck Process as a Model of Volatility

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

Chapter 15. Time Series: Descriptive Analyses, Models, and Forecasting

Time series model fitting via Kalman smoothing and EM estimation in TimeModels.jl

Why is Chinese Provincial Output Diverging? Joakim Westerlund, University of Gothenburg David Edgerton, Lund University Sonja Opper, Lund University

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

0.1 MAXIMUM LIKELIHOOD ESTIMATION EXPLAINED

Richard A. Davis Colorado State University Bojan Basrak Eurandom Thomas Mikosch University of Groningen

Macroeconomic Theory Ph.D. Qualifying Examination Fall 2005 ANSWER EACH PART IN A SEPARATE BLUE BOOK. PART ONE: ANSWER IN BOOK 1 WEIGHT 1/3

Solutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore

The Real Exchange Rate, Real Interest Rates, and the Risk Premium. Charles Engel University of Wisconsin

dy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page

Økonomisk Kandidateksamen 2005(II) Econometrics 2. Solution

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

Tourism forecasting using conditional volatility models

Transcription:

V Time Varying Covariance and Correlaion DEFINITION OF CONDITIONAL CORRELATIONS. ARE THEY TIME VARYING? WHY DO WE NEED THEM? FACTOR MODELS DYNAMIC CONDITIONAL CORRELATIONS Are correlaions/covariances ime varying? Yes Correlaions are ime varying Derivaive prices of correlaion sensiive producs imply changes. Derivaives on correlaion now are raded. Time series esimaes change. There are many varieies. 1

How do we define ime varying correlaions? Recall ha for a single asse, we defined he condiional variance as he variance of he unpredicable par. Tha is, if y hen he condiional 2 variance is given by h E F 1. Hence if here is predicabiliy in he mean, we don include ha variaion in he condiional variance, we remove i firs. The same is rue for he condiional variance covariance marix. We define he condiional variance covariance marix for he par of y ha is no predicable. The condiional variance covariance marix is given by Eεε F 1 where ε y μ If here are no dynamics in he mean hen y 2

If mean reurns are no predicable hen y = and a wo dimensional covariance marix looks like: 2 r 1, 1, 1,2, Eyy F 1 E r1, r2, F 1 2 r 2, 1,2, 2, Condiional Covariance: E r r 1, 2, F 1 Condiional Variances CONDITIONAL CORRELATIONS Condiional correlaions are hen given by 1,2, 1,2, 1, 2, E 1 r1, r2, 2 2 1, 2, E r E r 1 1 3

ESTIMATION HISTORICAL CORRELATIONS i, j, s r r i, m j, m m1 s s 2 2 ri, mrj, m m1 m1 Exponenial smooher: q r r 1 q i, j, i, 1 j, 1 i, j, 1 Then i, j, q q i, j, q ii,, j, j, Use he same value for for each series 4

.9 100 day hisorical correlaions beween AXP and GE.8.7.6.5.4.3.2.1.0 94 95 96 97 98 99 00 01 02 03 04 C100_AXP_GE MULTIVARIATE MODELS 5

SOME MODELS ONE FACTOR MODEL ONE FACTOR MODEL WITH GARCH IDEOSYNCRATIC ERRORS MANY FACTOR MODEL MULTIVARIATE GARCH DYNAMIC CONDITIONAL CORRELATION ONE FACTOR ARCH One facor model such as CAPM There is one marke facor wih fixed beas and consan variance idiosyncraic errors independen of he facor. The marke has 2 some ype of ARCH wih variance m,. r r e i, i m, i, 2 2 ii,, i m, i If he marke has asymmeric volailiy, hen individual socks will oo. 6

CORRELATIONS Beween sock i and sock j assuming idiosyncracies are uncorrelaed. 2 i, j, i j m, i, j, 2 i j m, 2 2 2 2 i im, j jm, Assuming beas are boh posiive, correlaions range from zero o one and increase wih marke volailiy. HOW TO ESTIMATE A ONE FACTOR MODEL FIT THE VOLATILITY OF THE MARKET ESTIMATE THE BETAS OF THE STOCKS AND THE VARIANCE OF THE IDIOSYNCRACIES 7

MARKET VOLATILITY.030.025.020.015.010.005.000 94 95 96 97 98 99 00 01 02 03 04 Condiional Sandard Deviaion CALCULATE DYNAMIC CORRELATIONS h 1 2 2 2 2 2 1h 1 2h 2 8

AXP AND GE AGAIN.9.8.7.6.5.4.3.2.1 94 95 96 97 98 99 00 01 02 03 04 C4_AXP_GE The mulivariae facor model Consider he case of n asses and le is he vecor of beas for he n socks. Then ' hm, Vn Where V n is he diagonal marix wih he idiosyncraic variances (v i s) on he diagonal. 9

The correlaions are given by: 1 ' 1 2 2 R diag diag FACTOR DOUBLE ARCH Idiosyncraic errors also follow a GARCH models 2 ' hm, V, V ~ diagonal garch sd where h 0 V h 0 h n 1 2 2 where i i iri1 hi 1 10

Mulifacor models Iniially, suppress he ime subscrip and consider our K facor model for asse i K ri i i 1F1i2F2... ikfk ei i ijfj ei j1 0 if i j Where Eee i j 2 I is convenien o wrie in marix noaion: r ' i i Bi F ei 1xK Kx1 vi if i=j If he facors are mean zero and he idiosyncraic errors are uncorrelaed wih he facors we ge ha he variance of asse i is: 2 Eri i EBF i eibf i ei EBF e F' B e i i i i E BFF B E BFe E e 2 ' 2 i i i i i BB v i i i EFFis he covariance marix of he facors. 0 11

The ideosyncraic errors are uncorrelaed across asses so ha he covariance of asse i wih asse j is E ri i rj i E BF i ei BjF e j E BF e F B e ' i i j j i ' j i j i j i j E BFF B B EFe E ef B E ee BB i j 0 0 0 Now le r denoe a vecor of n asse reurns. r α B F e nx1 nx1 nxk Kx1 nx1 E E F F E F F EB FFB ee BEFFB EeeBB V rαrα B eb e B e Be Where I is a nxn ideniy marix and V is a nxn marix wih diagonal elemen i equal o v i. B is a marix wih i h row given by B i 12

Now le s exend his o he ime varying case: E EFF E r α r α B B ee B B V Can be made ime varying as in: E 1 E 1FF E 1 r α r α B B ee B B V where E FF 1 If he facors are observable (like porfolios), we can apply a model for ime varying (co)variances o he facors. E FF 1 13

Principal Componens Can we idenify facors from he reurns daa raher han specifying hem? Recall ha he variance covariance marix of he K facor model reurns are given by: E ' r α r α B B V where E FF Wih he facors unknown, B is idenified only up o an orhogonal ransformaion. To see his, le G be ANY KxK marix such ha GG I hen B * =B G and F * =G F also yields he same variance covariance srucure since he marix *' * *' * E F F V GG GG B Β B' ' ' BV B' ΒV 14

This lack of idenificaion can be solved by imposing some condiion (if we don have names on he facors we don care which of he equivalen models we choose)! Typically, we impose ha he facors are uncorrelaed. This yields: B' DBV where D is a diagonal marix We can find a represenaion of he facors by finding he linear combinaion of he daa ha has he highes variance. In pracice, we solve Max xˆ x subjec o xx 1 x 1 1 1 1 1 where ˆ is he esimaed variance covariance marix of he reurns. 15

The soluion is ha x 1 will be he eigenvecor associaed wih he larges eigenvalue of ˆ. Nex we simply normalize x 1 o sum o one so ha i is a valid vecor of porfolio weighs. The second principal componen solves he problem Max xˆx subjec o xx 1 and xˆx 0 x 2 2 2 2 2 1 2 The new se of porfolio weighs correspond o he porfolio ha has he highes variance bu is uncorrelaed wih he firs porfolio. The soluion is ha x 2 will be he eigenvecor associaed wih he second larges eigenvalue of ˆ. Again, he weighs x 2 can be normalized o sum o 1. 16

This process can be repeaed K imes unil we have weighs w 1, w 2,, w K which are he normalized eigenvecors of K larges eigenvecors of he variance covariance marix of reurns. Once he weighs are known, we can consruc he porfolios and rea he facors as if hey are observed. How many facors? The larges eigenvalues should be associaed wih he common componens. The smaller eigenvalues will be associaed wih he ideosynraic variances saring wih he larges and hen becoming smaller. When n ges large, he ideosyncraic variances become an arbirarily small fracion of he overall variance. The facors do no go away. 17

Hence in large samples we can idenify he K facors up o an orhogonal ransformaion. Ofen we look a he scree plo of he eigenvalues from larges o smalles. Number of facors is deermined by he drop off. Looks like 4 facors here: Eigen value 18

Dynamic Condiional Correlaion o DCC is a new ype of mulivariae GARCH model ha is paricularly convenien for big sysems. See Engle(2002) or Engle(2005). o Jus las monh Engle s new ex came ou. DYNAMIC CONDITIONAL CORRELATION OR DCC 1. Esimae volailiies for each asse and compue he sandardized residuals or volailiy adjused reurns. 2. Esimae he ime varying covariances beween hese using a maximum likelihood crierion and one of several models for he correlaions. 3. Form he correlaion marix and covariance marix. They are guaraneed o be posiive definie. 19

HOW IT WORKS When wo asses move in he same direcion, he correlaion is increased slighly. This effec may be sronger in down markes (asymmery in correlaions). When hey move in he opposie direcion i is decreased. The correlaions ofen are assumed o only emporarily deviae from a long run mean CORRELATIONS UPDATE LIKE GARCH Approximaely, z z 1, 1 2, 1 1 1 1 z z 1, 1 2, 1 1 20

An Asymmeric model allows correlaions o increase more when boh prices move down ogeher (like our asymmeric GARCH models. z z z z ( I )( I ) 1, 1 2, 1 1, 1 2, 1 z 0 z 0 1 1 2.9 DCC Correlaions.8.7.6.5.4.3.2.1.0 95 96 97 98 99 00 01 02 03 04 C9_AXP_GE 21

.9.8.7.6.5.4.3.2.1.0 95 96 97 98 99 00 01 02 03 04 C100_AXP_GE C4_AXP_GE C9_AXP_GE A simple approach o esimaing and forecasing large covariance marices Engle(2002), (2006) H DRD, D is diagonal marix of condiional sandard deviaions R is correlaion marix Based on sandardized reurns z 1 D r, V z R z 1 22

The general, mulivariae version of he DCC model looks like: ' 1 1 2 2 R diag Q Q diag Q 1 1 1 1 Q S z z Q where and are scalars. Hence he dynamics are he same for all series. S is he uncondiional correlaion marix. Esimaion T 1 1 L log H r ' H r 2T 1 1 1 1 2T 2 T T T 2 1 2log D r ' D r log R z ' R z z 1 2T 1 T 1 The log likelihood can be wrien as he sum of a par ha depends on he variances and a par ha depends on he correlaions. I is addiive and can be separaed ino wo esimaion problems. ' z 23

Bivariae Esimaion of correlaion parameers 1 z z 2* z z L 2 1 qi, j, 1 ri, jzi, 1 zj, 1 qi, j, 1 qi, j, i, j, q q T 2 2 2 i, j, i, j, i, j, log(2 ) log 1 i, j, 2 1 i, j, ii,, ii,, o Engle proposes esimaing a model for all univariae pairs and hen averaging parameer esimaes (easy, bu no saisically sound). o An alernaive approach is o rea he pairs as a panel daase and simulaneously esimae DCC parameers (inference works here). Asymmeric Dynamic Correlaions of Global Equiy and Bond Reurns Lorenzo Capiello, Rober Engle and Kevin Sheppard Journal of Financial Economerics (2006) 24

Daa Weekly $ reurns Jan 1987 o Feb 2002 (785 observaions) 21 Counry Equiy Series from FTSE All-World Index 13 Daasream Benchmark Bond Indices wih 5 years average mauriy Europe AUSTRIA* BELGIUM* DENMARK* FRANCE* GERMANY* IRELAND* ITALY THE NETHERLANDS* SPAIN SWEDEN* SWITZERLAND* NORWAY UNITED KINGDOM* Ausralasia AUSTRALIA HONG KONG JAPAN* NEW ZEALAND SINGAPORE Americas CANADA* MEXICO UNITED STATES* *wih bond reurns 25

GARCH Models (asymmeric in orange) GARCH AVGARCH NGARCH EGARCH ZGARCH GJR-GARCH APARCH AGARCH NAGARCH 3EQ,8BOND 0 1BOND 6EQ,1BOND 8EQ,1BOND 3EQ,1BOND 0 1EQ,1BOND 0 26

AVERAGE EMU COUNTRY BOND RETURN CORRELATION 27

RESULTS Asymmeric Correlaions correlaions rise afer negaive reurns Shif in level of correlaions wih formaion of Euro Correlaions are rising no jus wihin EMU EMU Bond correlaions are especially high 28