拓端tecdat|R语言VAR模型的不同类型的脉冲响应分析

模型与数据

为了说明脉冲响应函数的概念，使用了Lütkepohl（2007）的示例。可以从教科书的网站上下载所需的数据集。它包含从1960年1季度到1982年4季度按季度和季节性调整的时间序列，这些序列是西德的固定投资，可支配收入和数十亿德国马克的消费支出。

# Download data
data <- read.table("http://www.jmulti.de/download/datasets/e1.dat", skip = 6, header = TRUE)

# Only use the first 76 observations so that there are 73 observations
# left for the estimated VAR(2) model after taking first differences.
data <- data[1:76, ]

# Convert to time series object
data <- ts(data, start = c(1960, 1), frequency = 4)

# Take logs and differences
data <- diff(log(data))

# Plot data
plot(data,  main = "Dataset E1 from Lütkepohl (2007)")

​

此数据用于估计具有常数项的VAR（2）模型。

估算值

可以使用vars软件包估算VAR模型：

# Look at summary statistics
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Σ=PP−1，其中PP是带有正对角线元素的下三角矩阵，通常通过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