Optimal policy in the New Keynesian Baseline model with Dynare

In this video I focus on Optimal Policy in the baseline New Keynesian model and discuss concepts like the Divine Coincidence, Indeterminacy, Taylor Principle, Optimal Policy under Commitment, Optimal Policy under Discretion, Simple Implementable Rules and how to compare different policy regimes.

Even though this video provides some insight on the theory behind these concepts, I try to not get caught up in the details, but focus mainly on the implementation aspects in Dynare.

Video

Timestamps

Optimal Policy in Baseline New Keynesian Model

Implementing Optimal Policy in Baseline New Keynesian Model with simple rules

Policy Trade-Offs and Breaking Divine Coincidence

Optimal Policy Under Commitment

Optimal Policy Under Discretion

Comparison of Optimal Policy

Simple Implementable Rules

Optimal Simple Implementable Rules

Outro & References

Slides

Presentation

Codes

NK_linear_common.inc

v ; v ; p ; m # # # # [ p [ y [ y [ w [ y [ a [ z [ r [ r [ y [ y @ [ i @ i @ i @ e s e a a a / O P L K n i n _ n = n n = n = n n _ n _ n _ n _ # n i # i # i # n h n r a z c y y y r r i p n w r e e r A B R R S V E T @ P P @ d C M S A A a e a g a a a = a c a R a = a n a r a n a g i a e e e d o d C C T i t ( I b M G F s _ _ _ _ i i e p p a L E H H I A P H # H H # e o E I M P m = m a m + m m ; m H m m a m e m a m a f m = l = l = n ; c ; o o h t h a t u E N o e n g n r e x s s m P T O O G R S E i I I e l m G _ B P e B e p e ( e S e e O e R e t e a e t e p e s s d k v v p p i e t t R U r e a a a e o _ _ e H A A Z M P I T f _ _ n ( p A Y D A = E = = = 1 = I = = A = H = = = l = = = = M = r e r e P i s a a y y s u i C t p t a a z t A A H L A P Y d l o = A A = ' T ' - ' - ' G ' ' * ' O ' - ' = ' P ' y O ' _ i _ H f ; r r r r n F y s W H G a < l e I O M I i i s ( = = L N A D 1 P A l M r T a P Z D S D i D S D - N I n f n I i i i d r o I A e h r N O E f n i 1 ( ( A e * y / r L a A e F ( r * e I e i e I e y P n a a _ e e g g s e e u d T N n c t s N e t - 1 1 M w p n S o P b * s P - e z f G f - f _ f _ O t t M t P p p h h r e r i H T e o t P a e A + - B i a I d H o c o 1 f ( i M i p i Y i n L e ; O I s s t t f s O A r p p O r L V T D K e m G u A r + u p ) e - n A n i n A n a r N + E _ _ r t S o t U B a y : L ) p P A H A e ( i M c ) V r r + r 1 i * i e i * i t = e P * z a ( ( e h o p r T I l $ $ $ $ $ $ $ $ $ $ ; a H R E * y + c A t * d A c o e e ) t P t ( t a t ; = s O P p C C e e f t i L / $ $ $ $ $ $ $ $ $ $ $ $ $ $ { { { { { { { { ! { { r A P T ( n 1 * i n e R e c p n i S i + i ; i t L H i ) ) t i b A I P f w { { { { { { { { { { { { { { \ \ \ \ \ \ \ \ = \ \ a ) H A S e ) I ( o ; m P e s c + o I o 1 o o 1 I e = = s t w o u N T u w a z c y y \ r r i \ n w \ \ a b r r s e t p p m / I ) I s + S i n a H c s _ e n _ n ) n n R = _ + 2 2 o e a n t Y Y b t w } } } } ^ t } } p } } v l e h h i a p h 1 h h e ( ) G i K i n I o s a e Y ; a = P P 1 1 0 0 f r r ) e l h . $ $ $ $ { i { $ $ i $ $ a a p t o o g r s e i i t 1 ( M a A C - f d * n ' ; s p n A r n o t I H ; ; 1 2 t m e d W o i e g n l n } r r h a _ _ m p i t _ _ e - ( 1 A n P u p u ' n s ] h s a * e a u e 2 E I 6 1 w s a A r c n a d a $ e e a } a { a h l a { { r A S - + P r i n ] ; t i _ t ( a t t * _ - a F n i R G u t e t p p } $ } z } i o } \ y s L I B ( P A v e c r f z u 1 l u p R p Y 2 W r o o y n R F L N . } } s s $ $ } $ } n $ p } P G E V h * e ( t a t ; r - r u u i * 0 i e f u A I i U o } y } i i } $ } i } H M T A i y ' + i i e a R i a t l e y 2 l : n l t N T c r $ } $ l l $ $ } $ A A A R l _ ] 1 o n r l H n l e + ; 1 l t d a h T N e G g $ o o } + * * P l g ) n t ' O t g : P i y h a t e Y E n e / n n $ A ( T H i a - ' ' ] r A e o a H J o e t e ; S s n l _ _ L 1 H I p p r ] ] a ) r u p E I o M u i r h S e e i a z P - E + s ; _ t * e t ' x _ h u G o o w r c } } H A T A n e a s p ] o Y a t c N n v p i F f a e $ $ A L A L C a + t u g * n s a U , e e t O o l n * P ) P u t o ( t e y n c n r h R r s E H / H r ) f 1 r ' n _ e h G e s t o P e P A T A v + - a ] o g s l r e i i h u A m u s S ) H ) e y i R t u a e e n t o a t o b / I + E / ' _ n H e s p P r d e h n t P r l > L V T ( ] g t O ' ; f i r e . e A e i . ( ( ( ( ( ( ( ( ( ( O A A 1 a e Z ] O e s a r i v R c ( ( ( ( ( ( ( ( ( ( ( ( ( ( l l l l l l l l l l N R * - p r ) n i t l t e T d l l l l l l l l l l l l l l o o o o o o o o o o ) P O A ( e * e f r v n I e L o o o o o o o o o o o o o o n n n n n n n n n n ; H M L + s z - e i P e w C t i n n n n n n n n n n n n n n g g g g g g g g g g I E P 1 t ; T r b u r i t U a c g g g g g g g g g g g g g g _ _ _ _ _ _ _ _ _ _ + G H ) ' o u b s l h L i e _ _ _ _ _ _ _ _ _ _ _ _ _ _ n n n n n n n n n n A A A ; ] - t l i l e A l n n n n n n n n n n n n n n n a a a a a a a a a a L ; ) O e i o R s s a a a a a a a a a a a a a a m m m m m m m m m m P ) n c n b i . e m m m m m m m m m m m m m m e e e e e e e e e e H ; e i e m P , e e e e e e e e e e e e e e = = = = = = = = = = A ' t L 3 p U = = = = = = = = = = = = = = ' ' ' ' ' ' ' ' ' ' ) ] i u l R ' ' ' ' ' ' ' ' ' ' ' ' ' ' o d a a i i D C i o ; a c o s i P t p c o n o n r n i h r t p n i u u n n i a n u n e f e e O e r o u a u a e o n o e e r e s t t v v x l f t d n f d S c e n t t t t a m f u a c e c o o e e i v l p / s t u E h f s p u p u l i l r l h f m o c c r r t o a u o e h l w . n e u u r u r n a s n e i u o o s s - t t r e , a o r m t a t a i a t w o r n n r r e e S p i a r l e p l l n l i w a l e u t r r t r o f m s L r S o n t ( g t o o g o n s e e E F i o n e o i a e g c i l o a i e i n r e g c f l l I r g b e d p c n e y e o o u p n r n k y e l a a a S i l a f d i u e t n g t t e t ( e ( a c t t ' s i b e b f b n y t s s p ( e s e l d l s s b t i i ) c t i e a y l s h h h ( d u l r t r o o h h o o o o h z l d c i e o e o o l e t o e e g ( g o o r r n n i b k s , f c c o v g s r s l c c ' e d t a h k k g ( t a t d o d k k s ) t p l e y c T e o s l d t e g e ' ' h e r a m ' k a d r p p d s o e r e r v v ) ) a c e s a ) y r r e ) g v a a d r h f t n T l b o o v ' t ( t s e s e n e i d a o y c c ) d s e l e s v s o r c y r e e s e s o ) ) i l e i e l s s s v ) ( g ( ' s ' n o n t l o R s s ) ' l l ) s ) g c y a r u ' s ) o d o ) p y e ' s l ( ( ) s g e g ' r ) t R e l l ) v ) o p p i u ' o o ' d d d r r c l ) g g ) e s e u o o i e v s v c c c t ' d d ) t e e y ) e e s ' s i s s ' v v s ) s o s s ) ) ) n ' ' s s ' ' ' ) ) s s ) ) ) ) ) ' ' ) )

NK_linear_optimal_rule1.mod

% % % @ @ S V T R R B A E s c % # # I A H H H E L P t h s E d i G R E O O T P S e e t x e n M P T Z A A H I a c o o f c A H A A L d k c g i l I O y ; h e n u N ; _ n e d s o e = = = = = = = = i u M m s O " 1 5 3 0 0 0 1 9 u N N ; ; / . . . / ; l o P K 4 5 9 9 4 ( n O _ ; ; ; 9 ; o e L l ; r - i d f = n e o e r r 1 a - r = o _ n c 1 e o , m i r m r u o f l n = e . 3 i 0 % % % n ) c ; "

NK_linear_optimal_rule2.mod

% % % @ @ S V T R R B A E P P s c s e P P e v d # # I A H H H E L P H H t h t s H H n a y o d i G R E O O T P S I I e e o t I I d r n p e n M P T Z A A H I _ _ a c c i _ _ ; o a t f c A H A A L P Y d k h m P Y b r i i l I O I y ; _ a I , s e m n u N E ; s t E _ a e d i e , 0 a s l e = = = = = = = = = = m d . e M u _ 1 1 n r O " 1 5 3 0 0 0 1 9 1 0 l p . 2 z s u N N ; ; / . . . / ; . . ( a 5 5 i l P K 4 5 9 9 4 5 1 o r , , t e O _ ; ; ; 9 ; ; 2 r a c i L l ; 5 d m 0 0 v w i ; e s , , i i = n r ; y t t e = 2 2 y h 2 a 1 ; ; ; r , y f _ i _ e c r n e o f a d m = t b m 3 a o 0 c n ) y k . ; _ i g t n a o c p " t a r r _ g n e a t t v a r r _ i r a e b a l l e s i % % % i p i e n w ; % % % d a o y l p n t t a e i r r o e n n _ a s s t _ e i . n v v s e a i r t a o i n b v d s i t q = y u i M r c _ e k . q e u w n i a d r y o e _ s t n o a v m a a e r s s o s ; b u s m e b l a o l c l k v a r i a b l e s a r e o b s e r v b a b l e :

NK_linear_costpush_common.inc

v ; v ; p ; m # [ p [ x [ y [ w [ y [ a [ z [ r [ r [ y [ y [ y [ u [ u e s % K V % p e a a a / P n i n n = n n = n = n n _ n _ n _ n _ n _ n n = n t A A i n r a z c y y y y x r r i p n w u r e e e r A B R R R S V E T K V d C S a e a = a a a = a c a R a = a e a r a e a g a g a = a R d e W P R A e d C C T i t ( I b M G F s _ _ _ _ _ i i e p p p a L E H H H I A P H A A e o I m = m m + m m ; m H m m = m e m = m a m a m m H ; a e P T l = ; o o h t h a t u E N o e n g e e r e x s s s m P T O O O G R S E P R l m _ e B e x e ( e S e e O e R e S e a e P e p e p e K e O d A H l 0 p p i e t t R U r e a a e o _ _ _ e H A A Z U M P I T P T ( p Y = E = ( = 1 = I = = A = H = I = l = S = = = = = A = U y n = E ; y y s u i C t p a a z u t A A H L A A H l o A ' T ' + ' - ' G ' ' * ' O ' G ' = ' I ' y ' x ' P ' * _ e ( T o y r r n F y s W H G a < l e I O E i s = N A D 1 P A l M r T a P Z D M D i D _ D - D + I P c u s e 1 A t _ i i i d r o I A e h r N T n i ( e * y ) r L a A e F ( r * e A e i e Y e y e ( m A o ( t d - = h g g g s e e u d T N n c t s A e t 1 w p n o P b * s P - e z f * f - f A f _ f y p * s - a T K e a h h r e r i H T e o t a e + i a - d H o c o 1 f ( i ( i p i * i n i _ l ( t 1 t t H A r p t t f s O A r p p r V K e m u A r + u p ) e - n y n i n a n a n e i y ) e h E P = r t S o t U B a y : ) p A e ( i 1 c ) V r r + r 1 i _ i e i ; i t i - c _ p + _ e T P v 0 ( ( e h o p r T I l $ $ $ $ $ $ $ $ $ $ $ ; a R y + c / t * d A c o e e ) t e t ( t t ; t y i e u e m s A A a ; C C e e f t i L / $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ { { { { { { { { { { { r P n 1 S i n e R e c p n i ( i + i i i _ t - s p o e ) r y ) ) t i b A I P f w { { { { { { { { { { { { { { { { { { \ \ \ \ \ \ \ \ \ \ \ a H e ) I I o ; m P e s c + o + o 1 o o o n y h s d E i _ s t w o u N T u w a z c y y \ y x r r i \ n w u \ \ \ a b r r r s e t k m I s + S G n a H c s _ e n 1 n ) n n n a d _ _ e p ( P a n 2 2 o e a n t Y Y b t w } } } } ^ t } ^ } } p } } } v l e h h h i a p h a a e ) i K M n I o s a e ) ; t e n p u l a 1 S b a 0 0 f r r ) e l h . $ $ $ $ { i { $ { $ $ i $ $ $ a a a p t o o o g r s e p r t / a A C A f d * n ' ; s p e - r e o l ) f a r ; ; r - I l t 1 2 t m e d W o i e g n l e e } r r r h a _ _ _ m p i t p t e ( n P u * u ' n s ] h s f y e f u i ; i t o a B L e = 6 1 w s a A r c n a d } } $ e e e a } a { { a h l a a h r S P r ( n ] ; t i _ f _ a f t n n ) c m E O s 0 - a F n i R G u t e } } p p p } $ } z u } i o } } e s I P A v i c r f z i e l i p k i ; e e T N ; 2 W r o o y n R F L N . } $ $ s s s $ $ } } $ } n $ $ t G h * e i t a t ; c ) c u i t s t A ; a y 0 i e f u A I i U o } y i i i } } $ } a M i x ' - i i e i + i i t n i s e * r = 2 l : n l t N T c r $ } l l l $ $ $ } A l ] p o n r e ( n e g o ' r T e 0 1 l t d a h T N e G g $ o o o $ * l + i n t ' n 1 t n g n ] s H ; i y h a t e Y E n e / n n n ( i e ' ' ] t - e t a v E z y J o e t e ; S s n l _ _ _ 1 p u ( ] ] R r p a o t T e _ o M u i r h S e e i a z u - s ; + i H e o ' r f o A r e h u G o o w r c } } } A 1 n O s u ] i ) o = a t c N n v p i F f a e $ $ $ L C ) t Z t t o n e / 0 n s a U , e e t O o l n P u - e ) p u a T i ; n c n r h R r s H r r r * r u s t a H n x e h G e s t o P e A v _ e z a t u l E = s l r e i i h u A m u s ) e e s ; t ' o r u T s 0 e e n t o a t o b / + ' ) t e ] u a a A t ; P r d e h n t P r l > V ] ; ' t l t * e i f i r e . e A e i . ( ( ( ( ( ( ( ( ( ( ( A r ] p e ( a i e s a r i v R c ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( l l l l l l l l l l l R a u o 1 d = i t l t e T d l l l l l l l l l l l l l l l l l l o o o o o o o o o o o P t t u l - y 0 f r v n I e L o o o o o o o o o o o o o o o o o o n n n n n n n n n n n H e t o A ; e i P e w C t i n n n n n n n n n n n n n n n n n n g g g g g g g g g g g I ' g p s L s r r b u r i t U a c g g g g g g g g g g g g g g g g g g _ _ _ _ _ _ _ _ _ _ _ + ] a u s P t _ u b s l h L i e _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ n n n n n n n n n n n A p t H a e t l i l e A l n n n n n n n n n n n n n n n n n n n a a a a a a a a a a a L s , f A t = e i o R s s a a a a a a a a a a a a a a a a a a m m m m m m m m m m m P ' u ) e 0 c n b i . e m m m m m m m m m m m m m m m m m m e e e e e e e e e e e H ] f n / ; i e m P , e e e e e e e e e e e e e e e e e e = = = = = = = = = = = A o c ( r t L 3 p U = = = = = = = = = = = = = = = = = = ' ' ' ' ' ' ' ' ' ' ' ) l t 1 _ i u l R ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' o d a a a i i D C C w ; l i - r a c o s i P t p c o n o e w e r n i h r c t p c n i u u u n n i a o e o o A e n e f e e O e r o u a u f e f e o n o e o e r o e s t t t v v x l m i w n L a d n f d S c e n t t t f l f a m f u a s c e s c o o o e e i v p g i ! P l / s t u E h f s p u p i f i l i l r l t h f t m o c c c r r t o o h n H = o e h l w . n e u u r u c a c n a s n e - i u o o o s s - s t g A 0 r e , a o r m t a t i r i i a t w p o r p n n r r r e e S p i + ; a r l e p l e e e n l i w a u l e u u t r r r t r t o f A c m s L r S o n t ( g n - n t o o g s o n s s e e e E F i o e f r L = o i a e g c i l o a t r t e i n r e h g c h f l l l I r g b o P 0 d p c n e y e o o u p e r n k y e l a a a a S i l a p x m H ; i u e t n g t o l i e t ( e ( p s a c t t t ' s i b a A n f b n y t s s p ( u e n s e l d l r s s h b t i i i ) c t i r i d * = y l s h h h ( d u l t v t t r o o o h h o o o o o o h z l a n e E 0 i e o e o o l e t o p a e e g ( g c o o c r r n n n i m f P ; s , f c c o v g u n r r s l e c c k ' e d t e l i S u h k k g ( t t e a t d o d s k k ' s ) t p c l e y t o n I = e o s l d s t e g e s ' ' ) h e r o a m ' e s i L 0 d r p p d s o e ( o t e r v v ) ) a c e s s a ) r s t O ; r r e ) g v l u a d ( r h f t t n i N a b o o v ' o t r ( t s e s l e n e - i d P f o ) = y c c ) d s g p a l e s v s o o r p c h u n * 0 e e s e s u t o ) ) g i l e u i e i n ( ; s s s v ) d t e g ( ' s ' n o n s t l l c o S z s s ) ' e l ) s ) d g c h y a l t f I = ' s ) v g ( d o ) e p y e ' s i i G 0 ( ( ) s a l e g ' v r p ) t p o u M ; l l ) s p o v ) o p p r i s n ' A w o o ' t g d s d r r o c ' ] + = g g ) e ( s e t u o o c i c ) ( 0 a l d s v e c c c e t u V ; d d d o e ) a t e e s y r A e e y g v ' s d i s s s ' v R v v ) s y o s s ' ) e P s d s ) n ' ' ) ' H s s t e s ' s ' ) ) ) I s s a v ) ) t ) + ) ) t ' a A ' ' e s ) t L ) ) ) t e P ' e ) H ) a ' A d ) ) y / ( s 1 t - a A t L e P ) H ' A ) ) ) ;

NK_linear_costpush_commitment.mod

@ s v v v e S V T R R B A E @ @ @ p r s e # h a a a n I A H H H E L P # # # l a t v i o r r r d G R E O O T P S i R % e R e a m o a n c ; M P T Z A A H I f H R l H n n s c l c k e e e A H A A L n O H s O d n e h u l s p p p I = O d U O e U i e y _ a u ; s s s = = 3 = = = = N e = U f r _ s t d _ _ _ 5 / = f 0 = = _ m i e e a u z 1 ; 4 0 0 0 1 ; 0 o o m _ ; ; . . . / 9 V . @ b d u p " = = = 5 9 9 4 ; A 8 { j e l l N ; ; 9 ; L ; V e l ( a K 0 1 0 ; U % A c ( o n _ ; ; ; E % L t i r n l R R U i n d e i % % H e R E v s e r n s s O s e R e t r _ e h h U p s H r = o a u u o p O p u 1 b r t t n o U i m , j _ s n } e e i e c o o e s ; ^ n r c o f f e 2 t f t s f f t s = i t o t + = 3 v p o V ( 0 e u t A i ) ; s r p R i h a e T ) x _ n r H , c s s E p p o i i T l i m t s A a e m o t * n o r e x n u n y n ^ e ; . t 2 r i c ; _ n o c d c s o i " t s s - t c p - o u p u s u n h s t h = s B h s E o h T c o A k c ) k ; u n u d n e d r e r c o c m o m m i m t i m t e m n e t n t

NK_linear_costpush_discretion.mod

@ s v v v e S V T R R B A E @ @ @ p d # h a a a n I A H H H E L P # # # l i i o r r r d G R E O O T P S i R % e R e a s n c ; M P T Z A A H I f H R l H n n c c k e e e A H A A L n O H s O d n r l s p p p I = O d U O e U i e e u ; s s s = = 3 = = = = N e = U f r t d _ _ _ 5 / = f 0 = = _ i e a u z 1 ; 4 0 0 0 1 ; 0 o o ; ; . . . / 9 V . @ b n " = = = 5 9 9 4 ; A 8 { j a N ; ; 9 ; L ; V e r K 0 1 0 ; U % A c y _ ; ; ; E % L t _ l R R U i p i % % H e R E v o n s s O s e R e l e h h U p s H i a u u o p O p c r t t n o U i y _ s n } e ( c o o e s ; ^ i o f f e 2 n s f f t s t o t + t p o V r u t A u s r p R m h a e T e _ n r H n c s s E t o i i T s m t s A = m o t * ( o r e x i n y n ^ i . t 2 ) i c ; , n o c i c s o r " t s f - t = p - 3 u p 0 s u , h s p h l s a h s n o h n c o e k c r k _ u d n u i d n s e d c r e o r u c n o c t m o = m m B i m E t i T m t A e m , n e d t n i t s c r e t i o n a r y _ t o l = 1 e - 1 2 ) x p i e u ;

NK_linear_costpush_osr.mod

v ; v ; p ; m # # [ p [ x [ y [ w [ y [ a [ z [ r [ r [ y [ y [ y [ u [ u [ i e s % p e s v v v e S V T R R R B A E P P o p y e o o P P e o a a a / P K n i n n = n n = n = n n _ n _ n _ n _ n _ n n = n i n t i n h a a a n I A H H H H E L P H H p i n s s H H n s r a z c y y y y x r r i p n w u r e e e r A B R R R S V E T P P d C S A a e a = a a a = a c a R a = a e a r a e a g a g a = a R a d e A e d o r r r d G R E O O O T P S I I t e 1 d r r I I d r C C T i t ( I b M G F s _ _ _ _ _ i i e p p p a L E H H H I A P H H H e o I P m = m m + m m ; m H m m = m e m = m a m a m m H m = ; a l = ; c ; M P T Z A U A H I _ _ i ; ; _ _ _ _ ; ( o o h t h a t u E N o e n g e e r e x s s s m P T O O O G R S E I I l m _ P e B e x e ( e S e e O e R e S e a e P e p e p e K e O e d l 0 k e e e A H A A L P Y m 1 p p P Y o o p p i e t t R U r e a a e o _ _ _ e H A A Z U M P I T _ _ ( p Y A = E = ( = 1 = I = = A = H = I = l = S = = = = = A = U = P y ; s p p p I = O I _ ; a a I , p s y y s u i C t p a a z u t A A H L A P Y l o A = ' T ' + ' - ' G ' ' * ' O ' G ' = ' I ' y ' x ' P ' * ' H _ o y ; s s s = = 3 = = = = = N E = w r r E t r r r n F y s W H G a < l e I O I i s = ( N A D 1 P A l M r T a P Z D M D i D _ D - D + I P c u I I s t _ _ _ _ 5 / = e a a , 0 _ . i i i d r o I A e h r N E n i ( 1 e * y ) r L a A e F ( r * e A e i e Y e y e ( m A o ( n _ t h g a u z 1 ; 4 0 0 0 0 1 = 0 i m m , a o g g s e e u d T N n c t s e t 1 - w p n o P b * s P - e z f * f - f A f _ f y p * s - t P a e a ; ; . . . . / 9 . g s s 0 l p h h r e r i H T e o t a e + T i a - d H o c o 1 f ( i ( i p i * i n i _ l ( t 1 e I t r p = = = 5 9 5 9 4 ; 1 1 h _ , 2 g t t t f s O A r p p r V H K e m u A r + u p ) e - n y n i n a n a n e i y ) r E e = ; ; ; 9 ; . 2 t P b ; o i r t S o t U B a y : ) p A E e ( i 1 c ) V r r + r 1 i _ i e i ; i t i - c _ p + e * _ v 0 0 1 0 ; 5 5 s H o 2 = m ( ( e h o p r T I l $ $ $ $ $ $ $ $ $ $ $ ; a R T y + c / t * d A c o e e ) t e t ( t t ; t y i e u e s p m a ; ; ; ; ; ; ; I u ; 9 _ C C e e f t i L / $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ { { { { { { { { { { { r P A n 1 S i n e R e c p n i ( i + i i i _ t - s p t i o r y _ n , p ) ) t i b A I P f w { { { { { { { { { { { { { { { { { { \ \ \ \ \ \ \ \ \ \ \ a H ) e ) I I o ; m P e s c + o + o 1 o o o n y h s e d i _ P d i a s t w o u N T u w a z c y y \ y x r r i \ n w u \ \ \ a b r r r s e t p p m I * s + S G n a H c s _ e n 1 n ) n n n a d _ _ R + e a n I s r r 2 2 o e a n t Y Y b t w } } } } ^ t } ^ } } p } } } v l e h h h i a p h h h e ) ( i K M n I o s a e ) ; t e n p u a P l b a E ; f a 0 0 f r r ) e l h . $ $ $ $ { i { $ { $ $ i $ $ $ a a a p t o o o g r s e i i t / 1 a A C A f d * n ' ; s p e - r e o l ) f a r ; t H ; l t = m 1 2 t m e d W o i e g n l e e } r r r h a _ _ _ m p i t _ _ e ( - n P u * u ' n s ] h s f y e f u i ; i t o e I e = P 3 s 6 1 w s a A r c n a d } } $ e e e a } a { { a h l a { { r S B P r ( n ] ; t i _ f _ a f t n n ) c _ s 0 H 0 - a F n i R G u t e } } p p p } $ } z u } i o } \ y s I E P A v i c r f z i e l i p k i ; e R Y ; I ) 2 W r o o y n R F L N . } $ $ s s s $ $ } } $ } n $ p } G T h * e i t a t ; c ) c u i t s u * a y _ 0 i e f u A I i U o } y i i i } } $ } i } M A i x ' - i i e i + i i t n i s l y r = Y u 2 l : n l t N T c r $ } l l l $ $ $ } $ A * l ] p o n r e ( n e g o ' e ; e 0 ; 1 l t d a h T N e G g $ o o o } * T l + i n t ' n 1 t n g n ] ' ; y i y h a t e Y E n e / n n n $ ( H i e ' ' ] t - e t a v ] z y J o e t e ; S s n l _ _ _ 1 E p u ( ] ] R r p a o e _ y o M u i r h S e e i a z u - T s ; + i H e o ' r f r e _ h u G o o w r c } } } A A 1 n O s u ] i o = n a t c N n v p i F f a e $ $ $ L ) C ) t Z t t o n 0 a n s a U , e e t O o l n P / u - e ) p u a i ; t n c n r h R r s H T r r r * r u s t n x e h G e s t o P e A H v _ e z a t u = y s l r e i i h u A m u s ) E e e s ; t ' o r s 0 _ e e n t o a t o b / + T ' ) t e ] u a t ; e P r d e h n t P r l > V A ] ; ' t l e i f i r e . e A e i . ( ( ( ( ( ( ( ( ( ( ( A * r ] p a i y e s a r i v R c ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( l l l l l l l l l l l R ( a u o d = _ i t l t e T d l l l l l l l l l l l l l l l l l l o o o o o o o o o o o P 1 t t u y 0 g f r v n I e L o o o o o o o o o o o o o o o o o o n n n n n n n n n n n H - e t ; a e i P e w C t i n n n n n n n n n n n n n n n n n n g g g g g g g g g g g I A ' g p s r p r b u r i t U a c g g g g g g g g g g g g g g g g g g _ _ _ _ _ _ _ _ _ _ _ + L ] a u t _ u b s l h L i e _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ n n n n n n n n n n n A P p t a e x t l i l e A l n n n n n n n n n n n n n n n n n n n a a a a a a a a a a a L H s , t = e i o R s s a a a a a a a a a a a a a a a a a a m m m m m m m m m m m P A ' e 0 p c n b i . e m m m m m m m m m m m m m m m m m m e e e e e e e e e e e H ) ] f ; i i e m P , e e e e e e e e e e e e e e e e e e = = = = = = = = = = = A / o r e t L 3 p U = = = = = = = = = = = = = = = = = = ' ' ' ' ' ' ' ' ' ' ' ) ( l _ ; i u l R ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' o d a a a i i D C i o ; 1 l r a c o s i P t p c o n o e w e r n i h r c t p c n i u u u n n i a n u - o e n e f e e O e r o u a u f e f e o n o e o e r o e s t t t v v x l f t A w a d n f d S c e n t t t f l f a m f u a s c e s c o o o e e i v l p L i l / s t u E h f s p u p i f i l i l r l t h f t m o c c c r r t o a u P n = o e h l w . n e u u r u c a c n a s n e - i u o o o s s - t t H g 0 r e , a o r m t a t i r i i a t w p o r p n n r r r e e S p i A ; a r l e p l e e e n l i w a u l e u u t r r r t r o f + f c m s L r S o n t ( g n - n t o o g s o n s s e e e E F i o n e A r = o i a e g c i l o a t r t e i n r e h g c h f l l l I r g b e L o 0 d p c n e y e o o u p e r n k y e l a a a a S i l a f d P m ; i u e t n g t o l i e t ( e ( p s a c t t t ' s i b e b H n f b n y t s s p ( u e n s e l d l r s s h b t i i i ) c t i e a A d = y l s h h h ( d u l t v t t r o o o h h o o o o o o h z l d c * e 0 i e o e o o l e t o p a e e g ( g c o o c r r n n n i b k E f ; s , f c c o v g u n r r s l e c c k ' e d t a P i u h k k g ( t t e a t d o d s k k ' s ) t p c l e y c T S n = e o s l d s t e g e s ' ' ) h e r o a m ' k a I i 0 d r p p d s o e ( o t e r v v ) ) a c e s s a ) y L t ; r r e ) g v l u a d ( r h f t t n T l O i a b o o v ' o t r ( t s e s l e n e - i d a o N o = y c c ) d s g p a l e s v s o o r p c y r ) n 0 e e s e s u t o ) ) g i l e u i e l * ; s s s v ) d t e g ( ' s ' n o n s t l o R ( o z s s ) ' e l ) s ) d g c h y a r u S f = ' s ) v g ( d o ) e p y e ' s l I 0 ( ( ) s a l e g ' v r p ) t R e G u ; l l ) s p o v ) o p p r i u ' M ' w o o ' t g d s d r r o c l ) A ] = g g ) e ( s e t u o o c i e + 0 a l d s v e c c c e t ' ( ; d d d o e ) a t e e s y ) V e e y g v ' s d i s s s ' A v v ) s y o s s ' ) R s d s ) n ' ' ) P s s t e s ' s ' ) ) H s s a v ) ) t ) I ) ) t ' a + ' ' e s ) t A ) ) ) t e L ' e ) P ) a ' H d ) A y ) / s ( t 1 a - t A e L ) P ' H ) A ) ) ;

ComparePermanentAndTransitory.m

c d i d i d i d i c i f s p h p a t l s p h p a t l s p h p a t l s p h p a t l s p h p a t l s p h p a t l l y r y r y r y r l r i u l o l x i e u l o l x i e u l o l x i e u l o l x i e u l o l x i e u l o l x i e e n f n f n f n f o f g b o l o i t g b o l o i t g b o l o i t g b o l o i t g b o l o i t g b o l o i t g a a s a s a s a s s _ u p t d t s l e p t d t s l e p t d t s l e p t d t s l e p t d t s l e p t d t s l e r r _ r _ r _ r _ e h r l ( ( e n l ( ( e n l ( ( e n l ( ( e n l ( ( e n l ( ( e n e d e d e c e c o e o 0 o 0 t ( d o 0 o 0 t ( d o 0 o 0 t ( d o 0 o 0 t ( d o 0 o 0 t ( d o 0 o 0 t ( d a i i o o a r t : n : i ' ( t : n : i ' ( t : n : i ' ( t : n : i ' ( t : n : i ' ( t : n : i ' ( l N s N s N m N m l i ( i i g O { ( i i g I { ( i i g C { ( i i g O { ( i i g I { ( i i g C { l K c K c K m K m l z 2 r r h u ' 2 r r h n ' 2 r r h o ' 2 r r h u ' 2 r r h n ' 2 r r h o ' ; _ r _ r _ i _ i ; o , f f t t D , f f t f D , f f t s D , f f t t D , f f t f D , f f t s D l e l e l t l t n 3 _ _ p i 3 _ _ l i 3 _ _ t i 3 _ _ p i 3 _ _ l i 3 _ _ t i c i t i t i m i m , h h u s , h h a s , h h - s , h h u s , h h a s , h h - s l n i n i n e n e = 1 o o t c 2 o o t c 3 o o P c 4 o o t c 5 o o t c 6 o o P c o e o e o e n e n ) r r r ) r r i r ) r r u r ) r r r ) r r i r ) r r u r s a n a n a t a t 1 i i g e i i o e i i s e i i g e i i o e i i s e e r _ r _ r _ r _ 5 z z a t z z n t z z h t z z a t z z n t z z h t _ t _ p _ t _ p ; o o p i o o i o o i o o p i o o i o o i a c r c e c r c e n n o n n ( o n n p o n n o n n ( o n n p o l o a o r o a o r - - ( n - - t n - - r n - - ( n - - p n - - r n l s n s s s n s s 1 1 t ' 1 1 r ' 1 1 o ' 1 1 p ' 1 1 e ' 1 1 o ' ; t s t i t s t i , , r , , a , , c , , e , , r , , c p i p s p i p s i i a ' i i n ' i i e ' i i r ' i i s ' i i e ' c u t u t u t u t r r n C r r s C r r s C r r s C r r i C r r s C l s o s e s o s e f f s o f f i o f f s o f f i o f f s o f f s o c h r h n h r h n s s i m s s t m s s m s s s m s s t m s s m ; _ y _ t _ y _ t _ _ t m _ _ o m _ _ ( m _ _ t m _ _ e m _ _ ( m d d c c d c o i d c r i d c t i d c e i d c n i d c p i i = i = o = o = i o r t i o y t i o r t i o n t i o t t i o e t s s m m s m y m s m m s m a m s m t m s m m s m r m c c m m c m e c m s e c m n e c m e c m s e c m s e r r i i r i s n r i h n r i s n r i s n r i h n r i i n e e t t e t h t e t o t e t i t e t h t e t o t e t s t t . t . m . m . t m o ' t m c ' t m t ' t m o ' t m c ' t m t ' i i i i e i e i i e c } i e k } i e o } i e c } i e k } i e e } o r o r n r n r o n k , o n ) , o n r , o n k , o n ) , o n n , n f n f t f t f n t ) ' n t ' ' n t y ' n t ) ' n t ' ' n t t ' s s s s _ _ ' L _ _ ) L _ _ ) L _ _ ' L _ _ ) L _ _ ) L - ; - ; - ; - ; t t ) o t t o t t ' o p t ) o p p o p p ' o D D D D r r c r r c r r ) c e r c e e c e e ) c V V V V a a a a a a a a a r a a r r a r r a A A A A n n t n n t n n t s n t s s t s s t L L L L s s i s s i s s i i s i i i i i i i U U U U i i o i i o i i o s i o s s o s s o E E E E t t n t t n t t n t t n t t n t t n R R R R o o ' o o ' o o ' e o ' e e ' e e ' H H H H r r , r r , r r , n r , n n , n n , O O O O y y ' y y ' y y ' t y ' t t ' t t ' U U U U . . S . . S . . S . . S . . S . . S = = = = x x o p p o u u o x x o p p o u u o 0 0 0 0 _ _ u i i u _ _ u _ _ u i i u _ _ u . . e e t e e t e e t e e t e e t e e t 8 8 p p h _ _ h p p h p p h _ _ h p p h s s E e e E s s E s s E e e E s s E _ _ a p p a _ _ a _ _ a p p a _ _ a u u s s s s u u s u u s s s s u u s ( ( t _ _ t ( ( t ( ( t _ _ t ( ( t 1 1 ' u u ' 1 1 ' 1 1 ' u u ' 1 1 ' : : ) ( ( ) : : ) : : ) ( ( ) : : ) i i 1 1 i i i i 1 1 i i r r : : r r r r : : r r f f i i f f f f i i f f _ _ r r _ _ _ _ r r _ _ h h f f h h h h f f h h o o _ _ o o o o _ _ o o r r h h r r r r h h r r i i o o i i i i o o i i z z r r z z z z r r z z o o i i o o o o i i o o n n z z n n n n z z n n ) ) o o ) ) ) ) o o ) ) , , n n , , , , n n , , ' ' ) ) ' ' ' ' ) ) ' ' , , , , x ' ' x x ' ' x ' ' ' ' ) x ) ) x ) ' ' ) )