Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- library(forecast)
- x <- AirPassengers
- mod_arima <- auto.arima(x, ic='aicc', stepwise=FALSE)
- mod_exponential <- ets(x, ic='aicc', restrict=FALSE)
- mod_neural <- nnetar(x, p=12, size=25)
- mod_tbats <- tbats(x, ic='aicc', seasonal.periods=12)
- par(mfrow=c(4, 1))
- plot(forecast(mod_arima, 12), include=36)
- plot(forecast(mod_exponential, 12), include=36)
- plot(forecast(mod_neural, 12), include=36)
- plot(forecast(mod_tbats, 12), include=36)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement