时间序列论文

  • 格式:doc
  • 大小:280.50 KB
  • 文档页数:13

时间序列分析

课程论文

基于ARMAX模型的财政收入与税收

的时间序列分析与预测

班级:13级应用统计学1班

学号:0

姓名:***

基于ARMAX模型的财政收入与税收

的时间序列分析与预测

摘要

财政收入,是指政府为履行其职能、实施公共政策和提供公共物品与服务需要而筹集的一切资金的总和,是衡量一国政府财力的重要指标;其中税收收入是国家财政收入的重要组成部分,一般占到财政收入的90%以上,是政府机器的经济基础;

本文利用应用时间序列分析 的知识通过sas 统计软件对1978-2012年中国财政收入与税收数据进行分析,通过单位根检验,发现两者都是非平稳时间序列,并且存在协整关系,所以拟合了ARIMAX模型;由于残差序列非白噪声,所以对残差序列又进行了进一步的拟合,最后对模型进行预测,做出预测图;

关键词: 财政收入与税收 ARIMAX模型 预测

一、引言

财政与税收关系到国家发展、民生大计;财政收入与税收对社会资源配置、收入分配、国民经济发展、企业经济活动、居民切身利益及政府决策行为都有重大影响;近年来,随着我国经济的持续高速发展和国家财政与税收的大幅度增长,以及我国经济体制改革的不断深化和国家对经济发展宏观调控力度的不断加大,国家也适时出台了一系列有关财政与税收管理的新规定、新政策和新的监管制度;可以看出两者地位越来越重要,作用越来越明显;通过本文的分析,旨在找出两者的关系,为我国财政与税收做出合理的解释,为以后的收入做出合理的预测;

二、数据分析

一、序列平稳性检验

1、时序图:

y0100002000030000400005000060000700008000090000100000110000120000t1970198019902000201020201978-2012年中国财政收入与税收数据

图 1 原数据时序图

图1中,红色为y财政收入序列书序图;黑色为x税收收入序列时序图;从时序图中可以看出x序列、y序列均显著非平稳;并且两者都有明显的增加趋势;

2、单位根检验:

表 1 序列x的单位根检验

The ARIMA Procedure

Augmented Dickey-Fuller Unit Root Tests

Type Lags Rho PrF

Zero Mean 0

1

Sing Mean 0

1

Trend 0

1

表 2 序列y的单位根检验 Augmented Dickey-Fuller Unit Root Tests

Type Lags Rho PrF

Zero Mean 0

1

Sing Mean 0

1

Trend 0

1

单位根检验的原假设H0:序列为非平稳序列,如果 P> ,则接受原假设,认为序列非平稳,否者序列为平稳序列;上面的X、Y序列单位根检验过后,P值均大于,所以认为俩序列均为非平稳序列;

时序图显示两个序列具有某种同变关系,所以考虑建立ARIMAX模型;

二、协整检验

多元非平稳序列之间能否建立回归模型,关键在于他们之间是否具有协整关系;所以建模前必须进行协整检验EG检验;

假设条件:

上述假设条件等价于

对回归残差序列进行平稳性检验:

图 2 残差的自相关图

图3 残差的偏相关图

表 3 残差序列单位根检验

Type Lags Rho Pr>Rho Tau PrF

Zero Mean 0

1

2

Single

Mean 0

1

2

Trend 0

1

2

由自相关图与偏相关图可以看出回归残差序列为平稳序列;另一方面,单位根检验的原假设H0:回归残差序列为非平稳序列,如果 P> ,则接受原假设,认为序列非平稳,否者序列为平稳序列;残差序列单位根检验图显示在第一种类型常数均值、无趋势项的情况下,有%即100%的把握断定残差序列平稳且不具有自相关性;

也就是说,可以以%的把握认为1978-2012年中国财政收入序列y和税收收入序列x之间存在协整关系;因此可以对两者进行建模;

三、ARIMAX模型建模

1、建立响应序列与输入序列之间的回归模型:

图 4 序列y与序列x之间的相关图

图中第一列为因变量序列的延迟阶数,第二列为延迟序列与输入序列之间的协方差,第三列为相关系数,后面为相关图;相关图显示序列y在延迟阶数为零时与序列x相关关系最强;因此可以将序列y与序列x同期建模;

tXˆY1t 1

2、模型的参数估计:

表 4 极大似然估计的模型的参数估计

Maximum Likelihood Estimation

Parameter Estimate Standard

Error T value Approx

pr>|t| Lag variable shift

MU 0 Y 0

NUM1 <.0001 0 x 0

参数估计是对模型中位置参数的值的估计,以确定模型的口径,并对模型进行显著性检验;的本题采用极大似然方法进行估计;

参数估计表从左到右分别为 参数名称MU为常数项,NUM1为θ1、参数估计值、估计值标准差、t检验值、t统计量p值、延迟阶数;参数估计结果显示,两个参数均显著p值小于;所以建立ARIMAX模型如下:

ttXY143.1 2

3、回归残差序列白噪声检验:

表 5

残差序列白噪声检验

To

Lag Chi-

Square DF Pr>ChiSq -------Autocorrelation-------

6 6 <.0001 12 12 18 18

24 24 <.0001

上表结果显示,延迟各阶的LB统计量的P值都小于,可以认为残差序列为非白噪声序列;说明残差序列中还有未提取到的因素,所以考虑对残差序列进行拟合;

4、回归残差序列的模型拟合:

表 6 残差序列参数估计

Conditional Least Squares Estimation

Parameter Estimate Standard Error t Value Approx Pr>|t|

Lag

MA1,1

1

MA1,2 2

AR1,1 1

AR1,2 2

AR1,3 3

由于残差序列为平稳序列,所以考虑通过ARMA模型对其进行的拟合;除去常数项,参数估计结果如下,在最小二乘估计下,AR1,1、MA1,1显著t统计量的p值小于;

5、残差序列模型的显著性检验:

表 7 残差自相关检验

Autocorrelation Check of Residuals

To

Lag Chi- DF Pr>ChiSq ------Autocorrelation----- Square

6 1 12 7 18 13 24 19

由上表结果显示,延迟各阶的LB统计量的p值均显著大于αα=,所以认为该拟合模型显著成立;

6、残差序列模型的拟合:

表 8 残差序列模型拟合

Autoregressive Factors

Factor 1:12+3

Moving Autoregressive Factors

Factor 1:1+2

输出形式等价于

132)58262.015117.035904.11(ttBBB 3

或记为

432158262.015117.035904.1ttttt 4

7、模型的最终形式:

432158262.015117.035904.1143.1tttttttXY 5

三、模型预测 Forecast for y-100000100002000030000400005000060000700008000090000100000110000120000t1970198019902000201020201978-2012年中国财政收入与税收数据

图 5 模型预测图

四、结论

一、时间序列分析中,尽管两组数据都是非平稳数据,但是只要他们之间存在协整关系,就可以拟合两者的关系ARIMAX模型,通过模型可以直接看出两者的函数关系;

二、经过上述分析,最后得到ARIMAX模型如下:

432158262.015117.035904.1143.1tttttttXY

根据模型,可以直观的看出财政收入与税收收入两者为正相关关系,中国财政收入呈现出快速递增趋势;表示序列的残差项;单位根检验显示,残差序列为平稳非白噪声序列,因此通过ARMA模型对残差序列又重新进行了拟合;第二个函数就是拟合出的残差函数; 三、通过本文的分析,找到了中国财政收入与税收收入的确切关系,由预测图可以观察到中国财政收入呈显著上升趋势;通过时间序列模型,预测出了近几年国家的财政收入情况;本次分析可以为政府提供出具政策的依据;