library(TSA)
y=qar.sim(n=100,const=0.0,phi0=3.97, phi1=-3.97,sigma=0,init=.377)
plot(x=1:100,y=y,type='l',ylab=expression(Y[t]),xlab='t')
y
cbind(time(y),y)
rbind(time(y),y)
plot(y=y,x=zlag(y,d=1))
plot(y=y,x=zlag(y,d=2))
plot(y=y,x=zlag(y,d=3))
plot(y=y,x=zlag(y,d=50))