SPSSAU-在线SPSS分析软件

var模型与vec模型区别(VAR模型Vectorautoregressive)(1)


VAR模型

在时间序列进行预测时, ARIMA可用于单一变量(比如GDP增长率)的预测,如果需要同时考虑几个变量的预测时(比如GDP增长率、失业率、储蓄率),此时可考虑分别针对研究变量进行,即多次重复进行。通常情况下同一系统的几个研究变量之间均有着相互依旧关系,因而为更好的利用各变量的此类关系,此时可以使用VAR模型(Vector autoregressive model)进行多变量预测。VAR模型的构建流程较为复杂,如下图所述:

var模型与vec模型区别(VAR模型Vectorautoregressive)(2)

通常情况下,VAR模型需要满足单位根检验,如果没有单位根则直接构建VAR模型即可,如果研究变量有单位根,则说明不适合进行VAR模型构建,但是如果有单位根且满足同阶单整,此时说明VAR模型构建是适合的,与此同时研究变量满足协整关系也是一种常见的前提条件。

VAR模型构建时,通常包括定阶这一步骤,即选择适合的滞后阶数。VAR模型构建完成后,接着还需要对模型的有效性进行分析,通常是针对AR特征根图进行分析。另外,理论上VAR模型的残差还满足满足正态性,并且通过自相关检验等,但通常对此类检验的关注度相对较少。

VAR模型构建之后,通常需要进行比如格兰杰困果检验,脉冲响应和方差分析,用于进一步分析研究变量之间的相互作用依存关系情况。

最后,可得到模型的预测数据,满足模型预测目的。

特别提示:

VAR模型案例

1 背景

当前有一项美国宏观联邦基金利率、通货膨胀率和失业率的数据,数据日期从1960年第1季度到2012年第1季度,单位为季度,共计209个数据。现希望建立该3项数据的VAR模型,通过VAR模型查看该3项的相互依存影响关系,并且进行模型预测等。部分数据如下图所示:

var模型与vec模型区别(VAR模型Vectorautoregressive)(3)

2 理论

VAR模型有着较为复杂的分析流程步骤,具体可参考下图和表格所述。下述分析步骤和流程仅供参考使用,实际研究中可能仅需要其中一部分的分析即可。比如很多时候并不需要关注残差自相关检验和残差正态性检验,也或者有时对格兰杰因果检验关注度较少等,具体以研究者实际研究情况为准即可。

var模型与vec模型区别(VAR模型Vectorautoregressive)(4)

1、关于单位根检验。

2、关于协整关系。

3、关于定阶。

4、关于AR特征根检验。

5、关于格兰杰因果。

6、关于脉冲响应。

7、关于方差分解。

8、关于残差自相关检验。

9、关于残差正态性检验。

10、关于预测。

3 操作

本例子操作截图如下,放入3个分析项(并且分析项需要注意放入顺序,便于后续分析使用),SPSSAU‘Trend类型’默认为常数,滞后阶数lags上默认是自动定阶,预测期数设置为40,即往后预期40期数据。

var模型与vec模型区别(VAR模型Vectorautoregressive)(5)

关于‘Trend类型’说明:本案例的3个数据并没有明显的趋势而且均值也不接近于0,因而选择默认的常数类型即可,绝大多数时候均应该使用常数类型。如果说研究变量的均值基本上接近于0,可以考虑使用‘None’类型,如果说研究变量有着较为明显的趋势项,此时可考虑使用‘线性趋势T’。具体更多说明可见本文档末尾的疑难解惑说明。

滞后阶数让系统自动定阶,当然研究者也可以自主设置阶数,具体设置原则建议结合专业知识,或者信息准则越小越好进行判断,多次构建模型对比不同滞后阶数时信息准则的大小,选择出信息准则最小对应的滞后阶数模型即可。

预测基数上,本研究设置为40,即向后预测40个季度的数据。

4 SPSSAU输出结果

SPSSAU分别8类表格,以及4种图形,分别如下说明。

var模型与vec模型区别(VAR模型Vectorautoregressive)(6)

5文字分析

通常情况下VAR模型构建之前需要查看研究变量的稳定性,可通过时序图或者单位根检验进行分析,如果存在单位根或者图示出来明显不平稳,则不能构建VAR模型。如果说差分后数据平稳(即同阶单整),也可以进行VAR模型构建,本次研究数据本身不平稳,但是1阶差分后的数据平稳即满足同阶单整前提。因而可以构建VAR模型,3个研究项进行1阶差分后的时序图如下图所示:

var模型与vec模型区别(VAR模型Vectorautoregressive)(7)

与此同时VAR模型之前一般还可进行协整检验,本次构建模型数据进行‘Johansen协整检验’后显示具有协整关系,即说明模型参数值具有长期稳定性,也即说明模型可靠,由于coint协整检验在SPSSAU中隶属于单独一个研究方法(其并不隶属于VAR模型),研究者可单独进行分析处理即可;除此之外,VAR模型分析后通常还需要研究格兰杰因果关系,其需要单独使用SPSSAU的格兰杰因果分析方法进行研究即可。

接下来逐一描述VAR模型输入的表格及图形结果等,如下:

var模型与vec模型区别(VAR模型Vectorautoregressive)(8)

上表格展示VAR模型参数设置情况,包括趋势类型、滞后阶数情况和预期期数;其仅为模型设置的参数展示,并不特别分析意义。本次模型时自动定阶为2阶,即构建VAR(2)模型。

var模型与vec模型区别(VAR模型Vectorautoregressive)(9)

上表格展示VAR模型参数结果值;通常情况下上表格中值的意味相对较小,如果需要分析,首先查看是否显著,如果呈现出显著性后再分析影响关系情况,通常情况下并不需要分析,本次也并不对上述参数进行分析,其仅作展示使用。另外,上表格中还提供llf似然值,以及AIC、SC和HQIC这3个信息准则值。

一般情况下,如果需要对比模型优劣,可通过对比3个信息准则值的大小(越小越优)来进行模型综合选择使用。

由于本次滞后阶数为2阶,因而本身有209个数据,但上表格中显示为207;

由于本次滞后阶数为2阶,因而参数中L1代表滞后1阶,L2代表滞后2阶。

var模型与vec模型区别(VAR模型Vectorautoregressive)(10)

上表格展示自动定阶的过程值,SPSSAU自动定阶的原则为信息准则值越小越优,共有提供4个指标值,分别是AIC、BIC、FPE和HQIC,如果按照AIC准则时应该以9阶为准, BIC准则时应该以2阶为准, FPE时应该以9阶为准, HQIC准则时应该以3阶为准。4个指标值中最小值为2阶,因而SPSSAU最终以2阶为准构建VAR模型。

一般情况下,滞后阶数越小时,可利用的样本信息越高,并且模型相对越简单,因而一般使用尽可能小的滞后阶数,建议直接使用SPSSAU的自动定阶建议值即可;

上述中共有4个判断指标,理论上研究者可使用任意其中1个作为标准,当然研究者也可分别设置不同的滞后阶数,然后对比不同模型信息准则值来判断哪个模型更优(根据信息准则越小越优标准)。

var模型与vec模型区别(VAR模型Vectorautoregressive)(11)

VAR模型构建后,可以通过AR根图判断VAR模型的稳定性;如果所有特征值均在单位圆内,即所有点均在圆内,此时说明模型具有稳定性;如果出现特征根在单位圆之外,意味着模型可能不具有长期持续稳定性。

AR根图是模型构建后模型质量的判断,一般情况需要查看此项,上图可以看到,所有特征根值均在单位圆之内,意味着构建的VAR模型稳定性较好,也即说明本次模型构建良好。

var模型与vec模型区别(VAR模型Vectorautoregressive)(12)

上图展示正交化脉冲响应图,并且下表格有具体数据。正交化脉冲响应图展示某变量冲击对该变量自身或其它变量产生的动态影响情况;如果大于0则为正向冲击,反之小于0则为负向冲击,绝对值越大意味着冲击越大,值越接近于0意味着冲击越小。

从上图可以看到:针对‘通货膨胀率’冲击变量,通货膨胀率首先对自身产生正向的冲击,然后冲击力度逐渐平稳减小,但一直均有着冲击且冲击力度较大。通货膨胀率对于失业率的冲击上,其先会产生负向抑制作用,前4期时此种负向冲击作用逐渐加强,但从第5期开始,负向冲击作用逐渐减弱并且趋于平稳,一直均为负向冲击。通货膨胀率对于联邦基金利率的冲击上,一直均有着微弱的正向冲击作用,但作用力度非常有限。

除此之外,关于失业率对于3个研究变量的脉冲响应分析表格或图,也或者联邦基金利率对于3个研究 变量的脉冲响应分析表格或图。可通过下拉选择进行切换展示分析。

var模型与vec模型区别(VAR模型Vectorautoregressive)(13)

var模型与vec模型区别(VAR模型Vectorautoregressive)(14)

上图展示预测误差方差分解情况,以及上表格展示具体数据;方差分解fevd目的在于分析变量影响比例情况;方差分解图的横坐标为阶数,纵坐标为方差分解值,其可用于反馈随着滞后阶数变化时方差分解变化情况;方差分解值越大意味着影响比例越大,方差分解值接近于0时意味着影响甚微;

从上图可以看到,通货膨胀率对于自身的方差分解率非常高,基本上均大于85%,也即说明通货膨胀率对于自身的解释力度非常强非常相对较为平稳。对于失业率而言,通货膨胀率对其解释力度有限,基本上均低于10%,但是解释力度呈现出逐渐加大的趋势性,也即说明通货膨胀率对于失业率的影响是具有连续稳定的解释作用,并且随着时间变化,解释作用会逐步提高。最后,通货膨胀率对于联邦基金利率的方差分解来看,其分解力度非常有限,基本上均接近于0,也即说明通货膨胀率对于联邦基金利率的变化解释力度很小。

var模型与vec模型区别(VAR模型Vectorautoregressive)(15)

var模型与vec模型区别(VAR模型Vectorautoregressive)(16)

上表格展示模型预测结果,本次案例设置为40期(篇幅限制,上表格中仅11期),下图可直观查看模型预期结果数据,均为向后预测数据。

var模型与vec模型区别(VAR模型Vectorautoregressive)(17)

var模型与vec模型区别(VAR模型Vectorautoregressive)(18)

上表格展示残差正态性检验结果;理论上VAR模型需要满足残差正态性要求,但多数情况下对其关注度相对较小。上表格展示为最终联合性正态检验(χ=1126.599,p值=0.000<0.05),即意味着拒绝原假定,即说明模型没有通过正态性检验。但本研究暂不关注此指标。

var模型与vec模型区别(VAR模型Vectorautoregressive)(19)

上表格展示Portmanteau自相关检验结果;理论上VAR模型需要满足残差无自相关性要求,但多数情况下对其关注度相对较小。上表格展示检验(χ=42.667,p值=0.000<0.05),即意味着拒绝原假定,即说明模型没有通过自相关性检验。但本研究暂不关注此指标。

6 剖析

涉及以下几个关键点,分别如下:

7疑难解惑

如果时间序列一阶单整和二阶单整,并且其它序列是二阶单整,是否可进行VAR模型构建?

如果时间序列经过一次差分变成平稳,则称原序列是1阶单整,记为I(1)。如果时间序列经过d次差分后变成平稳序列,则称原序列是d阶单整序列,记为I(d)。同阶单整是指两个或多个时间序列同为d阶单整。因而说明为2阶单整即I(2),此时如果原始序列通过协整检验,那么就可以进行VAR模型构建。

VAR模型时‘Trend类型’如何选择类型?

VAR模型时,SPSSAU共有三种类型可以选择,分别是截距(即常数),截距和趋势(即线性),无。其常见选择标准有两种,分别是‘主观判断法’和‘信息准则判断法’。

var模型与vec模型区别(VAR模型Vectorautoregressive)(20)

可结合SPSSAU输出的信息准则比如AIC值判断,比如分别使用3种类型得到3个不同的AIC值,然后对比找到AIC值最小时对应的类型使用即可。

VAR模型时,如果原始序列平稳此时不需要进行协整检验,直接进行VAR模型分析即可,并且可进行格兰杰因果检验。如果出现原始序列不平稳且满足同阶单整(且通过协整检验),此时是否进行格兰杰检验并没有固定要求,此时如果要进行格兰杰检验,通常应该以差分后的序列数据进行,而非原始不平稳序列数据。

,