-
2021-04-22 18:12:23
压缩包 : c2c6f64303164399995e96c534eb64.rar 列表
贝叶斯向量自回归MATLAB代码/var_bvar/ar.m
贝叶斯向量自回归MATLAB代码/var_bvar/arf.m
贝叶斯向量自回归MATLAB代码/var_bvar/arf_d.m
贝叶斯向量自回归MATLAB代码/var_bvar/becm.m
贝叶斯向量自回归MATLAB代码/var_bvar/becmf.m
贝叶斯向量自回归MATLAB代码/var_bvar/becmf_d.m
贝叶斯向量自回归MATLAB代码/var_bvar/becmf_g.m
贝叶斯向量自回归MATLAB代码/var_bvar/becmf_gd.m
贝叶斯向量自回归MATLAB代码/var_bvar/becm_d.m
贝叶斯向量自回归MATLAB代码/var_bvar/becm_g.m
贝叶斯向量自回归MATLAB代码/var_bvar/becm_gd.m
贝叶斯向量自回归MATLAB代码/var_bvar/bvar.m
贝叶斯向量自回归MATLAB代码/var_bvar/bvarf.m
贝叶斯向量自回归MATLAB代码/var_bvar/bvarf_d.m
贝叶斯向量自回归MATLAB代码/var_bvar/bvarf_g.m
贝叶斯向量自回归MATLAB代码/var_bvar/bvarf_gd.m
贝叶斯向量自回归MATLAB代码/var_bvar/bvar_d.m
贝叶斯向量自回归MATLAB代码/var_bvar/bvar_g.m
贝叶斯向量自回归MATLAB代码/var_bvar/bvar_gd.m
贝叶斯向量自回归MATLAB代码/var_bvar/commutation.m
贝叶斯向量自回归MATLAB代码/var_bvar/contents.html
贝叶斯向量自回归MATLAB代码/var_bvar/contents.m
贝叶斯向量自回归MATLAB代码/var_bvar/ecm.m
贝叶斯向量自回归MATLAB代码/var_bvar/ecmf.m
贝叶斯向量自回归MATLAB代码/var_bvar/ecmf_d.m
贝叶斯向量自回归MATLAB代码/var_bvar/ecm_d.m
贝叶斯向量自回归MATLAB代码/var_bvar/irf.m
贝叶斯向量自回归MATLAB代码/var_bvar/irf_d.m
贝叶斯向量自回归MATLAB代码/var_bvar/irf_d2.m
贝叶斯向量自回归MATLAB代码/var_bvar/lrratio.m
贝叶斯向量自回归MATLAB代码/var_bvar/lrratio_d.m
贝叶斯向量自回归MATLAB代码/var_bvar/make_html.m
贝叶斯向量自回归MATLAB代码/var_bvar/pftest.m
贝叶斯向量自回归MATLAB代码/var_bvar/pftest_d.m
贝叶斯向量自回归MATLAB代码/var_bvar/pgranger.m
贝叶斯向量自回归MATLAB代码/var_bvar/plt_var.m
贝叶斯向量自回归MATLAB代码/var_bvar/plt_varg.m
贝叶斯向量自回归MATLAB代码/var_bvar/prt_var.m
贝叶斯向量自回归MATLAB代码/var_bvar/prt_varg.m
贝叶斯向量自回归MATLAB代码/var_bvar/recm.m
贝叶斯向量自回归MATLAB代码/var_bvar/recm.out
贝叶斯向量自回归MATLAB代码/var_bvar/recmf.m
贝叶斯向量自回归MATLAB代码/var_bvar/recmf_d.m
贝叶斯向量自回归MATLAB代码/var_bvar/recmf_g.m
贝叶斯向量自回归MATLAB代码/var_bvar/recmf_gd.m
贝叶斯向量自回归MATLAB代码/var_bvar/recm_d.m
贝叶斯向量自回归MATLAB代码/var_bvar/recm_g.m
贝叶斯向量自回归MATLAB代码/var_bvar/recm_gd.m
贝叶斯向量自回归MATLAB代码/var_bvar/rvar.m
贝叶斯向量自回归MATLAB代码/var_bvar/rvarb.m
贝叶斯向量自回归MATLAB代码/var_bvar/rvarf.m
贝叶斯向量自回归MATLAB代码/var_bvar/rvarf_d.m
贝叶斯向量自回归MATLAB代码/var_bvar/rvarf_g.m
贝叶斯向量自回归MATLAB代码/var_bvar/rvarf_gd.m
贝叶斯向量自回归MATLAB代码/var_bvar/rvar_d.m
贝叶斯向量自回归MATLAB代码/var_bvar/rvar_g.m
贝叶斯向量自回归MATLAB代码/var_bvar/rvar_gd.m
贝叶斯向量自回归MATLAB代码/var_bvar/scstd.m
贝叶斯向量自回归MATLAB代码/var_bvar/svar.m
贝叶斯向量自回归MATLAB代码/var_bvar/svar.pdf
贝叶斯向量自回归MATLAB代码/var_bvar/svar.tex
贝叶斯向量自回归MATLAB代码/var_bvar/svar_d.m
贝叶斯向量自回归MATLAB代码/var_bvar/test.dat
贝叶斯向量自回归MATLAB代码/var_bvar/theilbf.m
贝叶斯向量自回归MATLAB代码/var_bvar/theilbv.m
贝叶斯向量自回归MATLAB代码/var_bvar/theil_g.m
贝叶斯向量自回归MATLAB代码/var_bvar/vare.m
贝叶斯向量自回归MATLAB代码/var_bvar/vare_d.m
贝叶斯向量自回归MATLAB代码/var_bvar/varf.m
贝叶斯向量自回归MATLAB代码/var_bvar/varf_d.m
贝叶斯向量自回归MATLAB代码/var_bvar/var_resid.m
贝叶斯向量自回归MATLAB代码/var_bvar/WS_FTP.LOG
贝叶斯向量自回归MATLAB代码/var_bvar
贝叶斯向量自回归MATLAB代码
更多相关内容 -
回归模型脉络(线性回归、自回归、向量自回归)
2020-08-20 17:16:13回归模型脉络(线性回归、自回归、向量自回归) -
贝叶斯向量自回归MATLAB代码,matlab贝叶斯算法,matlab
2021-09-10 17:10:32使用matlab实现贝叶斯向量自回归模型,可用于经济学中的预测 -
Python-VAR:Python中的向量自回归模型
2021-05-15 02:27:59Python变量 Python中的向量自回归模型 -
向量自回归模型建模步骤梳理(VAR、VMA、VARMA模型)
2020-08-20 17:12:37向量自回归模型建模步骤梳理(VAR、VMA、VARMA模型) -
向量自回归模型
2018-12-01 10:01:08向量自回归模型的原理和分析过程,包括其运用范围,原理讲解,应用实例 -
向量自回归模型(VAR)-Eviews实现.pptx
2021-09-22 09:08:04向量自回归模型(VAR)-Eviews实现.pptx -
向量自回归模型的理论方法及应用实例
2018-11-07 19:27:21详细讲解向量自回归的历史发展作用和实战应用,有着非常丰富的逻辑脉络和实战案例 -
matlab自相关代码-awesome-var:向量自回归资源的精选列表
2021-05-23 15:00:31向量自回归资源的精选列表。 目录: 的MATLAB 工具箱 :向量自回归模型 :收集Matlab例程以执行VAR分析(Ambrogio Cesa-Bianchi) :经验宏工具箱(F. Ferroni和F. Canova) :宏观经济建模工具箱 :贝叶斯估计,... -
向量自回归
2012-09-24 12:11:31这一章所讲述的向量自回归模型(Vector Auto regression, VAR)以及向量误差修正模型(Vector Error Correction, VEC)的估计与分析。同时也给出一些检验几个非稳定变量之间协整关系的工具 -
基于贝叶斯向量自回归的中国国债收益率预测.pdf
2020-07-01 15:21:17本文运用基于独立 Minnesota-Wishart 共轭先验分布的贝叶斯向量自回归模型( BVAR) ,并通过 Gibbs 抽样的马尔科夫链蒙特卡洛模拟方法预测中国银行间国债的收益率。此外,按照固定窗口的滚动预测规则, 采用统计性... -
大贝叶斯向量自回归-研究论文
2021-06-10 10:37:34贝叶斯向量自回归广泛用于宏观经济预测和结构分析。 然而,直到最近,由于参数扩散问题和计算限制,大多数实证工作只考虑了具有少量变量的小系统。 我们首先回顾了各种可用于解决大型贝叶斯 VAR 中参数增殖问题的... -
时变网络向量自回归模型-研究论文
2021-06-10 04:45:13在本文中,我们研究了时变网络向量自回归 (TV-NVAR) 模型。 构成一个超高维向量,节点被假定为 (a) 其先前值、(b) 其邻居先前值、(c) 外生变量和 (d) 独立噪声的线性组合。 我们采用非参数方法来获得时变和轮廓估计... -
【数学建模】基于matlab时变参数随机波动率向量自回归模型(TVP-VAR)【含Matlab源码 037期】.zip
2021-09-23 21:28:06【数学建模】基于matlab时变参数随机波动率向量自回归模型(TVP-VAR)【含Matlab源码 037期】.zip -
结构向量自回归模型操作步骤.pptx
2021-10-07 18:03:16结构向量自回归模型操作步骤.pptx -
将 VAR 参数转化为 MA 参数:将向量自回归模型参数转化为移动平均模型参数的函数。-matlab开发
2021-05-30 18:48:26此例程将向量自回归 (VAR) 的参数估计映射到相应移动平均 (MA) 模型的参数估计中。 此函数的输出可用于构建 VAR 模型的结构脉冲响应函数。 -
第十一章向量自回归(VAR)模型和向量误差修正(VEC)模型.pptx
2021-10-05 14:20:06第十一章向量自回归(VAR)模型和向量误差修正(VEC)模型.pptx -
Large Inverse Cholesky:计算结构化向量自回归模型的系数矩阵-matlab开发
2021-05-30 12:50:50此函数计算结构化自回归模型的系数矩阵,如下所示: L*x(n) = t + sum_{i=1}^KR(:,:,i)*x(ni) + w(n) n = [1,N];... 阿拉文德·克里希那穆提关于结构化向量自回归模型的系数矩阵计算,arXiv:1309.6290。 -
向量自回归模型VAR
2013-11-20 10:33:09向量自回归是一种很好的模型研究,要好好把握好,不然就很难做回归了 -
结构向量自回归(SVAR)模型操作步骤.pptx
2021-09-21 06:14:56结构向量自回归(SVAR)模型操作步骤.pptx -
数学建模之向量自回归模型
2020-07-16 13:16:242.向量自回归模型建模流程 1、变量平稳性检验——ADF检验: 原假设H0:序列不平稳 Matlab代码:(现阶段不知道啥意思) [h1,p1,adf,ljz]=adftest(x) h1 = 0 p1 = 0.9982 adf = 2.8923 ljz = -1.9489 2、如果非平稳...EverydayOneCat
非常快乐/十分快乐.gif
知识点
1.笔记
2.向量自回归模型建模流程
1、变量平稳性检验——ADF检验:
原假设H0:序列不平稳
Matlab代码:(现阶段不知道啥意思)
[h1,p1,adf,ljz]=adftest(x) h1 = 0 p1 = 0.9982 adf = 2.8923 ljz = -1.9489
2、如果非平稳,我们不忙着做差分,我们先来看看协整关系检验:简单的说就是变量自身非平稳,但其线性组合却是平稳的,这种组合反映了变量之间长期稳定的比例关系。
3、VAR§模型最优滞后阶数的确定
一般采用AIC、SBC等信息准则采用多种方法进行综合判断。
判别标准:在p的取值范围内,选择使AIC、SBS最小的p值为最优
从图中我们看出p=2为最优,但是我们做预测,一般选择简单的一阶,反正也相差不是很大。
4、估计VAR模型的参数
5、VAR模型的检验
3.案例分析
2010年上海世博会是首次在中国举办的世界博览会。从1851年伦敦开始,世博会正日益成为各国人民交流历史文化、展示科技成果、体现合作精神、展望未来发展的重要舞台。以1990-2008年的可能受世博会影响的我国8个指标数据定量分析2010年上海世博会的影响力。
SAS代码:
data ex; input x1 x2 x3 x4 x5 x6 x7 x8@@; y1=log(x1); y2=log(x2); y3=log(x3); y4=log(x4); y5=log(x5); y6=log(x6); y7=log(x7); y8=log(x8);/*对x1∽x8取对数,以消除趋势的变化和多重共线性*/ cards; 18667.82 1510.2 6767.2 174.73 1.8 17041 699.75 22.01 21781.5 1700.6 8542.5 263.3 1.8 17465 812.96 22.79 26923.48 2026.6 10317.8 355.6 3.2 17861 938.29 23.43 35333.92 2577.4 12093.1 438.6 3.2 18262 1051.5 24.58 48197.86 3496.2 13868.4 500.3 3.2 18653 1357.1 25.72 60793.73 4283 15643.7 588.67 2 19040 1702.4 26.86 71176.59 4838.9 17567.2 674.43 1.7 19922 2024.2 27.89 78973.04 5160.3 19194.2 742.8 1.4 20781 2208.2 28.29 84402.28 5425.1 22491.4 710.77 1.4 21616 2336.7 28.42 89677.05 5854 23732 843.23 1 22412 2475.2 28.32 99214.55 6280 26221.8 1016.04 1 23151 2694.7 28.44 109655.17 6859.6 30001.2 1122.64 1 23940 2945.7 28.61 120332.69 7702.8 35488.8 1343.95 0.7 24780 3184.9 28.72 135822.76 8472.2 45811.7 1140.29 0.7 25639 3558 29.32 159878.34 9421.6 59028.2 1693.25 0.7 26476 4163 30.72 183217.4 10493 75095.1 2025.51 0.7 27331 5153 31.96 211923.5 11759.5 93368.7 2221.03 0.7 28310 5828 33.35 257305.6 13785.8 117464.5 2610.97 0.7 29350 6769 34.47 300670 14306.9 148738.3 2432.53 0.7 30210 7983 35.63 ; /*VAR model unit root test原始序列*/ proc arima data=ex; identify var=y1 stationarity=(adf=3);/*平稳性检验-ADF检验,adf=3表示检验模型是滞后3阶*/ identify var=y2 stationarity=(adf=3); identify var=y3 stationarity=(adf=3); identify var=y4 stationarity=(adf=3); identify var=y5 stationarity=(adf=3); identify var=y6 stationarity=(adf=3); identify var=y7 stationarity=(adf=3); identify var=y8 stationarity=(adf=3); run; /*VAR model unit root test一阶差分*/ proc arima data=ex; identify var=y1(1) stationarity=(adf=3); identify var=y2(1) stationarity=(adf=3); identify var=y3(1) stationarity=(adf=3); identify var=y4(1) stationarity=(adf=3); identify var=y5(1) stationarity=(adf=3); identify var=y6(1) stationarity=(adf=3); identify var=y7(1) stationarity=(adf=3); identify var=y8(1) stationarity=(adf=3); run; /*VAR model unit root test二阶差分*/ proc arima data=ex; identify var=y1(1,1) stationarity=(adf=3); identify var=y2(1,1) stationarity=(adf=3); identify var=y3(1,1) stationarity=(adf=3); identify var=y4(1,1) stationarity=(adf=3); identify var=y5(1,1) stationarity=(adf=3); identify var=y6(1,1) stationarity=(adf=3); identify var=y7(1,1) stationarity=(adf=3); identify var=y8(1,1) stationarity=(adf=3); proc varmax data=ex; model y1-y8/minic=(type=aic p=2 q=0) /*根据AIC信息最小准则,在p=0,1时确定p的最优值。并在最优p值的条件下,估计VAR模型*/ print=(diagnose roots)/*显示模型的特征根其他检验,判断模型的稳定性和拟合效果*/ cointtest=(johansen=(type=trace));/*协整检验*/ output lead=5;/*向前5步预测*/ run;
VAR(1)模型参数估计结果:
稳定性检验:
单方程显著性检验:
预测:
作业
请结合下列数据,建立向量自回归模型,并预测2008-2010年GDP的值。
年份 GDP 亿元 进出口总额 城市基础设施投资 价格指数 就业率 城市人口比率 单位 亿美元 亿美元 亿元 1 % % 1990 781.66 74.3 47.20 108.20 98.5 67.4 1991 893.77 80.4 61.40 110.50 98.6 67.6 1992 1114.32 97.6 84.40 121.60 98.3 67.9 1993 1519.23 127.3 167.90 146.10 97.6 69.0 1994 1990.86 158.7 238.20 181.00 97.2 70.1 1995 2499.43 190.3 273.80 214.90 97.3 70.8 1996 2957.55 222.6 378.80 234.60 97.3 71.5 1997 3438.79 247.6 412.90 241.20 97.2 72.2 1998 3801.09 313.4 531.40 241.20 97.1 73.0 1999 4188.73 386.04 501.40 244.80 96.9 73.8 2000 4771.17 547.10 449.90 250.90 96.5 74.6 2001 5210.12 608.98 510.80 250.90 95.7 75.3 2002 5741.03 726.64 583.50 252.20 95.2 76.4 2003 6694.23 1123.97 604.62 252.5 95.1 77.6 2004 8072.83 1600.26 672.58 257.8 95.5 81.2 2005 9247.66 1863.65 885.74 260.3 95.6 84.5 2006 10572.2 2274.89 1125.54 263.4 95.6 85.8 2007 12494.0 2829.73 1466.33 271.7 95.7 86.8 1.平稳性检验
在平稳性检验中我们使用的是ADF 单位根检验法对6 组变量进行滞后3 期的平稳性检验。
SAS代码:
data ex; input x1-x6@@; y1=log(x1); y2=log(x2); y3=log(x3); y4=log(x4);y5=log(x5); y6=log(x6); cards; 781.66 74.3 47.20 108.20 98.5 67.4 893.77 80.4 61.40 110.50 98.6 67.6 1114.32 97.6 84.40 121.60 98.3 67.9 1519.23 127.3 167.90 146.10 97.6 69.0 1990.86 158.7 238.20 181.00 97.2 70.1 2499.43 190.3 273.80 214.90 97.3 70.8 2957.55 222.6 378.80 234.60 97.3 71.5 3438.79 247.6 412.90 241.20 97.2 72.2 3801.09 313.4 531.40 241.20 97.1 73.0 4188.73 386.04 501.40 244.80 96.9 73.8 4771.17 547.10 449.90 250.90 96.5 74.6 5210.12 608.98 510.80 250.90 95.7 75.3 5741.03 726.64 583.50 252.20 95.2 76.4 6694.23 1123.97 604.62 252.5 95.1 77.6 8072.83 1600.26 672.58 257.8 95.5 81.2 9247.66 1863.65 885.74 260.3 95.6 84.5 10572.2 2274.89 1125.54 263.4 95.6 85.8 12494.0 2829.73 1466.33 271.7 95.7 86.8 ; /*VAR model unit root test原始序列*/ proc arima data=ex; identify var=y1 stationarity=(adf=3); identify var=y2 stationarity=(adf=3); identify var=y3 stationarity=(adf=3); identify var=y4 stationarity=(adf=3); identify var=y5 stationarity=(adf=3); identify var=y6 stationarity=(adf=3); run; proc varmax data=ex; model y1-y6/minic=(type=aic p=2 q=0) print=(diagnose roots) cointtest=(johansen=(type=trace)); output lead=3; run;
检验的原假设为序列不是平稳序列,备择假设为序列是平稳序列。由上表可知x1 、x2、x4 、x5 中均存在Pr<Tau 的概率<0.05,因此拒绝原假设,认为x1 、x2、x4 、x5是平稳序列。但x3 、x6 中Pr<Tau 的概率均>0.05,所以接受原假设,认为序列不是平稳序列。因此我们还需要对序列进行协整关系检验。
2.协整关系检验
由于自变量本身非平稳,因此需要检验其线性组合之间是否平稳,这种组合反映了变量之间长期稳定的比例关系。我们采用特征根迹检验法,得到的检验表如下:
当r=0 时, trace=193.3870> 93.92; 当r=1 时, trace=101.0363> 68.68; 当r=2 时,trace=46.7559<47.21,说明变量间存在协整关系,个数是2 个;同理当r=3、4、5,时trace均小于5% Critical Value,说明变量间都存在协整关系,最大个数为5 个。因此,我们可以对变量使用向量自回归模型。
3.确定VAR§模型最优滞后阶数
采用AIC 信息准则采用多种方法进行综合判断,当信息保留值最大,即AIC 取值最小时p的取值即为最优滞后阶数。得到的AIC 取值表如下:
由AIC 取值表可以看出,当p=1 是,AIC 最小,所以我们确定的模型为VAR(1)。
4.参数估计
我们使用最小二乘法对VAR(1)模型参数进行估计,根据得到的参数写出VAR(1)的数学模型如下:
5.模型检验
首先对模型进行稳定性检验,得到的检验表如下:
由上检验表可知,模型的所有特征根都在单位圆之内,所以模型是稳定的。接着对单方程进行检验,得到检验表如下:
由单方程检验表可知,所有方程Pr>F 的概率均小于0.05,说明8 个方程的显著性均较好,该模型可以用来预测后三年的GDP 数值。预测所得的结果如下表:
转化为x做出如下表格:
结语
让你们见识一下什么叫做说了半天啥都没说:
有一说一,这种事情我见的多了,我只想说懂得都懂,不懂的我也不多说了,说了你也不明白,不如不说,细细品吧,你也别来问我怎么回事,这里面利益牵扯太大了,说了对你我都没有好处,你就当不知道就行了,其余的我只能说这里水很深,牵扯到很多东西,详细情况你们很难找到的,网上大部分都删干净了。
-
机器学习笔记:向量自回归模型VAR
2021-11-03 15:39:061 自回归模型 给定多元时间序列数据,对于任意第t个时间间隔,有: 换一个角度看,可以看成是个input 为N维,output为N维的fully-connectedlayer 2 自回归模型 最优解 我们令 则自回归模型可以改写为...1 向量自回归模型
时间序列分析从单一时间序列 拓展到了多元时间序列,在任意第t 个时刻,观测样本从 一维变成了N维
给定多元时间序列数据
,对于任意第t个时间间隔,有:
换一个角度看
,可以看成是个input 为N维,output为N维的fully-connected layer
2 自回归模型 最优解
我们令
则自回归模型可以改写为:
将向量拼成矩阵,有:
其中
对此采用最小二乘法,可以求得系数矩阵A的最优解
其中 第一行<——>第二行的推导可见
参考资料
-
【视频】向量自回归VAR数学原理及R软件经济数据脉冲响应分析实例
2022-04-01 22:17:36向量自回归 (VAR) 是一种用于多变量时间序列分析的统计模型,尤其是在变量具有相互影响关系的时间序列中,本视频中我们介绍了向量自回归并在R软件中进行实现。 视频:向量自回归VAR数学原理及R软件经济数据脉冲...原文链接:http://tecdat.cn/?p=9368
原文出处:拓端数据部落公众号
向量自回归 (VAR) 是一种用于多变量时间序列分析的统计模型,尤其是在变量具有相互影响关系的时间序列中,本视频中我们介绍了向量自回归并在R软件中进行实现。
视频:向量自回归VAR数学原理及R软件经济数据脉冲响应分析实例
【视频】向量自回归VAR数学原理及R语言软件经济数据脉冲响应分析实例
,时长12:01
为什么用向量自回归
为了能够理解几个变量之间的关系。允许动态变化。
为了能够得到更好的预测。
一组时间序列由多个单一序列组成。
我们在建立时间序列模型时说,简单的单变量ARMA 模型可以很好地进行预测。那么,为什么我们需要多个序列?
例子如:CPI反映的是通胀,CPI高了,通胀风险大,而抑制通胀最重要的手段就是加息,反之,当CPI很低,就说明经济不景气,那么就需要降息。降息之后刺激经济增长。因此,可能需要一个联合的动态模型来了解动态的相互关系 并可能做一个更好的预测工作。
在观察 ARMA 和 GARCH 模型时,您会立即注意到估计和预测是针对一个变量进行的。在现实生活中,这并不成立。实际上,还有许多其他变量可能会影响其他变量。市场参与者和经济学家总是对宏观经济变量与他们有兴趣购买的资产之间的动态关系感兴趣。此操作可以帮助他们预测市场上可能发生的潜在情况。
使用 VAR 模型的基本要求是:
具有至少两个变量的时间序列。
变量之间存在动态关系。
它被认为是一个自回归模型,因为模型所做的预测取决于过去的值,这意味着每个观测值都被建模为其滞后值的函数。
ARIMA 和 向量自回归 模型之间的基本区别在于,所有 ARIMA 模型都用于单变量时间序列,其中 向量自回归 模型适用于多变量时间序列。此外,ARIMA 模型是单向模型,这意味着因变量受其过去值或滞后值本身的影响,其中 向量自回归 是双向模型,这意味着因变量受其过去值或另一个变量值的影响或受这两件事的影响。
什么是向量自回归
向量自回归模型是统计分析中经常使用的模型,它探索了几个变量之间的相互关系。
在开始建模部分之前,让我们先了解一下模型背后的数学。
单变量时间序列的典型自回归模型 (AR(p)) 可以表示为
其中
y t -i 表示较早时期的变量值。
A是一个时不变的 ( k × k ) 矩阵。
e t是一个误差项。
c是模型的截距。
这里的阶数 p 的意思是,最多使用 y 的 p滞后。
众所周知,向量自回归模型处理的是多元时间序列,这意味着会有两个或多个变量相互影响。因此,向量自回归 模型方程随着时间序列中变量数量的增加而增加。
假设有两个时间序列变量 y1 和 y2,因此要计算 y1(t),向量自回归 模型将使用两个时间序列变量的滞后。
例如,具有两个时间序列变量(y1 和 y2)的 VAR(1) 模型的方程如下所示:
其中,Y{1,t-1} 是 y1 的第一个滞后值,Y{2,t-1} 是 y2 的第一个滞后值。
我们可以清楚地了解模型的方程将如何随着变量和滞后值的增加而增加。例如,具有 3 个时间序列变量的 VAR(3) 模型方程如下所示。
所以这就是 p 值将如何增加模型方程的长度,而变量的数量将增加方程的高度。
选择模型的滞后数
有两种主要方法可以选择模型的滞后数:
经验方法,我们使用信息标准
推理方法包括使用假设检验
我们只考虑信息标准。有 3 个流行的信息标准,即:
赤池信息准则 (AIC)
施瓦茨-贝叶斯 (BIC)
汉南-奎恩 (HQ)
实际上,最佳滞后数是信息标准最小的滞后数。然后,我们估计 p=0,...,pmax 的 VAR(p) 并选择最小化 AIC、BIC 或 HIQ 的值 p。
向量自回归模型的估计包括以下步骤:
选择最佳滞后长度
信息标准 (IC) 用于确定最佳滞后长度。最常用的是 Akaike IC、Hannah-Quinn 准则和 Schwarz 准则。
平稳性检验
下一步是估计变量的平稳性。一种广泛使用的估计平稳性的方法是增广迪基-富勒检验和菲利普斯-佩隆检验。如果变量是非平稳的,则应采用一阶差分并以相同的方式测试平稳性。
协整检验
变量可能是非平稳的,但具有相同阶数的积分。在这种情况下,可以使用矢量纠错模型 (VECM) 而不是 向量自回归 来分析它们。如果变量是协整的,则在以下分析中应用 VECM 而不是 向量自回归 模型。VECM 应用于非变换的非平稳序列,而 向量自回归 使用变换的或平稳的输入。
模型估计
使用选择的滞后数和具有标准误差的系数运行 向量自回归 模型,并计算相应的 t 统计量以评估统计显着性。
诊断测试
接下来,使用 Breusch-Godfrey 检验对模型进行序列相关性检验,使用 Breusch-Pagan 检验检验异方差性和稳定性。
脉冲响应函数 (IRF)
IRF 用于以图形方式表示 向量自回归 模型的结果,并预测变量对彼此的影响。
格兰杰因果检验
这些变量可能是相关的,但它们之间可能不存在因果关系,或者影响可能是双向的。Granger 检验表明变量之间的因果关系,并根据 向量自回归 系统中一对变量的当前值和过去值的交互作用显示因果关系的方向。
向量自回归 在以下几种情况下很有用
向量自回归 面临的一个批评是它们是非理论的。也就是说,它们不是建立在某些将理论结构强加于方程的经济理论之上。假设每个变量都会影响系统中的所有其他变量,这使得对估计系数的直接解释变得困难。尽管如此,向量自回归 在以下几种情况下很有用:
1. 在不需要明确解释的情况下预测相关变量的集合;
2. 测试一个变量是否对预测另一个变量有用(格兰杰因果检验的基础);
3. 脉冲响应分析,分析一个变量对另一个变量突然但暂时变化的响应;
4. 预测误差方差分解,其中每个变量的预测方差比例归因于其他变量的影响。
R语言用向量自回归(VAR)进行经济数据脉冲响应研究分析
自从Sims(1980)发表开创性的论文以来,向量自回归模型已经成为宏观经济研究中的关键工具。这篇文章介绍了VAR分析的基本概念,并指导了简单模型的估算过程。
单变量自回归
VAR代表_向量自回归_。为了理解这意味着什么,让我们首先来看一个简单的单变量(即仅一个因变量或内生变量)自回归(AR)模型,其形式为yt=a1yt−1+et。
平稳性
在估算此类模型之前,应始终检查所分析的时间序列是否稳定,即它们的均值和方差随时间变化是恒定的,并且不显示任何趋势行为。
有一系列统计检验,例如Dickey-Fuller,KPSS或Phillips-Perron检验,以检验序列是否稳定。另一种非常常见的做法是绘制序列并检查其是否围绕恒定的平均值(即水平线)移动。如果是这种情况,它很可能是稳定的。
自回归滞后模型
像AR(p)模型一样,仅凭其自身的滞后对宏观经济变量进行回归可能是一种限制性很大的方法。通常,更合适的假设是还有其他因素。通过包含因变量的滞后值以及其他(即,外生)变量的同期和滞后值的模型来实现这种想法。同样,这些外生变量应该是稳定的。对于内生变量yt和外生变量xt例如_自回归分布滞后_或ADL,模型可以写成
yt=a1yt−1+b0xt+b1xt−1+et.
这种ADL模型的预测性能可能会比简单的AR模型更好。但是,如果外生变量也依赖于内生变量的滞后值怎么办?这意味着xt也是内生的,还有进一步的空间可以改善我们的预测。
向量自回归模型
因此,如上所述,VAR模型可以重写为一系列单独的ADL模型。实际上,可以通过分别估计每个方程来估计VAR模型。
标准VAR模型的协方差矩阵是_对称的_,即,对角线右上角的元素(“上三角”)将对角线左下角的元素(“下三角”)镜像。这反映了这样一种想法,即内生变量之间的关系仅反映相关性,并且不允许做出因果关系的陈述,因为在每个方向上的影响都是相同的。
在所谓的_结构化_ VAR(SVAR)模型的背景下分析了同时因果关系,或更确切地说,是变量之间的结构关系,该模型对协方差矩阵施加了限制 。
在本文中,我考虑VAR(2)过程。
此示例的人工样本是在R中生成的
set.seed(123) # 由于可复制性的考虑,重置随机数发生器 # 生成样本 t <- 200 # 时间序列观察数 k <- 2 # 内生变量数 p <- 2 # 滞后阶数 # 生成系数矩阵 A.1 <- matrix(c(-.3, .6, -.4, .5), k) # 滞后系数矩阵1 A.2 <- matrix(c(-.1, -.2, .1, .05), k) # 滞后系数2 A <- cbind(A.1, A.2) # 系数矩阵 # 生成序列 series <- matrix(0, k, t + 2*p) # 带有0的原始序列 for (i in (p + 1):(t + 2*p)){ # 生成e ~ N(0,0.5)的序列 series[, i] <- A.1%*%series[, i-1] + A.2%*%series[, i-2] + rnorm(k, 0, .5) } series <- ts(t(series[, -(1:p)])) # 转换为时间序列格式 names <- c("V1", "V2") # 重命名变量 plot.ts(series) # 绘制序列
估算值
简单VAR模型的参数和协方差矩阵的估计很简单。
为了估计VAR模型,加载并指定数据(
y
)和 模型。比较
VAR分析中的一个中心问题是找到滞后的阶数,以产生最佳结果。模型比较通常基于信息标准,例如AIC,BIC或HQ。通常,由于是小样本预测,AIC优于其他标准。但是,BIC和HQ在大型样本中效果很好 。
可以计算标准信息标准以找到最佳模型。在此示例中,我们使用AIC:
通过查看,
summary
我们可以看到AIC建议使用2的阶数。summary(var.aic)
## VAR Estimation Results: ## ========================= ## Endogenous variables: Series.1, Series.2 ## Deterministic variables: none ## Sample size: 200 ## Log Likelihood: -266.065 ## Roots of the characteristic polynomial: ## 0.6611 0.6611 0.4473 0.03778 ## Call: ## VAR(y = series, type = "none", lag.max = 5, ic = "AIC") ## ## ## Estimation results for equation Series.1: ## ========================================= ## Series.1 = Series.1.l1 + Series.2.l1 + Series.1.l2 + Series.2.l2 ## ## Estimate Std. Error t value Pr(>|t|) ## Series.1.l1 -0.19750 0.06894 -2.865 0.00463 ** ## Series.2.l1 -0.32015 0.06601 -4.850 2.51e-06 *** ## Series.1.l2 -0.23210 0.07586 -3.060 0.00252 ** ## Series.2.l2 0.04687 0.06478 0.724 0.47018 ## --- ## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 ## ## ## Residual standard error: 0.4638 on 196 degrees of freedom ## Multiple R-Squared: 0.2791, Adjusted R-squared: 0.2644 ## F-statistic: 18.97 on 4 and 196 DF, p-value: 3.351e-13 ## ## ## Estimation results for equation Series.2: ## ========================================= ## Series.2 = Series.1.l1 + Series.2.l1 + Series.1.l2 + Series.2.l2 ## ## Estimate Std. Error t value Pr(>|t|) ## Series.1.l1 0.67381 0.07314 9.213 < 2e-16 *** ## Series.2.l1 0.34136 0.07004 4.874 2.25e-06 *** ## Series.1.l2 -0.18430 0.08048 -2.290 0.0231 * ## Series.2.l2 0.06903 0.06873 1.004 0.3164 ## --- ## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 ## ## ## Residual standard error: 0.4921 on 196 degrees of freedom ## Multiple R-Squared: 0.3574, Adjusted R-squared: 0.3443 ## F-statistic: 27.26 on 4 and 196 DF, p-value: < 2.2e-16 ## ## ## ## Covariance matrix of residuals: ## Series.1 Series.2 ## Series.1 0.21417 -0.03116 ## Series.2 -0.03116 0.24154 ## ## Correlation matrix of residuals: ## Series.1 Series.2 ## Series.1 1.000 -0.137 ## Series.2 -0.137 1.000` 仔细观察结果,我们可以将真实值 与模型的参数估计值进行比较:
真实值
A ## [,1] [,2] [,3] [,4] ## [1,] -0.3 -0.4 -0.1 0.10 ## [2,] 0.6 0.5 -0.2 0.05 # Extract coefficients, standard errors etc. from the object # produced by the VAR function est_coefs <- coef(var.aic) # 仅提取两个因变量的系数,并将它们组合为一个矩阵 # 输出四舍五入的估计值 round(est_coefs, 2) ## Series.1.l1 Series.2.l1 Series.1.l2 Series.2.l2 ## [1,] -0.20 -0.32 -0.23 0.05 ## [2,] 0.67 0.34 -0.18 0.07
所有估计值都有正确的符号,并且相对接近其真实值。
脉冲响应
一旦我们确定了最终的VAR模型,就必须解释其估计的参数值。由于VAR模型中的所有变量都相互依赖,因此单个参数值仅提供 有限信息。为了更好地了解模型的动态行为,使用了脉冲响应(IR)。可以绘制因变量的轨迹,产生在许多宏观论文中都可以找到的那些波浪曲线。
在下面的示例中,我们想知道受到冲击后序列2的行为。指定了我们想要脉冲响应的模型和变量后,我们将时间范围设置
n.ahead
为20。该图给出了序列2的响应。# 计算脉冲响应 # 绘制脉冲响应 plot(ir.1)
请注意,_正交_选项很重要,因为它说明了变量之间的关系。在我们的示例中,我们已经知道不存在这样的关系,因为真正的方差-协方差矩阵(或简称协方差矩阵)在非对角元素中是对角为零的对角线。但是,由于具有200个观测值的有限时间序列数据限制了参数估计的精度,因此协方差矩阵的非对角元素具有正值,这意味着 非零同时效应。为了在IR中排除这种情况,我们设置了
ortho = FALSE
。结果是,脉冲响应在周期0中从零开始。也可以尝试另一种方法并进行设置ortho = TRUE
,那么绘图从零开始。要了解这一点,还可以计算并绘制_累积_脉冲响应函数,以了解 总体长期影响:
# 计算脉冲响应 # 绘图 plot(ir.2)
我们看到,尽管序列2对序列1中的 反应在某些时期是负面的,但总体效果却是显着正面。
点击标题查阅往期内容
向量自回归VAR的迭代多元预测估计 GDP 增长率时间序列|数据分享
ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列
向量自回归(VAR)模型分析消费者价格指数 (CPI) 和失业率时间序列
Matlab创建向量自回归(VAR)模型分析消费者价格指数 (CPI) 和失业率时间序列
Stata广义矩量法GMM面板向量自回归 VAR模型选择、估计、Granger因果检验分析投资、收入和消费数据
R语言时变向量自回归(TV-VAR)模型分析时间序列和可视化
R语言arima,向量自回归(VAR),周期自回归(PAR)模型分析温度时间序列
R语言随机搜索变量选择SSVS估计贝叶斯向量自回归(BVAR)模型
R语言用LASSO,adaptive LASSO预测通货膨胀时间序列
Python中的ARIMA模型、SARIMA模型和SARIMAX模型对时间序列预测
-
第十一章向量自回归 ( VAR) 模型和向量误差修正 (VEC)模型
2022-03-23 12:26:32第十一章向量自回归 ( VAR) 模型和向量误差修正 (VEC)模型 -
VAR(向量自回归)模型
2020-06-14 10:03:34VAR(向量自回归)模型是基于数据统计性质建立起来的模型,它把系统中的每个内生变量作为系统里所有其它内生变量滞后值的函数进行构建模型,从而把单变量的自回归模型推广到了多元时间序列组成的向量自回归模型。... -
拓端tecdat:Matlab创建向量自回归(VAR)模型分析消费者价格指数 (CPI) 和失业率时间序列
2021-11-19 12:12:36var对象指定了p阶平稳的多变量向量自回归模型(VAR(p))模型的函数形式并存储了参数值。 varm对象的关键组成部分包括时间序列的数量和多元自回归多项式 ( p )的阶数,因为它们完全指定了模型结构。其他模型组件包括... -
区制转移向量自回归模型
2014-03-14 19:09:29OXMSVAR里面包括详细的ox使用资源 -
向量自回归模型(VAR)到底厉害在哪里?
2020-12-09 10:59:21下面我们主要阐述 VAR的性质,内生变量与外生变量,VAR的应用(厉害在哪里) 1)VAR的性质 向量自回归模型的数学表达如下,Zt是多元变量在t时刻的取值,通过时间推移算子可以将向量自回归的过程紧凑的表达出来;... -
结构向量自回归(SVAR)模型(二):操作步骤与结果解读
2021-04-23 23:21:07目录第一部分 经典论文解读第二部分 操作步骤余结果解读1、单位根检验,季节调整2、最优滞后阶数得选择3、建立约束矩阵4、建模分析5、模型稳定性的检验6、进行脉冲响应7、方差分解在20世纪80年代,传统的联立方程... -
向量自回归VAR模型、结构向量自回归SVAR模型、VMA模型
2020-08-20 17:30:54