Estimating Macroeconomic Models: A Likelihood Approach

Similar documents
FEDERAL RESERVE BANK of ATLANTA

Bayesian Estimation of DSGE Models: Lessons from Second-order Approximations

Why Nonlinear/Non-Gaussian DSGE Models?

Filtering and Likelihood Inference

State Space Models and Filtering. Jesús Fernández-Villaverde University of Pennsylvania

Labor-Supply Shifts and Economic Fluctuations. Technical Appendix

Monetary and Exchange Rate Policy Under Remittance Fluctuations. Technical Appendix and Additional Results

DSGE-Models. Calibration and Introduction to Dynare. Institute of Econometrics and Economic Statistics

High-dimensional Problems in Finance and Economics. Thomas M. Mertens

Sequential Monte Carlo Methods (for DSGE Models)

Nonlinear and/or Non-normal Filtering. Jesús Fernández-Villaverde University of Pennsylvania

Combining Macroeconomic Models for Prediction

Point, Interval, and Density Forecast Evaluation of Linear versus Nonlinear DSGE Models

Economics Discussion Paper Series EDP Measuring monetary policy deviations from the Taylor rule

Estimation and Inference on Dynamic Panel Data Models with Stochastic Volatility

Solution Methods. Jesús Fernández-Villaverde. University of Pennsylvania. March 16, 2016

Bayesian Inference for DSGE Models. Lawrence J. Christiano

Estimation of moment-based models with latent variables

Ambiguous Business Cycles: Online Appendix

Inference. Jesús Fernández-Villaverde University of Pennsylvania

Sequential Monte Carlo Methods (for DSGE Models)

DSGE Methods. Estimation of DSGE models: Maximum Likelihood & Bayesian. Willi Mutschler, M.Sc.

Bayesian Computations for DSGE Models

A Discussion of Arouba, Cuba-Borda and Schorfheide: Macroeconomic Dynamics Near the ZLB: A Tale of Two Countries"

1 Bewley Economies with Aggregate Uncertainty

Non-nested model selection. in unstable environments

Financial Factors in Economic Fluctuations. Lawrence Christiano Roberto Motto Massimo Rostagno

BEAR 4.2. Introducing Stochastic Volatility, Time Varying Parameters and Time Varying Trends. A. Dieppe R. Legrand B. van Roye ECB.

Taylor Rules and Technology Shocks

Nonlinear DSGE model with Asymmetric Adjustment Costs under ZLB:

PANEL DISCUSSION: THE ROLE OF POTENTIAL OUTPUT IN POLICYMAKING

WORKING PAPER NO NON-STATIONARY HOURS IN A DSGE MODEL. Yongsung Chang Seoul National University. Taeyoung Doh University of Pennsylvania

New Notes on the Solow Growth Model

The Kalman filter, Nonlinear filtering, and Markov Chain Monte Carlo

Technical appendices: Business cycle accounting for the Japanese economy using the parameterized expectations algorithm

Chapter 6. Maximum Likelihood Analysis of Dynamic Stochastic General Equilibrium (DSGE) Models

Citation Working Paper Series, F-39:

A Modern Equilibrium Model. Jesús Fernández-Villaverde University of Pennsylvania

Online appendix to On the stability of the excess sensitivity of aggregate consumption growth in the US

New Information Response Functions

Inference Based on SVARs Identified with Sign and Zero Restrictions: Theory and Applications

FEDERAL RESERVE BANK of ATLANTA

Estimating and Accounting for the Output Gap with Large Bayesian Vector Autoregressions

The Metropolis-Hastings Algorithm. June 8, 2012

Federal Reserve Bank of Chicago

Can News be a Major Source of Aggregate Fluctuations?

Identifying the Monetary Policy Shock Christiano et al. (1999)

A Comparison of Business Cycle Regime Nowcasting Performance between Real-time and Revised Data. By Arabinda Basistha (West Virginia University)

S TICKY I NFORMATION Fabio Verona Bank of Finland, Monetary Policy and Research Department, Research Unit

Small Open Economy RBC Model Uribe, Chapter 4

Perceived productivity and the natural rate of interest

An Introduction to Perturbation Methods in Macroeconomics. Jesús Fernández-Villaverde University of Pennsylvania

Assessing Structural Convergence between Romanian Economy and Euro Area: A Bayesian Approach

Bonn Summer School Advances in Empirical Macroeconomics

Estimating a Nonlinear New Keynesian Model with the Zero Lower Bound for Japan

Piecewise Linear Continuous Approximations and Filtering for DSGE Models with Occasionally-Binding Constraints

Graduate Macro Theory II: Notes on Quantitative Analysis in DSGE Models

Signaling Effects of Monetary Policy

Macroeconomics II. Dynamic AD-AS model

Dynamic AD-AS model vs. AD-AS model Notes. Dynamic AD-AS model in a few words Notes. Notation to incorporate time-dimension Notes

Sequential Monte Carlo samplers for Bayesian DSGE models

DSGE models: problems and some personal solutions. Fabio Canova EUI and CEPR. March 2014

Matching DSGE models,vars, and state space models. Fabio Canova EUI and CEPR September 2012

Lecture 2: From Linear Regression to Kalman Filter and Beyond

Economic and VAR Shocks: What Can Go Wrong?

Perturbation Methods I: Basic Results

Stochastic Trends & Economic Fluctuations

Macroeconomics Theory II

Estimating Macroeconomic Models: A Likelihood Approach

Dynamic Factor Models and Factor Augmented Vector Autoregressions. Lawrence J. Christiano

DSGE Model Forecasting

The B.E. Journal of Macroeconomics

Macroeconomics Theory II

Formulating, Estimating, and Solving Dynamic Equilibrium Models: an Introduction. Jesús Fernández-Villaverde University of Pennsylvania

Lecture 2: From Linear Regression to Kalman Filter and Beyond

DSGE Models in a Liquidity Trap and Japan s Lost Decade

The 2001 recession displayed unique characteristics in comparison to other

Escaping the Great Recession 1

Dynamic stochastic general equilibrium models. December 4, 2007

Perturbation Methods for Markov-Switching DSGE Models

Bayesian inference based only on simulated likelihood: particle filter analysis of dynamic economic models

Potential Output, the Output Gap, and the Labor Wedge

1.2. Structural VARs

Time-Varying Vector Autoregressive Models with Structural Dynamic Factors

Adaptive Learning and Applications in Monetary Policy. Noah Williams

Latent variables and shocks contribution in DSGE models with occasionally binding constraints

The Natural Rate of Interest and its Usefulness for Monetary Policy

Dynamics and Monetary Policy in a Fair Wage Model of the Business Cycle

Research Division Federal Reserve Bank of St. Louis Working Paper Series

Dynamic Models with Serial Correlation: Particle Filter Based Estimation

Open Economy Macroeconomics: Theory, methods and applications

MA Advanced Macroeconomics: 7. The Real Business Cycle Model

Indeterminacy and Sunspots in Macroeconomics

Are Structural VARs Useful Guides for Developing Business Cycle Theories? by Larry Christiano

Chapter 1. Introduction. 1.1 Background

Sentiment Shocks as Drivers of Business Cycles

Housing and the Business Cycle

Resolving the Missing Deflation Puzzle. June 7, 2018

Asset pricing in DSGE models comparison of different approximation methods

Discussion Global Dynamics at the Zero Lower Bound

Transcription:

Estimating Macroeconomic Models: A Likelihood Approach Jesús Fernández-Villaverde University of Pennsylvania, NBER, and CEPR Juan Rubio-Ramírez Federal Reserve Bank of Atlanta Estimating Dynamic Macroeconomic Models p. 1/45

Outline We are interested in performing likelihood-based inference in nonlinear and/or non-normal DSGE models. Estimating Dynamic Macroeconomic Models p. 2/45

Outline We are interested in performing likelihood-based inference in nonlinear and/or non-normal DSGE models. We apply particle filtering to evaluate the likelihood of the model. Estimating Dynamic Macroeconomic Models p. 2/45

Outline We are interested in performing likelihood-based inference in nonlinear and/or non-normal DSGE models. We apply particle filtering to evaluate the likelihood of the model. We estimate a neoclassical business cycle model with investment-specific technological change and stochastic volatility. Estimating Dynamic Macroeconomic Models p. 2/45

What is the Particle Filter? The particle filter is a monte carlo algorithm that tracks the distribution of states conditional on a sequence of observables. Estimating Dynamic Macroeconomic Models p. 3/45

What is the Particle Filter? The particle filter is a monte carlo algorithm that tracks the distribution of states conditional on a sequence of observables. Original idea seems to be by Handschin and Mayne (1969). Estimating Dynamic Macroeconomic Models p. 3/45

What is the Particle Filter? The particle filter is a monte carlo algorithm that tracks the distribution of states conditional on a sequence of observables. Original idea seems to be by Handschin and Mayne (1969). Applied in financial econometrics by Shephard and coauthors. Estimating Dynamic Macroeconomic Models p. 3/45

What is the Particle Filter? The particle filter is a monte carlo algorithm that tracks the distribution of states conditional on a sequence of observables. Original idea seems to be by Handschin and Mayne (1969). Applied in financial econometrics by Shephard and coauthors. Reviewed on Doucet, De Freitas, and Gordon (2001). Estimating Dynamic Macroeconomic Models p. 3/45

What is the Particle Filter? The particle filter is a monte carlo algorithm that tracks the distribution of states conditional on a sequence of observables. Original idea seems to be by Handschin and Mayne (1969). Applied in financial econometrics by Shephard and coauthors. Reviewed on Doucet, De Freitas, and Gordon (2001). We modify the filter to be more flexible with shocks. Estimating Dynamic Macroeconomic Models p. 3/45

Why Are Nonlinearities Important? Estimating Dynamic Macroeconomic Models p. 4/45

Why Are Nonlinearities Important? Most DSGE models are nonlinear. Estimating Dynamic Macroeconomic Models p. 4/45

Why Are Nonlinearities Important? Most DSGE models are nonlinear. Common practice: estimate a linearized version. Estimating Dynamic Macroeconomic Models p. 4/45

Why Are Nonlinearities Important? Most DSGE models are nonlinear. Common practice: estimate a linearized version. Linearization eliminates asymmetries, threshold effects, precautionary behavior, big shocks,... Estimating Dynamic Macroeconomic Models p. 4/45

Why Are Nonlinearities Important? Most DSGE models are nonlinear. Common practice: estimate a linearized version. Linearization eliminates asymmetries, threshold effects, precautionary behavior, big shocks,... Moreover, linearization induces an approximation error. Estimating Dynamic Macroeconomic Models p. 4/45

Why Are Nonlinearities Important? Most DSGE models are nonlinear. Common practice: estimate a linearized version. Linearization eliminates asymmetries, threshold effects, precautionary behavior, big shocks,... Moreover, linearization induces an approximation error. This is worse than you may think: Estimating Dynamic Macroeconomic Models p. 4/45

Why Are Nonlinearities Important? Most DSGE models are nonlinear. Common practice: estimate a linearized version. Linearization eliminates asymmetries, threshold effects, precautionary behavior, big shocks,... Moreover, linearization induces an approximation error. This is worse than you may think: 1. Theoretical arguments. Estimating Dynamic Macroeconomic Models p. 4/45

Why Are Nonlinearities Important? Most DSGE models are nonlinear. Common practice: estimate a linearized version. Linearization eliminates asymmetries, threshold effects, precautionary behavior, big shocks,... Moreover, linearization induces an approximation error. This is worse than you may think: 1. Theoretical arguments. 2. Computational evidence. Estimating Dynamic Macroeconomic Models p. 4/45

Theoretical Arguments Estimating Dynamic Macroeconomic Models p. 5/45

Theoretical Arguments We develop these arguments in "Convergence Properties of the Likelihood of Computed Dynamic Models" (2006) and in some work-in-progress. Estimating Dynamic Macroeconomic Models p. 5/45

Theoretical Arguments We develop these arguments in "Convergence Properties of the Likelihood of Computed Dynamic Models" (2006) and in some work-in-progress. 1. Second-order errors in the approximated policy function induce first-order errors in the likelihood function. Estimating Dynamic Macroeconomic Models p. 5/45

Theoretical Arguments We develop these arguments in "Convergence Properties of the Likelihood of Computed Dynamic Models" (2006) and in some work-in-progress. 1. Second-order errors in the approximated policy function induce first-order errors in the likelihood function. 2. As the sample size grows, the error in the likelihood function also grows and we may have inconsistent point estimates. Estimating Dynamic Macroeconomic Models p. 5/45

Theoretical Arguments We develop these arguments in "Convergence Properties of the Likelihood of Computed Dynamic Models" (2006) and in some work-in-progress. 1. Second-order errors in the approximated policy function induce first-order errors in the likelihood function. 2. As the sample size grows, the error in the likelihood function also grows and we may have inconsistent point estimates. 3. Linearization complicates the identification of parameters. Estimating Dynamic Macroeconomic Models p. 5/45

Computational Evidence Estimating Dynamic Macroeconomic Models p. 6/45

Computational Evidence We build the computational evidence in several papers. Estimating Dynamic Macroeconomic Models p. 6/45

Computational Evidence We build the computational evidence in several papers. Several researchers (Ann, King, Schorfheide, Winschel) have gathered similar evidence after our paper first circulated. Estimating Dynamic Macroeconomic Models p. 6/45

Computational Evidence We build the computational evidence in several papers. Several researchers (Ann, King, Schorfheide, Winschel) have gathered similar evidence after our paper first circulated. 1. Big differences in the level of the likelihood. Estimating Dynamic Macroeconomic Models p. 6/45

Computational Evidence We build the computational evidence in several papers. Several researchers (Ann, King, Schorfheide, Winschel) have gathered similar evidence after our paper first circulated. 1. Big differences in the level of the likelihood. 2. Often, important differences in point estimates. Estimating Dynamic Macroeconomic Models p. 6/45

Computational Evidence We build the computational evidence in several papers. Several researchers (Ann, King, Schorfheide, Winschel) have gathered similar evidence after our paper first circulated. 1. Big differences in the level of the likelihood. 2. Often, important differences in point estimates. 3. Big differences in smoothed shocks and states. Estimating Dynamic Macroeconomic Models p. 6/45

Computational Evidence We build the computational evidence in several papers. Several researchers (Ann, King, Schorfheide, Winschel) have gathered similar evidence after our paper first circulated. 1. Big differences in the level of the likelihood. 2. Often, important differences in point estimates. 3. Big differences in smoothed shocks and states. 4. Better identification of parameters. Estimating Dynamic Macroeconomic Models p. 6/45

Why Are Non-normalities Important? Estimating Dynamic Macroeconomic Models p. 7/45

Why Are Non-normalities Important? Evidence of time-varying volatility in time series. Estimating Dynamic Macroeconomic Models p. 7/45

Why Are Non-normalities Important? Evidence of time-varying volatility in time series. Fundamental issue in financial econometrics. Estimating Dynamic Macroeconomic Models p. 7/45

Why Are Non-normalities Important? Evidence of time-varying volatility in time series. Fundamental issue in financial econometrics. However, macro has lagged despite Engle (1982). Estimating Dynamic Macroeconomic Models p. 7/45

Why Are Non-normalities Important? Evidence of time-varying volatility in time series. Fundamental issue in financial econometrics. However, macro has lagged despite Engle (1982). Evidence by Geweke (1993 and 1994). Estimating Dynamic Macroeconomic Models p. 7/45

Why Are Non-normalities Important? Evidence of time-varying volatility in time series. Fundamental issue in financial econometrics. However, macro has lagged despite Engle (1982). Evidence by Geweke (1993 and 1994). Recent discussion about the Great Moderation: Estimating Dynamic Macroeconomic Models p. 7/45

Why Are Non-normalities Important? Evidence of time-varying volatility in time series. Fundamental issue in financial econometrics. However, macro has lagged despite Engle (1982). Evidence by Geweke (1993 and 1994). Recent discussion about the Great Moderation: 1. Kim and Nelson (1999), McConell and Pérez-Quirós (2000), and Stock and Watson (2002). Estimating Dynamic Macroeconomic Models p. 7/45

Why Are Non-normalities Important? Evidence of time-varying volatility in time series. Fundamental issue in financial econometrics. However, macro has lagged despite Engle (1982). Evidence by Geweke (1993 and 1994). Recent discussion about the Great Moderation: 1. Kim and Nelson (1999), McConell and Pérez-Quirós (2000), and Stock and Watson (2002). 2. Good luck? Sims and Zha (2005). Estimating Dynamic Macroeconomic Models p. 7/45

Why Are Non-normalities Important? Evidence of time-varying volatility in time series. Fundamental issue in financial econometrics. However, macro has lagged despite Engle (1982). Evidence by Geweke (1993 and 1994). Recent discussion about the Great Moderation: 1. Kim and Nelson (1999), McConell and Pérez-Quirós (2000), and Stock and Watson (2002). 2. Good luck? Sims and Zha (2005). 3. Or good policy? Clarida, Galí, and Gertler (2000). Estimating Dynamic Macroeconomic Models p. 7/45

State Space Representation of the Model How do we evaluate the likelihood p ( y T ;γ ) of a DSGE model? Estimating Dynamic Macroeconomic Models p. 8/45

State Space Representation of the Model How do we evaluate the likelihood p ( y T ;γ ) of a DSGE model? Transition equation: S t = f (S t 1, W t ;γ) p (S t S t 1 ;γ) Estimating Dynamic Macroeconomic Models p. 8/45

State Space Representation of the Model How do we evaluate the likelihood p ( y T ;γ ) of a DSGE model? Transition equation: S t = f (S t 1, W t ;γ) p (S t S t 1 ;γ) Measurement equation: Y t = g (S t, V t ;γ) p (Y t S t ;γ) Estimating Dynamic Macroeconomic Models p. 8/45

State Space Representation of the Model How do we evaluate the likelihood p ( y T ;γ ) of a DSGE model? Transition equation: S t = f (S t 1, W t ;γ) p (S t S t 1 ;γ) Measurement equation: Y t = g (S t, V t ;γ) p (Y t S t ;γ) We want to track conditional density p ( S t y t 1 ;γ ). Estimating Dynamic Macroeconomic Models p. 8/45

Factorization of the Likelihood Why? p ( y T ;γ ) = = T p ( y t y t 1 ;γ ) t=1 T t=1 p (y t S t ;γ)p ( S t y t 1 ;γ ) ds t Estimating Dynamic Macroeconomic Models p. 9/45

Factorization of the Likelihood Why? p ( y T ;γ ) = = T p ( y t y t 1 ;γ ) t=1 T t=1 p (y t S t ;γ)p ( S t y t 1 ;γ ) ds t Knowledge of { p ( S t y t 1 ;γ )} T t=1 of the likelihood of the model. allows the evaluation Estimating Dynamic Macroeconomic Models p. 9/45

Tracking the Conditional Distribution of States S t Filtering problem: forecast and update. Estimating Dynamic Macroeconomic Models p. 10/45

Tracking the Conditional Distribution of States S t Filtering problem: forecast and update. 1. Forecast: Chapman-Kolmogorov equation p ( S t y t 1 ;γ ) = p (S t S t 1 ;γ)p ( S t 1 y t 1 ;γ ) ds t 1 Estimating Dynamic Macroeconomic Models p. 10/45

Tracking the Conditional Distribution of States S t Filtering problem: forecast and update. 1. Forecast: Chapman-Kolmogorov equation p ( S t y t 1 ;γ ) = p (S t S t 1 ;γ)p ( S t 1 y t 1 ;γ ) ds t 1 2. Update: Bayes theorem p ( S t y t ;γ ) = p (y t S t ;γ)p ( S t y t 1 ;γ ) p (y t y t 1 ;γ) where: p ( y t y t 1 ;γ ) = p (y t S t ;γ)p ( S t y t 1 ;γ ) ds t Estimating Dynamic Macroeconomic Models p. 10/45

Intuition Previous equations involve complicated integrals. Estimating Dynamic Macroeconomic Models p. 11/45

Intuition Previous equations involve complicated integrals. In a linear and normal world: mean and variance are sufficient statistics translation and spread. Estimating Dynamic Macroeconomic Models p. 11/45

Intuition Previous equations involve complicated integrals. In a linear and normal world: mean and variance are sufficient statistics translation and spread. This can be done by the Kalman filter. Estimating Dynamic Macroeconomic Models p. 11/45

Intuition Previous equations involve complicated integrals. In a linear and normal world: mean and variance are sufficient statistics translation and spread. This can be done by the Kalman filter. In a nonlinear and/or non-normal world, we need to carry the distribution deformation. Estimating Dynamic Macroeconomic Models p. 11/45

Intuition Previous equations involve complicated integrals. In a linear and normal world: mean and variance are sufficient statistics translation and spread. This can be done by the Kalman filter. In a nonlinear and/or non-normal world, we need to carry the distribution deformation. This can be done by the Particle filter. Estimating Dynamic Macroeconomic Models p. 11/45

Intuition Previous equations involve complicated integrals. In a linear and normal world: mean and variance are sufficient statistics translation and spread. This can be done by the Kalman filter. In a nonlinear and/or non-normal world, we need to carry the distribution deformation. This can be done by the Particle filter. Alternatives? Unscented Kalman filter, Grid Filter,... Estimating Dynamic Macroeconomic Models p. 11/45

A Law of Large Numbers T t=1 p (y t S t ;γ)p ( S t y t 1 ;γ ) ds t Estimating Dynamic Macroeconomic Models p. 12/45

A Law of Large Numbers T t=1 p (y t S t ;γ)p ( S t y t 1 ;γ ) ds t { { } } N T Suppose we have s i t t 1 i=1 t=1 { p ( S t y t 1 ;γ )} T t=1. Estimating Dynamic Macroeconomic Models p. 12/45

A Law of Large Numbers T t=1 p (y t S t ;γ)p ( S t y t 1 ;γ ) ds t { { } } N T Suppose we have s i t t 1 i=1 t=1 Then: p ( y T ;γ ) T t=1 1 N { p ( S t y t 1 ;γ )} T t=1. N ) p (y t s i t t 1 ;γ i=1 Estimating Dynamic Macroeconomic Models p. 12/45

A Law of Large Numbers T t=1 p (y t S t ;γ)p ( S t y t 1 ;γ ) ds t { { } } N T Suppose we have s i t t 1 i=1 t=1 Then: p ( y T ;γ ) T t=1 1 N { p ( S t y t 1 ;γ )} T t=1. N ) p (y t s i t t 1 ;γ Evaluating the likelihood function Drawing from density: i=1 { p ( St y t 1 ;γ )} T t=1 Estimating Dynamic Macroeconomic Models p. 12/45

Introducing Particles { } N s i t 1 t 1 N i.i.d. draws from p ( S t 1 y t 1 ;γ ). i=1 Estimating Dynamic Macroeconomic Models p. 13/45

Introducing Particles { } N s i t 1 t 1 N i.i.d. draws from p ( S t 1 y t 1 ;γ ). i=1 { } N Each s i t 1 t 1 is a particle and s i t 1 t 1 particles. i=1 a swarm of Estimating Dynamic Macroeconomic Models p. 13/45

Introducing Particles { } N s i t 1 t 1 N i.i.d. draws from p ( S t 1 y t 1 ;γ ). i=1 { } N Each s i t 1 t 1 is a particle and s i t 1 t 1 i=1 particles. { } N s i t t 1 N i.i.d. draws from p ( S t y t 1 ;γ ). i=1 a swarm of Estimating Dynamic Macroeconomic Models p. 13/45

Introducing Particles { } N s i t 1 t 1 N i.i.d. draws from p ( S t 1 y t 1 ;γ ). i=1 { } N Each s i t 1 t 1 is a particle and s i t 1 t 1 i=1 particles. { } N s i t t 1 N i.i.d. draws from p ( S t y t 1 ;γ ). i=1 { } N Each s i t t 1 is a proposed particle and s i t t 1 of proposed particles. a swarm of i=1 a swarm Estimating Dynamic Macroeconomic Models p. 13/45

Introducing Particles { } N s i t 1 t 1 N i.i.d. draws from p ( S t 1 y t 1 ;γ ). i=1 { } N Each s i t 1 t 1 is a particle and s i t 1 t 1 i=1 particles. { } N s i t t 1 N i.i.d. draws from p ( S t y t 1 ;γ ). i=1 { } N Each s i t t 1 is a proposed particle and s i t t 1 of proposed particles. Weight of each proposed particle: q i t = ) p (y t s i t t 1 ;γ ) N i=1 (y p t s i t t 1 ;γ a swarm of i=1 a swarm Estimating Dynamic Macroeconomic Models p. 13/45

A Proposition Let { } N s i t t 1 be a draw from p ( S t y t 1 ;γ ). i=1 Estimating Dynamic Macroeconomic Models p. 14/45

A Proposition Let { } N s i t t 1 be a draw from p ( S t y t 1 ;γ ). i=1 Let { s i } { } N N t i=1 be a draw with replacement from s i t t 1 i=1 and probabilities qt. i Estimating Dynamic Macroeconomic Models p. 14/45

A Proposition Let { } N s i t t 1 be a draw from p ( S t y t 1 ;γ ). i=1 Let { s i } { } N N t i=1 be a draw with replacement from s i t t 1 i=1 and probabilities qt. i Then, { s i t} N i=1 is a draw from p ( S t y t ;γ ) : { } N s i t t = { s i N t} i=1 i=1 Estimating Dynamic Macroeconomic Models p. 14/45

A Proposition Let { } N s i t t 1 be a draw from p ( S t y t 1 ;γ ). i=1 Let { s i } { } N N t i=1 be a draw with replacement from s i t t 1 i=1 and probabilities qt. i Then, { s i t} N i=1 is a draw from p ( S t y t ;γ ) : { } N s i t t = { s i N t} i=1 i=1 Proof: Importance sampling and Bayes theorem. Estimating Dynamic Macroeconomic Models p. 14/45

Importance of the Proposition Estimating Dynamic Macroeconomic Models p. 15/45

Importance of the Proposition { } N 1. Update: We can use a draw s i t t 1 from p ( S t y t 1 ;γ ) i=1 { } N to get a draw s i t t from p ( S t y t ;γ ). i=1 Estimating Dynamic Macroeconomic Models p. 15/45

Importance of the Proposition { } N 1. Update: We can use a draw s i t t 1 from p ( S t y t 1 ;γ ) i=1 { } N to get a draw s i t t from p ( S t y t ;γ ). i=1 { } N 2. Forecast: We can use a draw s i t t from p ( S t y t ;γ ), a i=1 draw from p (W t+1 ;γ), and S t+1 = f (S t, W t+1 ;γ) to get a { } N draw s i t+1 t i=1. Estimating Dynamic Macroeconomic Models p. 15/45

Particle Filtering I Estimating Dynamic Macroeconomic Models p. 16/45

Particle Filtering I { } N Step 0, Initialization: Sample N values s i 1 0 from i=1 p (S 1 ;γ). Go to step 2. Estimating Dynamic Macroeconomic Models p. 16/45

Particle Filtering I { } N Step 0, Initialization: Sample N values s i 1 0 from i=1 p (S 1 ;γ). Go to step 2. Step 1, Forecast: Sample N values p (W t ;γ) and S t = f (S t 1, W t ;γ). { } N s i t t 1 i=1 with Estimating Dynamic Macroeconomic Models p. 16/45

Particle Filtering I { } N Step 0, Initialization: Sample N values s i 1 0 from i=1 p (S 1 ;γ). Go to step 2. Step 1, Forecast: Sample N values p (W t ;γ) and S t = f (S t 1, W t ;γ). { } N s i t t 1 i=1 with Step 2, Weighting: Assign to each draw s i t t 1 the weight q i t. Estimating Dynamic Macroeconomic Models p. 16/45

Particle Filtering I { } N Step 0, Initialization: Sample N values s i 1 0 from i=1 p (S 1 ;γ). Go to step 2. Step 1, Forecast: Sample N values p (W t ;γ) and S t = f (S t 1, W t ;γ). { } N s i t t 1 i=1 with Step 2, Weighting: Assign to each draw s i t t 1 the weight q i t. { } N Step 3, Update: Draw s i t t with replacement from i=1 { } N s i t t 1 with probabilities { q i } N t. If t < T set i=1 i=1 t t + 1 and go to step 1. Otherwise stop. Estimating Dynamic Macroeconomic Models p. 16/45

Particle Filtering II Use { { } } N T s i t t 1 i=1 t=1 to compute: p ( y T ;γ ) T t=1 1 N N ) p (y t s i t t 1 ;γ i=1 We can filter, forecast, and smooth Estimating Dynamic Macroeconomic Models p. 17/45

An Application: a Business Cycle Model Estimating Dynamic Macroeconomic Models p. 18/45

An Application: a Business Cycle Model I want to show the power of the particle filter. Estimating Dynamic Macroeconomic Models p. 18/45

An Application: a Business Cycle Model I want to show the power of the particle filter. And learn something about the evolution of the U.S. aggregate fluctuations over the last decades. Estimating Dynamic Macroeconomic Models p. 18/45

An Application: a Business Cycle Model I want to show the power of the particle filter. And learn something about the evolution of the U.S. aggregate fluctuations over the last decades. A business cycle model with: Estimating Dynamic Macroeconomic Models p. 18/45

An Application: a Business Cycle Model I want to show the power of the particle filter. And learn something about the evolution of the U.S. aggregate fluctuations over the last decades. A business cycle model with: 1. Investment-specific technological change. Greenwood, Herkowitz, and Krusell (1997 and 2000) Estimating Dynamic Macroeconomic Models p. 18/45

0.03 A Look at the Data 0.02 0.01 0-0.01-0.02-0.03 Per Capita Output Growth Q/Q 0.08 0.06 0.04 0.02 0-0.02-0.04-0.06-0.08 Per Capita Investment Growth Q/Q -0.04-0.1-1 Per Capita Log Hours Worked 0.04 Relative Price of Investment Growth Q/Q -1.05 0.03 0.02-1.1 0.01 0-1.15-0.01-1.2-0.02-0.03-1.25-0.04 Estimating Dynamic Macroeconomic Models p. 19/45

An Application: a Business Cycle Model I want to show the power of the particle filter. And learn something about the evolution of the U.S. aggregate fluctuations over the last decades. A business cycle model with: 1. Investment-specific technological change. Greenwood, Herkowitz, and Krusell (1997 and 2000) 2. Stochastic volatility. Estimating Dynamic Macroeconomic Models p. 20/45

2.5 Another x 10-3 look at the Data: The Variance Output Investment Labor Price 2 1.5 1 0.5 0 1960 1965 1970 1975 1980 1985 1990 1995 Estimating Dynamic Macroeconomic Models p. 21/45

Households Representative household with utility function: E 0 t=0 ( ) β t e d t log C t + ψ log (1 L t ) Law of motion of d t, the preference shock: d t = ρd t 1 + σ dt ε dt, ε dt N (0,1) We will explain later the law of motion of σ dt. Estimating Dynamic Macroeconomic Models p. 22/45

Technology Final Good: C t + X t = A t K α t L 1 α t Law of motion of capital: K t+1 = (1 δ) K t + V t X t Shocks: log A t = ζ + log A t 1 + σ at ε at, ζ 0 and ε at N (0,1) log V t = υ + log V t 1 + σ vt ε υt, υ 0 and ε υt N (0,1) Estimating Dynamic Macroeconomic Models p. 23/45

Stochastic Volatility We follow a standard specification: log σ dt = (1 λ d )log σ d + λ d log σ dt 1 + τ d η dt and η dt N (0,1) log σ at = (1 λ a ) log σ a + λ a log σ at 1 + τ a η at and η at N (0,1) log σ υt = (1 λ υ ) log σ υ + λ υ log σ υt 1 + τ υ η υt and η υt N (0,1) Estimating Dynamic Macroeconomic Models p. 24/45

Analysis of the Model It can be shown that along a balanced growth path the following variables are stationary: Y t /Z t, C t /Z t, X t /Z t, K t+1 /(Z t V t ), (Y t /L t )/Z t, and L t where Z t = A 1/(1 α) t V α/(1 α) t. Estimating Dynamic Macroeconomic Models p. 25/45

Analysis of the Model It can be shown that along a balanced growth path the following variables are stationary: Y t /Z t, C t /Z t, X t /Z t, K t+1 /(Z t V t ), (Y t /L t )/Z t, and L t where Z t = A 1/(1 α) t V α/(1 α) t. The growth rates for the exogenous shocks are: log A t log A t 1 log V t log V t 1 = ζ + σ at ε at = υ + σ vt ε υt Estimating Dynamic Macroeconomic Models p. 25/45

Analysis of the Model It can be shown that along a balanced growth path the following variables are stationary: Y t /Z t, C t /Z t, X t /Z t, K t+1 /(Z t V t ), (Y t /L t )/Z t, and L t where Z t = A 1/(1 α) t V α/(1 α) t. The growth rates for the exogenous shocks are: log A t log A t 1 log V t log V t 1 = ζ + σ at ε at = υ + σ vt ε υt Therefore, Y t, C t, X t, and Y t /L t grow at rate (ζ + αυ)/(1 α), K t grows at rate (ζ + υ)/(1 α), and L t is stationary. Estimating Dynamic Macroeconomic Models p. 25/45

Transforming the Model The model is nonstationary because of the presence of two unit roots, one in each technological process. We need to transform the model into a stationary problem. We scale variables as C t = C t Z t, X t = X t Z t, and K t = Then: K t Z t V t 1. E 0 t=0 ( β t e d t log C ) t + ψ log(1 L t ) C t + e γ+αυ+σ at ε at +ασ υt ε υt 1 α Kt+1 = e γ+σ atε at Kα t L 1 α t + (1 δ) e υ σ υtε υt Kt Estimating Dynamic Macroeconomic Models p. 26/45

Equilibrium Conditions 1. Euler equation: e d t e γ+αυ+ε at +ασ υt ε υt 1 α = C t βe t e d t+1 C t+1 ( αe γ+σ at+1ε at+1 Kα t+1 L 1 α t+1 + (1 δ) e υ σ υt+1ε υt+1 ) 2. A labor supply condition: 3. The resource constraint: ψ ed t Ct = (1 α) e γ+σ atε at Kα 1 L t L α t t C t +e γ+αυ+σ at ε at +ασ υt ε υt 1 α Kt+1 = e γ+σ atε at Kα t L 1 α t +(1 δ) e υ σ υtε υt Kt Estimating Dynamic Macroeconomic Models p. 27/45

The Steady State We can find the steady state of the transformed model. We have a cointegration relation between output and investment in nominal terms: ) (e ζ+αυ 1 α (1 δ) e υ X ss Ỹ ss = α exp( ζ+αυ 1 α ) β (1 δ) exp ( υ) Let log K t+1 = log K t+1 K ss and log X t = log X t X ss Estimating Dynamic Macroeconomic Models p. 28/45

Solution We solve the model using 2nd order perturbation. Other non-linear solution methods are possible. Aruoba, Fernández-Villaverde, and Rubio-Ramírez (2005). Structure of a 2nd order approximation: log K t+1 = Ψ k1 s t + 1 2 s tψ k2 s t log X t = Ψ x1 s t + 1 2 s tψ x2 s t log L t = Ψ l1 s t + 1 2 s tψ l2 s t where: s t = (1,log K t, d t, σ at ε at, σ υt ε υt,log σ dt,log σ at,log σ υt ) Estimating Dynamic Macroeconomic Models p. 29/45

State Space Representation I: Transition Equation f 1 (S t, W t ) = 1 f 2 (S t, W t ) = Ψ k1 s t + 1 2 s tψ k2 s t f 3 (S t, W t ) = ((1 λ a ) log σ a + λ a log σ at + τ a η at+1 ) ε at+1 f 4 (S t, W t ) = ((1 λ υ ) log σ υ + λ υ log σ υt + τ υ η υt+1 ) ε υt+1 f 5 (S t, W t ) = ρd t + e (1 λ d) log σ d +λ d log σ dt +τ d η dt+1 ε dt+1 f 6 (S t, W t ) = (1 λ a )log σ a + λ a log σ at + τ a η at+1 f 7 (S t, W t ) = (1 λ υ )log σ υ + λ υ log σ υt + τ υ η υt+1 f 8 (S t, W t ) = (1 λ d ) log σ d + λ d log σ dt + τ d η dt+1 f 9 16 (S t, W t ) = s t where S t = (s t, s t 1 ). Estimating Dynamic Macroeconomic Models p. 30/45

State Space Representation II: Measurement Equation log P t log Y t log X t log L t = υ γ+αυ 1 α γ+αυ 1 α log L ss + Ψ l3 ε υt + Ψ y1 (s t s t 1 ) + 1 ( 2 s t Ψ y2 s t s t 1 Ψ ) y2s t 1 + σ at 1 ε at 1 +ασ υt 1 ε υt 1 1 α Ψ x1 (s t s t 1 ) + 1 ( 2 s t Ψ x2 s t 1 2 s t 1 Ψ ) x2s t 1 + σ at 1 ε at 1 +ασ υt 1 ε υt 1 1 α Ψ l1 s t + 1 2 s tψ l2 s t + 0 ǫ 1t ǫ 2t ǫ 3t Estimating Dynamic Macroeconomic Models p. 31/45

Performing Likelihood-Based Inference Time series: 1. Relative price of capital, output, investment, and hours. 2. Sample: 1955:Q1 to 2000:Q4. Vector of parameters γ is: (ρ, β, ψ, α, δ, υ, ζ, τ d, τ a, τ υ,σ d,σ a,σ υ, λ a, λ υ, λ d, σ ǫ 1, σǫ 2, σǫ 3 ) Use a Random-walk Metropolis-Hastings to explore the likelihood: Classical and Bayesian. Estimating Dynamic Macroeconomic Models p. 32/45

MLE Table 5.1: Maximum Likelihood Estimates Parameter PointEstimate StandardError(x10 3 ) 0.967 0.999 2.343 8.960E-003 3.594E-005. 3.743 0.460 6.825 0.828 2.254 a d a d a d 1² 2² 3² 7.120E-002 1.589 7.772E-003 2.940 5.653E-002 2.034 4.008E-004 0.692 8.523E-003 0.101 5.016E-003 2.344 4.460E-002 6.788 0.998 8.248 0.998 2.302 1.031E-005 0.424 1.024E-004 0.495 1.110E-005 0.082 Estimating Dynamic Macroeconomic Models p. 33/45

Figure 6.1: Model versus Data 0.03 Fit of the Data 0.02 Output 0.1 0.05 Investment 0.01 0 0-0.01-0.05-0.02-0.03-0.04 Model Data -0.1-1 Hours 0.04 Relative Price of Investment -1.05 0.03 0.02-1.1 0.01 0-1.15-0.01-1.2-0.02-0.03-1.25-0.04 Estimating Dynamic Macroeconomic Models p. 34/45

Figure 6.2: Smoothed Capital and Shocks 0.7 Shocks 0.6 0.5 0.4 0.3 0.2 0.1 Mean of (+/- std) k t 1.5 x Mean of (+/- std) 10-3 a,t a,t 1 0.5 0-0.5-1 0-1.5 Mean of (+/- std),t,t Mean of (+/- std) d t 0.03 0.02 0.02 0.01 0.01 0 0-0.01-0.02-0.01-0.02-0.03-0.03-0.04-0.04-0.05-0.05 Estimating Dynamic Macroeconomic Models p. 35/45

Figure 6.3: Smoothed Volatilities -7.5 Volatility I Mean of (+/- std) a,t -4.65 Mean of (+/- std),t -7.6-4.7-7.7-4.75-7.8-7.9-4.8-8 -4.85-4.9-4.6 Mean of (+/- std) d,t -4.8-5 -5.2-5.4-5.6-5.8-6 -6.2 Estimating Dynamic Macroeconomic Models p. 36/45

Volatility II 10 x 10-3 Mean of (+/- std) Output 9 Figure 6.4: Instantaneous Standard Deviation Mean of (+/- std) Investment 0.045 0.04 8 0.035 7 6 0.03 5 0.025 4 0.02 3 0.015 0.06 0.055 0.05 0.045 0.04 0.035 0.03 0.025 Mean of (+/- std) Hours 10 x 10-3 Mean of (+/- std) Relative Price of Investment 9.5 9 8.5 8 7.5 7 0.02 6.5 Estimating Dynamic Macroeconomic Models p. 37/45

Counterfactuals I 10 x 10-3 Output 8 6 4 Actual Counterfactual Figure 6.5: Counterfactual Exercise 1 10 x 10-3 Output 8 6 4 10 x 10-3 Output 8 6 4 2 2 2 0.05 Investment 0.05 Investment 0.05 Investment 0.04 0.04 0.04 0.03 0.03 0.03 0.02 0.02 0.02 0.01 0.01 0.01 0.045 Hours 0.045 Hours 0.045 Hours 0.04 0.04 0.04 0.035 0.035 0.035 0.03 0.03 0.03 10 x 10-3 Relative Price of Investment 10 x 10-3 Relative Price of Investment 10 x 10-3 Relative Price of Investment 9 9 9 8 8 8 7 7 7 Estimating Dynamic Macroeconomic Models p. 38/45

Counterfactuals II 10 x 10-3 Output 8 6 4 Actual Counterfactual Figure 6.6: Counterfactual Exercise 2 10 x 10-3 Output 8 6 4 10 x 10-3 Output 8 6 4 2 2 2 0.05 Investment 0.05 Investment 0.05 Investment 0.04 0.04 0.04 0.03 0.03 0.03 0.02 0.02 0.02 0.01 0.01 0.01 0.045 Hours 0.045 Hours 0.045 Hours 0.04 0.04 0.04 0.035 0.035 0.035 0.03 0.03 0.03 10 x 10-3 Relative Price of Investment 10 x 10-3 Relative Price of Investment 10 x 10-3 Relative Price of Investment 9 9 9 8 8 8 7 7 7 Estimating Dynamic Macroeconomic Models p. 39/45

Counterfactuals III 10 x 10-3 Output 8 6 4 Actual Counterfactual Figure 6.7: Counterfactual Exercise 3 10 x 10-3 Output 8 6 4 10 x 10-3 Output 8 6 4 2 2 2 0.05 Investment 0.05 Investment 0.05 Investment 0.04 0.04 0.04 0.03 0.03 0.03 0.02 0.02 0.02 0.01 0.01 0.01 0.045 Hours 0.045 Hours 0.045 Hours 0.04 0.04 0.04 0.035 0.035 0.035 0.03 0.03 0.03 10 x 10-3 Relative Price of Investment 10 x 10-3 Relative Price of Investment 10 x 10-3 Relative Price of Investment 9 9 9 8 8 8 7 7 7 Estimating Dynamic Macroeconomic Models p. 40/45

Are Nonlinearities and Non-normalities Important? Estimating Dynamic Macroeconomic Models p. 41/45

Are Nonlinearities and Non-normalities Important? We estimate four version of the model: Table 7.1: Versions of the Model Solution No Stochastic Volatility Stochastic Volatility Linear Version 1 Version 2 Quadratic Version 3 Benchmark Estimating Dynamic Macroeconomic Models p. 41/45

Are Nonlinearities and Non-normalities Important? We estimate four version of the model: Table 7.1: Versions of the Model Solution No Stochastic Volatility Stochastic Volatility Linear Version 1 Version 2 Quadratic Version 3 Benchmark We use Likelihood Ratio tests to compare models, Rivers and Vuong (2002). Estimating Dynamic Macroeconomic Models p. 41/45

Are Nonlinearities and Non-normalities Important? We estimate four version of the model: Table 7.1: Versions of the Model Solution No Stochastic Volatility Stochastic Volatility Linear Version 1 Version 2 Quadratic Version 3 Benchmark We use Likelihood Ratio tests to compare models, Rivers and Vuong (2002). Loglike benchmark: 2350.6, loglike version 2: 2230.4 Estimating Dynamic Macroeconomic Models p. 41/45

Figure 7.1: Comparison of Smoothed Capital and Shocks 0.6 Comparison with Linear Model linear I 0.5 0.4 0.3 0.2 0.1 0-0.1 Mean of (+/- std) k t quadratic 1 x Mean of (+/- std) 10-3 a,t a,t 0.5 0-0.5-1 -1.5-2 -2.5-0.2-3 Mean of (+/- std),t,t Mean of (+/- std) d t 0.03 0.05 0.02 0.01 0 0-0.05-0.01-0.02-0.1-0.03-0.04-0.15-0.05-0.2 Estimating Dynamic Macroeconomic Models p. 42/45

Figure 7.2: Comparison of Smoothed Volatilities -6 Comparison with Linear Model linear II -6.5-7 -7.5 Mean of (+/- std) a,t quadratic -4.66-4.68-4.7-4.72-4.74-4.76-4.78-4.8-4.82 Mean of (+/- std),t -8-4.84 Mean of (+/- std) d,t -3.5-4 -4.5-5 -5.5-6 Estimating Dynamic Macroeconomic Models p. 43/45

What are We Doing Now? We are estimating a richer DSGE model with: 1. Nominal and real rigidities. 2. Monetary and fiscal policy. 3. Stochastic volatility. 4. Parameter drifting. Estimating Dynamic Macroeconomic Models p. 44/45

What are We Doing Now? We are estimating a richer DSGE model with: 1. Nominal and real rigidities. 2. Monetary and fiscal policy. 3. Stochastic volatility. 4. Parameter drifting. We are working on a model with micro heterogeneity. Estimating Dynamic Macroeconomic Models p. 44/45

What are We Doing Now? We are estimating a richer DSGE model with: 1. Nominal and real rigidities. 2. Monetary and fiscal policy. 3. Stochastic volatility. 4. Parameter drifting. We are working on a model with micro heterogeneity. We are exploring the semi-nonparametric estimation of DSGE models. Estimating Dynamic Macroeconomic Models p. 44/45

Conclusions 1. Particle filtering is a general purpose and efficient method to estimate DSGE models. 2. We learned about the importance of stochastic volatility to account for U.S. Business Cycle. 3. Much exciting work to do in the next few years! Estimating Dynamic Macroeconomic Models p. 45/45