# set dimensions
p <- 3; n <- 4; T <- 10
# set model parameters
SigmaE <- diag(p)/4
A <- createA(p, "chain")
# generate data
Y <- dataVAR1(n, T, A, SigmaE)
# fit VAR(1) model
VAR1hat <- ridgeVAR1(Y, 1, 1)
# impulse response analysis
impulseResponseVAR1(VAR1hat$A, 10, figure=TRUE)Run the code above in your browser using DataLab