Extract monetary policy shocks via Kalman filter

Extracts structural shocks from reduced-form residuals using the Kalman filter formula. Given the impact matrix Psi and the covariance matrices of residuals on event and control days, the function recovers the latent shocks for each observation. Optionally rescales shocks to unit variance.

Usage

kfpredict(Sig, SigR, Psi, et, tol = sqrt(.Machine$double.eps), scale = TRUE)

Arguments

Sig

Numeric matrix (N x N). Covariance matrix of reduced-form residuals on policy event days (e.g. et from hetiv()).

SigR

Numeric matrix (N x N). Covariance matrix of reduced-form residuals on control (non-event) days. Used to back out the implied shock variances when scale = TRUE.

Psi

Numeric matrix (N x E). Impact matrix, i.e. the contemporaneous responses of all N variables to the E structural shocks (e.g. Psi from hetiv()).

et

Numeric matrix or data frame (T x N). Reduced-form residuals on event days (e.g. et from hetiv()). Rows correspond to time periods, columns to variables.

tol

Numeric scalar. Tolerance for the generalized inverse (MASS::ginv()). Clipped from below at sqrt(.Machine$double.eps).

scale

Logical. If TRUE (default), shocks are rescaled to unit variance using the implied shock variances recovered from Sig and SigR. If FALSE, the raw Kalman filter projection is returned.

Value

A numeric matrix (T x E) of extracted structural shocks, with the same number of rows as et and one column per shock dimension.

References

Burri, M. and Kaufmann, D. (2026). Measuring monetary policy shocks. IRENE Working Paper 24-03, IRENE Institute of Economic Research, University of Neuchâtel.