1、季節調整中國化與NBS-SA軟件研發項目研究報告課題負責人：國家統計局馬建堂；南 開大學 張曉峭一、為什么對經濟序列進行季節調整(i)帶有季節性變化因素的經濟變量不能直接用來計算環比增長率在月度或季度的經濟序列中常包含有季節性周期。盡管這是真實的觀測數據，但是，卻不能用這種帶有季節性的數據直接計算環比增長率。(2)用季節調整數據才能真正揭示經濟變量的環比增長變化用季節調整方法對季度(月度)經濟序列進行季節調整，從序列中剔除季節因素、日歷 效應等因素的影響， 使數據處于一個“平等”的水平，這種情形下，計算出的環比增長率才能真正反映經濟變量的變化。(3)用季節調整方法從序列中分離出的季節成分含有重

2、要的經濟信息用季節調整方法從序列中分離出的季節成分序列含有重要的經濟信息。這種變化特征為商業部門和物資供應部門提供了有價值信息。按季節變化規律準備充足的商品和物資供不同月份的消費，對保障供給、滿足需求、繁榮經濟有重要意義。NBS-SA與X-13A-S的季節調整原理(1)時間序列的四種成分一個經濟時間序列(以后簡稱時間序列)通常受多種因素影響，一般地，可以把這些因 素分解為趨勢成分、循環成分、季節成分和不規則成分。經典的時間序列模型有兩種：(1)加法模型：Yt =1+ Ct+ St +It, At = Tt + Ct +It(2)乘法模型：Yt =下 x Ct x Stx It, A = TtX

3、Ctx ItX-13A-S程序的基本流程X-13A-S程序處理過程可以分為“建模”，“季節調整”和“診斷”三個階段， 如圖3-1中3個虛線框所示。圖3-1 X-13A-S 季節調整過程三階段第一階段，即利用 RegARIMAT法對原始序列建模階段。本階段的功能是首先從原始序 列中剔除各種日歷效應和假日效應以及離群值的影響并對誤差序列建立乘積季節模型，并對該序列進行前向、后向預測。第二階段是對序列的季節調整階段，即對上一階段(RegARIMA建模)輸出的序列進行季節調整。其中包括兩種調整方法。一種是“ X11 方法”；一種是“ SEATS方法。第三階段是對季節調整結果的診斷階段，即通過一系列統計

4、量的值檢驗，判斷季節調整結果是否符合要求。若用X11方法進行季節調整，對調整結果的診斷包括11個M統計量，2個Q統計量檢驗，季節調整后序列和不規則成分的譜分析，以及平移區間檢驗和修正歷史檢驗等。若用 SEATS方法進行季節調整，檢驗則包括4個信號成分的一階矩和二階矩檢驗，季節調整后序列和不規則成分的譜分析，以及平移區間檢驗和修正歷史檢驗等。下面對三個階段分別予以介紹。RegARIMA建模原理對時間序列yt寫出一個線性回歸模型Yt = 0+iXit + 2X3 + k Xkt + ut(3-1)其中yt是待調整時間序列，即前面所說的原始序列。Xt, i = 0, 1,，k,是解釋變量，其中包括各

5、種離群值變量，日歷效應變量，以及假日效應變量等。6, i = 0, 1,，k是回歸系數。ut是隨機誤差序列，應該滿足回歸模型設定的假定條件。對于季節性經濟時間序列來說，回歸后得到的殘差通常是自相關的。所謂RegARIMAt模方法就是在回歸式 (3-1 )的基礎上同時對誤差項 ut建立季節時間序列模型。表達式如下，p(L) P(Ls) dsDut = q(L) q(Ls) vt(3-2)其中ut即是來自式(3-1 )中的ut。是一階差分算子，d表示差分次數。s是s期的季節 差分算子，D表示季節差分次數。dsDut表示對ut進行一階差分d次，季節差分D次后為平穩 序列。vt是白噪聲序列。上式稱作(

6、p,d,q) (P,DQ)s階季節時間序列模型或乘積季節模型。注意：如果d = D = 0 (即無需差分)，通常將式(3-1)中的ut改用其對均值的離差代 替，即用ut - w代替ut,其中w = E(ut) 0將式(3-1 )寫成，k Ut = (Yt - Bo-iXit ),i 1并代入式(3-2 ),可以寫成，ks d Dsp(L) p(L) s(Y 3 0-iXit)= q(L)L) Vt(3-3)i 1RegARIMA莫型(3-3)可以理解為(1)在回歸模型(3-1 )的基礎上，令誤差項 ut服從 SARIMA莫型(3-2 ); (2)從yt中減去回歸因子的影響后，建立SARIMA莫

7、型。把式(3-3 )寫成，dsD(Ytk0-ii Xit )=1sQq(L) Fq(L ) -vtp(L)牛(L )(3-4)s e)q(L) Fq(L ) -vt p(L) Ap(L )w被假定為平穩序列，并表達為p( L) p( Ls) w =q( L) U Ls) Vt。模型(3-4)式的另一種寫法是k. d . s D - _d . s D(1-L) (1- L) Yt = ( 3o + iXit) (1 -L) (1- L) + w(3-5)i 1上式強調的是 RegARIMA莫型中的回歸變量 Xt, i = 0, 1,，k,同Y一樣，用ARIMA模 型的差分算子(1 -L )d(1

8、 - L s) D進行了差分。式(3-1 )中的回歸因子(解釋變量)Xt主要包括各種離群值(outlier )變量與日歷相關的影響變量以及假日效應變量等。X-13A-S軟件將這些功能集中到了第一階段 RegARIMA莫塊中。在 X-13A-S的X11季節 調整程序中，仍保留了在不規則成分中對這些效應的處理功能。但X-13A-S軟件說明書建議在第一階段處理離群值效應、日歷效應和假日效應等，對原始序列進行預調整。在X-13A-S程序中，離群值被分為5類，即加性離群值(AQ,水平漂移(LS),暫時變化(TC),斜線變動(RP)和季節性離群值 SQ分述如下，下面介紹7種可能在季節性經濟序列中產生影響的

9、效應，它們都屬于日歷效應。 分述如下。(1)固定季節效應。(2)閏年效應。(3)月份長度效應。(4)季度長度效應。(5)交易日效應。(6)工作日效應。(7)移動假日效應。X11在默認狀態下的計算原理以月度序列、加法模型為例，假設序列Y不含離群值、已經通過前向預測、后向預測對序列進行了擴展，所以序列兩端不需要使用修正公式。省略極端值( extreme value )調 整過程。由于在 X-11體系的季節調整程序中，對趨勢成分和循環成分不再細分，視為一種 成分，所以這里按照加法分解模型將Y分解為趨勢循環(C)、季節(S)和不規則(It)三種成分：Y = C + St + It,季節調整序列用 A表

10、示，則Y = A + St。X11季節調整計算過程也分為 3個階段。(1) 初始估計A, (2)再次估計 A, (3)估計最終Henderson趨勢和不規則成分。第1階段：初始估計A.使用“中心化12項”移動平均(2X12)估計趨勢循環成分Co.估計季節-不規則成分：(S+ I ) 1)= Yt - Ct(1).對每個月份應用 3X 3移動平均估計預備季節成分：S(1)= M3x3 (S I )(1).估計季節調整后序列A(1)。A(1) = ( C+ I ) t(1)= Y - S這是初次估計的季節調整后序列，所包含的季節性因素已經很少了。接下來X11方法將基于At序列繼續進行移動平均。第2

11、階段：再次估計A.用13項Henderson移動平均估計趨勢循環成分C。C (2)= H3 ( A(1).估計季節-不規則成分。(S+ I )(2)= Yt - Ct.對每個月份應用 3X5移動平均估計最終季節成分。g(2)= M3y (S I)(2).估計季節調整后序列At(2)。At(2)= ( C+ I ) (2)=Y - S(2)3階段：估計最終 Henderson趨勢和不規則成分.估計最終趨勢。C(3) = Hzh+i ( A(2).估計最終不規則成分。It(3)= At(2) - Ct最終得到分解序列的加法模型：Yt = Ct+ St+ It三、中國假日效應的處理中國固定假日效應的

12、處理1990:11995:51999:10中國每周工作天數體制轉換、固定假日、調休、黃金周等實施期間的沿革示意圖見圖4-1。每周6天工作制( 1995年5月以前)每周5天工作制(1995年5月始)元旦法定假日1天，2000年1月始伴以調休(休 3天)五一節法定假日1天 (2000年5月以前)五一勞動節法定假日3天(2000年5月始)五一節法定假日改回1天(2008年5月始)無五一黃金周(2000年5月以前)五一黃金周，伴以調休(2000: 52007: 5)取消五一黃金周 (2008年5月始取消)國慶節法定假日1天 (1999年10月以前)國慶節法定假日3天(1999年10月始)無H一黃金周(

13、1999年10月以前)十一黃金周，伴以調休 (1999年10月始)2008:52030:12圖4-1中國工作體制轉換、固定假日、調休、黃金周等因素變化示意圖關于中國月度交易日效應變量的定義與計算下面以2009年9月為例介紹當出現跨月調休時，中國交易日效應變量的計算過程。第1步：以公元紀年的日歷為基礎生成新的日歷。由于2009年9月涉及到黃金周調休，所以要先調整9月份的日歷。“國務院辦公廳關于2009年部分節假日安排的通知”（見本研究報告后附錄3）規定，將2009年9月27日（星期日）的公休日調至 10月7日（星期三），而把10月7日（星期三）調至 9月27日按 工作日處理。所以，對于 2009

14、年9月27日來說，其交易日的性質由周日變成周三。這導致 2009年9月的周日減少1天，周三增加1天。對于2009年10月7日來說，其交易日的性 質由周三變成了周日，導致2009年10月的周日增加1天，周三減少1天。見圖4-2。四大1-4-二十四中元節卜立十七89101112二一毛津.氧師節廿三廿四16117址廿廿九三+I舊22232425初四秋分初亢和七期晨2307 白114 g2128乳子一K公元攔9 ,月a13 廿五20選今天6519 忙2C 邕有廿五 如 初三27 貌工S25)9 后由F-四5六123十四料麻7 ,3g10十九囂廿一廿二1415n17廿六tH:廿N耽21222224就四幀

15、硒喊言-282930314-十二十二一四今天*十六 11廿三九月2tm農田己丑（牛，瑋圖4-2按國務院放假通知生成2009年9月和10月份新定義的工作日和休息日（新日歷）示意圖第2步：計算中國月度交易日效應變量。以新生成的日歷為基礎，則在計算2009年9月的交易日效應變量時，周日的天數就由原來的4天變為現在的3天，周三的天數由原來的5天變為現在的6天，而其他交易日的天數不變。若用 IT1 2009,9 , , 1162009,9 , 1772009,9分別表示2009年9月份調整后周一、周二、周日的天數，則，T1 2009,9=4T22009,9=5T32009,9=6T42009,9=4T5

16、2009,9=4T62009,9=4T7 2009,9=3則中國交易日效應變量 T1t, T2t,，T6 （分別表示調整后的周一、周二、周六效應變 量）在2009年9月的取值依次是T1 2009,9=T12009,9 772009,9=4 -3 :=1T22009,9 =二 T22009,9 772009,9 =5 -3 =:2T32009,9=T32009,9772009,9 =6 -3 =:3T42009,9=T42009,9772009,9 :=4 -3 ;=1T52009,9 =T52009,9 772009,9 =:4 -3 =：1T62009,9 =二 T62009,9 77200

17、9,9 =:4 -3 =：11990年1月至2030年12月中國交易日效應變量T1t, T2t,，T6t的賦值結果見本研究報告后附錄4。關于中國月度工作日效應變量的定義與計算中國月度工作日效應變量 WD勺定義方法是，以公元紀年日歷為基礎，根據我國固定假 日，國務院通知的黃金周調休方案重新計算每月中工作日和休息日的天數。由于1990年1月至1995年4月我國實行的是6天工作制，所以，中國月度工作日效應 變量WD勺值按下式計算WD t 二（調整后的每月工作日天數）-6 （調整后的每月周日天數）由于1995年5月至2030年12月我國實行的是 5天工作制，所以，中國月度工作日效 應變量WD勺值按下式

18、計算WD t 二（調整后的每月工作日天數）-（5/2）（調整后的每月周六、周日天數）仍以2009年9月月份為例，說明具體計算過程。根據公元紀年的日歷，2009年9月共有22個工作日，8個休息日。根據國務院調休方案，9月27日（周日）調到10月7日位置形成黃金周，同時把 10月7日（周三）的工作日調到 9月27日位置按工作日處理。于是， 經調整后，2009年9月份的工作日是23個（比原日歷增加一天），休息日7個（比原日歷 減少一天）。則2009年9月份 WD。2的值是WD09,9 = 23 - （5/2） 7 =1990 年1月至2030年12月中國月度工作日效應變量 WD勺值見研究報告后附錄

19、4的最 后一列。中國移動假日效應的處理中國移動假日春節、中秋節和端午節在陽歷中的移動范圍示意圖見圖4-10。1月21日2月20日端午節5月28日6月24日中秋節9月8日10月6日1990-1-12030-12-31 圖4-10 中國移動假日，春節、端午節、中秋節的移動范圍示意圖春節對經濟序列的影響方式是不一樣的。我們為春節效應變量設計了4種賦權方式，并對春節效應變量采用多段定義方式。4種賦權方式分別是在春節影響到的每天里（1）按等值賦權，（2）按”字型賦權，（3）按 7 字型賦權和（4）按“MT字型賦權。下面以月度序列為例，按字型賦權方式，仍設定春節對節前的影響天數為20天，節中節后的影響天數

20、為 15天，介紹計算春節效應變量的方法。字型賦權春節效應變量應該定義兩個，即節前效應變量spr1 t,和節后效應變量spr2t。如果按字型賦權，那么春節前20天至春節前1天，每天對經濟變量的影響權數應該依次是1/20 , 2/20 ,，20/20。春節向前大于 20天，則影響權數為零。春節日至春 節后15天，每天對經濟變量的影響權數應該依次是15/15 ,14/15,，1/15。春節向后大于15天，影響權數為零。如圖 4-15所示，圖中t=0處的垂線代表春節日，16-24-20 -16 -12-6-44 S 1 Z 16圖4-15春節對經濟序列的影響每天按字型賦權sprl2010,17ii 1

21、spr12010,2i 820 i 12020i20i 182/210 0.8666以2010年為例，春節在 2月14日。因為設定春節對節前的影響天數為20天，所以,前20天落在1月的天數為7天，落在2月的天數是13天，因此sprl t在1月的值是20i 28/210 0.1333i 1在2月的值是,因為春節節前效應不涉及對3至12月的影響，所以spr1 t在2010年3至12月的值應該都是零。則spr1 t在2010年1至12月的取值是spr1 2010,12010,12 = , , 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)因為春節節后效應設定的是 15天(從2月14至2月

22、28日)，全部落在2月份，所以 對應 2010年 2 月，spr2 2010,2 =1 ；對應 2010 年其它月份，spr2 2010.=0 , t = 1 , 3,，12。 則spr2 t在2010年1至12月的取值是spr2 2010,12010,12 = (0, 1,0, 0, 0, 0, 0, 0, 0, 0, 0, 0)同理，依據日歷，可計算出設定春節對節前的影響天數為20天，節中節后的影響天數為15天，按字型方式賦權的春節節前效應變量spr1 t和春節節后效應變量spr2 t從1990年1月至2030年12月的值。對序列 spr1 t和spr2 t中的每一相同月份的值分別做中心化

23、處 理，X 1t, j = spr1 t j - spr1.,j , j =1 , 2,12X 2t, j = spr2 t j - spr2.,j , j =1 , 2,，121990年1月至2030年12月的X1t, j 后附表5。X2t, j序列圖見圖4-16和圖4-17。數值見本報告199019952000200520102015202020252030圖5-3 NBS-SA軟件中的假日效應選項功能圖5-1 NBS-SA軟件主視窗圖4-16 ”字型權的春節節前效應變量Xlt，j圖4- 17字型權的春節節后效應變量X2t, j仿照上述定義方式，就可以計算出“ V字型權春節節前、節后效應變

24、量（2個）和“MT型賦權春節節前、節后效應變量（4個）。以月度序列為例，設定春節對節前的影響天數為20天，節中節后的影響天數為15天，按 7 字型和“ M字型賦權方式計算的6個春節效應變量（XVIt, j、X2t, j、XMt, j、XMt, j XM3t,j、XM1, j）的值見本報告后附表 5。關于處理中秋節效應，我們定義了節前、節中后 2個變量。節前最大影響天數設定為10天，節中和節后效應合并，最大影響天數設定為10天。也就是說，季節調整軟件 NBS-SA對中秋節節前效應白處理可以在110天中任意選擇。對中秋節節中后效應的處理可以在110天中任意選擇。對中秋節效應變量的賦權方式設計了等權

25、、V字型權和”字型權3種。關于端午節對經濟變量的影響，我們定義了節前效應變量。并設定端午節對節前的最大影響天數是10天。也就是說，季節調整軟件 NBS-SA中對端午節節前效應的影響可以在110天中任意選擇。對端午節效應變量的賦權方式設計了等值賦權一種。四、研制季節調整軟件NBS-SA國家統計局版季節調整軟件NBS-SA的研制是本項目的創新之一。季節調整軟件 NBS-SA研制成功為國家統計局確立季節調整數據和經濟指標的環比增長率發布體系奠定了堅實的 基礎。季節調整軟件 NBS-SA是在美國普查局季節調整軟件X-13A-S基礎上開發的。X-13A-S核心計算部分未作改動，只是在中國工作日效應和交易

26、日效應的處理方面嵌入了結合中國國 情所特有的一些節假日效應。所有X-13A-S軟件原有的計算功能依然可以運用程序語句在NBS-SA軟件中實現運算。NBS-SA與X-13A-S相比，我們主要做了如下創新：1. 季節調整軟件 NBS-SA在X-13A-S基礎上對視窗做了全面漢化。NBS-SA主視窗見圖5-1.在使用NBS-SA調整經濟序列之前， 增加了顯示經濟序列折線圖的畫圖功能。示意圖見圖5-2 。圖5-2 NBS-SA軟件中的畫圖功能.在保留X-13A-S軟件的交易日效應和工作日效應的同時，我們增加了中國交易日效 應和中國工作日效應兩個選項。關于移動假日效應的處理，增加了調整中國移動假日效應的

27、選項功能。選項視窗見圖5-3上部。關于春節的賦權方式有四種，即“均勻”、“型”、“ 型”和“M型”。 TOC o 1-5 h z 注意：NBS-SA軟件中，春節前的最大選擇天數是45天，即節前影響天數是在145天之間任意選擇。節后最大選擇天數是30天，即節后影響天數是在130天之間任意選擇。教提|咫喝整|覆調整史件|應明理期 回（變量|至節調生|審舂可審覽名汽推廣”柜桃pi-iF中取書G均加廣汽瑾rkF -I -r朗中節 &均為這是一個存量交易日I爻舄H熏匝-中怩恒縣白的應工也日r工作日效匝中國工作日刃（友R案的其他交量r皿檢利比如，春節對 GDPT值序列的影響應該是 V字型的。春節對鐵路客流

28、量的影響應該 是“MT字型的。關于中秋節效應，我們設計了三種賦權類型，即“均勻”、“型”、“V 型”。端午節定義了 1個效應變量，即節前效應變量。端午節效應變量是按等權設計的。怎樣選擇端午節對應的 2個天數選擇框，參見本研究報告第7章，國家統計局版季節調整軟件NBS-SA使用說明書。在圖5-3的下部“其它變量”的選擇框處可以輸入如季節變量、離群值變量、閏年效應變量、季度長度變量、月度長度變量等，用以消除這些因素的影響。.能輸出一個中文的“季節調整原始登記表”。這對于建立完整的季節調整檔案，保 留季節調整過程的原始信息有重要意義。.調整結束后能夠輸出一個 EXCEL格式的“序列與增長率分析表”。

29、其中包括原始序列、季節調整序列、季節因子分量、趨勢循環分量、不規則分量、同比增長率（%、環比增長率（為、折年率（%計算結果以及樣本外一年的原始序列和季節因子分量的預測值（即 月度數據預測12個月的值，季度數據預測 4個季度的值），見圖5-5。r-Pl, D % 4 , =nront- Mtnraoft Lncd- l 貨、而&0汨注式啕日冀13X_AB I CDEFGdI j rtffli聰如序列量節調IE序列年節Ff子片國司七一洶-卷分雪1和/1法司比Ci 產計中 葉王府店IQS ZOD9H119T90.SS630L117av (H,3B9版 1而5艮”修41 0. 9SEJ1851 20.

30、BO398 1.0S5B19 13 厘也小IC 9iOTSS1. It 176Se：436EB. 69S343 L. 00117i3 IS. 0C32： 0. 0ETS64 1. lTi64511 ：| 200M0110T5G.6 9T41.1. 054=76246637S5. ：-50SSl ,沔泥帕 1, 4的的 1. 052K.1 12.如7白 1111 2OD5S029323, ；Q 942. Ou93T5Ol 9935906009G3B. 394 2B 0. 9R2722111.55194b -2. 56063 -26.749112 2O0(r10S9317.6 9810fcS05?

31、6S0.加蚯附乳BB14. 074OSB I. RJ也07 14. TO9V7 3. 41 60.2170311 S 2O0&HX9343.2982,64989氏 9S5WF9*I849984. 4S5772 Q. BM81B1 14. 7&31321.89013114 25蠲05IQQE乩 4 10153. 471526 906251Z3?10163. 37654 09 9W635& 12 2ZW51 1, FT 39 23,45S4B115 2OQMQ6后 LQ31SL 595M4 96Z331473LOS工工麗$ 1. 0W434S l&gELS 1. B713S7 2i.9S6I二口

32、”口回l肝9S36.51鞏 bisasTTza1C5Z. .33E Q.把箍T3 1E. ltbOL 1. bE7M:G 2Z. ZJ631iit hwMoe10115. C LOTOS. 512&2il 943031=：107lC. 9752 L OOD3424 15.57347 1-7313&7 22.874031TB 2C0M0910912. S LuaG8.0391. 003M9717108S. 3912B 0. 9M9T32 15. 52215 L 5091/ 19. &911911B SODOiriO11T17.fi 11025. 674 38L 060M71DSU020.ESBE

33、B L 00046SB tfi. 2140 1.44dm 18. BUSS3O0W11ns-i n”已。1日境J. 014733799lrMCi9225 15, &12S 1.82729&12112filC HSi7,747OT1. 112&325T311302. 6410& 1.003990S IT, 5374 1. 572&4B 2cL 卬打口21Z2 2)010110112512109- 3I 210*02112619亂2124 MUMOa11252時,9123 ZO1M04112T193.5Fig BOIOBOB1盤制鮑鵬127 301 OtM1期3輒312B 201 MOT1207294. 4129 301 MOS1233194.6ILW 3C1WIT