《季节调整我国化和NBS--SA软件研发项目设计报告》由会员分享，可在线阅读，更多相关《季节调整我国化和NBS--SA软件研发项目设计报告（9页珍藏版）》请在金锄头文库上搜索。
1、 .wd.季节调整中国化与NBS-SA软件研发工程研究报告课题负责人：国家统计局马建堂；南开大学张晓峒一、 为什么对经济序列进展季节调整1 带有季节性变化因素的经济变量不能直接用来计算环比增长率在月度或季度的经济序列中常包含有季节性周期。尽管这是真实的观测数据，但是，却不能用这种带有季节性的数据直接计算环比增长率。2 用季节调整数据才能真正提醒经济变量的环比增长变化用季节调整方法对季度月度经济序列进展季节调整，从序列中剔除季节因素、日历效应等因素的影响，使数据处于一个“平等的水平，这种情形下，计算出的环比增长率才能真正反映经济变量的变化。3 用季节调整方法从序列中别离出的季节成分含有重要的经济
2、信息用季节调整方法从序列中别离出的季节成分序列含有重要的经济信息。这种变化特征为商业部门和物资供应部门提供了有价值信息。按季节变化规律准备充足的商品和物资供不同月份的消费，对保障供应、满足需求、繁荣经济有重要意义。二、 NBS-SA与X-13A-S的季节调整原理1 时间序列的四种成分一个经济时间序列以后简称时间序列通常受多种因素影响，一般地，可以把这些因素分解为趋势成分、循环成分、季节成分和不规那么成分。经典的时间序列模型有两种：1加法模型： Yt = Tt + Ct + St + It， At = Tt + Ct + It2乘法模型： Yt= Tt Ct St It， At = Tt Ct
3、It 2 X-13A-S程序的 基本流程X-13A-S程序处理过程可以分为“建模，“季节调整和“诊断三个阶段，如图3-1中3个虚线框所示。模型诊断M1M11，Q检验统计量，季节调整序列和不规那么分量的谱分析，平移区间、修正历史检验等残差序列进入X-11季节调整阶段建模与模型比较诊断regARIMA模型日历效应、移动假日、离群值回归，前向预测，后向预测，预调整对提取的4个信号分量的1阶矩、2阶矩检验。SA序列和不规那么分量的谱分析。平移区间、修正历史检验等残差序列进入SEATS季节调整阶段图3-1 X-13A-S季节调整过程三阶段第一阶段，即利用RegARIMA方法对原始序列建模阶段。本阶段的功
4、能是首先从原始序列中剔除各种日历效应和假日效应以及离群值的影响并对误差序列建设乘积季节模型，并对该序列进展前向、后向预测。第二阶段是对序列的季节调整阶段，即对上一阶段RegARIMA建模输出的序列进展季节调整。其中包括两种调整方法。一种是“X11方法；一种是“SEATS方法。第三阶段是对季节调整结果的诊断阶段，即通过一系列统计量的值检验，判断季节调整结果是否符合要求。假设用X11方法进展季节调整，对调整结果的诊断包括11个M统计量，2个Q统计量检验，季节调整后序列和不规那么成分的谱分析，以及平移区间检验和修正历史检验等。假设用SEATS方法进展季节调整，检验那么包括4个信号成分的一阶矩和二阶矩
5、检验，季节调整后序列和不规那么成分的谱分析，以及平移区间检验和修正历史检验等。下面对三个阶段分别予以介绍。3 RegARIMA建模原理对时间序列yt写出一个线性回归模型Yt= b0 +b1X1t+ b2X2t+ bkXkt + ut (3-1)其中yt是待调整时间序列，即前面所说的原始序列。Xit,i = 0, 1, , k，是解释变量，其中包括各种离群值变量，日历效应变量，以及假日效应变量等。i,i = 0, 1, , k是回归系数。ut是随机误差序列，应该满足回归模型设定的假定条件。 对于季节性经济时间序列来说，回归后得到的残差通常是自相关的。所谓RegARIMA建模方法就是在回归式3-1
6、的根基上同时对误差项ut建设季节时间序列模型。表达式如下，Fp(L) AP(Ls) DdDsDut = Qq(L) GQ(Ls) vt (3-2) 其中ut即是来自式3-1中的ut。D是一阶差分算子，d表示差分次数。Ds是s期的季节差分算子，D表示季节差分次数。DdDsDut表示对ut进展一阶差分d次，季节差分D次后为平稳序列。vt是白噪声序列。上式称作 (p, d, q) (P, D, Q)s 阶季节时间序列模型或乘积季节模型。 注意：如果d = D = 0即无需差分，通常将式3-1中的ut改用其对均值的离差代替，即用ut- 代替ut，其中 = E(ut)。将式3-1写成，ut = (Yt
7、0-)，并代入式3-2，可以写成，Fp(L) AP(Ls) DdDsD(Yt 0-) = Qq(L) GQ(Ls) vt (3-3)RegARIMA模型3-3可以理解为1在回归模型3-1的根基上，令误差项ut服从SARIMA模型3-2；2从yt中减去回归因子的影响后，建设SARIMA模型。把式3-3写成，DdDsD(Yt 0-) = (3-4)令wt =，wt被假定为平稳序列，并表达为Fp(L) AP(Ls)wt =Qq(L) GQ(Ls) vt。模型3-4式的另一种写法是(1-L)d(1- Ls)D Yt = (0 +) (1-L)d(1- Ls)D + wt (3-5)上式强调的是RegA
8、RIMA模型中的回归变量Xit,i = 0, 1, , k，同Yt一样，用ARIMA模型的差分算子(1-L)d(1- Ls)D进展了差分。式3-1中的回归因子解释变量Xit主要包括各种离群值outlier变量与日历相关的影响变量以及假日效应变量等。X-13A-S软件将这些功能集中到了第一阶段RegARIMA模块中。在X-13A-S的X11季节调整程序中，仍保存了在不规那么成分中对这些效应的处理功能。但X-13A-S软件说明书建议在第一阶段处理离群值效应、日历效应和假日效应等，对原始序列进展预调整。在X-13A-S程序中，离群值被分为5类，即加性离群值AO，水平漂移LS，暂时变化TC，斜线变动R
9、P和季节性离群值SO。分述如下，下面介绍7种可能在季节性经济序列中产生影响的效应，它们都属于日历效应。分述如下。1固定季节效应。2闰年效应。3月份长度效应。4季度长度效应。5交易日效应。6工作日效应。7移动假日效应。4 X11在默认状态下的计算原理以月度序列、加法模型为例，假设序列Yt不含离群值、已经通过前向预测、后向预测对序列进展了扩展，所以序列两端不需要使用修正公式。省略极端值extreme value调整过程。由于在X-11体系的季节调整程序中，对趋势成分和循环成分不再细分，视为一种成分，所以这里按照加法分解模型将Yt分解为趋势循环Ct、季节St和不规那么It三种成分：Yt = Ct +
10、 St + It ， 季节调整序列用At表示，那么Yt = At + St。X11季节调整计算过程也分为3个阶段。1初始估计At，2再次估计At，3估计最终Henderson趋势和不规那么成分。第1阶段：初始估计At1使用“中心化12项移动平均212估计趋势循环成分Ct。2估计季节不规那么成分： ( S+ I)= Yt -3对每个月份应用33移动平均估计预备季节成分：= M 334估计季节调整后序列A。A= ( C+ I)= Yt - 这是初次估计的季节调整后序列，所包含的季节性因素已经很少了。接下来X11方法将基于A序列继续进展移动平均。第2阶段：再次估计At1用13项Henderson移动
11、平均估计趋势循环成分Ct。C= H13 (A) 2估计季节不规那么成分。 ( S+ I)= Yt -3对每个月份应用35移动平均估计最终季节成分。= M 354估计季节调整后序列。= ( C+ I)=Yt -第3阶段：估计最终Henderson趋势和不规那么成分1估计最终趋势。C= H2H+1 (A)2估计最终不规那么成分。= - 最终得到分解序列的加法模型： Yt = +三、中国假日效应的处理1中国固定假日效应的处理 中国每周工作天数体制转换、固定假日、调休、黄金周等实施期间的沿革示意图见图4-1。每周6天工作制1995年5月以前每周5天工作制1995年5月始元旦法定假日1天，2000年1月
12、始伴以调休休3天五一节法定假日1天2000年5月以前五一劳动节法定假日3天2000年5月始五一节法定假日改回1天2008年5月始无五一黄金周2000年5月以前五一黄金周，伴以调休2000：52007：5取消五一黄金周2008年5月始取消国庆节法定假日1天1999年10月以前国庆节法定假日3天1999年10月始无十一黄金周1999年10月以前十一黄金周，伴以调休1999年10月始 1990:1 1995:5 1999:10 2008:5 2030:12图4-1 中国工作体制转换、固定假日、调休、黄金周等因素变化示意图关于中国月度交易日效应变量的定义与计算 下面以2009年9月为例介绍当出现跨月调
13、休时，中国交易日效应变量的计算过程。第1步：以公元纪年的日历为根基生成新的日历。由于2009年9月涉及到黄金周调休，所以要先调整9月份的日历。“国务院办公厅关于2009年局部节假日安排的通知见本研究报告后附录3规定，将2009年9月27日星期日的公休日调至10月7日星期三，而把10月7日星期三调至9月27日按工作日处理。所以，对于2009年9月27日来说，其交易日的性质由周日变成周三。这导致2009年9月的周日减少1天，周三增加1天。对于2009年10月7日来说，其交易日的性质由周三变成了周日，导致2009年10月的周日增加1天，周三减少1天。见图4-2。图4-2 按国务院放假通知生成2009
14、年9月和10月份新定义的工作日和休息日新日历示意图。第2步：计算中国月度交易日效应变量。以新生成的日历为根基，那么在计算2009年9月的交易日效应变量时，周日的天数就由原来的4天变为现在的3天，周三的天数由原来的5天变为现在的6天，而其他交易日的天数不变。假设用m12009,9，m62009,9，m72009,9分别表示2009年9月份调整后周一、周二、周日的天数，那么， m12009,9= 4 m22009,9= 5 m32009,9= 6 m42009,9= 4 m52009,9= 4 m62009,9= 4 m72009,9= 3那么中国交易日效应变量T1t，T2t，,T6t分别表示调整后的周一、周二、周六效应变量在2009年9月的取值依次是T12009,9 =
