模型与数据

估算值

预测误差脉冲响应

识别问题

正交脉冲响应

结构脉冲反应

广义脉冲响应

参考文献

最近我们被客户要求撰写关于脉冲响应分析的研究报告，包括一些图形和统计输出。脉冲响应分析是采用向量自回归模型的计量经济学分析中的重要一步。它们的主要目的是描述模型变量对一个或多个变量的冲击的演化。因此使它们成为评估经济时非常有用的工具。这篇文章介绍了VAR文献中常用的脉冲响应函数的概念和解释。

视频：向量自回归VAR数学原理及R软件经济数据脉冲响应分析实例

【视频】向量自回归VAR数学原理及R语言软件经济数据脉冲响应分析实例

，时长12:01

模型与数据

为了说明脉冲响应函数的概念，使用了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

最受欢迎的见解

R语言VAR模型的不同类型的脉冲响应分析相关推荐

R语言VAR模型方差分解代码
好的,您可以使用以下代码来进行 VAR 模型的方差分解: library(vars)# 建立 VAR 模型 model <- VAR(data, p = 2)# 进行方差分解 irf <- ...
R语言 VAR 模型的方差分解
可以使用以下代码来进行 VAR 模型的方差分解: library(vars)# 建立 VAR 模型 model <- VAR(data, p = 2)# 进行方差分解 irf <- irf ...
R语言VaR市场风险计算方法与回测、用LOGIT逻辑回归、PROBIT模型信用风险与分类模型...
全文链接:http://tecdat.cn/?p=27530 市场风险指的是由金融市场中资产的价格下跌或价格波动增加所导致的可能损失. 相关视频市场风险包含两种类型:相对风险和绝对风险.绝对风险关 ...
R语言︱机器学习模型评价指标+（转）模型出错的四大原因及如何纠错
笔者寄语:机器学习中交叉验证的方式是主要的模型评价方法,交叉验证中用到了哪些指标呢? 交叉验证将数据分为训练数据集.测试数据集,然后通过训练数据集进行训练,通过测试数据集进行测试,验证集进行验证. 模 ...
R语言分类模型：逻辑回归模型LR、决策树DT、推理决策树CDT、随机森林RF、支持向量机SVM、Rattle可视化界面数据挖掘、分类模型评估指标（准确度、敏感度、特异度、PPV、NPV）
R语言分类模型:逻辑回归模型LR.决策树DT.推理决策树CDT.随机森林RF.支持向量机SVM.Rattle可视化界面数据挖掘.分类模型评估指标(准确度.敏感度.特异度.PPV.NPV) 目录
R语言回归模型构建、回归模型基本假设（正态性、线性、独立性、方差齐性）、回归模型诊断、car包诊断回归模型、特殊观察样本分析、数据变换、模型比较、特征筛选、交叉验证、预测变量相对重要度
R语言回归模型构建.回归模型基本假设(正态性.线性.独立性.方差齐性).回归模型诊断.car包诊断回归模型.特殊观察样本分析.数据变换.模型比较.特征筛选.交叉验证.预测变量相对重要度目录
R语言KNN模型数据分类实战
R语言KNN模型数据分类实战目录 R语言KNN模型数据分类实战 #KNN基本原理 #案例分析 #数据集说明
R语言xgboost模型构建:基于prima糖尿病数据集
R语言xgboost模型构建:基于prima糖尿病数据集目录 R语言xgboost模型构建:基于prima糖尿病数据集
R语言回归模型协方差分析（Analysis of Covariance）
R语言回归模型协方差分析(Analysis of Covariance) 目录 R语言回归模型协方差分析(Analysis of Covariance) 输入数据 ANCOVA分析

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

视频：向量自回归VAR数学原理及R软件经济数据脉冲响应分析实例

模型与数据

估算值

预测误差脉冲响应

识别问题

正交脉冲响应

结构脉冲反应

广义脉冲响应

参考文献

R语言VAR模型的不同类型的脉冲响应分析相关推荐

最新文章

热门文章