require(graphics)
work <- diff(WWWusage)
par(mfrow = c(2, 1)); plot(WWWusage); plot(work)
## Not run:
# require(stats)
# aics <- matrix(, 6, 6, dimnames = list(p = 0:5, q = 0:5))
# for(q in 1:5) aics[1, 1+q] <- arima(WWWusage, c(0, 1, q),
# optim.control = list(maxit = 500))$aic
# for(p in 1:5)
# for(q in 0:5) aics[1+p, 1+q] <- arima(WWWusage, c(p, 1, q),
# optim.control = list(maxit = 500))$aic
# round(aics - min(aics, na.rm = TRUE), 2)
# ## End(Not run)
Run the code above in your browser using DataCamp Workspace