时间序列分析VAR模型实验
时间序列分析V A R模型
实验
The following text is amended on 12 November 2020.
基于VAR模型的我国房地产市场与汇率波动
的因果关系
————VAR模型实验
第一部分实验分析目的及方法
现选取人民币对美元汇率以及商品房房价作为变量构建VAR模型。对于不满足单位根检验的序列采取对数化或差分处理,使其成为平稳序列再进行模型的拟合。对于商品房房价这一变量,由于全国各省市差异较大,故此处采用全国房地产开发业综合景气指数这一变量。此外,为了消除春节假期不固定因素带来的影响,增强数据的可比性,按照国家统计制度,从2012年起,不单独对1月份统计数据进行调查,1-2月份数据一起调查,一起发布。所以国房景气指数p这一序列缺少每年一月份的相关数据,属于非随机、不可忽略缺失,在此采用平均值填充的方法,补足数据。
第二部分实验样本
数据来源
数据来源于中经网统计数据库。具体数据见附录表。
所选数据变量
由于我国于2005年7月实行第二次汇改,此次汇改以市场供求为基础、参考一篮子货币进行调节、有管理的浮动汇率制度取代了过去人民币汇率长达10年的紧盯美元的固定汇率体制。故本实验拟选取2005年07月到2014年10月我国以月为单位的数据。,用以上两个变量来构建VAR 模型,并利用该模型进行分析预测。
第四部分模型构建
判断序列的平稳性
首先绘制出E的折线图,结果如下图:
图汇率E的曲线图
从图中可以看出,汇率E序列较强的趋势性,由此可以初步判断该序列是非平稳的。为了减少m的变动趋势以及异方差性,先对m进行对数化处理,记为lm,其时序图如下:
图的曲线图
对数化后的趋势性减弱,但仍存在一定的趋势性,下面对lm进行一阶差分处理,去除趋势性,得到新变量dlm,观察dlm的曲线图。
图的曲线图
从图中可以看出,dle序列的趋势性基本已经消除,且新变量dle基本围绕0上下波
动,因此选择形式为y
t =y
t-1
+u
t
进行单位根检验:表单位根输出结果
NullHypothesis:DLEhasaunitroot
Exogenous:Constant
LagLength:2(Automatic-basedonSIC,maxlag=12)
t-Statistic Prob.*
AugmentedDickey-Fullerteststatistic
Testcriticalvalues:1%level
5%level
10%level
*MacKinnon(1996)one-sidedp-values.
AugmentedDickey-FullerTestEquation
DependentVariable:D(DLE)
Method:LeastSquares
Date:11/15/14Time:20:20
Sample(adjusted):2005M112014M10
Includedobservations:108afteradjustments
Variable Coefficient t-Statistic Prob.
DLE(-1)
D(DLE(-1))
D(DLE(-2))
C
R-squared M eandependentvar
AdjustedR-squared
单位根统计量ADF=小于临界值,且P为,因此该序列不是单位根过程,即该序列是平稳序列。
国房景气指数P序列
首先作出P序列的时序图:
图的曲线图
由于每年一月份的数据缺失,故取相邻两项进行平均补全数据,得到新序列的时序图如下:
图的曲线图(补全)
由上图可知,该序列P可能存在一定的趋势性和季节性,先进行单位根检验,确定改序列是否平稳。
由于序列
表单位根输出结果
NullHypothesis:Phasaunitroot
Exogenous:Constant,LinearTrend
LagLength:3(Automatic-basedonSIC,maxlag=12)
t-Statistic Prob.*
AugmentedDickey-Fullerteststatistic
Testcriticalvalues:1%level
5%level
10%level
*MacKinnon(1996)one-sidedp-values.
由单位根检验结果可知,T值小于临界值,且P=,在5%的置信水平下,该序列不存在单位根过程。
由于汇率E序列为一阶单整序列,并进行了一阶差分处理,因此样本数量减少,在下面的操作中,所有的样本序列调整为2005-08至2014-10。
模型参数识别
先进行VAR模型的拟合,初步选定滞后阶数为3:
表拟合输出结果
VectorAutoregressionEstimates
Date:11/22/14Time:22:20
Sample(adjusted):2005M112014M10
Includedobservations:108afteradjustments
Standarderrorsin()&t-statisticsin[]
DLE P
DLE(-1)
[][]
DLE(-2)
[][]
DLE(-3)
[][]
P(-1)
[][]
P(-2)
[][]
P(-3)
[][]
C
[][]
R-squared
再进行滞后阶数的确定:
表最优滞后阶数的判断
VARLagOrderSelectionCriteria
Endogenousvariables:DLEP
Exogenousvariables:C
Date:11/22/14Time:22:22
Sample:2005M072014M10
Includedobservations:99
Lag LogL LR FPE AIC SC HQ
0NA
1
2** 3***
4
5
6
7
8
9
10
11
12
*indicateslagorderselectedbythecriterion
LR:sequentialmodifiedLRteststatistic(eachtestat5%level)
FPE:Finalpredictionerror
AIC:Akaikeinformationcriterion
SC:Schwarzinformationcriterion
HQ:Hannan-Quinninformationcriterion
由上边可知,根据信息准则,采取少数服从多数原则,取滞后阶数为3,此外取滞后阶数为2(SC为或取滞后阶数为3(SC为时,两者SC值相差不是很大。
模型参数估计
选取了最优滞后阶数3,进行模型的拟合。拟合结果如下:
表(3)模型估计结果
VectorAutoregressionEstimates
Date:11/22/14Time:22:23
Sample(adjusted):2005M112014M10
Includedobservations:108afteradjustments
Standarderrorsin()&t-statisticsin[]
DLE P
DLE(-1)
[][]
DLE(-2)
[][]
DLE(-3)
[][]
P(-1)
[][]
P(-2)
[][]
P(-3)
[][]
C
[][]
R-squared
由回归结果可知,VAR模型的参数估计一部分显着。估计的方程为:DLE=*DLE(-1)+*DLE(-2)+*DLE(-3)-*P(-1)-*P(-2)+*P(-3)+
P=-*DLE(-1)+*DLE(-2)+*DLE(-3)+*P(-1)-*P(-2)-*P(-3)+
模型检验
首先对模型进行平稳性检验
表模型平稳性检验的表格显示
RootsofCharacteristicPolynomial
Endogenousvariables:DLEP
Exogenousvariables:C
Lagspecification:13
Date:11/22/14Time:22:27
Root Modulus
图模型平稳性检验的图形显示
由上表和上图可知,VAR模型的特征方程的根均在单位园内,因此VAR模型是平稳的。下面进行残差的自相关性的检验,检验结果如下:
图模型各方程残差项的自相关图
由上图可知,VAR模型允许不同方程的残差之间存在交叉相关性,但是残差自身不存在自相关性,因此,观察残差自身的自相关图,可以看出自相关系数均位于置信区间内,说明残差不存在自相关性。
第五部分模型应用
格兰杰因果检验
接下来做两两变量之间的格兰杰因果检验。序列P与序列DLE:
表序列P与序列DLE格兰杰因果检验表
PairwiseGrangerCausalityTests
Date:11/21/14Time:23:32
Sample:2005M072014M10
Lags:3
NullHypothesis:Obs F-Statistic Prob.
PdoesnotGrangerCauseDLE108
DLEdoesnotGrangerCauseP
由上述结果可知,在5%的置信水平下,P是dle的格兰杰原因,即全国房地产开发业综合景气指数是人民币对美元汇率变动幅度的格兰杰原因。
脉冲响应
由于脉冲响应函数收到变量顺序的影响,因此其结果与分析的主观因素有关,对于这三个变量:DLE、R、P,按照中国市场目前现状,认为DLE外生性最强,p其次最后为r。故选取顺序为DLE、P、R。
图脉冲响应图
方差分解
表方差分解结果
Varianc
eDecomp
osition
ofDLE:
Period.DLE P
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
Varianc
eDecomp
osition
ofP:
Period.DLE P
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
附录具体数据
1、在进行单位根检验时,对于阶数的确定根据信息准则,那么在单位根检验之前是需要先进行回归,确定阶数后在进行单位根检验么
2、当数据不全时如何补全数据。
3、初步选择滞后阶数时,是否可以参考单位根检验的滞后阶数。
4、分析脉冲响应函数图时,位于主对角线之外的图是否应该满足同向影响的关系,即分析A、B两个变量的脉冲响应函数,若A对B是负向影响,那B对A也应该是负向影响。