第五章 實證研究
第二節 模型分析
國
立 政 治 大 學
‧
N a tio na
l C h engchi U ni ve rs it y
第二節 模型分析
先針對上一節提到的切割問題作討論。表 5-1(a) 使用未經切割的資料，以向前一期 預測的方式所得到的結果，再將預測結果進行切割。表 5-1(b) 為切割時間序列後的結 果。比較雙季節指數平滑模型 (ds) 的 MAPE 、RMSE 與 MAE ，原始結果明顯優於切割 後結果，Holt-Winters 模型 (hw) 則是部分時段稍微優於切割後結果。
由於每個時段的時間序列都是以日為單位，不用考慮一日內的週期性。從表 5-1(b) 切割後結果可以看出，雙季節指數平滑模型的效果比 Holt-Winters 模型差，在表 5-1(a) 原始結果中，雙季節指數平滑模型的效果明顯優於 Holt-Winters 模型。
表 5-1(a) 原始結果
hwMAE dsMAE hwRMSE dsRMSE hwMAPE dsMAPE hr12 40.7785 12.8239 53.6892 18.3876 75.5364 23.0775 hr13 17.0110 16.5195 22.4340 22.7794 19.3047 18.4896 hr14 21.2681 18.1779 28.4254 24.6962 21.2290 17.6267 hr15 20.9838 17.3916 27.9896 23.9497 18.6574 14.1633 hr16 19.3669 16.9492 26.6916 22.8884 15.7448 12.9095 hr17 19.7582 17.0975 26.9127 24.3946 14.2363 12.6561 hr18 23.8890 18.7278 30.8040 26.1964 17.4071 14.2650
表 5-1(b) 切割後結果
hwMAE dsMAE hwRMSE dsRMSE hwMAPE dsMAPE hr12 12.6006 13.7750 17.9290 18.5899 21.5247 24.5132 hr13 18.4042 24.8268 25.9697 34.4116 19.6669 27.0415 hr14 21.9140 25.2312 31.8695 35.5935 19.4884 23.7716 hr15 23.9387 29.9521 34.9250 41.6261 18.7113 25.0338 hr16 25.3270 32.2892 37.0469 44.0801 18.4891 23.9149 hr17 27.0031 31.2908 39.6480 44.3243 18.2696 21.8242
‧ 國
立 政 治 大 學
‧
N a tio na
l C h engchi U ni ve rs it y
從表 5-2 可以看出，ETS(A,N,N) 跟 ETS(A,A,N) 的預測結果差距不大，因為 ETS(A,N,N) 與 ETS(A,A,N) 兩種模型的差別在於有沒有趨勢影響，剛好這筆來店人數資料沒有明顯 的趨勢，可以預期會有這樣的預測結果。ETS(A,N,N)+DOW 中的 DOW 代表外生變數
“day-of-week” (星期一至星期五設為 0，六日設為 1)，這個模型我們使用上一節所推導 的時間更新 (time update) 與測量更新 (measurement update) 來進行迭代更新。他比前兩 個模型多了一個外生變數 DOW，在預測上有略為改進。以上這三個模型都沒有考慮季 節性。
而 ETS(A,A,A) 考慮了季節性，預測結果明顯改善，再看到雙季節指數平滑模型的 三個預測指標又明顯下降，可見該筆資料的確具有雙季節影響，適合使用雙季節指數 平滑模型來預測。
表 5-2 實驗結果
RMSE MAE MAPE ETS(A,N,N) 47.03177 31.91049 38.11577 ETS(A,A,N) 47.04344 31.90937 38.12456 ETS(A,N,N)+DOW 46.66465 31.59696 37.26808 ETS(A,A,A) 32.43201 23.29364 26.01653 DSHW(7, 49) 23.43904 16.81247 16.16967
在圖5-2 雙季節指數平滑模型部分，除了少部分出現高估之外，在季節循環上有不 錯的預測結果。在 Holt-Winters 模型預測上，出現不少低估的情況。在一開始虛線包住 的區段為平日週三，卻出現假日的人潮，經查詢之後發現該日為 7 月 4 日美國獨立紀念 日。
看到圖5-3，黑色為真實來店人數，在 12 月 30 日附近，因為接近跨年的關係，店家 的人潮快速上升。Holt-Winters 模型在低點經常出現低估的情況，還出現小於 0 的預測 值。雖然雙季節指數平滑模型在高點低點都有高估低估的狀況，但與真實值相距不遠。
‧ 國
立 政 治 大 學
‧
N a tio na
l C h engchi U ni ve rs it y
圖 5-2 Holt-Winters vs. DSHW 1-250
圖 5-3 Holt-Winters vs. DSHW 900-1274
‧ 國
立 政 治 大 學
‧
N a tio na
l C h engchi U ni ve rs it y
圖 5-4 縮小區間範圍至 1200-1274 ，後期出現震幅放大的情形，Holt-Winters 模型在 這段時間的走勢中有晚一期的現象，人數下降時還保持上升的趨勢，雙季節指數平滑 模型在這方面表現穩定，幾乎能抓到上升下降的時機。
圖 5-4 Holt-Winters vs. DSHW 1200-1274
‧ 國
立 政 治 大 學
‧
N a tio na
l C h engchi U ni ve rs it y
圖 5-5 為兩模型的殘差圖，先看到 Holt-Winters 殘差圖，每隔一段時間殘差會異常 變大，在這裡推測可能是因為第二個季節因素導致預測結果不理想，在頭尾有較大 的變動，稍後會與雙季節指數平滑模型一起討論。再看到雙季節指數平滑模型殘差 圖，其他時間的殘差明顯比 Holt-Winters 模型小且穩定，頭尾部分一樣有不規則變動。
Holt-Winters 模型與雙季節指數平滑模型兩種方法的殘差在t ∈ (300, 800)之間有較小變 異，但在頭 (1, 250)、尾 (900, 1274) 部分，兩者的殘差就有明顯的起伏。在雙季節指數 平滑模型殘差圖第 237 的殘差突然增加且為正值，表示在該點的預測出現高估的情況，
推測該商店可能發生什麼事。
圖 5-5 殘差圖
‧ 國
立 政 治 大 學
‧
N a tio na
l C h engchi U ni ve rs it y
從圖 5-6 Holt-Winters 殘差的 ACF 圖可以看出預測結果還有一個季節影響，週期約 為 49，符合上面討論的長週期長度。DSHW 殘差的 ACF 則落在標準內，但仍能看出具 有些微週期性，或許是平日假日的季節影響。
圖 5-6 殘差 ACF
‧ 國
立 政 治 大 學
‧
N a tio na
l C h engchi U ni ve rs it y
第 六 章
結論與建議
在三種傳統的指數平滑模型比較中，Holt-Winters 模型明顯優於簡單指數平滑模型 與 Holt’s linear 模型。接著與 Taylor 的雙季節指數平滑模型作比較，雙季節指數平滑更 優於 Holt-Winters 模型，對於年末突然增長的時段依舊保有不錯的預測結果。
本論文將卡爾曼濾波器的推導邏輯，應用到觀測與狀態方程式中，誤差項具有線性 關係的創新狀態空間模型，推導得到一組新的更新方程式。因為創新狀態空間模型是 指數平滑法的一種模型表示法，目前文獻中，對於創新狀態的更新步驟，基本上仍沿 用指數平滑法的操作步驟，在處理外生變數的情況較為棘手。有了類似卡爾曼濾波器 的更新方法後，處理外生變數就容易許多。未來只要擴展狀態向量的維度，根據指數 平滑法的公式做調整，將相對應的參數放入矩陣 Ht ，就能套用到其他的指數平滑法。
若能套用至雙季節指數平滑模型，加上國定假日或競爭商家的資訊作為外生變數，預 期能得到更好的預測結果。
由於研究資料的時間長度只有一年，無法比較歷年的季節變化，另外從圖 6-1 可以 看出該筆資料可能存在兩個以上季節循環。設定星期一至星期日為一個區段，圖 6-1 中 以虛線作為分段點，可以看出每個區段包含七個尖點，前四個尖點約為 100 人，也就 是星期一至星期四的人數約為 100 人，第五至第七尖點明顯上升，因此一週七天又能 再分成兩個小季節循環。星期六的人數明顯高於其他天，或許能將星期六單獨視為一 個季節循環。
‧ 國
立 政 治 大 學
‧
N a tio na
l C h engchi U ni ve rs it y
圖 6-1 多季節循環
‧ 國
立 政 治 大 學
‧
N a tio na
l C h engchi U ni ve rs it y
參考文獻
[1] S. Tom Au, Guang-Qin Ma, and Shu-Ngai Yeung. Automatic forecasting of double sea-sonal time series with applications on mobility network traffic prediction. 2011.
[2] Robert Goodell Brown. Statistical forecasting for inventory control. McGraw-Hill, New York, 1959.
[3] Phillip G Gould, Anne B Koehler, J Keith Ord, Ralph D Snyder, Rob J Hyndman, and Farshid Vahid-Araghi. Forecasting time series with multiple seasonal patterns. European Journal of Operational Research, 2008.
[4] Charles C Holt. Forecasting seasonals and trends by exponentially weighted moving av-erages. Carnegie Institute of Technology, Graduate school of Industrial Administration, 1957.
[5] Rob J Hyndman and Yeasmin Khandakar. Automatic time series forecasting: the forecast package for R. Journal of Statistical Software, 26(3):1–22, 2008.
[6] Rob J. Hyndman, Anne B. Koehler, J. Keith Ord, and Ralph D. Snyder. Forecasting with Exponential Smoothing:The State Space Approach. Springer, 2008.
[7] Mohamed A. Ismail, Alyaa R. Zahran, and Eman M. Abd El-Metaal. Forecasting hourly electricity demand in egypt using double seasonal autoregressive integrated moving aver-age model. 2015.
[8] R. E. Kalman. A new approach to linear filtering and prediction problems. Journal of
‧ 國
立 政 治 大 學
‧
N a tio na
l C h engchi U ni ve rs it y
[9] R. E. Kalman and R. S. Bucy. New results in linear filtering and new results in linear filtering and prediction theory. Journal of Basic Engineering, 83:95–108, 1961.
[10] Reinaldo Castro Souza, Mônica Barros, and Cristina Vidigal C. de Miranda. Short term load forecasting using double seasonal exponential smoothing and interventions to account for holidays and temperature effects. 2007.
[11] Ivan Svetunkov. smooth: Forecasting Using Smoothing Functions, 2018.
[12] J. W. Taylor. Short-term electricity demand forecasting using double seasonal exponential smoothing. The Journal of the Operational Research Society, 2013.
[13] Greg Welch and Gary Bishop. An introduction to the kalman filter. Technical report, Department of Computer Science University of North Carolina at Chapel Hill, 2006.
[14] Peter R Winters. Forecasting sales by exponentially weighted moving averages. Manage-ment Science, 1960.
[15] 施佩吟. 指數平滑模型應用於來店人數預測之研究. 碩, 國立政治大學, 2015.