原文链接：http://tecdat.cn/?p=9384
原文出处：拓端数据部落公众号
最近我们被客户要求撰写关于VAR模型的研究报告，包括一些图形和统计输出。
目录
脉冲响应分析是采用[向量自回归模型的]()计量经济学分析中的重要一步。它们的主要目的是描述模型变量对一个或多个变量的冲击的演化。因此使它们成为评估经济时非常有用的工具。这篇文章介绍了VAR文献中常用的脉冲响应函数的概念和解释。
模型与数据
为了说明脉冲响应函数的概念，使用了Lütkepohl（2007）的示例。可以从教科书的[网站上]()下载所需的数据集。它包含从1960年1季度到1982年4季度按季度和季节性调整的时间序列，这些序列是西德的固定投资，可支配收入和数十亿德国马克的消费支出。
# 下载数据
data <- read.table("e1.dat", skip = 6, header = TRUE)
# 仅使用前76个观测值，因此有73个观测值
# 取一阶差分后，留给估计的VAR(2)模型。
data <- data[1:76, ]
# 转换为时间序列对象
data <- ts(data, start = c(1960, 1), frequency = 4)
# 取对数和差值
data <- diff(log(data))
# 绘图数据
plot(data, main = "Dataset E1 from Lütkepohl (2007)")
此数据用于估计具有常数项的VAR（2）模型。
估算值
可以使用vars软件包估算VAR模型：
# 查看摘要统计信息
summary(model)代码的结果应与Lütkepohl（2007）的3.2.3节中的结果相同。
预测误差脉冲响应
由于VAR模型中的所有变量都相互依赖，因此单独的系数估计仅提供有关反应的有限信息。为了更好地了解模型的动态行为，使用了脉冲响应（IR）。线性VAR模型的每个脉冲响应函数的出发点都是其移动平均值（MA）表示，这也是预测误差脉冲响应（FEIR）函数。
在R 中，程序包可用于获取预测误差脉冲响应。
识别问题
从上图可以看出，在第一期间FEIR为零。对于使用的数据集，估计为
## invest income cons
## invest 2.129629e-03 7.161667e-05 1.232404e-04
## income 7.161667e-05 1.373377e-04 6.145867e-05
## cons 1.232404e-04 6.145867e-05 8.920351e-05由于估计方差-协方差矩阵的非对角线元素不为零，因此我们可以假设VAR模型中的变量之间存在同期相关性。这由与Σ相对应的相关矩阵确认：
## invest income cons
## invest 1.0000000 0.1324242 0.2827548
## income 0.1324242 1.0000000 0.5552611
## cons 0.2827548 0.5552611 1.0000000但是，这些矩阵仅描述了误差之间的相关性，但不清楚因果关系的方向。识别这些因果关系是任何VAR分析的主要挑战之一。
正交脉冲响应
识别VAR模型的冲击的常用方法是使用正交脉冲响应（OIR）。基本思想是分解方差-协方差矩阵，使∑ = PP− 1，其中P是带有正对角线元素的下三角矩阵，通常通过Choleski分解获得。给定估计方差-协方差矩阵PP，可以通过以下方法获得分解
## invest income cons
## invest 0.046147903 0.000000000 0.000000000
## income 0.001551894 0.011615909 0.000000000
## cons 0.002670552 0.004934117 0.007597773从这个矩阵可以看出，收入冲击对消费具有同时性的影响，反之则不然。
在R 中，vars可以通过设置参数来使用包的功能来获得OIR：
plot(oir)
请注意，Choleski分解的输出是一个较低的三角矩阵，因此第一行中的变量永远不会对任何其他变量的同时冲击敏感，而系统中的最后一个变量将对所有其他变量的冲击敏感。因此，OIR的结果可能对变量的顺序很敏感，建议用不同的顺序估计上述VAR模型，以查看所产生的OIR受此影响的程度。
结构脉冲反应
在VAR模型的估计过程中，结构脉冲响应（SIR）已经考虑了识别问题。
广义脉冲响应
正交和结构响应都可以通过找到变量的正确顺序或通过识别估计的结构参数来约束。Koop等（1998）提出了一种不同类型的响应函数，即所谓的广义脉冲响应（GIR）。它们独立于变量顺序，因为它们将其他冲击的影响整合到响应之外。
对于难以识别结构关系的大型系统，GIR非常有用。
参考文献
Koop, G., Pesaran, M. H., Potter, S. M. (1996). Impulse response analysis in nonlinear multivariate models. Journal of Econometrics 74, 119-147. doi:10.1016/0304-4076(95)01753-4
最受欢迎的见解
1.在python中使用lstm和pytorch进行时间序列预测
2.python中利用长短期记忆模型lstm进行时间序列预测分析