Produção mensal média por vaca nos EUA, 01/1994 - 12/2005. 

library(tseries)
data(milk, package='TSA')
par(mfrow=c(1,1), mar=c(3,3,3,1), mgp=c(1.6,.6,0), pch=19)
plot(milk,xlab='Anos',main='Produção mensal média por vaca nos EUA entre 1994 e 2005',
  ylab='Mililitro (mL)',type='l')
points(milk, cex = 0.5)
grid()

Correlograma 

par(mfrow=c(1,1), mar=c(4,4,3,1), mgp=c(1.6,.6,0), pch=19)
plot(acf(c(milk),plot=F)[1:22],xlab='Defasagem',main=' ',ylab='Autocorrelação')
grid()

Observamos a evidente interferencia da tendência crescente desta série.

Considerando a tendência. 

par(mfrow=c(1,1), mar=c(3,3,3,1), mgp=c(1.6,.6,0), pch=19)
plot(milk,xlab='Anos',main='Produção mensal média por vaca nos EUA entre 1994 e 2005',
     ylab='Mililitro (mL)',type='l')
points(milk, cex = 0.5)
grid()
tempo=time(milk)
tendencia=ts(fitted(loess(milk~tempo)), start=c(1994,1),frequency=12)
lines(tendencia, col='turquoise', lwd=2)

Série sem tendência 

leite=milk-tendencia
par(mfrow=c(1,1), mar=c(3,3,3,1), mgp=c(1.6,.6,0), pch=19)
plot(leite,xlab='Anos',main='Produção ajustada mensal média por vaca nos EUA entre 1994 e 2005',
     ylab='Observado-Tendência (mL)',type='l')
points(leite, cex = 0.5)
abline(h = mean(leite), col = "red", lwd = 4)
grid()

Correlograma 

par(mfrow=c(1,1), mar=c(3,3,3,1), mgp=c(1.6,.6,0), pch=19)
plot(acf(c(leite),plot=F)[1:22],xlab='Defasagem',main=' ',ylab='Autocorrelação')
grid()