Affine Term Structure Models: An Introduction

Size: px
Start display at page:

Download "Affine Term Structure Models: An Introduction"

Transcription

1 / 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 I. European Universiy Insiue Florence, Sepember 7, 25 The views expressed here are solely he responsibiliy of he auhor and should no be inerpreed as reflecing he views of he Federal Reserve Bank of San Francisco or he Board of Governors of he Federal Reserve Sysem.

2 Overview of Course Day : Term Srucure Modeling in Normal Times Affine erm srucure models and heir esimaion Arbirage-free Nelson-Siegel models Sochasic volailiy Finie-sample bias and oher pracical issues Day 2: The Lower Bound Problem Models ha respec lower bounds for yields How do hey compare o affine models? Which challenges are solved? And which remain? Day 3: Applicaions o Policy Quesions How does QE work? Wha are he poenial coss or risks of QE? How o exrac marke-based inflaion expecaions? 2 / 53

3 Day : Term Srucure Modeling in Normal Times. Lecure: Affine Term Srucure Models: An Inroducion 2. Lecure: The Affine Arbirage-Free Class of Nelson-Siegel Term Srucure Models 3. Lecure: An Arbirage-Free Generalized Nelson-Siegel Term Srucure Model 4. Lecure: Can Spanned Term Srucure Facors Drive Sochasic Volailiy? 5. Lecure: How Efficien is he Kalman Filer a Esimaing Affine Term Srucure Models? 3 / 53

4 Hisory and Background The arbirage-free erm srucure lieraure sared wih founding papers like Vasiček (977) and Cox, Ingersoll, and Ross (985). Affine erm srucure models were hen and remain he workhorse model classes hanks o heir richness and racabiliy. Afer wo decades of sudying one- and wo-facor models, i was clear by he 99s ha more is needed 3 facors a leas. Alhough muli-dimensional affine models come in many varieies, hey can sill be caegorized as shown by Dai and Singleon (2). Wih muli-dimensional models he focus shifs from fi o forecasing and risk premiums. Duffee (22) and Cheridio, Filipović, and Kimmel (27) give us essenially and exended affine risk premiums Maximum affine flexibiliy is achieved! Wih flexibiliy comes esimaion problems. Aenion moves from fi and risk premiums o economeric issues and idenificaion. Mos recenly, lower bounds of yields pose ye anoher challenge. 4 / 53

5 The Ulimae Quesion Can we solve challenges and keep all benefis? Good fi Good forecasing Robus and racable esimaion Respec of lower bounds Goal of course is o provide a posiive answer: Yes! 5 / 53

6 Ouline of Presenaion Inroducion of Affine Models Duffie and Kan (996) Marke Prices of Risk Characerizaion of Canonical Affine Models Dai and Singleon (2) Risk Premiums in Affine Models Compleely affine Dai and Singleon (2) Essenially affine Duffee (22) Exended affine Cheridio, Filipović, and Kimmel (27) Kalman Filer Esimaion of Affine Models Oher Esimaion Schemes for Gaussian Affine Models Joslin, Singleon, and Zhu (2) Hamilon and Wu (22) Conclusion 6 / 53

7 The Class of Affine Term Srucure Models () Define a filered probabiliy space (Ω,F,(F ), Q), where he filraion (F ) = {F : } saisfies he usual condiions (Williams, 997). The sae variables X are assumed o be a Markov process defined on a se M R n ha solves he following sochasic differenial equaion where dx = K Q ()[θ Q () X ]d +Σ()D(X, )dw Q, W Q is a Brownian moion in R n conained in (F ), θ Q : [, T] R n is a bounded, coninuous funcion, K Q : [, T] R n n is a bounded, coninuous funcion, Σ : [, T] R n n is a bounded, coninuous funcion, while D : M [, T] R n n has a diagonal srucure... 7 / 53

8 The Class of Affine Term Srucure Models (2)... he sochasic par of he volailiy marix has a diagonal srucure: γ ()+δ ()X... D(X, ) =.....,... γn ()+δ n ()X where γ () δ ()... δ n() γ() =. and δ() =..... γ n () δ n()... δn n() wih γ : [, T] R n is a bounded, coninuous funcion, δ : [, T] R n n is a bounded, coninuous funcion, δ j () will be used o denoe he jh row of he δ()-marix. 8 / 53

9 The Class of Affine Term Srucure Models (3) Given his srucure we have an affine diffusion process dx = K Q ()[θ Q () X ]d +Σ()D(X, )dw Q. For asse pricing, we will also need he insananeous risk-free rae, which is assumed affine in X : r = ρ ()+ρ () X, where ρ : [, T] R is a bounded, coninuous funcion, ρ : [, T] R n is a bounded, coninuous funcion. Finally, for X o be well defined and admissible, i mus be he case ha for all X M R n. γ j ()+δ j ()X, j n, 9 / 53

10 The Class of Affine Term Srucure Models (4) Noe ha he possible ime dependence of ρ, ρ, K Q, θ Q, Σ, γ, and δ is suppressed in his noaion. / 53 For admissible X, Duffie and Kan (996) show ha arbirage-free bond prices will be exponenial-affine P(, T) = E Q [ ( T exp r u du )] = exp ( B(, T) X + A(, T) ), where B(, T) and A(, T) are soluions o a sysem of ordinary differenial equaions (ODEs) db(, T) = ρ +(K Q ) B(, T) n (Σ B(, T)B(, T) Σ) j,j (δ j ), d 2 da(, T) d = ρ B(, T) K Q θ Q 2 j= n (Σ B(, T)B(, T) Σ) j,j γ j, wih boundary condiions A(T, T) = and B(T, T) =. j=

11 The Class of Affine Term Srucure Models (5) This resul implies ha zero-coupon bond yields are affine: y(, T) = ln P(, T) = T T B(, T) X A(, T). T Key poin: Affine models only characerize he Q-dynamics. X is an affine diffusion process. r is affine in X. Bond yields are affine in X. A(, T) and B(, T) are sraighforward o solve. Noe: The resul is silen abou he P-dynamics of X. For pricing and calibraion exercises, only Q-dynamics are needed. However, for forecasing and erm premium decomposiions, we will need he marke prices of risk. / 53

12 Marke Prices of Risk () 2 / 53 Assume ha X and r are well-defined sochasic processes, bu now formulaed under he objecive P probabiliy measure dx = K P [θ P X ]d +ΣD(X )dw P. A sandard assumpion in finance is ha here exiss a sochasic discoun facor M wih dynamics given by: dm M = r d Λ dw P. M is also known as he pricing kernel or sae price deflaor. Noe ha he process Λ is an n-dimensional vecor.

13 Marke Prices of Risk (2) The pricing of he zero-coupon bond is given by [ P(, T) = E P MT ]. M By Io s lemma, Z = ln M has dynamics dz = r d Λ dw P 2 Λ Λ d. This implies ha ln M T = ln M or, equivalenly, M T M = exp( T T T r s ds r s ds T T T Λ s dw P s 2 Λ sdw P s 2 T T Λ s Λ sds Λ sλ s ds). 3 / 53

14 Marke Prices of Risk (3) Now, consider he process Provided E P [exp( 2 V = exp( T Λ s dw P s 2 Λ s Λ sds). Λ sλ s ds)] < (Novikov s condiion), i is he case ha {V } T is a P-maringale w.r.. F and dq = V T dp dq dp = V T (Radon-Nikodym derivaive) defines a probabiliy measure for which (Girsanov s heorem) dw Q = Λ d + dw P is a Brownian moion under his new Q-probabiliy measure. This has wo implicaions... 4 / 53

15 Marke Prices of Risk (4) Firs, bond prices can be wrien as [ P(, T) = E P MT ] M T = dp(ω) M Ω M T T = exp( r s ds Λ sdws P T Λ Ω 2 sλ s ds)dp(ω) T T = exp( r s ds) exp( Λ sdws P T Λ Ω 2 sλ s ds)dp(ω) T [ T ] = exp( r s ds)dq(ω) = E Q exp( r s ds) Ω as if invesors were risk-neural and jus discoun wih r. Second, under his risk-neural measure, he dynamics are dx dx dx = = K P [θ P X ]d +ΣD(X )dw P = K P [θ P X ]d +ΣD(X )[dw Q Λ d] ) (K P [θ P X ] ΣD(X )Λ d +ΣD(X )dw Q Thus, drif of X is adjused wih ΣD(X )Λ he price of risk! 5 / 53.

16 Risk Premiums in Affine Models To esablish he connecion beween he P- and Q-dynamics, we need he funcional form for he marke prices of risk Λ. Three specificaions dominae he lieraure: Compleely affine risk premiums Dai and Singleon (2). Essenially affine risk premiums Duffee (22). Exended affine risk premiums Cheridio, Filipović, and Kimmel (27). Common rai: P- and Q-dynamics remain affine in X. Noe: Wih risk premiums specified, i is also possible o characerize so-called canonical represenaions of affine models as demonsraed by Dai and Singleon (2). 6 / 53

17 Compleely Affine Risk Premiums () Early on, compleely affine risk premiums were popular where λ is an n vecor. The measure change is dw Q Λ = D(X )λ, = dw P Λ d. Thus, he erm o deduc from he drif of he P-dynamics is ΣD(X )Λ = ΣD(X ) 2 λ γ +δ X... ΣD(X )Λ = Σ γ n +δ n X λ. λ n. Noe: The covariance marix of Λ is affine and X is affine under boh P and Q Model is compleely affine. 7 / 53

18 Compleely Affine Risk Premiums (2) 8 / 53 As emphasized by Duffee (22), compleely affine risk premiums are proporional o yield volailiy. In Gaussian models, his becomes paricularly resricive: Consan volailiy Consan risk premiums K P = K Q. Thus, facor ineracions have o be he same under boh probabiliy measures. Also, he expecaions hypohesis mus hold!

19 Canonical Characerizaion of Affine Models () Dai and Singleon (2) classify affine erm srucure models and describe he canonical represenaion for each model class using compleely affine risk premiums Λ = D(X )λ. They inroduce he noaion of A m (n) model classes: n is he number of facors. Subscrip m is he number of square-roo processes. Problem: Invarian ransformaions imply ha he same model (idenical implicaions for he disribuion of bond yields) can be formulaed in an infinie number of ways. Approach: Find a sandard way of represening A m (n) models. Deermine he maximally flexible admissible specificaion ha is economerically idenifiable. This is he canonical represenaion of A m (n) models. 9 / 53

20 Canonical Characerizaion of Affine Models (2) DS consider admissible affine processes wih Q-dynamics dx = K Q [θ Q X ]d +ΣD(X )dw Q, r = ρ +ρ X. They use compleely affine risk premiums Λ = D(X )λ. This implies ha he P-dynamics are where dx = K Q [θ Q X ]d +ΣD(X )Λ d +ΣD(X )dw P = K P [θ P X ]d +ΣD(X )dw P, K P θ P = K Q Σdiag(λ )δ, = (K P ) [K Q θ Q +Σdiag(λ )γ]. 2 / 53

21 Canonical Characerizaion of Affine Models (3) The canonical represenaion of he A m (n) class is characerized by he following specificaions. ( ) K K P CC m m m (n m) = for m >. K UC (n m) m K UU (n m) (n m) For he Gaussian case (m = ), DS specify K P as riangular, while Singleon (26) and subsequen lieraure impose he riangular propery on K Q. ( ) (θ θ P = P ) C m. (n m) Means of unresriced processes are unidenified and fixed a zero. ( ) ( m m m m (n m) Σ = I, γ =, δ = (n m) B CU (n m) m (n m) (n m) Resricions on volailiy srucure ensure boh admissibiliy and idenificaion. 2 / 53 ).

22 Canonical Characerizaion of Affine Models (3) DS impose addiional parameric resricions. ρ (i), m+ i n. Unresriced processes mus have nonnegaive loadings in he shor rae oherwise heir facor values can swich sign. K P i, θp >, i m. Square-roo processes mus have posiive drif a zero oherwise zero is absorbing sae. K P i,j, i m, i j. Square-roo processes mus be posiively correlaed. θ P i, i m. Square-roo processes mus have nonnegaive drif. δ i,j, m+ i n, j m. Volailiy sensiiviies mus be nonnegaive. 22 / 53

23 Canonical Characerizaion of Affine Models (4) dx dx 2 Example: The canonical A (2) model has P-dynamics ) ( κ P = κ P 2 κ P 22 )[( θ P ) ( X X 2 )] ( X d+ +δ 2 X The compleely affine marke prices of risk are ( )( X ΣD(X )Λ = λ +δ 2X λ 2 Thus, he Q-dynamics are ). dx = K P [θ P X ]d +ΣD(X )dw Q ΣD(X )Λ d ( dx dx 2 ) = + [( κ P )( θ P κ P 2 κ P 22 λ 2/κ P 2 ( X +δ 2X ) ) ( dw,q dw 2,Q ( κ P +λ κ P 2 +δλ 2 2 κ P 22 ). ) ( dw,p dw 2,P )( X X 2 )] d ) Noe he limied flexibiliy beween P- and Q-dynamics! 23 / 53

24 Essenially Affine Risk Premiums () Essenially affine risk premiums were inroduced in Duffee (22). To describe hem, define he diagonal marix: { Dess (γ (X j +δ ) = j X ) /2, if infγ j ()+δ j ()X > ;, oherwise. Noe ha he denominaors in D (X ) already are well-defined processes No exra parameer resricions are required! Now, Duffee (22) considers he following risk premiums where λ 2 is an n n marix. Λ = D(X, )λ + D ess(x )λ 2 X, Noe: For he rows wih in D ess(x ), we will no be able o idenify he corresponding rows in λ 2 They are zero! 24 / 53

25 Essenially Affine Risk Premiums (2) The wo naural exremes are illusraive. A one exreme, we have he Gaussian models. For hose classes of models, Dess(X ) is given by (γ ) /2... Dess(X ) = (γ n ) /2 This means ha D ess(x )λ 2 X adds n n parameers o Λ ha all load on X. By implicaion, here is complee separaion beween he specificaions of K P and K Q in Gaussian models wih essenially affine risk premiums. 25 / 53

26 Essenially Affine Risk Premiums (3) 26 / 53 A he oher exreme, we have he A n (n) models wih square-roo processes only. For hose classes of models, Dess (X ) is given by... Dess (X ) =..... =!... Thus, compleely and essenially affine risk premiums are idenical for A n (n) models!

27 Essenially Affine Risk Premiums (4) dx dx 2 ( Example: The canonical A (2) model has P-dynamics ) ( κ P = κ P 2 κ P 22 )[( θ P ) ( X X 2 )] ( X d+ Since ( Dess (X ) = (+δ 2X ) /2 he essenially affine marke prices of risk are ( X ΣD(X )Λ = +δ 2X Thus, he Q-dynamics are )( λ λ 2 ) + +δ 2 X ), ( )( λ 2 2 λ 2 22 ) ( dw,p dw 2,P )( X X 2 dx = K P [θ P X ]d +ΣD(X )dw Q ΣD(X )Λ d ( dx dx 2 κ P +λ κ P 2 +δλ 2 2 +λ 2 2 κ P 22 +λ 2 22 ) ( κ P = κ P 2 κ P 22 )( X X 2 )( ) d+ θ P λ 2/κ P 2 ( X ) d +δ 2 X ) ( dw,q dw 2,Q ). ). 27 / 53 )

28 Exended Affine Risk Premiums () Cheridio, Filipović, and Kimmel (27) inroduce he exended affine risk premiums. To describe hem, define he following addiional marix: D ex (X ) = { (X j ) /2, if X j > ;, oherwise. Now, CFK consider he following exension of he essenially affine risk premiums Λ = D(X, )λ + Dess (X )λ 2 X + D ex (X )(λ 3 +λ 4 X ), where λ 3 is an n vecor { λ 3 λ 3 = j, if X j > ;, oherwise, while λ 4 is an n n diagonal marix { λ 4 λ 4 = ij, if i j and X i, X j > ;, oherwise. 28 / 53

29 Exended Affine Risk Premiums (2) While for Gaussian A (n) models here is no difference beween essenially and exended affine risk premiums, he exension is paricularly powerful for he A n (n) models: λ 3 = λ 3.. λ 3 n λ λ 4 n and λ 4 λ λ 4 2n = λ 4 n... λ 4 n(n ) Thus, here is a oal of n n exra free parameers in he canonical A n (n) model wih exended affine risk premiums. In principle, his implies ha he P- and Q-dynamics are fully flexible relaive o each oher dx dx = K P [θ P X ]d +ΣD(X )dw P = K Q [θ Q X ]d +ΣD(X )dw Q. However, his is neglecing he sign resricions in λ / 53

30 Exended Affine Risk Premiums (3) Essenially affine risk premiums ruly nes compleely affine risk premiums. CFK claim (p. 29, l. 34) ha he exended affine risk premiums always nes boh essenially and compleely affine risk premiums. A face value, his is rue (more parameers are indeed indispuably allowed!), bu i does no come for free! Under he exended affine risk premiums, he resriced processes have o saisfy Feller condiions under boh probabiliy measures in addiion o all he usual resricions for exisence and admissibiliy: (K P θ P ) j >.5σ 2 jj for j m, (K Q θ Q ) j >.5σ 2 jj for j m. 3 / 53

31 Exended Affine Risk Premiums (4) Problem: We need E P [ exp( T If so, exp( T Λ sdws P 2 Q = exp( Λ s dw P s 2 T T T ] Λ s Λ sds) =. Λ sλ s ds) is a maringale and Λ s dw P s 2 T Λ s Λ sds)p defines an equivalen probabiliy measure and absence of arbirage is ensured. Quesion: If X is a square-roo process, is E P [exp( T Xs dw P s 2 T X s ds)] <? This is abou properies of sochasic exponenials. Answer: Yes provided Feller condiions are saisfied! 3 / 53

32 Exended Affine Risk Premiums (5) 32 / 53 dx dx 2 Example: The canonical A (2) model has P-dynamics ) ( κ P = κ P 2 κ P 22 )[( θ P ) ( X X 2 )] ( X d+ +δ 2 X ) ( dw,p dw 2,P ) Since D ess (X ) = ( (+δ 2 X ) /2 ) ( and D (X ex (X ) = ) /2 ), he exended affine marke prices of risk are ΣD(X )Λ = + ( X +δ 2X ( (X ) /2 )( λ λ 2 )[( λ 3 ) ( + ) + ( )( )]. λ 2 2 λ 2 22 )( X X 2 )

33 Exended Affine Risk Premiums (6) Thus, he Q-dynamics are dx = K P [θ P X ]d +ΣD(X )dw Q ΣD(X )Λ d, which is equivalen o ( ) ( )( ) dx κ P dx 2 = θ P λ 3 κ P 2 κ P /κp 22 λ d 2 /κp 2 ( κ P +λ κ P 2 +δ2 λ 2 +λ2 2 κ P 22 +λ2 22 ( )( X + dw,q +δ 2 X dw 2,Q )( X X 2 ). ) d Noe: The Q-dynamics of X are now compleely deached from is P-dynamics, bu here are Feller condiions κ P θp > 2 σ2 and κ P θp λ 3 > 2 σ2. 33 / 53

34 Esimaion of Affine Models Dai and Singleon (2) use simulaed mehod of momens. Duffee (22) esimaes affine models wih QML and assume n yields observed wihou error. Cheridio, Filipović, and Kimmel (27) use approximae maximum likelihood esimaion and n yields observed wihou error. More recenly, Joslin, Singleon, and Zhu (2) and Hamilon and Wu (22) have offered efficien esimaion algorihms for Gaussian affine models. Throughou he course, he Kalman filer is used for model esimaion, so deails of his esimaion procedure are briefly provided in he following. 34 / 53

35 Kalman Filer Esimaion of Affine Models () 35 / 53 For affine Gaussian models, in general, he condiional mean vecor and he condiional covariance marix are E P [X T F ] = (I exp( K P ))θ P + exp( K P )X, V P [X T F ] = e K P s ΣΣ e (K P ) s ds, where = T is he ime beween observaions. We compue condiional momens of discree observaions and obain he sae ransiion equaion X = (I exp( K P ))θ P + exp( K P )X +ξ.

36 Kalman Filer Esimaion of Affine Models (2) In he sandard Kalman filer, he measuremen equaion is y = A+BX +ε. The assumed error srucure is ( ) [( ξ N ε ) ( Q, H )], where he marix H is assumed diagonal, while he marix Q has he following srucure: Q = e K P s ΣΣ e (K P ) s ds. In addiion, he ransiion and measuremen errors are assumed orhogonal o he iniial sae. Now, Kalman filering is used o evaluae he likelihood funcion. 36 / 53

37 Kalman Filer Esimaion of Affine Models (3) Under a saionariy assumpion he filer can be iniialized a he uncondiional mean and covariance marix: X = θ P, Σ = e K P s ΣΣ e (K P ) s ds. Le Y = (y, y 2,...,y ) be he informaion available a ime, and denoe model parameers by ψ. Consider period and suppose ha he sae updae X and is mean square error marix Σ have been obained. The predicion sep is where Φ X, X = E P [X Y ] = Φ X, (ψ)+φ X, (ψ)x, Σ = Φ X, (ψ)σ Φ X, (ψ) + Q (ψ), (ψ) = (I exp( K P ))θ P, Φ X, (ψ) = exp( K P ), Q (ψ) = e K P s ΣΣ e (K P ) s ds. 37 / 53

38 Kalman Filer Esimaion of Affine Models (4) In he ime- updae sep, X is improved by using he addiional informaion conained in Y : where X = E[X Y ] = X +Σ B(ψ) F v, Σ = Σ Σ B(ψ) F B(ψ)Σ, v = y E[y Y ] = y A(ψ) B(ψ)X, F = cov(v ) = B(ψ)Σ B(ψ) + H(ψ), H(ψ) = diag(σ 2 ε(τ ),...,σ 2 ε(τ N )). A his poin, he Kalman filer has delivered all ingrediens needed o evaluae he Gaussian log likelihood, he predicion-error decomposiion of which is T ( log l(y,...,y T ;ψ) = N 2 log(2π) 2 log F 2 v F v ), = where N is he number of observed yields. 38 / 53

39 Kalman Filer Esimaion of Affine Models (5) I ypically maximize he likelihood wih respec o ψ using he Nelder-Mead simplex algorihm. Upon convergence, we obain sandard errors from he esimaed covariance marix Ω( ψ) = [ T log l ( ψ) log l ( ψ) ], T T ψ ψ = where ψ denoes he esimaed model parameers. Noe: For non-gaussian affine models, he only change is ha Q and Σ are calculaed wih T exp( K P (T s))σd(e P [X s X ])D(E P [X s X ]) Σ exp( (K P ) (T s))ds using formulas from Fisher and Gilles (996). 39 / 53

40 Problems wih Canonical Models 4 / 53 Empirical problems wih canonical models: Since X are laen facors, hey may roae during he esimaion. This leaves muliple maxima wih close o idenical likelihood values bu differen yield decomposiions. Consequence: Two researchers may come up wih very differen esimaion resuls despie he fac ha hey use he same daa, he same model, and he same esimaion mehod. Kim and Orphanides (22) describe hese problems for he specific case of he A (3) model. Duffee (2) also conains an elaborae discussion. Wha o do?

41 Joslin, Singleon, and Zhu (2) () Joslin, Singleon, and Zhu (2) consider a generic discree-ime represenaion of Gaussian affine models X = K P + K P X +Σ X ε P, X = K Q + K Q X +Σ X ε Q, r = ρ +ρ X. Canonical: maximally flexible subjec only o resricions required for economeric idenificaion. Le X be n-dimensional and assume ha we observe J > n bond yields denoed y. For any full-rank marix W R n J, P = Wy represens n porfolios of yields. 4 / 53

42 Joslin, Singleon, and Zhu (2) (2) Assuming ha he n porfolios are priced perfecly, JSZ use invarian affine ransformaions o demonsrae ha any canonical Gaussian affine model has a unique observaionally equivalen represenaion P = K P P + K P P P +Σ P ε P, P = K Q P + K Q P P +Σ P ε Q, r = ρ P +ρ P P, where K Q P, K Q P, ρ P, and ρ P are funcions of (λ Q, k Q,Σ P): λ Q are he eigenvalues of K Q. Σ P Σ P is he covariance marix of innovaions o he porfolios. k Q is long-run mean parameer. 42 / 53

43 Joslin, Singleon, and Zhu (2) (3) Key observaion: Since P are observed facors, he P-condiional likelihood funcion of he observed yields is f(y y ;ψ) = f(y P ;λ Q, k Q,Σ P,Σ ε ) f(p P ; K P P, K P P,Σ P). Criical resul: KP P and K P P can be esimaed by OLS independen of Σ P. This delivers K P P and K P P. In he final sep, λ Q, k Q,Σ P,Σ ε are esimaed by maximizing he likelihood funcion using K P P and K P P as inpu and Σ P Σ P as iniial guess for Σ P Σ P : T f(y y ;ψ) = (2π) (J n) Σ ε = exp ( 2 Σ ε (y o A(λ Q, k Q,Σ P Σ P) B(λ Q, k Q ) P ) 2) f(p P ; K P P, K P P,Σ P). 43 / 53

44 Joslin, Singleon, and Zhu (2) (4) The pricing facors P are linear combinaions of yields. They focus on he firs N principal componens: Their P-dynamics follow an unconsrained VAR. Therefore, no-arbirage assumpions do no improve forecas performance by hemselves. Two conribuions: Compuaionally efficien esimaion of Gaussian affine models. Forecas gains from using a Gaussian affine model relaive o an unconsrained VAR mus come from addiional resricions on he P-dynamics and no from imposing absence of arbirage. Noe: Malab code is available on Ken Singleon s websie. 44 / 53

45 Hamilon and Wu (22) Overview Hamilon and Wu (22, HW) describe anoher novel way o idenify and esimae Gaussian affine models. Unlike oher papers ha emphasize how affine models should be represened, HW s focus is on how o esimae any given Gaussian affine model wihou prejudice. In he paper, HW only consider he popular class of Gaussian affine models where i can be assumed ha M yields are observed wihou error, and M is idenical o he number of pricing facors, i.e., he facors become observable. Implicaion: Their so-called reduced-form represenaion is a resriced vecor auoregression ha can be esimaed wih OLS. Model esimaion involves solving muliple equaions wih as many unknowns, bu presens a more manageable problem han brue force full join ML esimaion of all model parameers. 45 / 53

46 Hamilon and Wu (22) Canonical Model The general Gaussian affine model in discree ime is given by (Eq. (), (2), (3)) X = c +ρx +Σu, X = c Q +ρ Q X +Σu Q, r = δ +δ X. In he case of M = 3, here are 37 free parameers. HW use he following normalizaion: Σ = I M (9 resricions); ρ Q is (upper) riangular (3 resricions); c = (3 resricions); δ. Wih hese 5 resricions he model is canonical wih 22 parameers ha are all economerically idenified. 46 / 53

47 Hamilon and Wu (22) HW Base Model The saring poin for HW s analysis is he model X = ρx + u, wih a shor rae given by r = δ +δ X, while zero-coupon bond yields are given by y = B(ρ Q,δ ) X + A(c Q,δ,ρ Q,δ ), where c Q is an M vecor; ρ Q is an upper riangular M M marix; δ is an M vecor; δ is a consan. 47 / 53

48 Hamilon and Wu (22) Reduced-Form Model () Since M linear combinaions of yields are priced wihou error, i is possible o presen he model in he following reduced-form represenaion ( Y Y 2 ) = ( A A 2 wih u e i.i.d.n(, I N M ). ) ( B + B 2 No Pricing error X becomes observable: ) ( X + Σ e ) u e, Y = A + B X X = B (Y A ). This defines an invarian affine ransformaion A + B X = A + B ρb B X + B u = A B ρb A + B ρb (A A + B X )+B u. This is equivalen o Y = A +φ Y + u. 48 / 53

49 Hamilon and Wu (22) Reduced-Form Model (2) The reduced-from model is where Y = A +φ Y + u, A φ u and u i.i.d.n(, B B ). = A B ρb A, = B ρb, = B u Poin: Reduced-form model is a VAR() esimaed wih OLS. This gives us Â, φ, and he covariance marix of u : Ω = T (Y T Â φ Y )(Y Â φ Y ) = = B B. 49 / 53

50 Hamilon and Wu (22) Reduced-Form Model (3) 5 / 53 The equaions for he N M yields wih errors ake he form Y 2 = A 2 +φ 2Y + u 2. This equaion can also be esimaed wih OLS o yield  2, φ 2, and he covariance marix of u 2 Ω 2 = T T = (Y 2  2 φ 2Y )(Y 2  2 φ 2Y ). This was he easy par. The remaining seps are a lile more compuaionally involved o ge back o (δ,δ,ρ, c Q,ρ Q ).

51 Hamilon and Wu (22) B(ρ Q,δ ) Parameers HW consider he following se of equaions: φ 2 Ω = B 2 B B B = B 2B, Ω = B B. The firs line conains M equaions, while he second line conains M(M + )/2 equaions. Combined hey mach he number of unknowns: (ρ Q,δ ). Guess (ρ Q,δ ), calculae B(ρ Q,δ ), and define π = (vec( φ 2 Ω ), vech( Ω )), g(ρ Q,δ ) = (vec(b 2 B ), vech(b B )). Keep adjusing (ρ Q,δ ) o solve ( ρ Q, δ ) = argmin ρ Q,δ ( π g(ρq,δ )) ( π g(ρ Q,δ )). Noe: If N M >, here are more equaions han unknowns - how o weigh deviaions? For N M =, problem is exacly idenified. 5 / 53

52 Hamilon and Wu (22) ρ and A(c Q,δ,ρ Q,δ ) Parameers Wih ( ρ Q, δ ) in hand, we can now recall ha φ = B ρ B ρ = B φ B. Finally, we exrac he las M + parameers, δ and c Q, from: Â = (I B ρ B )Â, Â 2 = Â2 B 2 B Â. Noe: If N M >, here are again more equaions han unknowns - weighs? For N M =, problem is exacly idenified. Procedure: Guess (c Q,δ ), calculae A(c Q,δ, ρ Q, δ ), and compue sum of squared differences in he equaions. Sop once minimized! This is minimum-chi-square esimaion. 52 / 53

53 Limiaions of JSZ and HW 53 / 53 A shor lis of some limiaions: JSZ and HW only apply o Gaussian affine models. They rely on observable sae variables. Hence, missing observaions or irregular daa frequencies are problemaic o handle. No easy exension o nonlinear measuremen equaions. No exensions o non-gaussian dynamics. Forunaely, here exiss an alernaive approach ha is able o overcome hese limiaions...

Affine term structure models

Affine term structure models 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,

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

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

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

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

Institute for Mathematical Methods in Economics. University of Technology Vienna. Singapore, May Manfred Deistler MULTIVARIATE TIME SERIES ANALYSIS AND FORECASTING Manfred Deisler E O S Economerics and Sysems Theory Insiue for Mahemaical Mehods in Economics Universiy of Technology Vienna Singapore, May 2004 Inroducion

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

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

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

An Introduction to Backward Stochastic Differential Equations (BSDEs) PIMS Summer School 2016 in Mathematical Finance.

An Introduction to Backward Stochastic Differential Equations (BSDEs) PIMS Summer School 2016 in Mathematical Finance. 1 An Inroducion o Backward Sochasic Differenial Equaions (BSDEs) PIMS Summer School 2016 in Mahemaical Finance June 25, 2016 Chrisoph Frei cfrei@ualbera.ca This inroducion is based on Touzi [14], Bouchard

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

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

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

f(s)dw Solution 1. Approximate f by piece-wise constant left-continuous non-random functions f n such that (f(s) f n (s)) 2 ds 0.

f(s)dw Solution 1. Approximate f by piece-wise constant left-continuous non-random functions f n such that (f(s) f n (s)) 2 ds 0. Advanced Financial Models Example shee 3 - Michaelmas 217 Michael Tehranchi Problem 1. Le f : [, R be a coninuous (non-random funcion and W a Brownian moion, and le σ 2 = f(s 2 ds and assume σ 2

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

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

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

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

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

Cash Flow Valuation Mode Lin Discrete Time

Cash Flow Valuation Mode Lin Discrete Time IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728,p-ISSN: 2319-765X, 6, Issue 6 (May. - Jun. 2013), PP 35-41 Cash Flow Valuaion Mode Lin Discree Time Olayiwola. M. A. and Oni, N. O. Deparmen of Mahemaics

More information

6. Stochastic calculus with jump processes

6. Stochastic calculus with jump processes A) Trading sraegies (1/3) Marke wih d asses S = (S 1,, S d ) A rading sraegy can be modelled wih a vecor φ describing he quaniies invesed in each asse a each insan : φ = (φ 1,, φ d ) The value a of a porfolio

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

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

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

What Ties Return Volatilities to Price Valuations and Fundamentals? On-Line Appendix Wha Ties Reurn Volailiies o Price Valuaions and Fundamenals? On-Line Appendix Alexander David Haskayne School of Business, Universiy of Calgary Piero Veronesi Universiy of Chicago Booh School of Business,

More information

An Introduction to Malliavin calculus and its applications

An Introduction to Malliavin calculus and its applications An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214

More information

LECTURE 1: GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS

LECTURE 1: GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS LECTURE : GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS We will work wih a coninuous ime reversible Markov chain X on a finie conneced sae space, wih generaor Lf(x = y q x,yf(y. (Recall ha q

More information

An Arbitrage-Free Generalized Nelson-Siegel Term Structure Model

An Arbitrage-Free Generalized Nelson-Siegel Term Structure Model FEDERAL RESERVE BANK OF SAN FRANCISCO WORKING PAPER SERIES An Arbirage-Free Generalized Nelson-Siegel Term Srucure Model Jens H E Chrisensen Federal Reserve Bank of San Francisco Francis X Diebold Universiy

More information

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

Georey E. Hinton. University oftoronto.   Technical Report CRG-TR February 22, Abstract Parameer Esimaion for Linear Dynamical Sysems Zoubin Ghahramani Georey E. Hinon Deparmen of Compuer Science Universiy oftorono 6 King's College Road Torono, Canada M5S A4 Email: zoubin@cs.orono.edu Technical

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

Lecture Notes 2. The Hilbert Space Approach to Time Series

Lecture Notes 2. The Hilbert Space Approach to Time Series Time Series Seven N. Durlauf Universiy of Wisconsin. Basic ideas Lecure Noes. The Hilber Space Approach o Time Series The Hilber space framework provides a very powerful language for discussing he relaionship

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

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

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

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

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

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

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal

More information

On a Fractional Stochastic Landau-Ginzburg Equation

On a Fractional Stochastic Landau-Ginzburg Equation Applied Mahemaical Sciences, Vol. 4, 1, no. 7, 317-35 On a Fracional Sochasic Landau-Ginzburg Equaion Nguyen Tien Dung Deparmen of Mahemaics, FPT Universiy 15B Pham Hung Sree, Hanoi, Vienam dungn@fp.edu.vn

More information

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix

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

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

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

DEPARTMENT OF ECONOMICS AND FINANCE COLLEGE OF BUSINESS AND ECONOMICS UNIVERSITY OF CANTERBURY CHRISTCHURCH, NEW ZEALAND DEPARTMENT OF ECONOMICS AND FINANCE COLLEGE OF BUSINESS AND ECONOMICS UNIVERSITY OF CANTERBURY CHRISTCHURCH, NEW ZEALAND Asymmery and Leverage in Condiional Volailiy Models Michael McAleer WORKING PAPER

More information

Math 334 Fall 2011 Homework 11 Solutions

Math 334 Fall 2011 Homework 11 Solutions Dec. 2, 2 Mah 334 Fall 2 Homework Soluions Basic Problem. Transform he following iniial value problem ino an iniial value problem for a sysem: u + p()u + q() u g(), u() u, u () v. () Soluion. Le v u. Then

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

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

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

Stochastic Modelling in Finance - Solutions to sheet 8

Stochastic Modelling in Finance - Solutions to sheet 8 Sochasic Modelling in Finance - Soluions o shee 8 8.1 The price of a defaulable asse can be modeled as ds S = µ d + σ dw dn where µ, σ are consans, (W ) is a sandard Brownian moion and (N ) is a one jump

More information

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation: M ah 5 7 Fall 9 L ecure O c. 4, 9 ) Hamilon- J acobi Equaion: Weak S oluion We coninue he sudy of he Hamilon-Jacobi equaion: We have shown ha u + H D u) = R n, ) ; u = g R n { = }. ). In general we canno

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

Excel-Based Solution Method For The Optimal Policy Of The Hadley And Whittin s Exact Model With Arma Demand

Excel-Based Solution Method For The Optimal Policy Of The Hadley And Whittin s Exact Model With Arma Demand Excel-Based Soluion Mehod For The Opimal Policy Of The Hadley And Whiin s Exac Model Wih Arma Demand Kal Nami School of Business and Economics Winson Salem Sae Universiy Winson Salem, NC 27110 Phone: (336)750-2338

More information

Optimal Investment under Dynamic Risk Constraints and Partial Information

Optimal Investment under Dynamic Risk Constraints and Partial Information Opimal Invesmen under Dynamic Risk Consrains and Parial Informaion Wolfgang Puschögl Johann Radon Insiue for Compuaional and Applied Mahemaics (RICAM) Ausrian Academy of Sciences www.ricam.oeaw.ac.a 2

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

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

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

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

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution Physics 7b: Saisical Mechanics Fokker-Planck Equaion The Langevin equaion approach o he evoluion of he velociy disribuion for he Brownian paricle migh leave you uncomforable. A more formal reamen of his

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

Backward stochastic dynamics on a filtered probability space

Backward stochastic dynamics on a filtered probability space Backward sochasic dynamics on a filered probabiliy space Gechun Liang Oxford-Man Insiue, Universiy of Oxford based on join work wih Terry Lyons and Zhongmin Qian Page 1 of 15 gliang@oxford-man.ox.ac.uk

More information

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

Chapter 2. Models, Censoring, and Likelihood for Failure-Time Data Chaper 2 Models, Censoring, and Likelihood for Failure-Time Daa William Q. Meeker and Luis A. Escobar Iowa Sae Universiy and Louisiana Sae Universiy Copyrigh 1998-2008 W. Q. Meeker and L. A. Escobar. Based

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

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

) were both constant and we brought them from under the integral.

) were both constant and we brought them from under the integral. YIELD-PER-RECRUIT (coninued The yield-per-recrui model applies o a cohor, bu we saw in he Age Disribuions lecure ha he properies of a cohor do no apply in general o a collecion of cohors, which is wha

More information

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]

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

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

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION SUPPLEMENTARY INFORMATION DOI: 0.038/NCLIMATE893 Temporal resoluion and DICE * Supplemenal Informaion Alex L. Maren and Sephen C. Newbold Naional Cener for Environmenal Economics, US Environmenal Proecion

More information

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

Recursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems 8 Froniers in Signal Processing, Vol. 1, No. 1, July 217 hps://dx.doi.org/1.2266/fsp.217.112 Recursive Leas-Squares Fixed-Inerval Smooher Using Covariance Informaion based on Innovaion Approach in Linear

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

Speaker Adaptation Techniques For Continuous Speech Using Medium and Small Adaptation Data Sets. Constantinos Boulis

Speaker Adaptation Techniques For Continuous Speech Using Medium and Small Adaptation Data Sets. Constantinos Boulis Speaker Adapaion Techniques For Coninuous Speech Using Medium and Small Adapaion Daa Ses Consaninos Boulis Ouline of he Presenaion Inroducion o he speaker adapaion problem Maximum Likelihood Sochasic Transformaions

More information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

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

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

Modeling Economic Time Series with Stochastic Linear Difference Equations

Modeling Economic Time Series with Stochastic Linear Difference Equations A. Thiemer, SLDG.mcd, 6..6 FH-Kiel Universiy of Applied Sciences Prof. Dr. Andreas Thiemer e-mail: andreas.hiemer@fh-kiel.de Modeling Economic Time Series wih Sochasic Linear Difference Equaions Summary:

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

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

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

Ensamble methods: Bagging and Boosting

Ensamble methods: Bagging and Boosting Lecure 21 Ensamble mehods: Bagging and Boosing Milos Hauskrech milos@cs.pi.edu 5329 Senno Square Ensemble mehods Mixure of expers Muliple base models (classifiers, regressors), each covers a differen par

More information

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec

More information

Financial Econometrics Introduction to Realized Variance

Financial Econometrics Introduction to Realized Variance Financial Economerics Inroducion o Realized Variance Eric Zivo May 16, 2011 Ouline Inroducion Realized Variance Defined Quadraic Variaion and Realized Variance Asympoic Disribuion Theory for Realized Variance

More information

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

EMS SCM joint meeting. On stochastic partial differential equations of parabolic type EMS SCM join meeing Barcelona, May 28-30, 2015 On sochasic parial differenial equaions of parabolic ype Isván Gyöngy School of Mahemaics and Maxwell Insiue Edinburgh Universiy 1 I. Filering problem II.

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

The Optimal Stopping Time for Selling an Asset When It Is Uncertain Whether the Price Process Is Increasing or Decreasing When the Horizon Is Infinite

The Optimal Stopping Time for Selling an Asset When It Is Uncertain Whether the Price Process Is Increasing or Decreasing When the Horizon Is Infinite American Journal of Operaions Research, 08, 8, 8-9 hp://wwwscirporg/journal/ajor ISSN Online: 60-8849 ISSN Prin: 60-8830 The Opimal Sopping Time for Selling an Asse When I Is Uncerain Wheher he Price Process

More information

Robert Kollmann. 6 September 2017

Robert Kollmann. 6 September 2017 Appendix: Supplemenary maerial for Tracable Likelihood-Based Esimaion of Non- Linear DSGE Models Economics Leers (available online 6 Sepember 207) hp://dx.doi.org/0.06/j.econle.207.08.027 Rober Kollmann

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

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

2.160 System Identification, Estimation, and Learning. Lecture Notes No. 8. March 6, 2006 2.160 Sysem Idenificaion, Esimaion, and Learning Lecure Noes No. 8 March 6, 2006 4.9 Eended Kalman Filer In many pracical problems, he process dynamics are nonlinear. w Process Dynamics v y u Model (Linearized)

More information

THE MYSTERY OF STOCHASTIC MECHANICS. Edward Nelson Department of Mathematics Princeton University

THE MYSTERY OF STOCHASTIC MECHANICS. Edward Nelson Department of Mathematics Princeton University THE MYSTERY OF STOCHASTIC MECHANICS Edward Nelson Deparmen of Mahemaics Princeon Universiy 1 Classical Hamilon-Jacobi heory N paricles of various masses on a Euclidean space. Incorporae he masses in he

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

Estimation for Parameters in Partially Observed Linear Stochastic System

Estimation for Parameters in Partially Observed Linear Stochastic System Esimaion for Parameers in Parially Observed Linear Sochasic Sysem Chao Wei Absrac This paper is concerned wih he esimaion problem for parially observed linear sochasic sysem. The sae esimaion equaion is

More information

DEPARTMENT OF STATISTICS

DEPARTMENT OF STATISTICS A Tes for Mulivariae ARCH Effecs R. Sco Hacker and Abdulnasser Haemi-J 004: DEPARTMENT OF STATISTICS S-0 07 LUND SWEDEN A Tes for Mulivariae ARCH Effecs R. Sco Hacker Jönköping Inernaional Business School

More information

Estimation of Poses with Particle Filters

Estimation of Poses with Particle Filters Esimaion of Poses wih Paricle Filers Dr.-Ing. Bernd Ludwig Chair for Arificial Inelligence Deparmen of Compuer Science Friedrich-Alexander-Universiä Erlangen-Nürnberg 12/05/2008 Dr.-Ing. Bernd Ludwig (FAU

More information

Lecture 20: Riccati Equations and Least Squares Feedback Control

Lecture 20: Riccati Equations and Least Squares Feedback Control 34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he

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

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

EXERCISES FOR SECTION 1.5

EXERCISES FOR SECTION 1.5 1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler

More information

Sequential Importance Resampling (SIR) Particle Filter

Sequential Importance Resampling (SIR) Particle Filter Paricle Filers++ Pieer Abbeel UC Berkeley EECS Many slides adaped from Thrun, Burgard and Fox, Probabilisic Roboics 1. Algorihm paricle_filer( S -1, u, z ): 2. Sequenial Imporance Resampling (SIR) Paricle

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

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details! MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his

More information

Chapter 3 Boundary Value Problem

Chapter 3 Boundary Value Problem Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le

More information

Empirical Process Theory

Empirical Process Theory Empirical Process heory 4.384 ime Series Analysis, Fall 27 Reciaion by Paul Schrimpf Supplemenary o lecures given by Anna Mikusheva Ocober 7, 28 Reciaion 7 Empirical Process heory Le x be a real-valued

More information

Probabilistic Robotics

Probabilistic Robotics Probabilisic Roboics Bayes Filer Implemenaions Gaussian filers Bayes Filer Reminder Predicion bel p u bel d Correcion bel η p z bel Gaussians : ~ π e p N p - Univariae / / : ~ μ μ μ e p Ν p d π Mulivariae

More information

CHAPTER 2 Signals And Spectra

CHAPTER 2 Signals And Spectra CHAPER Signals And Specra Properies of Signals and Noise In communicaion sysems he received waveform is usually caegorized ino he desired par conaining he informaion, and he undesired par. he desired par

More information