导入vars包

library('vars')

导入数据，这里假设有两个变量，y和x

data <- read.csv('data.csv', header = TRUE)

将数据转化为时间序列对象

ts_data <- ts(data[, c('y', 'x')], start = c(2000, 1), frequency = 12)

估计VAR模型

var_model <- VAR(ts_data, p = 2)

估计SVAR模型

svar_model <- SVAR(var_model, Amat = matrix(c(0.5, 0.5, 0, 1), ncol = 2))

打印SVAR模型结果

summary(svar_model)

进行SVAR冲击响应分析

irf <- irf(svar_model, impulse = 'y', response = c('y', 'x'), boot = TRUE, runs = 1000)

绘制SVAR冲击响应图

plot(irf)