ts()
auto.arima()
forecast()
ets()
accuracy()
Unlike for an exponential smoothing model, missing values are ok when fitting a seasonal ARIMA model
fulldat <- window(chinookts, c(1990,1), c(1999,12)) fit <- forecast::auto.arima(fulldat) fr <- forecast::forecast(fit, h=12) plot(fr)