【时间序列】最完整的时间序列分析和预测（含实例及代码）
時間序列
在生產和科學研究中，對某一個或者一組變量 ?進行觀察測量，將在一系列時刻所得到的離散數字組成的序列集合，稱之為時間序列。
pandas生成時間序列
過濾數據
重采樣
插值
滑窗
數據平穩性與差分法
pandas生成時間序列
時間戳（timestamp）
固定周期（period）
時間間隔（interval）
truncate過濾
數據重采樣
時間數據由一個頻率轉換到另一個頻率
降采樣
升采樣
插值方法：
ffill 空值取前面的值
bfill 空值取后面的值
interpolate 線性取值
Pandas滑動窗口：
滑動窗口就是能夠根據指定的單位長度來框住時間序列，從而計算框內的統計指標。
相當于一個長度指定的滑塊在刻度尺上面滑動，每滑動一個單位即可反饋滑塊內的數據。
滑動窗口可以使數據更加平穩，浮動范圍會比較小，具有代表性，單獨拿出一個數據可能或多或少會離群，有差異或者錯誤，使用滑動窗口會更規范一些。
數據平穩性與差分法：
基本模型：自回歸移動平均模型(ARMA(p，q))是時間序列中最為重要的模型之一。
它主要由兩部分組成：AR代表p階自回歸過程，MA代表q階移動平均過程。
平穩性
要求經由時間序列所得到的的擬合曲線在未來一段時間內仍能順著現有形態‘慣性’延續下去
即均值和方差不發生明顯變化
ARIMA 模型對時間序列的要求是平穩型。
因此，當你得到一個非平穩的時間序列時，首先要做的即是做時間序列的差分，直到得到一個平穩時間序列。
如果你對時間序列做d次差分才能得到一個平穩序列，那么可以使用ARIMA(p,d,q)模型，其中d是差分次數
ARIMA(p,d,q)
當數據差異特別大時，為了使數據變得平穩些，可以使用差分法
即時間序列在t與t-1時刻的差值
二階差分是指在一階差分基礎上再做一階差分。
自回歸 AR
用自身變量的歷史時間對自己預測
自回歸模型必須滿足平穩性（可以使用差分）
p階自回歸過程公式：? y = u + 求和a*y(t-i) + e
y 是當前值， u是常數項， e 是誤差項（服從獨立同分布） y(t-i)當前預測的值與前P天相關 ,a是自相關系數
自回歸模型限制
用自身來預測
平穩性
自相關性 ?判斷自相關系數！！
只適用于預測與自身前期相關的現象
移動平均模型（MA）
關注自回歸模型中的誤差項的累加
q階自回歸過程的 定義: ?y = u + e + b*e(t-i)
移動平均能有效消除預測中的隨機波動
ARIMA
I表示差分項，1是一階，0是不用做，一般做1階就夠了
原理:將非平穩時間序列轉化為平穩時間序列 ，然后將隱變量僅對它的滯后值以及隨機誤差項的現值和滯后值進行回歸所建立的模型。（滯后指階數）
自相關函數ACF
有序的隨機變量與其自身相比較
ACF反映了同一序列在不同時序的取值之間的相關性
ACF(k) = cov（y(t),y(t-k)）/var(y(t)) ? ? ?[-1,1]
如何確定 pq參數？
利用ACF 和 PCAF
實例操作
主要分為4部分
用pandas處理時序數據
檢驗序數據的穩定性
處理時序數據變成穩定數據
時序數據的預測
1 用pandas導入和處理時序數據
數據集是：航空乘客數量預測例子數據集international-airline-passengers.csv
網上一大推：下載地址：https://github.com/sunlei-1997/ML-DL-datasets/blob/master/international-airline-passengers.csv
import?numpy?as?np import?pandas?as?pd from?datetime?import?datetime import?matplotlib.pylab?as?plt import?tqdm import?statsmodels from?statsmodels.graphics.tsaplots?import?plot_acf,?plot_pacf from?statsmodels.tsa.arima_model?import?ARIMA import?warnings warnings.filterwarnings('ignore')#?讀取數據，pd.read_csv默認生成DataFrame對象，需將其轉換成Series對象 df?=?pd.read_csv('international-airline-passengers.csv',?encoding='utf-8',?index_col='Month') df.index?=?pd.to_datetime(df.index)??#?將字符串索引轉換成時間索引 ts?=?df['Passengers']??#?生成pd.Series對象 ts?=?ts.astype('float') ts.head()Month 1949-01-01 112.0 1949-02-01 118.0 1949-03-01 132.0 1949-04-01 129.0 1949-05-01 121.0 Name: Passengers, dtype: float64ts.indexDatetimeIndex(['1949-01-01', '1949-02-01', '1949-03-01', '1949-04-01','1949-05-01', '1949-06-01', '1949-07-01', '1949-08-01','1949-09-01', '1949-10-01',...'1960-03-01', '1960-04-01', '1960-05-01', '1960-06-01','1960-07-01', '1960-08-01', '1960-09-01', '1960-10-01','1960-11-01', '1960-12-01'],dtype='datetime64[ns]', name='Month', length=144, freq=None)ts['1949-01-01']112.0ts[datetime(1949,1,1)]112.0ts['1949-1'?:?'1949-6']Month 1949-01-01 112.0 1949-02-01 118.0 1949-03-01 132.0 1949-04-01 129.0 1949-05-01 121.0 1949-06-01 135.0 Name: Passengers, dtype: float642 檢驗序數據的穩定性
因為ARIMA模型要求數據是穩定的，所以這一步至關重要。
2.1 判斷數據是穩定的?；趯τ跁r間是常量的幾個統計量：
常量的均值
常量的方差
與時間獨立的自協方差
2.2 python判斷時序數據穩定
平穩性檢驗一般采用觀察法和單位根檢驗法。
觀察法：需計算每個時間段內的平均的數據均值和標準差。
單位根檢驗法：通過Dickey-Fuller Test 進行判斷，大致意思就是在一定置信水平下，對于時序數據假設 Null hypothesis: 非穩定。這是一種常用的單位根檢驗方法，它的原假設為序列具有單位根，即非平穩，對于一個平穩的時序數據，就需要在給定的置信水平上顯著，拒絕原假設。
#?移動平均圖 def?draw_trend(timeseries,?size):f?=?plt.figure(facecolor='white')#?對size個數據進行移動平均rol_mean?=?timeseries.rolling(window=size).mean()#?對size個數據移動平均的方差rol_std?=?timeseries.rolling(window=size).std()timeseries.plot(color='blue',?label='Original')rol_mean.plot(color='red',?label='Rolling?Mean')rol_std.plot(color='black',?label='Rolling?standard?deviation')plt.legend(loc='best')plt.title('Rolling?Mean?&?Standard?Deviation')plt.show()def?draw_ts(timeseries):f?=?plt.figure(facecolor='white')timeseries.plot(color='blue')plt.show()#Dickey-Fuller?test: def?teststationarity(ts,max_lag?=?None):dftest?=?statsmodels.tsa.stattools.adfuller(ts,maxlag=?max_lag)#?對上述函數求得的值進行語義描述dfoutput?=?pd.Series(dftest[0:4],?index=['Test?Statistic','p-value','#Lags?Used','Number?of?Observations?Used'])for?key,value?in?dftest[4].items():dfoutput['Critical?Value?(%s)'%key]?=?valuereturn?dfoutput#查看原始數據的均值和方差 draw_trend(ts,12)通過上圖，我們可以發現數據的移動平均值/標準差有越來越大的趨勢，是不穩定的。接下來我們再看Dickey-Fuller的結果
此時p值為0.991880，說明并不能拒絕原假設。通過DF的數據可以明確的看出，在任何置信度下，數據都不是穩定的。
3 處理時序數據變成穩定數據
數據不穩定的原因主要有以下兩點：
趨勢（trend）-數據隨著時間變化。比如說升高或者降低。
季節性(seasonality)-數據在特定的時間段內變動。比如說節假日，或者活動導致數據的異常。
3.1 對數變換
對數變換主要是為了減小數據的振動幅度，使其線性規律更加明顯，同時保留其他信息。這里強調一下，變換的序列需要滿足大于0，小于0的數據不存在對數變換。
ts_log?=?np.log(ts) draw_trend(ts_log,12)可以看出經過對數變換后，數據值域范圍縮小了，振幅也沒那么大了。
3.2 平滑法
根據平滑技術的不同，平滑法具體分為移動平均法和指數平均法。
移動平均即利用一定時間間隔內的平均值作為某一期的估計值，而指數平均則是用變權的方法來計算均值。
移動平均：
def?draw_moving(timeSeries,?size):f?=?plt.figure(facecolor='white')#?對size個數據進行移動平均rol_mean?=?timeSeries.rolling(window=size).mean()#?對size個數據進行加權移動平均rol_weighted_mean?=?pd.Series.ewm(timeSeries,?span=size)rol_weighted_mean=timeSeries.ewm(halflife=size,min_periods=0,adjust=True,ignore_na=False).mean()timeSeries.plot(color='blue',?label='Original')rol_mean.plot(color='red',?label='Rolling?Mean')rol_weighted_mean.plot(color='black',?label='Weighted?Rolling?Mean')plt.legend(loc='best')plt.title('Rolling?Mean')plt.show() draw_moving(ts_log,12)從上圖可以發現窗口為12的移動平均能較好的剔除年周期性因素，
而指數平均法是對周期內的數據進行了加權，能在一定程度上減小年周期因素，但并不能完全剔除，如要完全剔除可以進一步進行差分操作。
3.3 差分
時間序列最常用來剔除周期性因素的方法當屬差分了，它主要是對等周期間隔的數據進行線性求減。
ARIMA模型相對ARMA模型，僅多了差分操作，ARIMA模型幾乎是所有時間序列軟件都支持的，差分的實現與還原都非常方便。
從上面的統計檢驗結果可以看出，經過12階滑動平均和1階差分后，該序列滿足平穩性的要求了。
3.4 分解
所謂分解就是將時序數據分離成不同的成分。
statsmodels使用的X-11分解過程，它主要將時序數據分離成長期趨勢、季節趨勢和隨機成分。
與其它統計軟件一樣，statsmodels也支持兩類分解模型，加法模型和乘法模型，這里我只實現加法，乘法只需將model的參數設置為"multiplicative"即可。
如圖所示，數據的均值和方差趨于常數，幾乎無波動(看上去比之前的陡峭，但是要注意他的值域只有[-0.05,0.05]之間)
所以直觀上可以認為是穩定的數據。另外DFtest的結果顯示，Statistic值原小于1%時的Critical value，所以在99%的置信度下，數據是穩定的。
4 時序數據的預測
在前面的分析可知，該序列具有明顯的年周期與長期成分。
對于年周期成分我們使用窗口為12的移動平進行處理，對于長期趨勢成分我們采用1階差分來進行處理。
觀察其統計量發現該序列在置信水平為95%的區間下并不顯著，我們對其進行再次一階差分。
再次差分后的序列其自相關具有快速衰減的特點，t統計量在99%的置信水平下是顯著的，這里我不再做詳細說明。
數據平穩后，需要對模型定階，即確定p、q的階數。先畫出ACF,PACF的圖像,代碼如下：
觀察上圖，發現自相關和偏相系數都存在拖尾的特點，并且他們都具有明顯的一階相關性，所以我們設定p=1, q=1。
下面就可以使用ARMA模型進行數據擬合了。（Ps.PACF是判定AR模型階數的，也就是p。ACF是判斷MA階數的，也就是q）
模型擬合完后，我們就可以對其進行預測了。由于ARMA擬合的是經過相關預處理后的數據，故其預測值需要通過相關逆變換進行還原。
我們使用均方根誤差（RMSE）來評估模型樣本內擬合的好壞。利用該準則進行判別時，需要剔除“非預測”數據的影響。
總結
以上是生活随笔為你收集整理的【时间序列】最完整的时间序列分析和预测（含实例及代码）的全部內容，希望文章能夠幫你解決所遇到的問題。
- 上一篇： 【竞赛相关】Kaggle知识点：入门到进
- 下一篇： JDK线程池CompletionServ