5.12 Forecast using a seasonal model

Forecasting works the same using the forecast() function.

fr <- forecast::forecast(fit, h = 12)
plot(fr)
points(testdat)