Gaussian Mixture Approximations of Impulse Responses and the Non-Linear Effects of Monetary Shocks

Transcription

1 Gaussian Mixture Approximations of Impulse Responses and the Non-Linear Effects of Monetary Shocks Regis Barnichon (CREI, Universitat Pompeu Fabra) Christian Matthes (Richmond Fed)

2 Effects of monetary policy on the real economy Broad consensus on the average effect of monetary policy on output A monetary contraction (expansion) leads to a decline (increase) in output How about non-linear effects of monetary policy? 1. Sign-dependence? Pushing vs. pulling on a string 2. State-dependence? Is pushing more effective in recessions than in expansions?

3 Non-linear effects of monetary policy No consensus in the literature Not surprising: standard approach to identify effect of money shock is based on linear models SVAR with recursive identification scheme (Christiano, Eichenbaum and Evans, 1999) SVAR with external instrument (Gertler-Karadi, 214)

4 This paper How an economy in a given state responds to a shock of a given sign? New method: Does not assume existence of VAR representation Approximate the impulse responses with Gaussian basis functions Directly parametrically estimate the structural MA representation of the system (great for prior elicitation) Easily amenable to estimation of non-linear effects

5 Model Y t a vector of stationary macro variables Structural vector MA representation of the economy with K finite or infinite. Y t = K Ψ k ε t k (1) k= With {ε t } the structural shocks affecting the economy (E ε t = and E ε t ε t = I) We assume throughout that Ψ is lower trianglular If Ψ(L) invertible, (1) can be estimated from a VAR on Y t

6 Approximating IRFs with mixtures of Gaussians Intuition: IRFs often monotonic or hump-shaped Idea: Use Gaussian basis functions to approximate (parametrize) the impulse response functions

7 Approximating IRFs with mixtures of Gaussians Approximate IRF ψ(k) by a sum of Gaussian basis functions: ψ(k) N a n e n=1 k bn ( ) cn 2 With N Gaussians, we have a GMA(N) In practice, only a very small number of Gaussian basis functions are needed

8 Interest rate Price level Unemployment Example/Motivation Standard SVAR with (U, d ln P, r) over IRFs to a monetary contraction

9 GMA(1): a one-gaussian approximation Approximate the IRFs with one Gaussian function { k b ψ(k) = ae ( c ) 2, k > ψ() unconstrained

10 Interest rate Price level Unemployment Fit of a one-gaussian approximation Approximation with a GMA(1) VAR GMA Quarters

11 GMA(1): a one-gaussian approximation Approximate the IRFs with one Gaussian function { k b ψ(k) = ae ( c ) 2, k > ψ() unconstrained Only 4 free parameters to capture one IRF {ψ(k)} k=1 For a trivariate system, nb of free parameters: GMA(1): 39 VAR(12): 117

12 Advantage of a one-gaussian parametrization Parsimonious and yet flexible k b ψ(k) = ae ( c ) 2, k > Can summarize IRF with 3 parameters: The peak effect of the shock: a The time needed to reach the peak effect: b The "persistence" of the effect of the shock: c

13 Using one Gaussian: GMA(1) t b ψ(t) = ae ( c ) 2 t

14 Using one Gaussian: GMA(1) t b ψ(t) = ae ( c ) 2 a: peak effect of the shock a t

15 Using one Gaussian: GMA(1) b t b ψ(t) = ae ( c ) 2 b: time to peak effect t

16 Using one Gaussian: GMA(1) t b ψ(t) = ae ( c ) 2 c: persistence in effect of shock a 2 c ln 2 t

17 Using one Gaussian: GMA(1) t b ψ(t) = ae ( c ) 2 b a c ln 2 a 2 t

18 Using two Gaussians: GMA(2) t

19 Using two Gaussians: GMA(2) ( ψ(t) = a 1 e t b1 c 1 ) 2 + a 2 e ( t b2 c 2 ) 2 t

20 Interest rate Interest rate Price level Price level Unemployment Unemployment GMA(2): a two-gaussian approximation Approximation with a GMA(1) Approximation with a GMA(2) VAR GMA.8.4 VAR GMA Quarters Quarters

21 Interest rate Inflation Unemployment The Gaussian basis functions Gaussian basis functions GB 1.4 GB Quarters

22 Introducing non-linearities Non-linear vector MA representation of the economy Y t = Ψ k (ε t k, Z t k )ε t k k= With {ε t } the structural shocks affecting the economy And Z t a function of past endogenous variables { } Y t j j> or a function of exogenous variables

23 Asymmetry (1) Asymmetric model Y t = ] [Ψ + k 1 ε t k >ε t k + Ψ k 1 ε t k <ε t k k=

24 Asymmetry (2) Asymmetric GMA(N) ( ) ψ + (k) = N a n + e k b + 2 n c n + n=1 ψ (k) similarly

25 Asymmetry and state-dependence GMA(N) with asymmetry and state dependence ψ + (k) = (1 + γ + z t k ) N n=1 ( ) a n + e k b n + 2 c n + Parameter γ + captures state dependence at time of shock The state of the cycle allowed to stretch/contract the IRF, but the overall shape is fixed Analogy with varying-coeffi cient model

26 Bayesian estimation Similar to linear case Taking into account Jacobian of mapping between (Gaussian) structural shocks and (non-gaussian) forecast errors

27 Bayesian estimation Construct conditional likelihood assuming Normal shocks Estimation routine: Multiple-block Metropolis-Hastings algorithm Initialize chain with GMA parameters chosen to fit VAR-based IRFs

28 Constructing the likelihood Assume Normal shocks ε t Parameter vector θ Decompose the likelihood using p(y 1,..., Y T θ) = p(y T Y 1,..., Y T 1, θ)... p(y 2 Y 1, θ)p(y 1 θ) Proceed recursively with p(y t+1 Y 1,..., Y t, θ) = p(ψ ε t+1 Y 1,..., Y t, θ) And obtain shock t + 1 from Ψ ε t+1 = Y t+1 K k= Ψ k+1 ε t k. Need Ψ to be invertible, guaranteed by struct. identifying assumption (Ψ lower triangular) Initialize recursion by setting first K shocks to zero

29 Other technical issues Loose priors centered on VAR-based IRFs Choosing the order of the GMA (ie the number of Gaussians): => compare marginal densities of GMAs of increasing order N Dealing with non-stationary data First-difference Can allow for deterministic trend

30 Back to monetary policy Are there non-linearities in the effects of monetary shocks? Small-scale model similar to Primiceri (25) (U, Π,r) over Unemployment rate, PCE price index, fed funds rate

31 Expansionary shock Contractionary shock Asymmetry: IRFs from GMA(2) Unemployment Price level Interest rate VAR ( ) Unemployment ( ) Price level ( ) Interest rate VAR

32 Robustness Same results when: Using output growth Using detrended output

33 Asymmetry and state-dependence GMA(1) Similarly for ψ (k) z t : Unemployment rate ψ + (k) = (1 + γ + z t k )a + e ( k b + c + ) 2

34

35 Unemployment Distribution of monetary shocks over cycle 12 1 UR Monetary shocks

36 Shocks Distribution Economic slack and effects of monetary shocks.2.1 Peak effect on U Contractionary shock Expansionary shock VAR GMA VAR GMA Peak effect on Π Unemployment rate Unemployment rate

37 Robustness Same results with alternative business cycle indicator z t : Detrended unemployment Output growth

38 Take away: monetary policy has highly non-linear effects Asymmetry: large effect of contractionary shocks small effect of expansionary shocks State dependence: In tight labor market, expansionary shock only generates inflation

39 Conclusion New method to identify the non-linear effects of shocks strong non-linearities in the effects of monetary policy shocks Can apply method to non-linear effects fiscal policy, credit supply shocks, etc.. Not limited to just-identified models, can generalize to under-identified models (sign-restrictions)

40 Theories behind Non-linearities Asymmetry Credit constraints Binding vs. non-binding Downward wage (price) rigidity State dependence Less inflationary pressure in periods of slack -> CB has more room to stimulate Y before P increases

41 Literature 1. Get non-linear effects from single equation regression (distributed lags model or Jorda s Local Projection) Asymmetry: Cover (1992), DeLong and Summers (1988), Morgan (1993). State-dependence: Thwaites and Tenreyro (213), Santoro et al. (214)

42 Literature 1. Get non-linear effects from single equation regression (distributed lags model or Jorda s Local Projection) Asymmetry: Cover (1992), DeLong and Summers (1988), Morgan (1993). State-dependence: Thwaites and Tenreyro (213), Santoro et al. (214) Problem: non-parametric nature comes at effi ciency cost

43 Literature 1. Get non-linear effects from single equation regression (distributed lags model or Jorda s Local Projection) Asymmetry: Cover (1992), DeLong and Summers (1988), Morgan (1993). State-dependence: Thwaites and Tenreyro (213), Santoro et al. (214) Problem: non-parametric nature comes at effi ciency cost 2. Regime-switching models: allows the economy to respond differently depending on its state Beaudry and Koop (1993), Thoma (1994), Potter (1995), Kandil (1995), Koop, Pesaran and Potter (1996), Koop and Potter, (1998), Ravn and Sola (1996, 24), Weise (1999), Lo and Piger (22)

44 Literature 1. Get non-linear effects from single equation regression (distributed lags model or Jorda s Local Projection) Asymmetry: Cover (1992), DeLong and Summers (1988), Morgan (1993). State-dependence: Thwaites and Tenreyro (213), Santoro et al. (214) Problem: non-parametric nature comes at effi ciency cost 2. Regime-switching models: allows the economy to respond differently depending on its state Beaudry and Koop (1993), Thoma (1994), Potter (1995), Kandil (1995), Koop, Pesaran and Potter (1996), Koop and Potter, (1998), Ravn and Sola (1996, 24), Weise (1999), Lo and Piger (22) Problem: more effi cient but restrictive and not flexible -> switching variable is not the current shock, but a function of past shocks -> only a small number of states

45 Literature 1. Get non-linear effects from single equation regression (distributed lags model or Jorda s Local Projection) Asymmetry: Cover (1992), DeLong and Summers (1988), Morgan (1993). State-dependence: Thwaites and Tenreyro (213), Santoro et al. (214) Problem: non-parametric nature comes at effi ciency cost 2. Regime-switching models: allows the economy to respond differently depending on its state Beaudry and Koop (1993), Thoma (1994), Potter (1995), Kandil (1995), Koop, Pesaran and Potter (1996), Koop and Potter, (1998), Ravn and Sola (1996, 24), Weise (1999), Lo and Piger (22) Problem: more effi cient but restrictive and not flexible -> switching variable is not the current shock, but a function of past shocks -> only a small number of states 3. Angrist et al. (213) semi-parametric method based on propensity score weighting

46 Approximating IRFs with mixtures of Gaussians (2) Theorem Let f be a continuous function over [, [ that satisfies f (x)2 dx <. Then, there exists a function f N defined by f N (x) = N a n e n=1 x bn ( ) cn 2 with a n, b n, c n R for n N, such that f N f on every interval of R.

47 Monte Carlo simulations Illustrate performance of GMA in finite sample with 3 simulations: 1. Linear 2. Asymmetry 3. Asymmetry and state dependence

48 Lessons from Monte-Carlo simulations GMA performs as well as VAR in linear models Successfully detects asymmetry or state-dependence Delivers good estimates of magnitude of non-linearities

49 Simulations - Estimate a VAR model - Invert VAR to get ˆΨ k - Modify ˆΨ k to introduce non-linearity (asym, state dep) - Simulate data using Normal shocks

50 Simulations - Estimate a VAR model - Invert VAR to get ˆΨ k - Modify ˆΨ k to introduce non-linearity (asym, state dep) - Simulate data using Normal shocks In all sims, GMA is misspecified: - Conservative - Consistent with "GMA meant to approximate true DGP, and yet can recover non-linearities"

51 DGP #1: Linear case Use trivariate VAR with (U,Π,r) over Simulate data with T = 2 over 5 replications Compare MSE, avg length and coverage rate of GMA vs. VAR

52 Summary statistics of simulation #1 Table 1: Summary statistics for Monte Carlo simulation with linear model U FFR VAR GMA VAR GMA VAR GMA MSE Avg length (at peak effect) Coverage rate (at peak effect) Note: Summary statistics over 5 Monte Carlo replications. MSE is the mean squared error of the estimated impulse response function over horizons 1 to 25. Avg length is the average distance between the lower (2.5%) and upper (97.5%) confidence bands at the time of peak effect of the monetary shock. The coverage rate is the frequency with which the true value lays within 95 percent of the posterior distribution. The VAR estimates and confidence bands are obtained from a Bayesian VAR with Normal Whishart priors. U, Π and FFR denote respectively unemployment, inflation and the fed funds rate.

53 Interest rate Inflation Output DGP #2: Asymmetry (lny,dlnp,r) in response to monetary shock Linear Contractionary shock Expansionary shock

54 Summary statistics of simulation #2 Table 2: Summary statistics for Monte Carlo exercise with asymmetric model a + a y dp ffr Mean (true value).82 ( 1.).5 (.6).3 (.) Std dev Frequency of rejection of zero coefficient Coverage rate Note: Summary statistics over 5 Monte Carlo replications. For each coefficient of interest, "Frequency of rejection of zero coefficient" is the frequency that lies outside 9 percent of the posterior distribution, and "Coverage rate" is the frequency with which the true value lies within 9 percent of the posterior distribution. y, dp and ffr denote respectively output, inflation and the fed funds rate.

55 DGP #3: Asymmetry and state dep. idem as DGP #1 but in addition γ + y > with z t the US unemployment rate A positive monetary shock has a larger effect on output in recessions

56 Summary statistics of simulation #3 Table 3: Summary statistics for Monte Carlo exercise with asymmetry and state dependence γ + γ γ + γ α + α y dp y dp y dp y dp Mean (true value).96 (1.).2 (.).71 (1.). (.).21 (.). (.).78 ( 1.).48 (.6) Std dev Frequency of rejection of zero coefficient Coverage rate Note: Summary statistics over 5 Monte Carlo replications. For each coefficient of interest, "Frequency of rejection of zero coefficient" is the frequency that lies outside 9 percent of the posterior distribution, and "Coverage rate" is the frequency with which the true value lies within 9 percent of the posterior distribution. y and dp denote respectively output, inflation and the fed funds rate.

57 Difference in IRFs Testing for asymmetry Posterior distribution of difference between IRFs to positive and negative shocks.3 Unemployment.1 Inflation 1.2 Fed funds rate

58 Expansionaryshock Contractionary shock GMA(1) to GMA(3) Unemployment Price level Fed funds rate ( ) Unemployment ( ) Price level ( ) Fed funds rate GMA 1 GMA 2 GMA

59 Using Romer and Romer shocks Using Romer and Romer (24) shocks ε MP over Consider Y = (ε MP, U, Π, r) with recursive ordering

60 Expansionary shock Contractionary shock Using Romer and Romer shocks.3.2 Unemployment.5 Price level Interest rate ( ) Unemployment.5 ( ) Price level ( ) Interest rate

61 Expansionary shock Contractionary shock Peak effect and state of business cycle Output growth Peak effect on Unemployment.5 Peak effect on Inflation.2.1 VAR GMA VAR GMA Output growth Output growth

62 Expansionary shock Contractionary shock Peak effect and state of business cycle detrended UR Peak effect on Unemployment.1 Peak effect on Inflation VAR GMA VAR GMA (detrended) UR (detrended) UR

63 Expansionary shock Contractionary shock Peak effect on FFR and state of business cycle 1 Cumulative change in FFR VAR Output gap (in %)