时间序列分析VAR模型实验.docx

1、时间序列分析VAR模型实验基于VAR模型的我国房地产市场与汇率波动的因果关系 VAR模型实验第一部分实验分析目的及方法现选取人民币对美元汇率以及商品房房价作为变量构建 VAR模型。对于不满足单位根检验的序 列采取对数化或差分处理， 使其成为平稳序列再进行模型的拟合。 对于商品房房价这一变量， 由于 全国各省市差异较大， 故此处采用全国房地产开发业综合景气指数这一变量。 此外， 为了消除春节 假期不固定因素带来的影响，增强数据的可比性，按照国家统计制度，从 2012 年起，不单独对 1月份统计数据进行调查， 1-2 月份数据一起调查，一起发布。所以国房景气指数 p 这一序列缺少每 年一月份的相关

2、数据，属于非随机、不可忽略缺失，在此采用平均值填充的方法，补足数据。第二部分实验样本数据来源数据来源于中经网统计数据库。具体数据见附录表。所选数据变量由于我国于 2005 年 7 月实行第二次汇改，此次汇改以市场供求为基础、参考一篮子货币进行 调节、有管理的浮动汇率制度取代了过去人民币汇率长达 10 年的紧盯美元的固定汇率体制。故本 实验拟选取2005年07月到2014年10月我国以月为单位的数据。，用以上两个变量来构建 VAR模 型，并利用该模型进行分析预测。第四部分模型构建判断序列的平稳性首先绘制出 E 的折线图，结果如下图：图汇率E的曲线图从图中可以看出， 汇率 E 序列较强的趋势性，由

3、此可以初步判断该序列是非平稳的。为了减少m的变动趋势以及异方差性，先对 m进行对数化处理，记为Im，其时序图如下:图的曲线图对数化后的趋势性减弱，但仍存在一定的趋势性，下面对Im进行一阶差分处理, 去除趋势性，得到新变量dim，观察dim的曲线图。图的曲线图从图中可以看出，die序列的趋势性基本已经消除，且新变量 die基本围绕0上下波动，因此选择形式为yt=yt-i +ut进行单位根检验：表单位根输出结果NullHypothesis:DLEhasaunitrootExogenous:ConstantLagLength:2(Automatic-basedonSIC,maxlag=12)t-St

4、atisticProb.AugmentedDickey-FullerteststatisticTestcriticalvalues: 1%level5%level10%levelIncludedobservations:108afteradjustmentsVariable Coefficient t-Statistic Prob.DLE(-1)D(DLE(-1)D(DLE(-2)C单位根统计量ADF=J、于临界值，且P为，因此该序列不是单位根过程，即该序 列是平稳序列。国房景气指数P序列首先作出P序列的时序图：图的曲线图由于每年一月份的数据缺失，故取相邻两项进行平均补全数据，得到新序列的时序

5、图如下：图的曲线图(补全)由上图可知，该序列P可能存在一定的趋势性和季节性， 先进行单位根检验，确定改序列是否平稳。由于序列表单位根输出结果NullHypothesis:PhasaunitrootExogenous:Constant, LinearTrendLagLength:3(Automatic-basedonSIC,maxlag=12)AugmentedDickey-FullerteststatisticTestcriticalvalues:1%level5%level10%level*MacKinnon(1996)one-sidedp-values.由单位根检验结果可知，T值小于临界值

6、，且 P=,在5%勺置信水平下，该序列不存在单位根过程。由于汇率E序列为一阶单整序列，并进行了一阶差分处理， 因此样本数量减少，在下面的操作中，所有的样本序列调整为 2005-08至2014-10。模型参数识别先进行VAR模型的拟合，初步选定滞后阶数为 3：表拟合输出结果VectorAutoregressionEstimatesDate:11/22/14Time:22:20Sample(adjusted):2005M112014M10lncludedobservations:108afteradjustmentsStandarderrorsin( )& t-statisticsinDLE PD

7、LE(-1) DLE(-2) DLE(-3)P(-1)P(-2)P(-3)CR-squared再进行滞后阶数的确定:表最优滞后阶数的判断VARLagOrderSelectionCriteriaEndogenousvariables:DLEPExogenousvariables:CDate:11/22/14Time:22:22Sample:2005M072014M10lncludedobservations:99Lag LogLLR FPE AIC SC HQ01NA2* 34-k k k56789101112SC:SchwarzinformationcriterionHQ:Hannan-Qui

8、nninformationcriterion2(SC由上边可知，根据信息准则，采取少数服从多数原则，取滞后阶数为3,此外取滞后阶数为为或取滞后阶数为3(SC为时，两者SC值相差不是很大。模型参数估计选取了最优滞后阶数3,进行模型的拟合。拟合结果如下: 表(3)模型估计结果VectorAutoregressionEstimatesDate:11/22/14Time:22:23Sample(adjusted):2005M112014M10lncludedobservations:108afteradjustmentsStandarderrorsin( )& t-statisticsinDLE PR

9、-squaredP(-2)P(-3)C由回归结果可知，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:13Date:11/22/14Time:22:

10、27RootModulusInverse Roots of AR Characteristic Polynomial1_5 1 .o - 0.5 - *0.0 - + * *-0.5 - 40 - * * H=4-51 | | | | | -1.5 -10 -0.5 00 0.5 1.0 1.5图模型平稳性检验的图形显示由上表和上图可知，VAR模型的特征方程的根均在单位园内，因此 VAR模型是平稳的。下面进行残差的自相关性的检验，检验结果如下：图模型各方程残差项的自相关图由上图可知，VAR模型允许不同方程的残差之间存在交叉相关性，但是残差自 身不存在自相关性，因此，观察残差自身的自相关图，可以

11、看出自相关系数均位于 置信区间内，说明残差不存在自相关性。第五部分模型应用格兰杰因果检验接下来做两两变量之间的格兰杰因果检验。序列 P与序列DLE:表序列P与序列DLE格兰杰因果检验表PairwiseGrangerCausalityTestsDate:11/21/14Time:23:32Sample:2005M072014M10Lags:3NullHypothesis:Obs F-StatisticProb.PdoesnotGrangerCauseDLE108DLEdoesnotGrangerCauseP由上述结果可知，在5%的置信水平下，P是dle的格兰杰原因，即全国房地产开发业综合景气指

12、数是人民币对美元汇率变动幅度的格兰杰原因。脉冲响应由于脉冲响应函数收到变量顺序的影响，因此其结果与分析的主观因素有关，对于这三个变量：DLE R、P,按照中国市场目前现状，认为 DLE外生性最强，p其次最后为r。故选取顺序为DLE P、Ro图脉冲响应图方差分解表方差分解结果VarianceDecomp osition ofDLE:PeriodDLE P123456789101112131415161718192021222324252627282930313233343536VarianceDecompositionofP:PeriodDLE附录具体数据国房景气指数_当人民币对美元期末指标汇率

13、地区全国全国频度月月123456789101112131415161718192021222324252627282930313233343536单位- 人民币/美元2005-072005-082005-092005-102005-112005-122006-012006-022006-032006-042006-052006-0682006-072006-082006-092006-102006-112006-122007-012007-022007-032007-042007-052007-062007-071042007-082007-092007-102007-112007-12200

14、8-012008-022008-032008-0472008-052008-062008-072008-082008-092008-102008-112008-122009-012009-022009-032009-042009-052009-062009-072009-082009-092009-102009-112009-122010-012010-022010-032010-042010-052010-062010-072010-082010-092010-102010-112010-122011-012011-022011-032011-042011-052011-062011-072

15、011-082011-092011-102011-112011-122012-012012-022012-032012-042012-052012-062012-072012-082012-092012-102012-112012-122013-012013-022013-032013-042013-052013-062013-072013-082013-092013-102013-112013-122014-012014-022014-032014-042014-052014-062014-072014-082014-092014-10存在问题1、 在进行单位根检验时，对于阶数的确定根据信息准则，那么在单位根检验之前是需要先进行回归，确定阶数后在进行单位根检验么2、 当数据不全时如何补全数据。3、 初步选择滞后阶数时，是否可以参考单位根检验的滞后阶数。4、 分析脉冲响应函数图时，位于主对角线之外的图是否应该满足同向影响的关系，即分析 A、B两个变量的脉冲响应函数，若 A对B是负向影响,那B对A也应该是负向影响。