# NOT RUN {
# set dimensions
p <- 3
n <- 4
T <- 10
# set model parameters
SigmaE <- diag(p)/4
A <- matrix(c(-0.1, -0.3, 0.6, 0.5, -0.4, 0, 0.3, -0.5, -0.2),
byrow=TRUE, ncol=3)
# generate data
Y <- dataVAR1(n, T, A, SigmaE)
# fit VAR(1) model
VAR1hat <- ridgeVAR1(Y, 1, 1)
# impulse response analysis
impulseResponseVAR1(VAR1hat$A, 10)
# }
Run the code above in your browser using DataLab