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

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

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

3、E的折线图，结果如下图：图 汇率E的曲线图从图中可以看出，汇率E序列较强的趋势性，由此可以初步判断该序列是非平稳的。为了减少m的变动趋势以及异方差性，先对m进行对数化处理，记为lm，其时序图如下：图 lm的曲线图对数化后的趋势性减弱，但仍存在一定的趋势性，下面对lm进行一阶差分处理，去除趋势性，得到新变量dlm，观察dlm的曲线图。图 DLE的曲线图从图中可以看出，dle序列的趋势性基本已经消除，且新变量dle基本围绕0上下波动，因此选择形式为yt=yt-1+ut 进行单位根检验：表 单位根输出结果Null Hypothesis: DLE has a unit rootExogenous: C

4、onstantLag Length: 2 (Automatic - based on SIC, maxlag=12)t-StatisticProb.*Augmented Dickey-Fuller test statisticTest critical values:1% level5% level10% level*MacKinnon (1996) one-sided p-values.Augmented Dickey-Fuller Test EquationDependent Variable: D(DLE)Method: Least SquaresDate: 11/15/14 Time:

5、 20:20Sample (adjusted): 2005M11 2014M10Included observations: 108 after adjustmentsVariableCoefficientStd. Errort-StatisticProb.DLE(-1)D(DLE(-1)D(DLE(-2)CR-squaredMean dependent varAdjusted R-squared. dependent var. of regressionAkaike info criterionSum squared residSchwarz criterionLog likelihoodH

6、annan-Quinn criter.F-statisticDurbin-Watson statProb(F-statistic)单位根统计量ADF=小于临界值，且P为，因此该序列不是单位根过程，即该序列是平稳序列。国房景气指数P序列首先作出P序列的时序图：图 P的曲线图由于每年一月份的数据缺失，故取相邻两项进行平均补全数据，得到新序列的时序图如下：图 P的曲线图(补全)由上图可知，该序列P可能存在一定的趋势性和季节性，先进行单位根检验，确定改序列是否平稳。由于序列表 单位根输出结果Null Hypothesis: P has a unit rootExogenous: Constant, L

7、inear TrendLag Length: 3 (Automatic - based on SIC, maxlag=12)t-StatisticProb.*Augmented Dickey-Fuller test statisticTest critical values:1% level5% level10% level*MacKinnon (1996) one-sided p-values.由单位根检验结果可知，T值小于临界值，且P=，在5%的置信水平下，该序列不存在单位根过程。由于汇率E序列为一阶单整序列，并进行了一阶差分处理，因此样本数量减少，在下面的操作中，所有的样本序列调整为20

8、05-08至2014-10。模型参数识别先进行VAR模型的拟合，初步选定滞后阶数为3：表 拟合输出结果Vector Autoregression EstimatesDate: 11/22/14 Time: 22:20Sample (adjusted): 2005M11 2014M10Included observations: 108 after adjustmentsStandard errors in ( ) & t-statistics in DLEPDLE(-1) DLE(-2) DLE(-3) P(-1) P(-2)P(-3) C R-squaredAdj. R-squaredSum

9、 sq. resids. equationF-statisticLog likelihoodAkaike AICSchwarz SCMean dependent. dependentDeterminant resid covariance (dof adj.)Determinant resid covarianceLog likelihoodAkaike information criterionSchwarz criterion再进行滞后阶数的确定：表 最优滞后阶数的判断VAR Lag Order Selection CriteriaEndogenous variables: DLE P?E

10、xogenous variables: C?Date: 11/22/14 Time: 22:22Sample: 2005M07 2014M10Included observations: 99LagLogLLRFPEAICSCHQ0NA?12*3*456789101112* indicates lag order selected by the criterionLR: sequential modified LR test statistic (each test at 5% level)FPE: Final prediction errorAIC: Akaike information c

11、riterionSC: Schwarz information criterionHQ: Hannan-Quinn information criterion由上边可知，根据信息准则，采取少数服从多数原则，取滞后阶数为3，此外取滞后阶数为2(SC为 或取滞后阶数为3(SC为时，两者SC值相差不是很大。模型参数估计选取了最优滞后阶数3，进行模型的拟合。拟合结果如下：表 VAR(3)模型估计结果Vector Autoregression EstimatesDate: 11/22/14 Time: 22:23Sample (adjusted): 2005M11 2014M10Included obs

12、ervations: 108 after adjustmentsStandard errors in ( ) & t-statistics in DLEPDLE(-1) DLE(-2) DLE(-3) P(-1) P(-2)P(-3) C R-squaredAdj. R-squaredSum sq. resids. equationF-statisticLog likelihoodAkaike AICSchwarz SCMean dependent. dependentDeterminant resid covariance (dof adj.)Determinant resid covari

13、anceLog likelihoodAkaike information criterionSchwarz criterion由回归结果可知，VAR模型的参数估计一部分显着。估计的方程为：模型检验首先对模型进行平稳性检验表 VAR模型平稳性检验的表格显示Roots of Characteristic PolynomialEndogenous variables: DLE P?Exogenous variables: C?Lag specification: 1 3Date: 11/22/14 Time: 22:27RootModulus - + - + No root lies outside

14、 the unit circle.VAR satisfies the stability condition.图 VAR模型平稳性检验的图形显示由上表和上图可知，VAR模型的特征方程的根均在单位园内，因此VAR模型是平稳的。下面进行残差的自相关性的检验，检验结果如下：图 VAR模型各方程残差项的自相关图由上图可知，VAR模型允许不同方程的残差之间存在交叉相关性，但是残差自身不存在自相关性，因此，观察残差自身的自相关图，可以看出自相关系数均位于置信区间内，说明残差不存在自相关性。第五部分 模型应用格兰杰因果检验接下来做两两变量之间的格兰杰因果检验。序列P与序列DLE:表 序列P与序列DLE格兰杰

15、因果检验表Pairwise Granger Causality TestsDate: 11/21/14 Time: 23:32Sample: 2005M07 2014M10Lags: 3Null Hypothesis:ObsF-StatisticProb.P does not Granger Cause DLE108DLE does not Granger Cause P由上述结果可知，在5%的置信水平下，P是dle的格兰杰原因，即全国房地产开发业综合景气指数是人民币对美元汇率变动幅度的格兰杰原因。脉冲响应由于脉冲响应函数收到变量顺序的影响，因此其结果与分析的主观因素有关，对于这三个变量：DL

16、E、R、P,按照中国市场目前现状，认为DLE外生性最强，p其次最后为r。故选取顺序为DLE、P、R。图 脉冲响应图方差分解 表 方差分解结果Variance Decomposition of DLE:Period.DLEP123456789101112131415161718192021222324252627282930313233343536Variance Decomposition of P:Period.DLEP123456789101112131415161718192021222324252627282930313233343536附录具体数据指标国房景气指数_当月人民币对美元期

17、末汇率地区全国全国频度月月单位-人民币/美元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-122008-012008-022008-032008-0472008-052008-062008-072008-082008-09

18、2008-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-072011-082011-092011-102011-112011-122012-012012-022012-032012-0

19、42012-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也应该是负向影响。