使用季節差分達到定態
對於具有季節性的資料,比起連續期間之間取差分,更常見的是在相鄰年份的同一季之間取差分。舉例來說,若是季資料,就會用某一年的 Q1 減去前一年的 Q1。這稱為季節差分。
有時你需要在同一條序列上同時套用季節差分與 lag-1 差分,也就是在已差分的序列上再做一次差分。
在本練習中,你會同時使用差分與轉換,讓時間序列看起來更接近定態。這裡的資料集是 h02,包含澳洲皮質類固醇藥品每月銷售量,共 17 年。它已載入到你的工作環境中。
本練習屬於課程
R 的時間序列預測
練習說明
- 繪製資料以觀察趨勢與季節性。
- 對
h02先取log(),再使用diff()並搭配合適的lag值進行季節差分。將結果指定給difflogh02。 - 繪製取對數且已差分後的資料。
- 由於
difflogh02看起來仍非定態,請再對它做一次 lag-1 差分(對自身再套用一次diff()),並將結果存為ddifflogh02。繪製此序列。 - 使用適當的函式繪製最終
ddifflogh02序列的 ACF。
動手互動練習
試著完成這個範例程式碼,體驗一下這個練習。
# Plot the data
___
# Take logs and seasonal differences of h02
difflogh02 <- diff(log(___), lag = ___)
# Plot difflogh02
___
# Take another difference and plot
ddifflogh02 <- ___
___
# Plot ACF of ddifflogh02
___