季節(jié)性時間序列模型

1、上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,1,非平穩(wěn)和季節(jié)時間序列模型分析方法,在第四章中，我們介紹了非平穩(wěn)時間序列模型，但是在前面的討論中，對于時間序列的特性分析，以及模型的統(tǒng)計分析都集中于平穩(wěn)時間序列問題上。本章將介紹幾個非平穩(wěn)時間序列的建模方法，并且分析不同的非平穩(wěn)時間序列模型的動態(tài)性質(zhì),上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,2,8.1 arima模型的分析方法,8.1.1 arima模型的結(jié)構(gòu) 具有如下結(jié)構(gòu)的模型稱為求和自回歸移動平均(autoregressive integrated moving average)，簡記為arima(p,d,q)模型： (8.1) 式中,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,3

2、,式(8.1)可以簡記為： 式中， 為零均值白噪聲序列。 由式(8.2)顯而易見，arima模型的實質(zhì)就是差分運算與arma模型的組合。這一關(guān)系意義重大，這說明任何非平穩(wěn)序列只要通過適當(dāng)階數(shù)的差分運算實現(xiàn)差分后平穩(wěn)，就可以對差分后序列進(jìn)行arma模型擬合了。而arma模型的分析方法非常成熟，這意味著對差分平穩(wěn)序列的分析也將是非常簡單、非?？煽康牧?8.2,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,4,例如，設(shè)arima(1,1,1)模型 圖8.1是給出的arima(1,1,1)模型一個模擬數(shù)據(jù)，樣本容量為200，可以看出時間趨勢是非常明顯的。圖8.2是經(jīng)過一階差分得到的數(shù)據(jù)。經(jīng)過一階差分我們看到下降的時間

3、趨勢被去掉，新的序列看起來是平穩(wěn)的,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,5,圖8.1 arima(1,1,1)模型一個模擬數(shù)據(jù) 圖8.2 模擬數(shù)據(jù)的一階差分?jǐn)?shù)據(jù),上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,6,求和自回歸移動平均模型這個名字的由來是因為階差分后序列可以表示為： 式中 ，即差分后序列等于原序 列的若干序列值的加權(quán)和，而對它又可以擬合自回歸移動平均(arma)模型，所以稱它為求和自回歸移動平均模型,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,7,特別地， 當(dāng)d=0時，arima(p,d,q)模型實際上就是arma(p,q)模型； 當(dāng)p=0時，arima(o,d,q)模型可以簡記為iam(d,q)模型； 當(dāng)q=0時，

4、arima(p,d,0)模型可以簡記為ari(p,d)模型. 當(dāng)d=1,p=q=0時，arima(0,1,0)模型為： (8.3) 該模型被稱為隨機(jī)游走(random walk)模型,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,8,隨機(jī)游走模型的產(chǎn)生有一個有趣的典故。它最早于1905年7月由卡爾皮爾遜(karl pearson)在自然雜志上作為一個問題提出：假如有一個醉漢醉得非常嚴(yán)重，完全喪失方向感，把他放在荒郊野外，一段時間之后再去找他，在什么地方找到他的概率最大呢？ 考慮到他完全喪失方向感，那么他第步的位置將是他第步的位置再加一個完全隨機(jī)的位移。用數(shù)學(xué)模型來描述任意時刻這個醉漢可能的位置，即為一個隨即游

5、走模型(8.3,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,9,1905年8月，雷利爵士(lord rayleigh)對卡爾皮爾遜的這個問題作出了解答。他算出這個醉漢離初始點的距離為至的概率為： 且當(dāng)n很大時，該醉漢離初始點的距離服從零均值正態(tài)分布。這意味著，假如有人想去尋找醉漢的話，最好是去初始點附近找他，該地點是醉漢未來位置的無偏估計值。 作為一個最簡單的arima模型，隨機(jī)游走模型目前廣泛應(yīng)用于計量經(jīng)濟(jì)學(xué)領(lǐng)域。傳統(tǒng)的經(jīng)濟(jì)學(xué)家普遍認(rèn)為投機(jī)價格的走勢類似于隨機(jī)游走模型，隨機(jī)游走模型也是有效市場理論(efficient market theory)的核心,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,10,8.1.2 ar

6、ima模型的性質(zhì),一、平穩(wěn)性 假如服從arima(p,d,q)模型： 式中： 記 ， 被稱為廣義自回歸系數(shù)多項式。顯然arima模型的平穩(wěn)性完全由 的根的性質(zhì)決定,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,11,因為階差分后平穩(wěn)，服從arma(p,q)模型，所以不妨設(shè) 則 (8.4) 由式(8.4)容易判斷，arima(p,d,q)模型的廣義自回歸系數(shù)多項式共有p+d個特征根，其中p個在單位圓內(nèi)，d個在單位圓上。因為有d個特征根在單位圓上而非單位圓內(nèi)，所以當(dāng) 時，arima(p,d,q)模型不平穩(wěn),上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,12,二、方差齊性 對于arima(p,d,q)模型，當(dāng) 時，不僅均值非平穩(wěn)，

7、序列方差也非平穩(wěn)。以最簡單的隨機(jī)游走模型arima(0,1,0)為例： 則 這是一個時間的遞增函數(shù)，隨著時間趨向無窮，序列 的方差也趨向無窮。 但1階差分之后， 差分后序列方差齊性,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,13,8.1.3 arima模型建模,在掌握了arma模型建模的方法之后，嘗試使用arima模型對觀察序列建模是一件比較簡單的事情。它遵循如下的操作流程，如下圖所示,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,14,圖8.3 arima模型建模流程,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,15,8.1.4 arima模型預(yù)測,在最小均方誤差預(yù)測原理下，arima模型的預(yù)測和arma模型的預(yù)測方法非常類似。 a

8、rima(p,d,q)模型的一般表示方法為： 和arma模型一樣，也可以用歷史觀測值的線性函數(shù)表示它： 式中 的值由如下等式確定,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,16,如果把 記為廣義自相關(guān)函數(shù)，有 容易驗證 的值滿足如下遞推公式： 式中 那么， 的真實值為,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,17,由于 的不可獲得性，所以 的估計值只能為： 真實值與預(yù)報值之間的均方誤差為： 要使均方誤差最小，當(dāng)且僅當(dāng),上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,18,所以，在均方誤差最小的原則下，期預(yù)報值為： 期預(yù)報誤差為： 真實值等于預(yù)報值加上預(yù)報誤差： 期預(yù)報的方差為,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,19,例8.1 對1950

9、年2005年我國進(jìn)出口貿(mào)易總額數(shù)據(jù)（單位：億元人民幣）序列建立arima模型（數(shù)據(jù)見附錄1.15）1. 對原序列（nx）的分析 (1) 做出1950年2005年我國進(jìn)出口貿(mào)易總額數(shù)據(jù)（nx）的時序圖及自相關(guān)圖，如圖8.4，圖8.5,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,20,圖8.4 圖8.5,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,21,2) 對該序列做單位根檢驗，原假設(shè)：；備擇假設(shè)：，檢驗結(jié)果如圖8.4,圖8.6 根據(jù)圖8.6的檢驗結(jié)果，我們可以認(rèn)為這一序列非平穩(wěn),上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,22,2. 對原序列取對數(shù)并分析 由于這一序列有著非常明顯的指數(shù)趨勢，因此我們對它進(jìn)行取對數(shù)的運算，以消除指數(shù)趨勢

10、的影響，將取對數(shù)后的序列命名為 ，即 。 作出序列 的時序圖與自相關(guān)圖分別如圖8.7，8.8,圖8.7 圖8.8,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,23,依然對序列 做單位根檢驗，檢驗結(jié)果如圖8.9。 圖8.9 根據(jù)這一檢驗結(jié)果，我們看到這一序列依然沒有平穩(wěn)，結(jié)合圖8.7和圖8.8，我們看到在序列 中有著明顯的增長趨勢，因此我們還需要對其進(jìn)行差分處理,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,24,圖8.10 圖8.11 根據(jù)上述結(jié)果，可以認(rèn)為這一序列已經(jīng)平穩(wěn)，接下來，可以針對該序列做進(jìn)一步的建模擬合,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,25,4. 針對平穩(wěn)序列 的建立arma模型 (1) 畫出序列 的自相關(guān)圖，如

11、圖。根據(jù)該圖，我們可以初步判斷該序列的偏自相關(guān)圖一階截尾，而針對自相關(guān)圖并不能馬上做出判斷,圖8.12,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,26,2) 針對序列 我們嘗試幾種不同的模型擬合，比如arma（1，1），arma（1，2），arma（1，3）等。經(jīng)過不斷的嘗試，我們最終選擇了arma（1，6）模型，并且該模型中移動平均部分的系數(shù)只有ma（6）的系數(shù)是顯著的，這樣我們就把1-5階的系數(shù)全部放棄，最終的估計結(jié)果如圖8.13。 圖8.13 通過圖8.11，我們可以看到最終選擇的模型的整體檢驗效果還是良好 的,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,27,5) 對擬合模型后的殘差序列做純隨機(jī)性檢驗，檢驗結(jié)果

12、如圖8.14。 圖8.14 通過這一檢驗，我們看到殘差序列已經(jīng)可以認(rèn)為是一個純白噪聲的序列，說明我們的模型已經(jīng)將有用信息充分提取了。 這一模型的整體擬合效果見圖8.15,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,28,圖8.15 綜合上述分析過程，實際上我們是針對原序列（nx）：1950年2005年我國進(jìn)出口貿(mào)易總額數(shù)據(jù)序列，建立了一個arima（1，1，6）模型進(jìn)行擬合，模型機(jī)構(gòu)如下,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,29,8.2 季節(jié)時間序列模型的分析方法,8.2.1季節(jié)時間序列的重要特征 一、季節(jié)時間序列表示 許多商業(yè)和經(jīng)濟(jì)時間序列都包含季節(jié)現(xiàn)象，例如，冰淇淋的銷量的季 度序列在夏季最高，序列在每年都會

13、重復(fù)這一現(xiàn)象。相應(yīng)的周期為4。 類似地，在美國汽車的月度銷售量和銷售額數(shù)據(jù)在每年的7月和8月也 趨于下降，因為每年這時汽車廠家將會推出新的產(chǎn)品；在西方，玩具 的銷售量在每年12月份會增加，主要是因為圣誕節(jié)的緣故；在中國， 每年農(nóng)歷5月份糯米的銷售量大大地增加，這是因為中國的端午節(jié)有 吃粽子的習(xí)慣。以上三種情況的季節(jié)周期都是12個月。由上面的例子 可以看到，很多的實際問題中，時間序列會顯示出周期變化的規(guī)律， 這種周期性是由于季節(jié)變化或其他物理因素所致，我們稱這類序列為 季節(jié)性序列。單變量的時間序列為了分析方便，可以編制成一個二維 的表格，其中一維表示周期，另一維表示某個周期的一個觀測值，如表8.

14、1所示,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,30,表8.1 單變量時間序列觀測數(shù)據(jù)表 例如，19932000年各月中國社會消費品零售總額序列，是一個月度資料，其周期s=12，起點為1993年1月，具體數(shù)據(jù)見附錄,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,31,二、季節(jié)時間序列的重要特征 季節(jié)性時間序列的重要特征表現(xiàn)為周期性。在一個序列中，如果經(jīng)過s個時間間隔后觀測點呈現(xiàn)出相似性，比如同處于波峰或波谷，我們就說該序列具有以s為周期的周期特性。具有周期特性的序列稱為季節(jié)時間序列，s為周期的長度，不同的季節(jié)時間序列會表現(xiàn)出不同的周期，季度資料的一個周期表現(xiàn)為一年的四個季度，月度資料的周期表現(xiàn)為一年的12各月，周資料表

15、現(xiàn)為一周的7天或5天。 例如，圖8.16的數(shù)據(jù)是1993年1月到2000年12月的中國社會消費品月銷售總額,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,32,圖8.16 1993年1月2000年12月的中國社會消費品月銷售總額 當(dāng)然影響一個季節(jié)性時間序列的因素除了季節(jié)因素外，還存在趨勢變動和不規(guī)則變動等。我們研究季節(jié)性時間序列的目的就是分解影響經(jīng)濟(jì)指標(biāo)變量的季節(jié)因素、趨勢因素和不規(guī)則因素，據(jù)以了解它們對經(jīng)濟(jì)的影響,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,33,8.2.2 季節(jié)時間序列模型,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,34,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,35,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,36,上海財經(jīng)大學(xué) 統(tǒng)計與管

16、理學(xué)院,37,三、常見的隨機(jī)季節(jié)模型,為了讀者學(xué)習(xí)起來方便，這里列舉幾個常見的隨機(jī)季節(jié)模型，并簡介其生成的過程。 在實際問題中，季節(jié)性時間序列所含有的成分不同，記憶性長度各異，因而模型形式也是多種多樣的。這里以季節(jié)周期s=12為例，介紹幾種常見的季節(jié)模型,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,38,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,39,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,40,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,41,8.2.3 季節(jié)性檢驗和季節(jié)模型的建立,檢驗一個時間序列是否具有季節(jié)性是十分必要的，如果一個時間序列季節(jié)性顯著，那么擬合適應(yīng)的季節(jié)時間序列模型是合理的，否則會有欠擬合之嫌。如果不是一個具有顯著季節(jié)性的

17、時間序列，即使是一個月度數(shù)據(jù)資料，也不應(yīng)該擬合季節(jié)性時間序列模型。下面我們討論如何識別一個時間序列的季節(jié)性。 一、季節(jié)性時間序列自相關(guān)函數(shù)和偏自相關(guān)函數(shù)的檢驗 根據(jù)box-jenkins的建模方法，自相關(guān)函數(shù)和偏自相關(guān)函數(shù)的特征是識別非季節(jié)性時間序列的工具。從第七章第二節(jié)的討論已經(jīng)看到季節(jié)性時間序列模型實際上是一種特殊的arima模型，不同的是它的系數(shù)是稀疏的，即部分系數(shù)為零，所以對于乘積季節(jié)模型的階數(shù)識別，基本上可以采用box-jenkins的方法，考察序列樣本自相關(guān)函數(shù)和偏自相關(guān)函數(shù)，從而對季節(jié)性進(jìn)行檢驗,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,42,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,43,上海財經(jīng)大學(xué)

18、 統(tǒng)計與管理學(xué)院,44,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,45,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,46,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,47,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,48,例8.3 繪制1993年1月至2000年12月中國社會消費品零售總額序列的自相關(guān)和偏自相關(guān)圖（圖8.17）。 圖8.17 圖8.17顯示中國社會消費品零售總額月度時間序列的自相關(guān)函數(shù)緩慢下降，且在滯后期為周期倍數(shù)時出現(xiàn)峰值，滯后期為12的自相關(guān)函數(shù)為0.645，滯后期為24的自相關(guān)函數(shù)為0.318，說明該時間序列是一個典型的既有趨勢又有季節(jié)變動的序列，由于該序列不是一個平穩(wěn)的時間序列，所以我們不能由其偏自相關(guān)函數(shù)簡單建立一個自回

19、歸模型，該序列建模必須將序列進(jìn)行差分變化，使其平穩(wěn)化,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,49,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,50,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,51,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,52,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,53,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,54,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,55,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,56,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,57,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,58,eviews軟件介紹(,一、x-12季節(jié)調(diào)整方法簡介 x-12-arima方法最早由美國普查局findley等人在20世紀(jì)90年代左右提出，現(xiàn)已成為對重要時間序列進(jìn)行深入處理和分析的工具

20、，也是處理最常用經(jīng)濟(jì)類指標(biāo)的工具，在美國和加拿大被廣泛使用。其在歐洲統(tǒng)計界也得到推薦，并在包括歐洲中央銀行在內(nèi)的歐洲內(nèi)外的許多中央銀行、統(tǒng)計部門和其他經(jīng)濟(jì)機(jī)構(gòu)被廣泛應(yīng)用。 x-12-arima方法提供了四個方面的改進(jìn)和提高，（1）可選擇季節(jié)、交易日及假日進(jìn)行調(diào)整，包括調(diào)整用戶定義的回歸自變量估計結(jié)果，選擇輔助季節(jié)和趨勢過濾器，以及選擇季節(jié)、趨勢和不規(guī)則因素的分解形式；（2）對各種選項條件下調(diào)整的質(zhì)量和穩(wěn)定性做出新診斷；（3） 對具有arima誤差及可選擇穩(wěn)健估計系數(shù)的線性回歸模型，進(jìn)行廣泛的時間序列建模和模型選擇能力分析；（4）提供一個新的易于分批處理大量時間序列能力的用戶界面。 x-12-a

21、rima方法現(xiàn)已廣泛應(yīng)用于世界各國的中央銀行、統(tǒng)計部門和其他經(jīng)濟(jì)機(jī)構(gòu)，并且已成為對重要時間序列進(jìn)行深入處理和分析的工具,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,59,二、案例：1993-2000年中國社會消費品零售總額月度序列（單位：億元） 通過1993-2000年中國社會消費品零售總額月度序列的時序圖（圖8.16），我們可以觀察到該序列有著很強(qiáng)的季節(jié)特征。通過該序列的自相關(guān)函數(shù)圖（圖8.17）及單位根檢驗結(jié)果（圖8.19）的進(jìn)一步判斷，認(rèn)為該序列非平穩(wěn)，并且有著很強(qiáng)的季節(jié)特征,圖8.19,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,60,上海財經(jīng)大學(xué) 統(tǒng)計與管理學(xué)院,61,首先顯示的是seasonal adjustment（季節(jié)調(diào)整）模塊（圖8.19），該模塊共有5個選項區(qū)。 在x11 method（x11方法）選項區(qū)選multiplicative（乘法模型）。 在seasonal filter（季節(jié)濾子）選項區(qū)選auto（自動）。 在trend filter（趨勢濾子）選項區(qū)選auto（自動）。 在component series to save（保存分量）選項區(qū)選季節(jié)調(diào)整序列（_sa）、季節(jié)因子序列（_sf）、趨勢循環(huán)序列（_tc）、不規(guī)則序列（_ir）四個分量（通過在小方格內(nèi)勾選保存在工作文件