如果只是上图可用如下过程
#导入数据集为a
b<-ts(a[,1], frequency=12, start=c(1985,1))
c<-diff(b,12)
c
plot(b)
plot(c)