1ts.plot(lynx) 2title("Canadian lynx trapping (1821-1934)") 3pause() 4cat("Data are now log-transformed\n") 5log.lynx <- log(lynx) 6sm.ts.pdf(log.lynx, lags=4:5) 7pause() 8sm.autoregression(log.lynx, maxlag=5, se=TRUE) 9pause() 10sm.autoregression(log.lynx, lags=cbind(4,5)) 11