covariance”) #保存协方差
2、同时绘制两组数据的时序图
d=read.csv(“double.csv”,header=F)
double=ts(d,start=1964,freq=1)
plot...(double, plot.type = “multiple”) #两组数据两个图
plot(double, plot.type = “single”) #两组数据一个图
plot(double, plot.type...acf.3(x) #同时绘制3个相关图,acf函数的扩展
ur.df.01(x) #进行单位根检验,得到更加舒服的结果
tsdiag2(x) #返回x的
arma.choose(x,ari=3,mai...=3) #选择合适的AR和MA,基于包tseries的arma函数
#########################附属自编函数
#…
acf.3=function(x,lag.max=10,…){...type=”covariance”)
acf(x,lag.max= lag.max,type=”partial”)
par(ol)
}
#…
#…类似于tsgiag函数的扩展
tsdiag2=function