Affine term structure models

Size: px
Start display at page:

Download "Affine term structure models"

Transcription

1 Affine erm srucure models A. Inro o Gaussian affine erm srucure models B. Esimaion by minimum chi square (Hamilon and Wu) C. Esimaion by OLS (Adrian, Moench, and Crump) D. Dynamic Nelson-Siegel model (Chrisensen, Diebold, and Rudebusch) E. Small-sample bias (Bauer, Rudebusch, and Wu) 1

2 A. Inro o Gaussian affine erm srucure models P n price a of pure-discoun n-period bond P n E M 1 P n1,1 One approach: specify M 1 and derive bond prices. Today: reverse engineer sar wih convenien empirical model of risk and hen figure ou wha M 1 his requires. Laer: will give a parial equilibrium example which would imply his M 1. 2

3 Suppose here is an r 1 vecor of possibly unobserved facors ha summarize everyhing ha maers for deermining ineres raes. Suppose log of any bond price is affine funcion of hese facors: p n n n e.g., r 3 (level, slope, curvaure). 3

4 Conjecure ha facors follow a firsorder homoskedasic Gaussian VAR: 1 c u 1 u 1 i.i.d. N0, I r summarizes unpredicabiliy of facors and risk premia should be funcions of. 4

5 Thus for he informaion se a, 1 N, E 1 c 5

6 Consider asse ha pays i,1 dollars nex period. How much would you pay for each of hese asses i 1,..,r oday? If risk neural, price would be e r i. 6

7 If risk averse, maybe I would only pay e r i Q Q r1 r1 rrr1 i price of facor i risk If i 0, ac as if column i of 0 (uncerainy abou facor i does no affec price of any securiy). 7

8 True disribuion of facors (someimes called "hisorical disribuion" or "P measure") 1 P N, 8

9 Risk-averse invesors behave he same way as a risk-neural invesor would if ha person believed he disribuion was insead he "Q measure" or "risk-neural disribuion" Q 1 N Q, Q 9

10 Quesion: wha pricing kernel M 1 would imply his? e r Q E M 1 1 M ;, d 1 We would obain desired answer if M 1 1 ;, e r 1 ; Q, 10

11 1 ; Q, 1 2 r/2 exp 1 Q 1 1 Q 2 1 Q 1 1 Q

12 Or since 1 c u u 1 1 ; Q, 1 ;, exp1/2 u 1 12

13 Conclusion: 1 ; Q, 1 ;, exp

14 Summary: M 1 1 ;, e r 1 ; Q, calls for specifying M 1 expr 1/2 u 1 14

15 Suppose we furher conjecure ha price of risk is also affine funcion: r1 r1 rr r1 15

16 Then Q c c Q Q c Q c Q 16

17 P-measure dynamics: 1 c u 1 u 1 P N0, I r Q-measure dynamics: c Q Q Q u 1 1 Q u 1 Q N0, Ir 17

18 Invesors ac he way a risk-neural invesor would who hough he facors follow he Q-measure disribuion, ha is, P n E Q e r P n1,1 E M 1 P n1,1 18

19 Recall ha if z N, 2 hen Ee z exp 2 /2 19

20 Thus if p n n n, we require e p n E Q e r e p n1,1 exp n n exp r n1 n1 Q 1/2 n1 n1 20

21 exp n n exp r n1 n1 Q 1/2 n1 n1 Or since r p Q c Q Q we require n n 1 1 n1 n1 c Q Q 1/2 n1 n1 21

22 n n 1 1 n1 n1 c Q Q 1/2 n1 n1 n n1 Q 1 n n1 n1 c Q 1/2 n1 n1 1 22

23 Given c Q, Q, 1, 1, we can calculae he log of he price of any bond p n. If c Q c and Q his would correspond o he expecaions hypohesis of he erm srucure. 23

24 Gives us a way of summarizing dynamics of yield curve in erms of separae conribuions of risk premia and expecaions E r j. 24

25 B. Esimaion by minimum chi square (Hamilon and Wu) Model implies y n a n b n for yield on any mauriy n and an r 1 vecor. If number of observed yields r, sysem is sochasically singular. 25

26 One soluion: assume any observed yield n differs from model predicion by measuremen or specificaion error: y n a n b n n 26

27 Collec sysem for observed yields in a vecor y y n1,, y n2,,..., y nn, N1 y a B c 1 v This is sae-space sysem observaion equaion: y sae equaion: parameers a and B are highly nonlinear funcions of c,,,,, 1, 1. 27

28 Alernaive popular approach: assume ha model holds exacly for r of he observed y n y 6 a 6 b 6 y 24, a 24 b 24 y 120, a 120 b 120 a n n /n b n n /n n n1 Q 1 n n1 n1 c Q 1/2 n1 n1 1 y 1 a 1 B 1 28

29 For oher yields, model holds wih error y 2 y 3, y 12,, y 36,, y 60,, y 84, y 2 a 2 B

30 y 1 a 1 B 1 c 1 v y a y 1,1 1 a 1 a 1 B 1 c B 1 B 1 1 a 1 11 B 1 B1 1 1 B 1 v y 1 follows VAR(1) 30

31 y 2 a 2 B 2 2 y 1 a 1 B 1 y a y 1 2 a 2 a 2 B 2 B 1 1 a 1 21 B 2 B1 1 31

32 y 1 a 1 11 y 1,1 1 y 2 a 2 21 y 1 2 y y 1,y 2 follows resriced VAR(1) whose coefficiens are nonlinear funcions of c,,,,, 1, 1. 32

33 Can firs esimae unresriced parameersa 1, 11,a 2, 21, by OLS, hen find srucural parameers ha make prediced values close o observed (minimize chi square saisic for es ha resricions are valid). Asympoically equivalen o MLE, bu simpler. 33

34 Can easily generalize above o suppose ha here are r linear combinaions of y for which model holds wihou error: y 1 r1 y 1 r1 H y rn N1 H a rnn1 H B rnnr y 1 a 1 11 y 1,1 1 e.g., y 1 firs r principal componens of y ( H firs r eigenvecors of T 1 T 1 y yy y ) 34

35 Noe he model as wrien is unidenified. If q and q, he model would be observaionally idenical Same if Q, QQ 1, b 1 b 1 Q 1 35

36 (1) Sample normalizaion: c 0 Q lower riangular Q ii Q jj for i j I r elemens of b 1 nonposiive 36

37 This normalizaion is inernally inconsisen. If observe Hy direcly and y a B, hen HaHB requiring Ha 0 and HB I r. Upside: in pracice Ha urns ou o be close o 0 and HB close o I r even wihou imposing. 37

38 (2) Joslin, Singleon Zhu normalizaion (Rev Financial Sudies, 2011) as implemened by Hamilon-Wu (J. Economerics, 2014): unknown parameers are c,,, 1, 2, 1 1, 2,..., r are eigenvalues of Q. 38

39 n x n 1 n j0 11 x j n1 1 n2 1 nn 1 K rn n1 2 n2 2 nn 2 n1 r n2 r nn r 1 0 V rr 0 r Q K H 1 V K H rr 1 K H 1 1 r for 1 r 1,1,...,1 39

40 Then HB I r. A relaed calculaion for he inerceps guaranees Ha 0. 40

41 Benefis: c,, esimaed by simple OLS: c 1 1 E 1 1, 2, 1 esimaed by MCS on y 2 a E

42 C. Esimaion by OLS Simpler approach (Adrian, Crump, and Moench, JFE 2013): Don impose any resricions, ge everyhing by OLS. firs r 5 principal componens of observed se of yields y n for n 6m,12m,18m,24m,30m, 36m, 42m, 48m, 54m, 60m,7y,10y 1 c v 1 can esimae by OLS 42

43 log price of n-period bond: p n n n excess reurn of n-period bond: x n1,1 p n1,1 p n r p n1,1 p n p 1 n1 n1 c u 1 n n

44 x n1,1 n1 n1 c v 1 n n 1 1 a n1 b n1 v 1 c n1 a n1 n1 n1 c n 1 b n1 c n1 n1 n1 n 1 44

45 b n1 c n1 n1 n1 n 1 Prediced coefficiens: affine model implies n n1 1 c n1 b n1 45

46 c n1 b n1 Proposal: esimae x n1,1 a n1 b n1 v 1 c n1 e n1,1 for v 1 1 ĉ by unresriced OLS separaely for each n 46

47 Then esimae o minimize rr sum of squared discrepancies beween ĉ n1 and b n1 across n n ĉ n1 b n1 n b n1b n1 1 47

48 Prediced inerceps: a n1 n1 n1 c n 1 Affine model implies n n1 n1 c 1 1/2 n1 n1 2 a n1 n1 1/2 n1 n1 2 48

49 a n1 n1 1/2 n1 n1 2 Le ã n1 â n1 1/2b n1 b n1 2 Esimae by minimizing difference beween ã n1 and b n1 n b n1b n1 1 n b n1ã n1 49

50 Summary: we have now esimaed c,,,,, 2 using only OLS regressions. Can esimae 1 and 1 by OLS regression r 1 1 e 1 And hen calculae n and n by recursion: n n1 1 n n1 n1 c 1 1/2 n1 n1 2 50

51 From hese we can calculae he prediced yield on any bond y n n 1 n n. 51

52 We can hen redo he recursions seing 0 and 0 o ge prediced yields if invesors were risk neural RF n RF n1 1 RF n RF n1 RF n1 c 1 n1 2 1/2 n1 RF n 1 RF n RF n y n 52

53 Can calculae he risk premium RF as he difference y n y n 53

54 Updaed daily a hps:// research/daa_indicaors/acmtermpremium.xls year yield and erm premium Apr Mar Feb Jan Dec Nov Oc Sep Aug Jul Jun May May Apr Mar Feb Jan Dec Nov Oc Sep Aug Jul Jun May Apr Mar Feb Jan-2016 erm premium yield

55 D. Dynamic Nelson-Siegel model (Chrisensen, Diebold, and Rudebusch) Recall ha we ofen summarize forward curve a dae using funcion such as f n 0 1 expn/ 1 2 n/ 1 expn/ 1 3 n/ 2 expn/ 2 Nelson-Siegel: Describe forward rae of any mauriy n as smooh funcion of hree magniudes a (level, slope, and curvaure). 55

56 Consider following special case of GATSM normalizaion and addiional resricions (Bauer, 2011, following Chrisensen, Diebold, and Rudebusch): Q has eigenvalues 1,, c Q 0 b 1 1, 1, 0 56

57 Implicaions: 1 c Q Q Q v 1 r f n b 1 E Q n b 1 Q n n n n n n 2 n n

58 f n 1 n 2 n n1 1 3 So forward rae f n loads on facor 1 wih weigh 1 facor 2 wih weigh n 2 facor 3 wih weigh n n1 1 58

59 Facor loadings for yields of differen mauriy (gamma = 0.98) Mauriy (in monhs) Facor 1 (level) Facor 2 (slope) Facor 3 (curvaure) 59

60 Dynamic Nelson-Siegel: can hen esimae P-measure dynamics for sae vecor as c 1 v which gives complee dynamic descripion of process for all yields. 60

61 E. Small-sample bias (Bauer, Rudebusch, Wu) c 1 v Cauion: for unresriced OLS esimaion of, eigenvalues biased downard. Implicaion: if 0, n 0 for large n 61

62 Empirical models wan o aribue mos of flucuaion in y n for large n o (changes in risk premium) no n (expecaions componen). Bauer, Rudebusch and Wu (JBES, 2012). Correc esimae of for small-sample bias. 62

63 4-year forward raes (black) and expeced 4-year-ahead shor raes wih and wihou bias correcion 63

64 4-year forward raes (black) and componen aribued o risk premium wih and wihou bias correcion 64

65 Big picure: (1) Mehods exis o decompose long yield ino expecaions componen and risk premium. (2) Idenificaion comes from fac ha predicable excess reurns aribued o risk premium. (3) Specific answer sensiive o assumed underlying forecasing model. 65

Affine term structure models

Affine term structure models /5/07 Affie erm srucure models A. Iro o Gaussia affie erm srucure models B. Esimaio by miimum chi square (Hamilo ad Wu) C. Esimaio by OLS (Adria, Moech, ad Crump) D. Dyamic Nelso-Siegel model (Chrisese,

More information

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

Dynamic models for largedimensional. Yields on U.S. Treasury securities (3 months to 10 years) y t Dynamic models for largedimensional vecor sysems A. Principal componens analysis Suppose we have a large number of variables observed a dae Goal: can we summarize mos of he feaures of he daa using jus

More information

Forward guidance. Fed funds target during /15/2017

Forward guidance. Fed funds target during /15/2017 Forward guidance Fed funds arge during 2004 A. A wo-dimensional characerizaion of moneary shocks (Gürkynak, Sack, and Swanson, 2005) B. Odyssean versus Delphic foreign guidance (Campbell e al., 2012) C.

More information

13.3 Term structure models

13.3 Term structure models 13.3 Term srucure models 13.3.1 Expecaions hypohesis model - Simples "model" a) shor rae b) expecaions o ge oher prices Resul: y () = 1 h +1 δ = φ( δ)+ε +1 f () = E (y +1) (1) =δ + φ( δ) f (3) = E (y +)

More information

C. Theoretical channels 1. Conditions for complete neutrality Suppose preferences are E t. Monetary policy at the zero lower bound: Theory 11/22/2017

C. Theoretical channels 1. Conditions for complete neutrality Suppose preferences are E t. Monetary policy at the zero lower bound: Theory 11/22/2017 //7 Moneary policy a he zero lower bound: Theory A. Theoreical channels. Condiions for complee neuraliy (Eggersson and Woodford, ). Marke fricions. Preferred habia and risk-bearing (Hamilon and Wu, ) B.

More information

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

Introduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate. Inroducion Gordon Model (1962): D P = r g r = consan discoun rae, g = consan dividend growh rae. If raional expecaions of fuure discoun raes and dividend growh vary over ime, so should he D/P raio. Since

More information

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

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance

More information

Dynamic Econometric Models: Y t = + 0 X t + 1 X t X t k X t-k + e t. A. Autoregressive Model:

Dynamic Econometric Models: Y t = + 0 X t + 1 X t X t k X t-k + e t. A. Autoregressive Model: Dynamic Economeric Models: A. Auoregressive Model: Y = + 0 X 1 Y -1 + 2 Y -2 + k Y -k + e (Wih lagged dependen variable(s) on he RHS) B. Disribued-lag Model: Y = + 0 X + 1 X -1 + 2 X -2 + + k X -k + e

More information

GMM - Generalized Method of Moments

GMM - Generalized Method of Moments GMM - Generalized Mehod of Momens Conens GMM esimaion, shor inroducion 2 GMM inuiion: Maching momens 2 3 General overview of GMM esimaion. 3 3. Weighing marix...........................................

More information

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

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8) I. Definiions and Problems A. Perfec Mulicollineariy Econ7 Applied Economerics Topic 7: Mulicollineariy (Sudenmund, Chaper 8) Definiion: Perfec mulicollineariy exiss in a following K-variable regression

More information

Linear Gaussian State Space Models

Linear Gaussian State Space Models Linear Gaussian Sae Space Models Srucural Time Series Models Level and Trend Models Basic Srucural Model (BSM Dynamic Linear Models Sae Space Model Represenaion Level, Trend, and Seasonal Models Time Varying

More information

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

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

Smoothing. Backward smoother: At any give T, replace the observation yt by a combination of observations at & before T Smoohing Consan process Separae signal & noise Smooh he daa: Backward smooher: A an give, replace he observaion b a combinaion of observaions a & before Simple smooher : replace he curren observaion wih

More information

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Kriging Models Predicing Arazine Concenraions in Surface Waer Draining Agriculural Waersheds Paul L. Mosquin, Jeremy Aldworh, Wenlin Chen Supplemenal Maerial Number

More information

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

L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS NA568 Mobile Roboics: Mehods & Algorihms Today s Topic Quick review on (Linear) Kalman Filer Kalman Filering for Non-Linear Sysems Exended Kalman Filer (EKF)

More information

OBJECTIVES OF TIME SERIES ANALYSIS

OBJECTIVES OF TIME SERIES ANALYSIS OBJECTIVES OF TIME SERIES ANALYSIS Undersanding he dynamic or imedependen srucure of he observaions of a single series (univariae analysis) Forecasing of fuure observaions Asceraining he leading, lagging

More information

Unit Root Time Series. Univariate random walk

Unit Root Time Series. Univariate random walk Uni Roo ime Series Univariae random walk Consider he regression y y where ~ iid N 0, he leas squares esimae of is: ˆ yy y y yy Now wha if = If y y hen le y 0 =0 so ha y j j If ~ iid N 0, hen y ~ N 0, he

More information

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

Comparing Means: t-tests for One Sample & Two Related Samples Comparing Means: -Tess for One Sample & Two Relaed Samples Using he z-tes: Assumpions -Tess for One Sample & Two Relaed Samples The z-es (of a sample mean agains a populaion mean) is based on he assumpion

More information

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

Testing the Random Walk Model. i.i.d. ( ) r he random walk heory saes: esing he Random Walk Model µ ε () np = + np + Momen Condiions where where ε ~ i.i.d he idea here is o es direcly he resricions imposed by momen condiions. lnp lnp µ ( lnp lnp

More information

Generalized Least Squares

Generalized Least Squares Generalized Leas Squares Augus 006 1 Modified Model Original assumpions: 1 Specificaion: y = Xβ + ε (1) Eε =0 3 EX 0 ε =0 4 Eεε 0 = σ I In his secion, we consider relaxing assumpion (4) Insead, assume

More information

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

Vectorautoregressive Model and Cointegration Analysis. Time Series Analysis Dr. Sevtap Kestel 1 Vecorauoregressive Model and Coinegraion Analysis Par V Time Series Analysis Dr. Sevap Kesel 1 Vecorauoregression Vecor auoregression (VAR) is an economeric model used o capure he evoluion and he inerdependencies

More information

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

Outline. lse-logo. Outline. Outline. 1 Wald Test. 2 The Likelihood Ratio Test. 3 Lagrange Multiplier Tests Ouline Ouline Hypohesis Tes wihin he Maximum Likelihood Framework There are hree main frequenis approaches o inference wihin he Maximum Likelihood framework: he Wald es, he Likelihood Raio es and he Lagrange

More information

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

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t Exercise 7 C P = α + β R P + u C = αp + βr + v (a) (b) C R = α P R + β + w (c) Assumpions abou he disurbances u, v, w : Classical assumions on he disurbance of one of he equaions, eg. on (b): E(v v s P,

More information

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

Two Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017 Two Popular Bayesian Esimaors: Paricle and Kalman Filers McGill COMP 765 Sep 14 h, 2017 1 1 1, dx x Bel x u x P x z P Recall: Bayes Filers,,,,,,, 1 1 1 1 u z u x P u z u x z P Bayes z = observaion u =

More information

Distribution of Estimates

Distribution of Estimates Disribuion of Esimaes From Economerics (40) Linear Regression Model Assume (y,x ) is iid and E(x e )0 Esimaion Consisency y α + βx + he esimaes approach he rue values as he sample size increases Esimaion

More information

Lecture 3: Exponential Smoothing

Lecture 3: Exponential Smoothing NATCOR: Forecasing & Predicive Analyics Lecure 3: Exponenial Smoohing John Boylan Lancaser Cenre for Forecasing Deparmen of Managemen Science Mehods and Models Forecasing Mehod A (numerical) procedure

More information

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates Biol. 356 Lab 8. Moraliy, Recruimen, and Migraion Raes (modified from Cox, 00, General Ecology Lab Manual, McGraw Hill) Las week we esimaed populaion size hrough several mehods. One assumpion of all hese

More information

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

Bias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé Bias in Condiional and Uncondiional Fixed Effecs Logi Esimaion: a Correcion * Tom Coupé Economics Educaion and Research Consorium, Naional Universiy of Kyiv Mohyla Academy Address: Vul Voloska 10, 04070

More information

ACE 562 Fall Lecture 4: Simple Linear Regression Model: Specification and Estimation. by Professor Scott H. Irwin

ACE 562 Fall Lecture 4: Simple Linear Regression Model: Specification and Estimation. by Professor Scott H. Irwin ACE 56 Fall 005 Lecure 4: Simple Linear Regression Model: Specificaion and Esimaion by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Simple Regression: Economic and Saisical Model

More information

Financial Econometrics Jeffrey R. Russell Midterm Winter 2009 SOLUTIONS

Financial Econometrics Jeffrey R. Russell Midterm Winter 2009 SOLUTIONS Name SOLUTIONS Financial Economerics Jeffrey R. Russell Miderm Winer 009 SOLUTIONS You have 80 minues o complee he exam. Use can use a calculaor and noes. Try o fi all your work in he space provided. If

More information

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

State-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter Sae-Space Models Iniializaion, Esimaion and Smoohing of he Kalman Filer Iniializaion of he Kalman Filer The Kalman filer shows how o updae pas predicors and he corresponding predicion error variances when

More information

Forecasting optimally

Forecasting optimally I) ile: Forecas Evaluaion II) Conens: Evaluaing forecass, properies of opimal forecass, esing properies of opimal forecass, saisical comparison of forecas accuracy III) Documenaion: - Diebold, Francis

More information

3.1 More on model selection

3.1 More on model selection 3. More on Model selecion 3. Comparing models AIC, BIC, Adjused R squared. 3. Over Fiing problem. 3.3 Sample spliing. 3. More on model selecion crieria Ofen afer model fiing you are lef wih a handful of

More information

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...

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... Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger

More information

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

More information

How to Deal with Structural Breaks in Practical Cointegration Analysis

How to Deal with Structural Breaks in Practical Cointegration Analysis How o Deal wih Srucural Breaks in Pracical Coinegraion Analysis Roselyne Joyeux * School of Economic and Financial Sudies Macquarie Universiy December 00 ABSTRACT In his noe we consider he reamen of srucural

More information

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

Testing for a Single Factor Model in the Multivariate State Space Framework esing for a Single Facor Model in he Mulivariae Sae Space Framework Chen C.-Y. M. Chiba and M. Kobayashi Inernaional Graduae School of Social Sciences Yokohama Naional Universiy Japan Faculy of Economics

More information

Ready for euro? Empirical study of the actual monetary policy independence in Poland VECM modelling

Ready for euro? Empirical study of the actual monetary policy independence in Poland VECM modelling Macroeconomerics Handou 2 Ready for euro? Empirical sudy of he acual moneary policy independence in Poland VECM modelling 1. Inroducion This classes are based on: Łukasz Goczek & Dagmara Mycielska, 2013.

More information

SZG Macro 2011 Lecture 3: Dynamic Programming. SZG macro 2011 lecture 3 1

SZG Macro 2011 Lecture 3: Dynamic Programming. SZG macro 2011 lecture 3 1 SZG Macro 2011 Lecure 3: Dynamic Programming SZG macro 2011 lecure 3 1 Background Our previous discussion of opimal consumpion over ime and of opimal capial accumulaion sugges sudying he general decision

More information

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

Volatility. Many economic series, and most financial series, display conditional volatility Volailiy Many economic series, and mos financial series, display condiional volailiy The condiional variance changes over ime There are periods of high volailiy When large changes frequenly occur And periods

More information

Linear Dynamic Models

Linear Dynamic Models Linear Dnamic Models and Forecasing Reference aricle: Ineracions beween he muliplier analsis and he principle of acceleraion Ouline. The sae space ssem as an approach o working wih ssems of difference

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and

More information

Financial Econometrics Kalman Filter: some applications to Finance University of Evry - Master 2

Financial Econometrics Kalman Filter: some applications to Finance University of Evry - Master 2 Financial Economerics Kalman Filer: some applicaions o Finance Universiy of Evry - Maser 2 Eric Bouyé January 27, 2009 Conens 1 Sae-space models 2 2 The Scalar Kalman Filer 2 21 Presenaion 2 22 Summary

More information

Y, where. 1 Estimate St.error

Y, where. 1 Estimate St.error 1 HG Feb 2014 ECON 5101 Exercises III - 24 Feb 2014 Exercise 1 In lecure noes 3 (LN3 page 11) we esimaed an ARMA(1,2) for daa) for he period, 1978q2-2013q2 Le Y ln BNP ln BNP (Norwegian Model: Y Y, where

More information

Notes on Kalman Filtering

Notes on Kalman Filtering Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

ACE 564 Spring Lecture 7. Extensions of The Multiple Regression Model: Dummy Independent Variables. by Professor Scott H.

ACE 564 Spring Lecture 7. Extensions of The Multiple Regression Model: Dummy Independent Variables. by Professor Scott H. ACE 564 Spring 2006 Lecure 7 Exensions of The Muliple Regression Model: Dumm Independen Variables b Professor Sco H. Irwin Readings: Griffihs, Hill and Judge. "Dumm Variables and Varing Coefficien Models

More information

BU Macro BU Macro Fall 2008, Lecture 4

BU Macro BU Macro Fall 2008, Lecture 4 Dynamic Programming BU Macro 2008 Lecure 4 1 Ouline 1. Cerainy opimizaion problem used o illusrae: a. Resricions on exogenous variables b. Value funcion c. Policy funcion d. The Bellman equaion and an

More information

Regression with Time Series Data

Regression with Time Series Data Regression wih Time Series Daa y = β 0 + β 1 x 1 +...+ β k x k + u Serial Correlaion and Heeroskedasiciy Time Series - Serial Correlaion and Heeroskedasiciy 1 Serially Correlaed Errors: Consequences Wih

More information

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

Let us start with a two dimensional case. We consider a vector ( x, Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our

More information

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

Robust estimation based on the first- and third-moment restrictions of the power transformation model h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,

More information

14 Autoregressive Moving Average Models

14 Autoregressive Moving Average Models 14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class

More information

Lecture 10 Estimating Nonlinear Regression Models

Lecture 10 Estimating Nonlinear Regression Models Lecure 0 Esimaing Nonlinear Regression Models References: Greene, Economeric Analysis, Chaper 0 Consider he following regression model: y = f(x, β) + ε =,, x is kx for each, β is an rxconsan vecor, ε is

More information

m = 41 members n = 27 (nonfounders), f = 14 (founders) 8 markers from chromosome 19

m = 41 members n = 27 (nonfounders), f = 14 (founders) 8 markers from chromosome 19 Sequenial Imporance Sampling (SIS) AKA Paricle Filering, Sequenial Impuaion (Kong, Liu, Wong, 994) For many problems, sampling direcly from he arge disribuion is difficul or impossible. One reason possible

More information

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

Linear Combinations of Volatility Forecasts for the WIG20 and Polish Exchange Rates Eliza Buszkowska Universiy of Poznań, Poland Linear Combinaions of Volailiy Forecass for he WIG0 and Polish Exchange Raes Absrak. As is known forecas combinaions may be beer forecass hen forecass obained

More information

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

Chapter 5. Heterocedastic Models. Introduction to time series (2008) 1 Chaper 5 Heerocedasic Models Inroducion o ime series (2008) 1 Chaper 5. Conens. 5.1. The ARCH model. 5.2. The GARCH model. 5.3. The exponenial GARCH model. 5.4. The CHARMA model. 5.5. Random coefficien

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

Affine Term Structure Models: An Introduction

Affine Term Structure Models: An Introduction / 53 Affine Term Srucure Models: An Inroducion Jens H. E. Chrisensen Federal Reserve Bank of San Francisco Term Srucure Modeling and he Lower Bound Problem Day : Term Srucure Modeling in Normal Times Lecure

More information

Vector autoregression VAR. Case 1

Vector autoregression VAR. Case 1 Vecor auoregression VAR So far we have focused mosl on models where deends onl on as. More generall we migh wan o consider oin models ha involve more han one variable. There are wo reasons: Firs, we migh

More information

Wednesday, November 7 Handout: Heteroskedasticity

Wednesday, November 7 Handout: Heteroskedasticity Amhers College Deparmen of Economics Economics 360 Fall 202 Wednesday, November 7 Handou: Heeroskedasiciy Preview Review o Regression Model o Sandard Ordinary Leas Squares (OLS) Premises o Esimaion Procedures

More information

STATE-SPACE MODELLING. A mass balance across the tank gives:

STATE-SPACE MODELLING. A mass balance across the tank gives: B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing

More information

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms

More information

THE UNIVERSITY OF TEXAS AT AUSTIN McCombs School of Business

THE UNIVERSITY OF TEXAS AT AUSTIN McCombs School of Business THE UNIVERITY OF TEXA AT AUTIN McCombs chool of Business TA 7.5 Tom hively CLAICAL EAONAL DECOMPOITION - MULTIPLICATIVE MODEL Examples of easonaliy 8000 Quarerly sales for Wal-Mar for quarers a l e s 6000

More information

USP. Surplus-Production Models

USP. Surplus-Production Models USP Surplus-Producion Models 2 Overview Purpose of slides: Inroducion o he producion model Overview of differen mehods of fiing Go over some criique of he mehod Source: Haddon 2001, Chaper 10 Hilborn and

More information

Pricing of Risk in Natural Gas Futures Contracts: A New Approach to Affine Term Structure Models Subject to Changes in Regime

Pricing of Risk in Natural Gas Futures Contracts: A New Approach to Affine Term Structure Models Subject to Changes in Regime Pricing of Risk in Naural Gas Fuures Conracs: A New Approach o Affine Term Srucure Models Subjec o Changes in Regime Irina Zhecheva Deparmen of Economics Universiy of California, San Diego Sepember 7,

More information

Økonomisk Kandidateksamen 2005(II) Econometrics 2. Solution

Økonomisk Kandidateksamen 2005(II) Econometrics 2. Solution Økonomisk Kandidaeksamen 2005(II) Economerics 2 Soluion his is he proposed soluion for he exam in Economerics 2. For compleeness he soluion gives formal answers o mos of he quesions alhough his is no always

More information

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

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H. ACE 56 Fall 005 Lecure 5: he Simple Linear Regression Model: Sampling Properies of he Leas Squares Esimaors by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Inference in he Simple

More information

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

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals

More information

Exponential Smoothing

Exponential Smoothing Exponenial moohing Inroducion A simple mehod for forecasing. Does no require long series. Enables o decompose he series ino a rend and seasonal effecs. Paricularly useful mehod when here is a need o forecas

More information

Exercise: Building an Error Correction Model of Private Consumption. Part II Testing for Cointegration 1

Exercise: Building an Error Correction Model of Private Consumption. Part II Testing for Cointegration 1 Bo Sjo 200--24 Exercise: Building an Error Correcion Model of Privae Consumpion. Par II Tesing for Coinegraion Learning objecives: This lab inroduces esing for he order of inegraion and coinegraion. The

More information

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

Advanced time-series analysis (University of Lund, Economic History Department) Advanced ime-series analysis (Universiy of Lund, Economic Hisory Deparmen) 30 Jan-3 February and 6-30 March 01 Lecure 9 Vecor Auoregression (VAR) echniques: moivaion and applicaions. Esimaion procedure.

More information

Article from. Predictive Analytics and Futurism. July 2016 Issue 13

Article from. Predictive Analytics and Futurism. July 2016 Issue 13 Aricle from Predicive Analyics and Fuurism July 6 Issue An Inroducion o Incremenal Learning By Qiang Wu and Dave Snell Machine learning provides useful ools for predicive analyics The ypical machine learning

More information

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

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid

More information

Distribution of Least Squares

Distribution of Least Squares Disribuion of Leas Squares In classic regression, if he errors are iid normal, and independen of he regressors, hen he leas squares esimaes have an exac normal disribuion, no jus asympoic his is no rue

More information

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II Roland Siegwar Margaria Chli Paul Furgale Marco Huer Marin Rufli Davide Scaramuzza ETH Maser Course: 151-0854-00L Auonomous Mobile Robos Localizaion II ACT and SEE For all do, (predicion updae / ACT),

More information

Lecture 33: November 29

Lecture 33: November 29 36-705: Inermediae Saisics Fall 2017 Lecurer: Siva Balakrishnan Lecure 33: November 29 Today we will coninue discussing he boosrap, and hen ry o undersand why i works in a simple case. In he las lecure

More information

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

Explaining Total Factor Productivity. Ulrich Kohli University of Geneva December 2015 Explaining Toal Facor Produciviy Ulrich Kohli Universiy of Geneva December 2015 Needed: A Theory of Toal Facor Produciviy Edward C. Presco (1998) 2 1. Inroducion Toal Facor Produciviy (TFP) has become

More information

Mathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation

Mathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation Mahcad Lecure #8 In-class Workshee Curve Fiing and Inerpolaion A he end of his lecure, you will be able o: explain he difference beween curve fiing and inerpolaion decide wheher curve fiing or inerpolaion

More information

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

References are appeared in the last slide. Last update: (1393/08/19) SYSEM IDEIFICAIO Ali Karimpour Associae Professor Ferdowsi Universi of Mashhad References are appeared in he las slide. Las updae: 0..204 393/08/9 Lecure 5 lecure 5 Parameer Esimaion Mehods opics o be

More information

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004 ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform

More information

Policy regimes Theory

Policy regimes Theory Advanced Moneary Theory and Policy EPOS 2012/13 Policy regimes Theory Giovanni Di Barolomeo giovanni.dibarolomeo@uniroma1.i The moneary policy regime The simple model: x = - s (i - p e ) + x e + e D p

More information

20. Applications of the Genetic-Drift Model

20. Applications of the Genetic-Drift Model 0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0

More information

Chapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws

Chapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws Chaper 5: Phenomena Phenomena: The reacion (aq) + B(aq) C(aq) was sudied a wo differen emperaures (98 K and 35 K). For each emperaure he reacion was sared by puing differen concenraions of he 3 species

More information

System of Linear Differential Equations

System of Linear Differential Equations Sysem of Linear Differenial Equaions In "Ordinary Differenial Equaions" we've learned how o solve a differenial equaion for a variable, such as: y'k5$e K2$x =0 solve DE yx = K 5 2 ek2 x C_C1 2$y''C7$y

More information

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence MATH 433/533, Fourier Analysis Secion 6, Proof of Fourier s Theorem for Poinwise Convergence Firs, some commens abou inegraing periodic funcions. If g is a periodic funcion, g(x + ) g(x) for all real x,

More information

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.

More information

Applications in Industry (Extended) Kalman Filter. Week Date Lecture Title

Applications in Industry (Extended) Kalman Filter. Week Date Lecture Title hp://elec34.com Applicaions in Indusry (Eended) Kalman Filer 26 School of Informaion echnology and Elecrical Engineering a he Universiy of Queensland Lecure Schedule: Week Dae Lecure ile 29-Feb Inroducion

More information

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17 EES 16A Designing Informaion Devices and Sysems I Spring 019 Lecure Noes Noe 17 17.1 apaciive ouchscreen In he las noe, we saw ha a capacior consiss of wo pieces on conducive maerial separaed by a nonconducive

More information

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check

More information

Review - Quiz # 1. 1 g(y) dy = f(x) dx. y x. = u, so that y = xu and dy. dx (Sometimes you may want to use the substitution x y

Review - Quiz # 1. 1 g(y) dy = f(x) dx. y x. = u, so that y = xu and dy. dx (Sometimes you may want to use the substitution x y Review - Quiz # 1 (1) Solving Special Tpes of Firs Order Equaions I. Separable Equaions (SE). d = f() g() Mehod of Soluion : 1 g() d = f() (The soluions ma be given implicil b he above formula. Remember,

More information

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

1. Consider a pure-exchange economy with stochastic endowments. The state of the economy Answer 4 of he following 5 quesions. 1. Consider a pure-exchange economy wih sochasic endowmens. The sae of he economy in period, 0,1,..., is he hisory of evens s ( s0, s1,..., s ). The iniial sae is given.

More information

Bond Risk Premia. November 6, Abstract

Bond Risk Premia. November 6, Abstract Bond Risk Premia John H. Cochrane and Monika Piazzesi November 6, 2001 Absrac This paper sudies risk premia in he erm srucure. We sar wih regressions of annual holding period reurns on forward raes. We

More information

Solutions: Wednesday, November 14

Solutions: Wednesday, November 14 Amhers College Deparmen of Economics Economics 360 Fall 2012 Soluions: Wednesday, November 14 Judicial Daa: Cross secion daa of judicial and economic saisics for he fify saes in 2000. JudExp CrimesAll

More information

Econ Autocorrelation. Sanjaya DeSilva

Econ Autocorrelation. Sanjaya DeSilva Econ 39 - Auocorrelaion Sanjaya DeSilva Ocober 3, 008 1 Definiion Auocorrelaion (or serial correlaion) occurs when he error erm of one observaion is correlaed wih he error erm of any oher observaion. This

More information

2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS

2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS Andrei Tokmakoff, MIT Deparmen of Chemisry, 2/22/2007 2-17 2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS The mahemaical formulaion of he dynamics of a quanum sysem is no unique. So far we have described

More information

5.2. The Natural Logarithm. Solution

5.2. The Natural Logarithm. Solution 5.2 The Naural Logarihm The number e is an irraional number, similar in naure o π. Is non-erminaing, non-repeaing value is e 2.718 281 828 59. Like π, e also occurs frequenly in naural phenomena. In fac,

More information

Has the Business Cycle Changed? Evidence and Explanations. Appendix

Has the Business Cycle Changed? Evidence and Explanations. Appendix Has he Business Ccle Changed? Evidence and Explanaions Appendix Augus 2003 James H. Sock Deparmen of Economics, Harvard Universi and he Naional Bureau of Economic Research and Mark W. Wason* Woodrow Wilson

More information

The Simple Linear Regression Model: Reporting the Results and Choosing the Functional Form

The Simple Linear Regression Model: Reporting the Results and Choosing the Functional Form Chaper 6 The Simple Linear Regression Model: Reporing he Resuls and Choosing he Funcional Form To complee he analysis of he simple linear regression model, in his chaper we will consider how o measure

More information

Online Appendix to Solution Methods for Models with Rare Disasters

Online Appendix to Solution Methods for Models with Rare Disasters Online Appendix o Soluion Mehods for Models wih Rare Disasers Jesús Fernández-Villaverde and Oren Levinal In his Online Appendix, we presen he Euler condiions of he model, we develop he pricing Calvo block,

More information