library(TSA) set.seed(12317) y=arima.sim(model=list(ar=.8),n=100) plot(y,type='l') acf(y) acf(y)$acf .69^(1:20) pacf(y)