forecasting using rsigni˝cant spike at lag 1 in acf suggests non-seasonal ma(1) component....
TRANSCRIPT
![Page 1: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/1.jpg)
Forecasting using R 1
Forecasting using R
Rob J Hyndman
2.5 Seasonal ARIMA models
![Page 2: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/2.jpg)
Outline
1 Backshift notation reviewed
2 Seasonal ARIMA models
3 ARIMA vs ETS
4 Lab session 12
Forecasting using R Backshift notation reviewed 2
![Page 3: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/3.jpg)
Backshift notationA very useful notational device is the backward shiftoperator, B, which is used as follows:
Byt = yt−1 .
In other words, B, operating on yt, has the effect ofshifting the data back one period. Two applications of Bto yt shifts the data back two periods:
B(Byt) = B2yt = yt−2 .
For monthly data, if we wish to shift attention to “thesame month last year,” then B12 is used, and the notationis B12yt = yt−12.
Forecasting using R Backshift notation reviewed 3
![Page 4: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/4.jpg)
Backshift notationA very useful notational device is the backward shiftoperator, B, which is used as follows:
Byt = yt−1 .
In other words, B, operating on yt, has the effect ofshifting the data back one period. Two applications of Bto yt shifts the data back two periods:
B(Byt) = B2yt = yt−2 .
For monthly data, if we wish to shift attention to “thesame month last year,” then B12 is used, and the notationis B12yt = yt−12.
Forecasting using R Backshift notation reviewed 3
![Page 5: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/5.jpg)
Backshift notationA very useful notational device is the backward shiftoperator, B, which is used as follows:
Byt = yt−1 .
In other words, B, operating on yt, has the effect ofshifting the data back one period. Two applications of Bto yt shifts the data back two periods:
B(Byt) = B2yt = yt−2 .
For monthly data, if we wish to shift attention to “thesame month last year,” then B12 is used, and the notationis B12yt = yt−12.
Forecasting using R Backshift notation reviewed 3
![Page 6: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/6.jpg)
Backshift notationA very useful notational device is the backward shiftoperator, B, which is used as follows:
Byt = yt−1 .
In other words, B, operating on yt, has the effect ofshifting the data back one period. Two applications of Bto yt shifts the data back two periods:
B(Byt) = B2yt = yt−2 .
For monthly data, if we wish to shift attention to “thesame month last year,” then B12 is used, and the notationis B12yt = yt−12.
Forecasting using R Backshift notation reviewed 3
![Page 7: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/7.jpg)
Backshift notation
First difference: 1− B.Double difference: (1− B)2.dth-order difference: (1− B)dyt.Seasonal difference: 1− Bm.Seasonal difference followed by a first difference:(1− B)(1− Bm).Multiply terms together together to see thecombined effect:
(1− B)(1− Bm)yt = (1− B− Bm + Bm+1)yt= yt − yt−1 − yt−m + yt−m−1.
Forecasting using R Backshift notation reviewed 4
![Page 8: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/8.jpg)
Backshift notation for ARIMAARMA model:yt = c + φ1yt−1 + · · · + φpyt−p + et + θ1et−1 + · · · + θqet−q
= c + φ1Byt + · · · + φpBpyt + et + θ1Bet + · · · + θqBqetφ(B)yt = c + θ(B)et
where φ(B) = 1− φ1B− · · · − φpBp
and θ(B) = 1 + θ1B + · · · + θqBq.
ARIMA(1,1,1) model:
(1− φ1B) (1− B)yt = c + (1 + θ1B)et↑ ↑ ↑
AR(1) First MA(1)difference
Forecasting using R Backshift notation reviewed 5
![Page 9: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/9.jpg)
Backshift notation for ARIMAARMA model:yt = c + φ1yt−1 + · · · + φpyt−p + et + θ1et−1 + · · · + θqet−q
= c + φ1Byt + · · · + φpBpyt + et + θ1Bet + · · · + θqBqetφ(B)yt = c + θ(B)et
where φ(B) = 1− φ1B− · · · − φpBp
and θ(B) = 1 + θ1B + · · · + θqBq.
ARIMA(1,1,1) model:
(1− φ1B) (1− B)yt = c + (1 + θ1B)et↑ ↑ ↑
AR(1) First MA(1)difference
Forecasting using R Backshift notation reviewed 5
![Page 10: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/10.jpg)
Backshift notation for ARIMA
ARIMA(p, d, q) model:
(1− φ1B− · · · − φpBp) (1− B)dyt = c + (1 + θ1B + · · · + θqBq)et↑ ↑ ↑
AR(p) d differences MA(q)
Forecasting using R Backshift notation reviewed 6
![Page 11: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/11.jpg)
Outline
1 Backshift notation reviewed
2 Seasonal ARIMA models
3 ARIMA vs ETS
4 Lab session 12
Forecasting using R Seasonal ARIMA models 7
![Page 12: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/12.jpg)
Seasonal ARIMA models
ARIMA (p, d, q)︸ ︷︷ ︸ (P,D,Q)m︸ ︷︷ ︸↑ ↑
Non-seasonal part Seasonal part ofof the model of the model
wherem = number of observations per year.
Forecasting using R Seasonal ARIMA models 8
![Page 13: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/13.jpg)
Seasonal ARIMA modelsE.g., ARIMA(1, 1, 1)(1, 1, 1)4 model (without constant)
(1−φ1B)(1−Φ1B4)(1−B)(1−B4)yt = (1+θ1B)(1+Θ1B4)et.
6 6 6 6 6 6(Non-seasonal
AR(1)
)(SeasonalAR(1)
)(Non-seasonaldifference
)(Seasonaldifference
)(Non-seasonal
MA(1)
)(SeasonalMA(1)
)
Forecasting using R Seasonal ARIMA models 9
![Page 14: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/14.jpg)
Seasonal ARIMA modelsE.g., ARIMA(1, 1, 1)(1, 1, 1)4 model (without constant)
(1−φ1B)(1−Φ1B4)(1−B)(1−B4)yt = (1+θ1B)(1+Θ1B4)et.
6 6 6 6 6 6(Non-seasonal
AR(1)
)(SeasonalAR(1)
)(Non-seasonaldifference
)(Seasonaldifference
)(Non-seasonal
MA(1)
)(SeasonalMA(1)
)
Forecasting using R Seasonal ARIMA models 9
![Page 15: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/15.jpg)
Seasonal ARIMA modelsE.g., ARIMA(1, 1, 1)(1, 1, 1)4 model (without constant)
(1−φ1B)(1−Φ1B4)(1−B)(1−B4)yt = (1+θ1B)(1+Θ1B4)et.
6 6 6 6 6 6(Non-seasonal
AR(1)
)(SeasonalAR(1)
)(Non-seasonaldifference
)(Seasonaldifference
)(Non-seasonal
MA(1)
)(SeasonalMA(1)
)
Forecasting using R Seasonal ARIMA models 9
![Page 16: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/16.jpg)
Seasonal ARIMA modelsE.g., ARIMA(1, 1, 1)(1, 1, 1)4 model (without constant)
(1−φ1B)(1−Φ1B4)(1−B)(1−B4)yt = (1+θ1B)(1+Θ1B4)et.
All the factors can be multiplied out and the general modelwritten as follows:
yt = (1 + φ1)yt−1 − φ1yt−2 + (1 + Φ1)yt−4− (1 + φ1 + Φ1 + φ1Φ1)yt−5 + (φ1 + φ1Φ1)yt−6−Φ1yt−8 + (Φ1 + φ1Φ1)yt−9 − φ1Φ1yt−10+ et + θ1et−1 + Θ1et−4 + θ1Θ1et−5.
Forecasting using R Seasonal ARIMA models 10
![Page 17: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/17.jpg)
Common ARIMA models
In the US Census Bureau uses the following models mostoften:
ARIMA(0,1,1)(0,1,1)m with log transformationARIMA(0,1,2)(0,1,1)m with log transformationARIMA(2,1,0)(0,1,1)m with log transformationARIMA(0,2,2)(0,1,1)m with log transformationARIMA(2,1,2)(0,1,1)m with no transformation
Forecasting using R Seasonal ARIMA models 11
![Page 18: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/18.jpg)
Seasonal ARIMA modelsThe seasonal part of an AR or MA model will be seen in theseasonal lags of the PACF and ACF.
ARIMA(0,0,0)(0,0,1)12 will show:
a spike at lag 12 in the ACF but no other significantspikes.The PACF will show exponential decay in the seasonallags; that is, at lags 12, 24, 36, . . . .
ARIMA(0,0,0)(1,0,0)12 will show:
exponential decay in the seasonal lags of the ACFa single significant spike at lag 12 in the PACF.
Forecasting using R Seasonal ARIMA models 12
![Page 19: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/19.jpg)
European quarterly retail trade
autoplot(euretail) + xlab("Year") + ylab("Retail index")
92
96
100
2000 2005 2010Year
Ret
ail i
ndex
Forecasting using R Seasonal ARIMA models 13
![Page 20: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/20.jpg)
European quarterly retail trade
ggtsdisplay(diff(euretail,4))
−2
0
2
2000 2005 2010Time
x
diff(euretail, 4)
0.0
0.5
1.0
4 8 12 16Lag
AC
F
0.0
0.5
1.0
4 8 12 16Lag
PAC
F
Forecasting using R Seasonal ARIMA models 14
![Page 21: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/21.jpg)
European quarterly retail trade
ggtsdisplay(diff(diff(euretail,4)))
−2
−1
0
1
2000 2005 2010Time
x
diff(diff(euretail, 4))
−0.6
−0.4
−0.2
0.0
0.2
4 8 12 16Lag
AC
F
−0.6
−0.4
−0.2
0.0
0.2
4 8 12 16Lag
PAC
F
Forecasting using R Seasonal ARIMA models 15
![Page 22: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/22.jpg)
European quarterly retail trade
d = 1 and D = 1 seems necessary.Significant spike at lag 1 in ACF suggests non-seasonalMA(1) component.Significant spike at lag 4 in ACF suggests seasonalMA(1) component.Initial candidate model: ARIMA(0,1,1)(0,1,1)4.We could also have started with ARIMA(1,1,0)(1,1,0)4.
Forecasting using R Seasonal ARIMA models 16
![Page 23: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/23.jpg)
European quarterly retail trade
fit <- Arima(euretail, order=c(0,1,1),seasonal=c(0,1,1))
ggtsdisplay(residuals(fit))
−1.0
−0.5
0.0
0.5
1.0
2000 2005 2010Time
x
residuals(fit)
−0.2
0.0
0.2
4 8 12 16Lag
AC
F
−0.2
0.0
0.2
4 8 12 16Lag
PAC
F
Forecasting using R Seasonal ARIMA models 17
![Page 24: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/24.jpg)
European quarterly retail trade
ACF and PACF of residuals show significant spikes atlag 2, and maybe lag 3.AICc of ARIMA(0,1,2)(0,1,1)4 model is 74.36.AICc of ARIMA(0,1,3)(0,1,1)4 model is 68.53.
fit <- Arima(euretail, order=c(0,1,3),seasonal=c(0,1,1))
ggtsdisplay(residuals(fit))
Forecasting using R Seasonal ARIMA models 18
![Page 25: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/25.jpg)
European quarterly retail trade
ACF and PACF of residuals show significant spikes atlag 2, and maybe lag 3.AICc of ARIMA(0,1,2)(0,1,1)4 model is 74.36.AICc of ARIMA(0,1,3)(0,1,1)4 model is 68.53.
fit <- Arima(euretail, order=c(0,1,3),seasonal=c(0,1,1))
ggtsdisplay(residuals(fit))
Forecasting using R Seasonal ARIMA models 18
![Page 26: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/26.jpg)
European quarterly retail trade
−1.0
−0.5
0.0
0.5
2000 2005 2010Time
x
residuals(fit)
−0.2
−0.1
0.0
0.1
0.2
4 8 12 16Lag
AC
F
−0.2
−0.1
0.0
0.1
0.2
4 8 12 16Lag
PAC
F
Forecasting using R Seasonal ARIMA models 19
![Page 27: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/27.jpg)
European quarterly retail trade
res <- residuals(fit)Box.test(res, lag=16, fitdf=4, type="Ljung")
#### Box-Ljung test#### data: res## X-squared = 7.0105, df = 12, p-value = 0.8569
Forecasting using R Seasonal ARIMA models 20
![Page 28: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/28.jpg)
European quarterly retail trade
autoplot(forecast(fit, h=12))
90
95
100
2000 2005 2010 2015Time
y
level
80
95
Forecasts from ARIMA(0,1,3)(0,1,1)[4]
Forecasting using R Seasonal ARIMA models 21
![Page 29: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/29.jpg)
European quarterly retail trade
auto.arima(euretail)
## Series: euretail## ARIMA(1,1,2)(0,1,1)[4]#### Coefficients:## ar1 ma1 ma2 sma1## 0.7345 -0.4655 0.2162 -0.8413## s.e. 0.2239 0.1995 0.2096 0.1869#### sigma^2 estimated as 0.1592: log likelihood=-29.69## AIC=69.37 AICc=70.51 BIC=79.76
Forecasting using R Seasonal ARIMA models 22
![Page 30: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/30.jpg)
European quarterly retail trade
auto.arima(euretail, stepwise=FALSE, approximation=FALSE)
## Series: euretail## ARIMA(0,1,3)(0,1,1)[4]#### Coefficients:## ma1 ma2 ma3 sma1## 0.2625 0.3697 0.4194 -0.6615## s.e. 0.1239 0.1260 0.1296 0.1555#### sigma^2 estimated as 0.1564: log likelihood=-28.7## AIC=67.4 AICc=68.53 BIC=77.78
Forecasting using R Seasonal ARIMA models 23
![Page 31: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/31.jpg)
Cortecosteroid drug sales
0.50
0.75
1.00
1.25
−0.8
−0.4
0.0
H02 sales (m
illion scripts)Log H
02 sales
1995 2000 2005Year
Forecasting using R Seasonal ARIMA models 24
![Page 32: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/32.jpg)
Cortecosteroid drug sales
0.0
0.2
0.4
1995 2000 2005Year
x
Seasonally differenced H02 scripts
−0.2
0.0
0.2
0.4
12 24 36Lag
AC
F
−0.2
0.0
0.2
0.4
12 24 36Lag
PAC
F
Forecasting using R Seasonal ARIMA models 25
![Page 33: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/33.jpg)
Cortecosteroid drug sales
Choose D = 1 and d = 0.Spikes in PACF at lags 12 and 24 suggest seasonalAR(2) term.Spikes in PACF suggests possible non-seasonal AR(3)term.Initial candidate model: ARIMA(3,0,0)(2,1,0)12.
Forecasting using R Seasonal ARIMA models 26
![Page 34: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/34.jpg)
Cortecosteroid drug sales
Model AICc
ARIMA(3,0,0)(2,1,0)12 −475.12ARIMA(3,0,1)(2,1,0)12 −476.31ARIMA(3,0,2)(2,1,0)12 −474.88ARIMA(3,0,1)(1,1,0)12 −463.40ARIMA(3,0,1)(0,1,1)12 −483.67ARIMA(3,0,1)(0,1,2)12 −485.48ARIMA(3,0,1)(1,1,1)12 −484.25
Forecasting using R Seasonal ARIMA models 27
![Page 35: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/35.jpg)
Cortecosteroid drug sales
(fit <- Arima(h02, order=c(3,0,1), seasonal=c(0,1,2),lambda=0))
## Series: h02## ARIMA(3,0,1)(0,1,2)[12]## Box Cox transformation: lambda= 0#### Coefficients:## ar1 ar2 ar3 ma1 sma1 sma2## -0.1603 0.5481 0.5678 0.3827 -0.5222 -0.1768## s.e. 0.1636 0.0878 0.0942 0.1895 0.0861 0.0872#### sigma^2 estimated as 0.004278: log likelihood=250.04## AIC=-486.08 AICc=-485.48 BIC=-463.28
Forecasting using R Seasonal ARIMA models 28
![Page 36: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/36.jpg)
Cortecosteroid drug sales
ggtsdisplay(residuals(fit))
−0.2
−0.1
0.0
0.1
0.2
1995 2000 2005Time
x
residuals(fit)
−0.1
0.0
0.1
0.2
12 24 36Lag
AC
F
−0.1
0.0
0.1
0.2
12 24 36Lag
PAC
F
Forecasting using R Seasonal ARIMA models 29
![Page 37: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/37.jpg)
Cortecosteroid drug sales
Box.test(residuals(fit), lag=36, fitdf=6,type="Ljung")
#### Box-Ljung test#### data: residuals(fit)## X-squared = 50.712, df = 30, p-value = 0.01045
Forecasting using R Seasonal ARIMA models 30
![Page 38: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/38.jpg)
Cortecosteroid drug salesfit <- auto.arima(h02, lambda=0, d=0, D=1, max.order=9,
stepwise=FALSE, approximation=FALSE)ggtsdisplay(residuals(fit))
−0.2
−0.1
0.0
0.1
1995 2000 2005Time
x
residuals(fit)
−0.1
0.0
0.1
0.2
12 24 36Lag
AC
F
−0.1
0.0
0.1
0.2
12 24 36Lag
PAC
F
Forecasting using R Seasonal ARIMA models 31
![Page 39: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/39.jpg)
Cortecosteroid drug sales
Box.test(residuals(fit), lag=36, fitdf=8,type="Ljung")
#### Box-Ljung test#### data: residuals(fit)## X-squared = 44.766, df = 28, p-value = 0.02329
Forecasting using R Seasonal ARIMA models 32
![Page 40: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/40.jpg)
Cortecosteroid drug salesModel RMSE
ARIMA(3,0,0)(2,1,0)[12] 0.0661ARIMA(3,0,1)(2,1,0)[12] 0.0646ARIMA(3,0,2)(2,1,0)[12] 0.0645ARIMA(3,0,1)(1,1,0)[12] 0.0679ARIMA(3,0,1)(0,1,1)[12] 0.0644ARIMA(3,0,1)(0,1,2)[12] 0.0622ARIMA(3,0,1)(1,1,1)[12] 0.0630ARIMA(4,0,3)(0,1,1)[12] 0.0648ARIMA(3,0,3)(0,1,1)[12] 0.0639ARIMA(4,0,2)(0,1,1)[12] 0.0648ARIMA(3,0,2)(0,1,1)[12] 0.0644ARIMA(2,1,3)(0,1,1)[12] 0.0634ARIMA(2,1,4)(0,1,1)[12] 0.0632ARIMA(2,1,5)(0,1,1)[12] 0.0640
Forecasting using R Seasonal ARIMA models 33
![Page 41: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/41.jpg)
Cortecosteroid drug sales
Models with lowest AICc values tend to give slightlybetter results than the other models.AICc comparisons must have the same orders ofdifferencing. But RMSE test set comparisons caninvolve any models.No model passes all the residual tests.Use the best model available, even if it does not passall tests.In this case, the ARIMA(3,0,1)(0,1,2)12 has the lowestRMSE value and the best AICc value for models withfewer than 6 parameters.
Forecasting using R Seasonal ARIMA models 34
![Page 42: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/42.jpg)
Cortecosteroid drug sales
fit <- Arima(h02, order=c(3,0,1), seasonal=c(0,1,2),lambda=0)
autoplot(forecast(fit)) +ylab("H02 sales (million scripts)") + xlab("Year")
0.5
1.0
1.5
1995 2000 2005 2010Year
H02
sal
es (
mill
ion
scrip
ts)
level
80
95
Forecasts from ARIMA(3,0,1)(0,1,2)[12]
Forecasting using R Seasonal ARIMA models 35
![Page 43: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/43.jpg)
Outline
1 Backshift notation reviewed
2 Seasonal ARIMA models
3 ARIMA vs ETS
4 Lab session 12
Forecasting using R ARIMA vs ETS 36
![Page 44: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/44.jpg)
ARIMA vs ETS
Myth that ARIMA models are more general thanexponential smoothing.Linear exponential smoothing models all special casesof ARIMA models.Non-linear exponential smoothing models have noequivalent ARIMA counterparts.Many ARIMA models have no exponential smoothingcounterparts.ETS models all non-stationary. Models withseasonality or non-damped trend (or both) have twounit roots; all other models have one unit root.
Forecasting using R ARIMA vs ETS 37
![Page 45: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/45.jpg)
Equivalences
ETS model ARIMA model Parameters
ETS(A,N,N) ARIMA(0,1,1) θ1 = α− 1ETS(A,A,N) ARIMA(0,2,2) θ1 = α + β − 2
θ2 = 1− αETS(A,A,N) ARIMA(1,1,2) φ1 = φ
θ1 = α + φβ − 1− φθ2 = (1− α)φ
ETS(A,N,A) ARIMA(0,0,m)(0,1,0)mETS(A,A,A) ARIMA(0,1,m + 1)(0,1,0)mETS(A,A,A) ARIMA(1,0,m + 1)(0,1,0)m
Forecasting using R ARIMA vs ETS 38
![Page 46: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/46.jpg)
Outline
1 Backshift notation reviewed
2 Seasonal ARIMA models
3 ARIMA vs ETS
4 Lab session 12
Forecasting using R Lab session 12 39
![Page 47: Forecasting using RSigni˝cant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Signi˝cant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model:](https://reader034.vdocument.in/reader034/viewer/2022050520/5fa4262e778d3a7f9c10d2ec/html5/thumbnails/47.jpg)
Lab Session 12
Forecasting using R Lab session 12 40