VAR impulse response

From lag equations to dynamic responses

A VAR$(2)$ models every variable using two lags of the full vector:

Its moving-average form writes $z_t=\sum_{h=0}^{\infty}\Psi_hu_{t-h}$. A structural impact matrix $P$ maps orthogonal shocks into reduced-form innovations (Lutkepohl, 2005).

A two-variable propagation table

Suppose a one-unit shock to the first variable produces these responses.

The second variable reacts on impact because the contemporaneous identification matrix is not diagonal.

Assumptions and implementation

The fitted VAR should be stable. The shock needs an identification rule. Cholesky identification is recursive, so variable order changes $P$ and the responses.

Packages differ in whether horizon zero is included and how cumulative responses are labeled. Confidence intervals can use asymptotic formulas or simulation.

A coefficient-recursion diagram connects the VAR parameters to aligned response paths for each shock.

Interpretation and reporting

Each response traces the effect of one identified shock through future horizons. In the table, $z_2$ moves by $0.40$ on impact, peaks at $0.50$, and then declines. Those numbers depend on shock normalization and the matrix $P$.

Report the shocked variable, responding variable, shock size, identification rule, horizon indexing, and interval method. State whether horizon zero appears. For cumulative responses, name the transformation before interpreting long-horizon levels.

A Cholesky response changes when variable order changes because the ordering imposes contemporaneous exclusions. Stability of the reduced-form VAR does not validate that structural ordering. Long-horizon estimates can also accumulate substantial sampling uncertainty.

Reproducible implementation

Estimate every VAR equation on the same aligned sample and declared lag order. Check companion-matrix roots before computing responses. Build the moving-average matrices recursively from $A_1$ and $A_2$.

Choose an impact matrix and verify that it reproduces the innovation covariance through $PP'=\widehat\Sigma_u$. Multiply each $\Psi_h$ by $P$, then label rows and columns with response and shock names.

Compare the first horizons with a direct recursion calculation. Generate intervals with the stated asymptotic or resampling method. Reorder variables or use another defensible identification as a sensitivity check. Store the exact normalization and random seed with simulated intervals (Lutkepohl, 2005).

Further reading

Sims develops vector autoregressions as multivariate dynamic models (Sims, 1980). Lütkepohl covers impulse-response computation and inference (Lutkepohl, 2005).