基于VAR模型的我国房地产市场与汇率波动的因果关系————VAR模型实验第一部分实验分析目的及方法现选取人民币对美元汇率以及商品房房价作为变量构建VAR模型。
对于不满足单位根检验的序列采取对数化或差分处理,使其成为平稳序列再进行模型的拟合。
对于商品房房价这一变量,由于全国各省市差异较大,故此处采用全国房地产开发业综合景气指数这一变量。
此外,为了消除春节假期不固定因素带来的影响,增强数据的可比性,按照国家统计制度,从2012年起,不单独对1月份统计数据进行调查,1-2月份数据一起调查,一起发布。
所以国房景气指数p这一序列缺少每年一月份的相关数据,属于非随机、不可忽略缺失,在此采用平均值填充的方法,补足数据。
第二部分实验样本2.1数据来源数据来源于中经网统计数据库。
具体数据见附录表。
2.2所选数据变量由于我国于2005年7月实行第二次汇改,此次汇改以市场供求为基础、参考一篮子货币进行调节、有管理的浮动汇率制度取代了过去人民币汇率长达10年的紧盯美元的固定汇率体制。
故本实验拟选取2005年07月到2014年10月我国以月为单位的数据。
,用以上两个变量来构建VAR模型,并利用该模型进行分析预测。
第四部分模型构建4.1判断序列的平稳性4.1.1汇率E序列首先绘制出E的折线图,结果如下图:图4.1 汇率E的曲线图从图中可以看出,汇率E序列较强的趋势性,由此可以初步判断该序列是非平稳的。
为了减少m的变动趋势以及异方差性,先对m进行对数化处理,记为lm,其时序图如下:图4.2 lm的曲线图对数化后的趋势性减弱,但仍存在一定的趋势性,下面对lm进行一阶差分处理,去除趋势性,得到新变量dlm,观察dlm的曲线图。
图4.3 DLE的曲线图从图中可以看出,dle序列的趋势性基本已经消除,且新变量dle基本围绕0上下波动,因此选择形式为y t=y t-1+u t进行单位根检验:表4.1 单位根输出结果Null Hypothesis: DLE has a unit rootExogenous: ConstantLag Length: 2 (Automatic - based on SIC, maxlag=12)t-Statistic Prob.*Augmented Dickey-Fuller test statistic -3.031673 0.0351Test critical values: 1% level -3.4919285% level -2.88841110% level -2.581176*MacKinnon (1996) one-sided p-values.Augmented Dickey-Fuller Test EquationDependent Variable: D(DLE)Method: Least SquaresDate: 11/15/14 Time: 20:20Sample (adjusted): 2005M11 2014M10Included observations: 108 after adjustmentsVariable Coefficient Std. Error t-Statistic Prob.DLE(-1) -0.353005 0.116439 -3.031673 0.0031D(DLE(-1)) -0.502730 0.115417 -4.355768 0.0000D(DLE(-2)) -0.311531 0.093265 -3.340258 0.0012C -0.000888 0.000470 -1.887592 0.0619R-squared 0.450240 Mean dependent var 1.15E-05Adjusted R-squared 0.434382 S.D. dependent var 0.005058S.E. of regression 0.003804 Akaike info criterion -8.269046Sum squared resid 0.001505 Schwarz criterion -8.169708Log likelihood 450.5285 Hannan-Quinn criter. -8.228768F-statistic 28.39119 Durbin-Watson stat 2.061613Prob(F-statistic) 0.000000单位根统计量ADF=-3.031673小于临界值,且P为 0.0351,因此该序列不是单位根过程,即该序列是平稳序列。
4.1.2国房景气指数P序列首先作出P序列的时序图:图4.4 P的曲线图由于每年一月份的数据缺失,故取相邻两项进行平均补全数据,得到新序列的时序图如下:图4.5 P的曲线图(补全)由上图可知,该序列P可能存在一定的趋势性和季节性,先进行单位根检验,确定改序列是否平稳。
由于序列表4.2 单位根输出结果Null Hypothesis: P has a unit rootExogenous: Constant, Linear TrendLag Length: 3 (Automatic - based on SIC, maxlag=12)t-Statistic Prob.*Augmented Dickey-Fuller test statistic -3.972457 0.0124Test critical values: 1% level -4.0452365% level -3.45195910% level -3.151440*MacKinnon (1996) one-sided p-values.由单位根检验结果可知,T值小于临界值,且P=0.0124,在5%的置信水平下,该序列不存在单位根过程。
由于汇率E序列为一阶单整序列,并进行了一阶差分处理,因此样本数量减少,在下面的操作中,所有的样本序列调整为2005-08至2014-10。
4.2模型参数识别先进行VAR模型的拟合,初步选定滞后阶数为3:表4.3 拟合输出结果Vector Autoregression EstimatesDate: 11/22/14 Time: 22:20Sample (adjusted): 2005M11 2014M10Included observations: 108 after adjustmentsStandard errors in ( ) & t-statistics in [ ]DLE PDLE(-1) 0.063183 -19.12274(0.09626) (14.1374)[ 0.65638] [-1.35263]DLE(-2) 0.116798 15.42129(0.09604) (14.1052)[ 1.21615] [ 1.09330]DLE(-3) 0.245260 16.39171(0.09617) (14.1243)[ 2.55030] [ 1.16053]P(-1) -9.04E-05 1.490708(0.00066) (0.09765)[-0.13593] [ 15.2656]P(-2) -0.000583 -0.355442(0.00118) (0.17380)[-0.49226] [-2.04508]P(-3) 0.000346 -0.160740(0.00067) (0.09872)[ 0.51479] [-1.62821]C 0.031328 2.571540(0.01274) (1.87084)[ 2.45943] [ 1.37454]R-squared 0.295033 0.979509Adj. R-squared 0.253154 0.978292Sum sq. resids 0.001390 29.99247S.E. equation 0.003710 0.544936F-statistic 7.044848 804.6767Log likelihood 454.8094 -84.06138Akaike AIC -8.292766 1.686322Schwarz SC -8.118924 1.860164Mean dependent -0.002527 100.2406S.D. dependent 0.004293 3.698585Determinant resid covariance (dof adj.) 4.08E-06Determinant resid covariance 3.57E-06Log likelihood 370.8871Akaike information criterion -6.609021Schwarz criterion -6.261337再进行滞后阶数的确定:表4.4 最优滞后阶数的判断VAR Lag Order Selection CriteriaEndogenous variables: DLE PExogenous variables: CDate: 11/22/14 Time: 22:22Sample: 2005M07 2014M10Included observations: 99Lag LogL LR FPE AIC SC HQ0 134.7784 NA 0.000234 -2.682392 -2.629965 -2.6611801 302.5627 325.3999 8.57e-06 -5.991165 -5.833886 -5.9275302 329.0230 50.24783 5.45e-06 -6.444909 -6.182775* -6.338849*3 334.3733 9.943949* 5.30e-06* -6.472187* -6.105200 -6.3237044 337.4531 5.599742 5.40e-06 -6.453598 -5.981758 -6.2626915 339.7589 4.099176 5.60e-06 -6.419372 -5.842679 -6.1860416 345.0498 9.192324 5.46e-06 -6.445451 -5.763905 -6.1696967 345.5484 0.846076 5.87e-06 -6.374716 -5.588316 -6.0565378 346.7369 1.968760 6.23e-06 -6.317917 -5.426663 -5.9573149 352.5801 9.443639 6.01e-06 -6.355154 -5.359047 -5.95212810 353.7714 1.877082 6.39e-06 -6.298411 -5.197451 -5.85296111 354.3649 0.911279 6.87e-06 -6.229594 -5.023780 -5.74172012 356.4617 3.134644 7.18e-06 -6.191146 -4.880479 -5.660848* indicates lag order selected by the criterionLR: sequential modified LR test statistic (each test at 5% level)FPE: Final prediction errorAIC: Akaike information criterionSC: Schwarz information criterionHQ: Hannan-Quinn information criterion由上边可知,根据信息准则,采取少数服从多数原则,取滞后阶数为3,此外取滞后阶数为2(SC为 -6.182775)或取滞后阶数为3(SC为-6.105200)时,两者SC值相差不是很大。