时间序列分析ARMA模型实验

基于ARMA模型的社会融资规模增长分析————ARMA模型实验-2-第一部分实验分析目的及方法一般说来,若时间序列满足平稳随机过程的性质,则可用经典的ARMA模型进行建模和预则。但是,由于金融时间序列随机波动较大,很少满足ARMA模型的适用条件,无法直接采用该模型进行处理。通过对数化及差分处理后,将原本非平稳的序列处理为近似平稳的序列,可以采用ARMA模型进行建模和分析。第二部分实验数据2.1数据来源数据来源于中经网统计数据库。具体数据见附录表5.1。2.2所选数据变量社会融资规模指一定时期内（每月、每季或每年）实体经济从金融体系获得的全部资金总额，为一增量概念，即期末余额减去期初余额的差额，或当期发行或发生额扣除当期兑付或偿还额的差额。社会融资规模作为重要的宏观监测指标，由实体经济需求所决定，反映金融体系对实体经济的资金量支持。本实验拟选取2005年11月到2014年9月我国以月为单位的社会融资规模的数据来构建ARMA模型，并利用该模型进行分析预测。第三部分ARMA模型构建3.1判断序列的平稳性首先绘制出M的折线图，结果如下图：-3-图3.1社会融资规模M曲线图从图中可以看出，社会融资规模M序列具有一定的趋势性，由此可以初步判断该序列是非平稳的。此外，m在每年同时期出现相同的变动趋势，表明m还存在季节特征。下面对m的平稳性和季节性·进行进一步检验。为了减少m的变动趋势以及异方差性，先对m进行对数化处理，记为lm，其时序图如下：图3.2lm曲线图-4-对数化后的趋势性减弱，但仍存在一定的趋势性，下面观察lm的自相关图表3.1lm的自相关图上表可以看出，该lm序列的PACF只在滞后一期、二期和三期是显著的，ACF随着滞后结束的增加慢慢衰减至0，由此可以看出该序列表现出一定的平稳性。进一步进行单位根检验，由于存在较弱的趋势性且均值不为零，选择存在趋势项的形式，并根据AIC自动选择之后结束，单位根检验结果如下：表3.2单位根输出结果NullHypothesis:LMhasaunitrootExogenous:Constant,LinearTrendLagLength:0(Automatic-basedonSIC,maxlag=12)t-StatisticProb.*-5-AugmentedDickey-Fullerteststatistic-8.6746460.0000Testcriticalvalues:1%level-4.0469255%level-3.45276410%level-3.151911*MacKinnon(1996)one-sidedp-values.单位根统计量ADF=-8.674646小于临界值，且P为0.0000，因此该序列不存在单位根，即该序列是平稳序列。由于趋势性会掩盖季节性，从lm图中可以看出，该序列有一定的季节性，为了分析季节性，对lm进行差分处理，进一步观察季节性：图3.3dlm曲线图观察dlm的自相关表：表3.3dlm的自相关图Date:11/02/14Time:22:35Sample:2005M112014M09Includedobservations:106AutocorrelationPartialCorrelationACPACQ-StatProb****|.|****|.|1-0.566-0.56634.9340.000.|*|**|.|20.113-0.30536.3410.000.|.|*|.|30.032-0.09336.4550.000*|.|*|.|4-0.084-0.11437.2440.000-6-.|*|.|.|50.1050.01538.4940.000*|.|*|.|6-0.182-0.18242.2960.000.|*|*|.|70.105-0.15643.5630.000.|.|*|.|8-0.058-0.17143.9540.000.|.|*|.|9-0.019-0.19643.9960.000.|*|.|.|100.110-0.04545.4290.000**|.|**|.|11-0.242-0.32952.5010.000.|***|.|.|120.3630.02368.5160.000*|.|.|.|13-0.2020.03273.5340.000.|*|.|*|140.1010.12574.8150.000.|.|.|*|150.0040.14174.8170.000*|.|*|.|16-0.161-0.08978.1100.000.|**|.|.|170.2190.03784.2520.000**|.|.|.|18-0.221-0.03690.6230.000.|*|.|.|190.089-0.04691.6620.000*|.|*|.|20-0.080-0.15892.5160.000.|.|.|.|210.067-0.03993.1150.000.|.|.|.|220.0680.05693.7490.000**|.|*|.|23-0.231-0.130101.080.000.|***|.|*|240.3590.116119.040.000*|.|.|*|25-0.1890.123124.090.000.|.|.|.|260.0320.034124.230.000.|.|.|.|270.0590.037124.740.000*|.|.|.|28-0.1260.044127.080.000.|*|*|.|290.087-0.079128.210.000.|.|.|*|30-0.0500.092128.580.000.|.|.|.|31-0.037-0.019128.790.000.|.|*|.|32-0.035-0.113128.970.000.|.|.|.|330.041-0.056129.240.000.|*|.|.|340.078-0.027130.210.000**|.|*|.|35-0.215-0.197137.640.000.|***|.|*|360.3800.130161.260.000由dlm的自相关图可知，dlm在滞后期为12、24、36等差的自相关系数均显著异于零。因此该序列为以12为周期呈现季节性，而且季节自相关系数并没有衰减至零，因此为了考虑这种季节性，进行季节性差分，得新变量sdlm：观察sdlm的自相关图：表3.4sdlm的自相关图Date:11/02/14Time:22:40Sample:2005M112014M09Includedobservations:94-7-AutocorrelationPartialCorrelationACPACQ-StatProb****|.|****|.|1-0.505-0.50524.7670.000.|.|***|.|2-0.057-0.41925.0820.000.|.|**|.|30.073-0.29225.6090.000.|*|.|.|40.1600.06728.1690.000**|.|.*|.|5-0.264-0.12535.2520.000.|*|.*|.|60.098-0.11036.2440.000.|*|.|.|70.0980.01937.2430.000.|.|.|*|8-0.0410.08237.4190.000.*|.|.|.|9-0.132-0.03839.2750.000.|*|.*|.|100.076-0.13939.9020.000.|**|.|**|110.2270.24745.4850.000***|.|**|.|12-0.459-0.25968.6470.000.|*|**|.|130.193-0.25172.7770.000.|*|.*|.|140.132-0.10174.7530.000.*|.|.*|.|15-0.142-0.18977.0560.000.|.|.|.|16-0.053-0.05677.3780.000.|**|.|*|170.2330.09183.7510.000**|.|.*|.|18-0.234-0.17990.2580.000.|*|.|.|190.1020.05491.5050.000.|.|.|.|20-0.052-0.03591.8410.000.|*|.|.|210.123-0.00993.7140.000.|.|.|*|22-0.0590.12094.1500.000.|.|.|**|23-0.0110.21594.1660.000.|.|.*|.|24-0.032-0.17094.3010.000.|*|.*|.|250.088-0.13795.3030.000.*|.|.|.|26-0.105-0.03496.7600.000.|*|.*|.|270.077-0.11697.5620.000.|.|.*|.|28-0.054-0.17897.9670.000.|.|.|.|290.0100.03297.9820.000.|*|.|.|300.1020.03999.4570.000.*|.|.*|.|31-0.179-0.099104.060.000.|.|.|.|320.071-0.058104.790.000.|.|.*|.|330.031-0.066104.930.000.*|.|.*|.|34-0.089-0.144106.130.000.|.|.|*|350.0360.082106.320.000.|*|.*|.|360.105-0.102108.050.000Sdlm在滞后期24之后的季节ACF和PACF已衰减至零，下面对sdlm建立SARMA模型。3.2模型参数识别-8-由表3.4sdlm的自相关图的自相关图可知，偏自相关系数在3阶后都落在两倍标准差的范围以内，即不显著异于零。自相关系数在1阶和12阶显著异于零。因此SARMA(p,q)模型中选择p、q均不超过3。此外，由于高阶移动平均模型估计较为困难而且自回归模型可以表示无穷阶的移动平均过程，因此Q尽可能取小。拟选择SARMA(1,0)（1,0）12、SARMA(1,0)（1,1）12、SARMA(1,1)（1,0）12、SARMA(1,1)（1,1）12、SARMA(2,0)（1,0）12、SARMA(2,0)（1,1）12、SARMA(3,0)（1,0）12、SARMA(3,0)（1,1）12八个模型来拟合sdlnm。3.3模型参数估计以SARMA(1,0)（1,0）12模型为例，分析该模型的估计及残差的检验，其他模型类似。回归结果为：表3.5SARMA(1,0)（1,0）12模型估计结果DependentVariable:SDLMMethod:LeastSquaresDate:11/02/14Time:22:50Sample(adjusted):2008M012014M09Includedobservations:81afteradjustmentsConvergenceachievedafter6iterationsVariableCoefficientStd.Errort-StatisticProb.C-0.0053050.023352-0.2271650.8209AR(1)-0.4908550.098580-4.9792560.0000SAR(12)-0.5485090.096987-5.6554710.0000R-squared0.448053Meandependentvar-0.004983AdjustedR-squared0.433901S.D.dependentvar0.644876S.E.ofregression0.485202Akaikeinfocriterion1.427829Sumsquaredresid18.36280Schwarzcriterion1.516512Loglikelihood-54.82707Hannan-Quinncriter.1.463410F-statistic31.65